Selection of long oligonucleotides for gene expression microarrays using weighted rank-sum strategy

Sterile protection against P. vivax malaria by repeated blood stage infection in the Aotus monkey model

The malaria parasite Plasmodium vivax remains a major global public health challenge, and no vaccine is approved for use in humans. Here, we assessed whether P. vivax strain-transcendent immunity can be achieved by repeated infection in Aotus monkeys. Sterile immunity was achieved after two homologous infections, whereas subsequent heterologous challenge provided only partial protection. IgG levels based on P. vivax lysate ELISA and protein microarray increased with repeated infections and correlated with the level of homologous protection. Parasite transcriptional profiles provided no evidence of major antigenic switching upon homologous or heterologous challenge. However, we observed significant sequence diversity and transcriptional differences in the P. vivax core gene repertoire between the two strains used in the study, suggesting that partial protection upon heterologous challenge is due to molecular differences between strains rather than immune evasion by antigenic switching. Our study demonstrates that sterile immunity against P. vivax can be achieved by repeated homologous blood stage infection in Aotus monkeys, thus providing a benchmark to test the efficacy of candidate blood stage P. vivax malaria vaccines.

Cyclical regression covariates remove the major confounding effect of cyclical developmental gene expression with strain-specific drug response in the malaria parasite Plasmodium falciparum

Background

The cyclical nature of gene expression in the intraerythrocytic development cycle (IDC) of the malaria parasite, Plasmodium falciparum, confounds the accurate detection of specific transcriptional differences, e.g. as provoked by the development of drug resistance. In lab-based studies, P. falciparum cultures are synchronized to remove this confounding factor, but the rapid detection of emerging resistance to artemisinin therapies requires rapid analysis of transcriptomes extracted directly from clinical samples. Here we propose the use of cyclical regression covariates (CRC) to eliminate the major confounding effect of developmentally driven transcriptional changes in clinical samples. We show that elimination of this confounding factor reduces both Type I and Type II errors and demonstrate the effectiveness of this approach using a published dataset of 1043 transcriptomes extracted directly from patient blood samples with different patient clearance times after treatment with artemisinin.

Results

We apply this method to two publicly available datasets and demonstrate its ability to reduce the confounding of differences in transcript levels due to misaligned intraerythrocytic development time. Adjusting the clinical 1043 transcriptomes dataset with CRC results in detection of fewer functional categories than previously reported from the same data set adjusted using other methods. We also detect mostly the same functional categories, but observe fewer genes within these categories. Finally, the CRC method identifies genes in a functional category that was absent from the results when the dataset was adjusted using other methods. Analysis of differential gene expression in the clinical data samples that vary broadly for developmental stage resulted in the detection of far fewer transcripts in fewer functional categories while, at the same time, identifying genes in two functional categories not present in the unadjusted data analysis. These differences are consistent with the expectation that CRC reduces both false positives and false negatives with the largest effect on datasets from samples with greater variance in developmental stage.

Conclusions

Cyclical regression covariates have immediate application to parasite transcriptome sequencing directly from clinical blood samples and to cost-constrained in vitro experiments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08281-y.

Simultaneous genome-wide gene expression and transcript isoform profiling in the human malaria parasite

Gene expression DNA microarrays have been vital for characterizing whole-genome transcriptional profiles. Nevertheless, their effectiveness relies heavily on the accuracy of genome sequences, the annotation of gene structures, and the sequence-dependent performance of individual probes. Currently available gene expression arrays for the malaria parasite Plasmodium falciparum rely on an average of 2 probes per gene, usually positioned near the 3′ end of genes; consequently, existing designs are prone to measurement bias and cannot capture complexities such as the occurrence of transcript isoforms arising from alternative splicing or alternative start/ stop sites. Here, we describe two novel gene expression arrays with exon-focused probes designed with an average of 12 and 20 probes spanning each gene. This high probe density minimizes signal noise inherent in probe-to-probe sequence-dependent hybridization intensity. We demonstrate that these exon arrays accurately profile genome-wide expression, comparing favorably to currently available arrays and RNA-seq profiling, and can detect alternatively spliced transcript isoforms as well as non-coding RNAs (ncRNAs). Of the 964 candidate alternate splicing events from published RNA-seq studies, 162 are confirmed using the exon array. Furthermore, the exon array predicted 330 previously unidentified alternate splicing events. Gene expression microarrays continue to offer a cost-effective alternative to RNA-seq for the simultaneous monitoring of gene expression and alternative splicing events. Microarrays may even be preferred in some cases due to their affordability and the rapid turn-around of results when hundreds of samples are required for fine-scale systems biology investigations, including the monitoring of the networks of gene co-expression in the emergence of drug resistance.

Histone 4 lysine 8 acetylation regulates proliferation and host-pathogen interaction in Plasmodium falciparum

Background

The dynamics of histone modifications in Plasmodium falciparum indicates the existence of unique mechanisms that link epigenetic factors with transcription. Here, we studied the impact of acetylated histone code on transcriptional regulation during the intraerythrocytic developmental cycle (IDC) of P. falciparum.

Results

Using a dominant-negative transgenic approach, we showed that acetylations of histone H4 play a direct role in transcription. Specifically, these histone modifications mediate an inverse transcriptional relationship between the factors of cell proliferation and host–parasite interaction. Out of the four H4 acetylations, H4K8ac is likely the rate-limiting, regulatory step, which modulates the overall dynamics of H4 posttranslational modifications. H4K8ac exhibits maximum responsiveness to HDAC inhibitors and has a highly dynamic distribution pattern along the genome of P. falciparum during the IDC. Moreover, H4K8ac functions mainly in the euchromatin where its occupancy shifts from intergenic regions located upstream of 5′ end of open reading frame into the protein coding regions. This shift is directly or indirectly associated with transcriptional activities at the corresponding genes. H4K8ac is also active in the heterochromatin where it stimulates expression of the main antigenic gene family (var) by its presence in the promoter region.

Conclusions

Overall, we demonstrate that H4K8ac is a potential major regulator of chromatin-linked transcriptional changes during P. falciparum life cycle which is associated not only with euchromatin but also with heterochromatin environment. This is potentially a highly significant finding that suggests a regulatory connection between growth and parasite–host interaction both of which play a major role in malaria parasite virulence.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-017-0147-z) contains supplementary material, which is available to authorized users.

Integrated analysis of the Plasmodium species transcriptome

The genome sequence available for different Plasmodium species is a valuable resource for understanding malaria parasite biology. However, comparative genomics on its own cannot fully explain all the species-specific differences which suggests that other genomic aspects such as regulation of gene expression play an important role in defining species-specific characteristics. Here, we developed a comprehensive approach to measure transcriptional changes of the evolutionary conserved syntenic orthologs during the intraerythrocytic developmental cycle across six Plasmodium species. We show significant transcriptional constraint at the mid-developmental stage of Plasmodium species while the earliest stages of parasite development display the greatest transcriptional variation associated with critical functional processes. Modeling of the evolutionary relationship based on changes in transcriptional profile reveal a phylogeny pattern of the Plasmodium species that strictly follows its mammalian hosts. In addition, the work shows that transcriptional conserved orthologs represent potential future targets for anti-malaria intervention as they would be expected to carry out key essential functions within the parasites. This work provides an integrated analysis of orthologous transcriptome, which aims to provide insights into the Plasmodium evolution thereby establishing a framework to explore complex pathways and drug discovery in Plasmodium species with broad host range.

•

Comparison of variations in mRNA abundance across six different Plasmodium species.

•

Transcriptional conservation and divergence of Plasmodium syntenic orthologs.

•

Pattern of Plasmodium transcriptome evolution are established.

•

Transcriptionally conserved orthologs represent attractive intervention targets.

Malaria remains a major public health concern despite global efforts in the fight against this disease. The intraerythrocytic stage of the malaria parasites is currently in the spotlight for anti-malarial intervention and vaccine targets. The primary goal of this study is to generate a comprehensive and directly comparable transcriptome dataset across multiple Plasmodium species originating from different hosts. We establish that specific pathways and intraerythrocytic stages are more transcriptionally diverged than others, reflecting transcriptional evolutionary diversity. We further propose a panel of transcriptionally conserved genes as potential drug targets.

DNA damage regulation and its role in drug-related phenotypes in the malaria parasites

DNA of malaria parasites, Plasmodium falciparum, is subjected to extraordinary high levels of genotoxic insults during its complex life cycle within both the mosquito and human host. Accordingly, most of the components of DNA repair machinery are conserved in the parasite genome. Here, we investigated the genome-wide responses of P. falciparum to DNA damaging agents and provided transcriptional evidence of the existence of the double strand break and excision repair system. We also showed that acetylation at H3K9, H4K8, and H3K56 play a role in the direct and indirect response to DNA damage induced by an alkylating agent, methyl methanesulphonate (MMS). Artemisinin, the first line antimalarial chemotherapeutics elicits a similar response compared to MMS which suggests its activity as a DNA damaging agent. Moreover, in contrast to the wild-type P. falciparum, two strains (Dd2 and W2) previously shown to exhibit a mutator phenotype, fail to induce their DNA repair upon MMS-induced DNA damage. Genome sequencing of the two mutator strains identified point mutations in 18 DNA repair genes which may contribute to this phenomenon.

Genome-wide transcriptome profiling reveals functional networks involving the Plasmodium falciparum drug resistance transporters PfCRT and PfMDR1

BACKGROUND

The acquisition of multidrug resistance by Plasmodium falciparum underscores the need to understand the underlying molecular mechanisms so as to counter their impact on malaria control. For the many antimalarials whose mode of action relates to inhibition of heme detoxification inside infected erythrocytes, the digestive vacuole transporters PfCRT and PfMDR1 constitute primary resistance determinants.

RESULTS

Using gene expression microarrays over the course of the parasite intra-erythrocytic developmental cycle, we compared the transcriptomic profiles between P. falciparum strains displaying mutant or wild-type pfcrt or varying in pfcrt or pfmdr1 expression levels. To account for differences in the time of sampling, we developed a computational method termed Hypergeometric Analysis of Time Series, which combines Fast Fourier Transform with a modified Gene Set Enrichment Analysis. Our analysis revealed coordinated changes in genes involved in protein catabolism, translation initiation and DNA/RNA metabolism. We also observed differential expression of genes with a role in transport or coding for components of the digestive vacuole. Interestingly, a global comparison of all profiled transcriptomes uncovered a tight correlation between the transcript levels of pfcrt and pfmdr1, extending to dozens of other genes, suggesting an intricate regulatory balance in order to maintain optimal physiological processes.

CONCLUSIONS

This study provides insight into the mechanisms by which P. falciparum adjusts to the acquisition of mutations or gene amplification in key transporter loci that mediate drug resistance. Our results implicate several biological pathways that may be differentially regulated to compensate for impaired transporter function and alterations in parasite vacuole physiology.

Heterochromatin protein 1 secures survival and transmission of malaria parasites

Clonally variant expression of surface antigens allows the malaria parasite Plasmodium falciparum to evade immune recognition during blood stage infection and secure malaria transmission. We demonstrate that heterochromatin protein 1 (HP1), an evolutionary conserved regulator of heritable gene silencing, controls expression of numerous P. falciparum virulence genes as well as differentiation into the sexual forms that transmit to mosquitoes. Conditional depletion of P. falciparum HP1 (PfHP1) prevents mitotic proliferation of blood stage parasites and disrupts mutually exclusive expression and antigenic variation of the major virulence factor PfEMP1. Additionally, PfHP1-dependent regulation of PfAP2-G, a transcription factor required for gametocyte conversion, controls the switch from asexual proliferation to sexual differentiation, providing insight into the epigenetic mechanisms underlying gametocyte commitment. These findings show that PfHP1 is centrally involved in clonally variant gene expression and sexual differentiation in P. falciparum and have major implications for developing antidisease and transmission-blocking interventions against malaria.

Hybridisation of N4-methylcytosine-containing amplicons on DNA microarrays

5'-Cy5-labelled PCR amplicons containing the analogue base, N(4)-methylcytosine, instead of cytosines were compared in microarray hybridisation experiments with the corresponding amplicons containing the canonical set of bases, with respect to the intensity of the fluorescence signal obtained, and cross hybridisation to non-corresponding probes. In general, higher hybridisation temperatures resulted in reduced signal intensities, particularly in the case of the N(4)-methylcytosine containing amplicons. At the lower hybridisation temperatures tested (40 °C, 30 °C), these modified amplicons gave about equal or stronger fluorescence signal than the corresponding regular amplicons. With the two GC-richest amplicons tested, in one instance the N(4)-methylated target gave a dramatically higher signal intensity than the unmodified amplicon, interpreted as reflecting the reduced formation of hairpin structures in the target sequence, due to the lower thermodynamic stability of the G:N(4)-methylC base pair, making the target more accessible, while in the other case no hybridisation was observed with either version of the amplicon, probably due to interference from a G-tetrad structure. Both for the regular and the N(4)-methylated amplicons, no significant cross hybridisation was seen in these experiments.

Genome-wide analysis in Plasmodium falciparum reveals early and late phases of RNA polymerase II occupancy during the infectious cycle

Background

Over the course of its intraerythrocytic developmental cycle (IDC), the malaria parasite Plasmodium falciparum tightly orchestrates the rise and fall of transcript levels for hundreds of genes. Considerable debate has focused on the relative importance of transcriptional versus post-transcriptional processes in the regulation of transcript levels. Enzymatically active forms of RNAPII in other organisms have been associated with phosphorylation on the serines at positions 2 and 5 of the heptad repeats within the C-terminal domain (CTD) of RNAPII. We reasoned that insight into the contribution of transcriptional mechanisms to gene expression in P. falciparum could be obtained by comparing the presence of enzymatically active forms of RNAPII at multiple genes with the abundance of their associated transcripts.

Results

We exploited the phosphorylation state of the CTD to detect enzymatically active forms of RNAPII at most P. falciparum genes across the IDC. We raised highly specific monoclonal antibodies against three forms of the parasite CTD, namely unphosphorylated, Ser5-P and Ser2/5-P, and used these in ChIP-on-chip type experiments to map the genome-wide occupancy of RNAPII. Our data reveal that the IDC is divided into early and late phases of RNAPII occupancy evident from simple bi-phasic RNAPII binding profiles. By comparison to mRNA abundance, we identified sub-sets of genes with high occupancy by enzymatically active forms of RNAPII and relatively low transcript levels and vice versa. We further show that the presence of active and repressive histone modifications correlates with RNAPII occupancy over the IDC.

Conclusions

The simple early/late occupancy by RNAPII cannot account for the complex dynamics of mRNA accumulation over the IDC, suggesting a major role for mechanisms acting downstream of RNAPII occupancy in the control of gene expression in this parasite.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-959) contains supplementary material, which is available to authorized users.

Role of calcium signaling in the transcriptional regulation of the apicoplast genome of Plasmodium falciparum

Calcium is a universal second messenger that plays an important role in regulatory processes in eukaryotic cells. To understand calcium-dependent signaling in malaria parasites, we analyzed transcriptional responses of Plasmodium falciparum to two calcium ionophores (A23187 and ionomycin) that cause redistribution of intracellular calcium within the cytoplasm. While ionomycin induced a specific transcriptional response defined by up- or downregulation of a narrow set of genes, A23187 caused a developmental arrest in the schizont stage. In addition, we observed a dramatic decrease of mRNA levels of the transcripts encoded by the apicoplast genome during the exposure of P. falciparum to both calcium ionophores. Neither of the ionophores caused any disruptions to the DNA replication or the overall apicoplast morphology. This suggests that the mRNA downregulation reflects direct inhibition of the apicoplast gene transcription. Next, we identify a nuclear encoded protein with a calcium binding domain (EF-hand) that is localized to the apicoplast. Overexpression of this protein (termed PfACBP1) in P. falciparum cells mediates an increased resistance to the ionophores which suggests its role in calcium-dependent signaling within the apicoplast. Our data indicate that the P. falciparum apicoplast requires calcium-dependent signaling that involves a novel protein PfACBP1.

Structural polymorphism in the promoter of pfmrp2 confers Plasmodium falciparum tolerance to quinoline drugs

Drug resistance in Plasmodium falciparum remains a challenge for the malaria eradication programmes around the world. With the emergence of artemisinin resistance, the efficacy of the partner drugs in the artemisinin combination therapies (ACT) that include quinoline-based drugs is becoming critical. So far only few resistance markers have been identified from which only two transmembrane transporters namely PfMDR1 (an ATP-binding cassette transporter) and PfCRT (a drug-metabolite transporter) have been experimentally verified. Another P. falciparum transporter, the ATP-binding cassette containing multidrug resistance-associated protein (PfMRP2) represents an additional possible factor of drug resistance in P. falciparum. In this study, we identified a parasite clone that is derived from the 3D7 P. falciparum strain and shows increased resistance to chloroquine, mefloquine and quinine through the trophozoite and schizont stages. We demonstrate that the resistance phenotype is caused by a 4.1 kb deletion in the 5′ upstream region of the pfmrp2 gene that leads to an alteration in the pfmrp2 transcription and thus increased level of PfMRP2 protein. These results also suggest the importance of putative promoter elements in regulation of gene expression during the P. falciparum intra-erythrocytic developmental cycle and the potential of genetic polymorphisms within these regions to underlie drug resistance.

Asexual populations of the human malaria parasite, Plasmodium falciparum, use a two-step genomic strategy to acquire accurate, beneficial DNA amplifications

Malaria drug resistance contributes to up to a million annual deaths. Judicious deployment of new antimalarials and vaccines could benefit from an understanding of early molecular events that promote the evolution of parasites. Continuous in vitro challenge of Plasmodium falciparum parasites with a novel dihydroorotate dehydrogenase (DHODH) inhibitor reproducibly selected for resistant parasites. Genome-wide analysis of independently-derived resistant clones revealed a two-step strategy to evolutionary success. Some haploid blood-stage parasites first survive antimalarial pressure through fortuitous DNA duplications that always included the DHODH gene. Independently-selected parasites had different sized amplification units but they were always flanked by distant A/T tracks. Higher level amplification and resistance was attained using a second, more efficient and more accurate, mechanism for head-to-tail expansion of the founder unit. This second homology-based process could faithfully tune DNA copy numbers in either direction, always retaining the unique DNA amplification sequence from the original A/T-mediated duplication for that parasite line. Pseudo-polyploidy at relevant genomic loci sets the stage for gaining additional mutations at the locus of interest. Overall, we reveal a population-based genomic strategy for mutagenesis that operates in human stages of P. falciparum to efficiently yield resistance-causing genetic changes at the correct locus in a successful parasite. Importantly, these founding events arise with precision; no other new amplifications are seen in the resistant haploid blood stage parasite. This minimizes the need for meiotic genetic cleansing that can only occur in sexual stage development of the parasite in mosquitoes.

Malaria parasites kill up to a million people around the world every year. Emergence of resistance to drugs remains a key obstacle against elimination of malaria. In the laboratory, parasites can efficiently acquire resistance to experimental antimalarials by changing DNA at the target locus. This happens efficiently even for an antimalarial that the parasite has never encountered in a clinical setting. In this study, we formally demonstrate how parasites achieve this feat: first, individual parasites in a population of millions randomly amplify large regions of DNA between short sequence repeats of adenines (A) or thymines (T) that are peppered throughout the malaria parasite genome. The rare lucky parasite that amplifies DNA coding for the target of the antimalarial, along with dozens of its neighboring genes, gains an evolutionary advantage and survives. In a second step, to withstand increasing drug pressure and to achieve higher levels of resistance, each parasite line makes additional copies of this region. This second expansion does not rely on the random A/T-based DNA rearrangements but, instead, a more precise amplification mechanism that retains the unique signature of co-amplified genes created earlier in each parasite. Generation of multiple copies of the target genes in the parasite genome may be the beginning of other beneficial changes for the parasite, including the future acquisition of mutations.

Dynamic epigenetic regulation of gene expression during the life cycle of malaria parasite Plasmodium falciparum

Epigenetic mechanisms are emerging as one of the major factors of the dynamics of gene expression in the human malaria parasite, Plasmodium falciparum. To elucidate the role of chromatin remodeling in transcriptional regulation associated with the progression of the P. falciparum intraerythrocytic development cycle (IDC), we mapped the temporal pattern of chromosomal association with histone H3 and H4 modifications using ChIP-on-chip. Here, we have generated a broad integrative epigenomic map of twelve histone modifications during the P. falciparum IDC including H4K5ac, H4K8ac, H4K12ac, H4K16ac, H3K9ac, H3K14ac, H3K56ac, H4K20me1, H4K20me3, H3K4me3, H3K79me3 and H4R3me2. While some modifications were found to be associated with the vast majority of the genome and their occupancy was constant, others showed more specific and highly dynamic distribution. Importantly, eight modifications displaying tight correlations with transcript levels showed differential affinity to distinct genomic regions with H4K8ac occupying predominantly promoter regions while others occurred at the 5′ ends of coding sequences. The promoter occupancy of H4K8ac remained unchanged when ectopically inserted at a different locus, indicating the presence of specific DNA elements that recruit histone modifying enzymes regardless of their broad chromatin environment. In addition, we showed the presence of multivalent domains on the genome carrying more than one histone mark, highlighting the importance of combinatorial effects on transcription. Overall, our work portrays a substantial association between chromosomal locations of various epigenetic markers, transcriptional activity and global stage-specific transitions in the epigenome.

Malaria is a devastating parasitic disease caused by the protozoan protist Plasmodium falciparum. The complex life cycle of P. falciparum comprises various morphological and functionally distinct forms and is completed in two different hosts. Various regulatory mechanisms are employed by these parasites to complete their life cycle and survive in human hosts. Epigenetic mechanisms, though not fully explored, have been implicated as one of the key players in gene regulation, morphological differentiation and antigenic variation. Here, we present a comprehensive epigenetic map of 12 histone post-translational modifications during the intraerythrocytic life cycle of P. falciparum. We have been able to identify at least eight histone modifications whose dynamic patterns correlate with the transcriptional regulation across the life cycle. In particular, we have shown that a set of euchromatic histone marks work in synergy, creating a dynamic unique histone code that is linked with gene expression during the progression of the Plasmodium intraerythrocytic developmental cycle. These findings enhance our knowledge of complex gene regulation and will help to identify novel targets for fighting malaria.

H2A.Z and H2B.Z double-variant nucleosomes define intergenic regions and dynamically occupy var gene promoters in the malaria parasite Plasmodium falciparum

Histone variants are important components of eukaryotic chromatin and can alter chromatin structure to confer specialized functions. H2B variant histones are rare in nature but have evolved independently in the phyla Apicomplexa and Trypanasomatida. Here, we investigate the apicomplexan-specific Plasmodium falciparum histone variant Pf H2B.Z and show that within nucleosomes Pf H2B.Z dimerizes with the H2A variant Pf H2A.Z and that Pf H2B.Z and Pf H2A.Z occupancy correlates in the subset of genes examined. These double-variant nucleosomes also carry common markers of euchromatin like H3K4me3 and histone acetylation. Pf H2B.Z levels are elevated in intergenic regions across the genome, except in the var multigene family, where Pf H2A.Z/Pf H2B.Z double-variant nucleosomes are only enriched in the promoter of the single active var copy and this enrichment is developmentally regulated. Importantly, this pattern seems to be specific for var genes and does not apply to other heterochromatic gene families involved in red blood cell invasion which are also subject to clonal expression. Thus, Pf H2A.Z/Pf H2B.Z double-variant nucleosomes appear to have a highly specific function in the regulation of P. falciparum virulence.

Saccharopolyspora erythraea's genome is organised in high-order transcriptional regions mediated by targeted degradation at the metabolic switch

BACKGROUND

Actinobacteria form a major bacterial phylum that includes numerous human pathogens. Actinobacteria are primary contributors to carbon cycling and also represent a primary source of industrial high value products such as antibiotics and biopesticides. Consistent with other members of the actinobacterial phylum, Saccharopolyspora erythraea undergo a transitional switch. This switch is characterized by numerous metabolic and morphological changes.

RESULTS

We performed RNA sequencing to analyze the transcriptional changes that occur during growth of Saccharopolyspora erythraea in batch culture. By sequencing RNA across the fermentation time course, at a mean coverage of 4000X, we found the vast majority of genes to be prominently expressed, showing that we attained close to saturating sequencing coverage of the transcriptome. During the metabolic switch, global changes in gene expression influence the metabolic machinery of Saccharopolyspora erythraea, resetting an entirely novel gene expression program. After the switch, global changes include the broad repression of half the genes regulated by complex transcriptional mechanisms. Paralogous transposon clusters, delineate these transcriptional programs. The new transcriptional program is orchestrated by a bottleneck event during which mRNA levels are severely restricted by targeted mRNA degradation.

CONCLUSIONS

Our results, which attained close to saturating sequencing coverage of the transcriptome, revealed unanticipated transcriptional complexity with almost one third of transcriptional content originating from un-annotated sequences. We showed that the metabolic switch is a sophisticated mechanism of transcriptional regulation capable of resetting and re-synchronizing gene expression programs at extraordinary speed and scale.

DNA microarray-based genome-wide analyses of Plasmodium parasites

DNA microarray is presently one of the most powerful and fastest growing technologies for genomic research of infectious diseases. Accordingly, DNA microarray-based global analyses of Plasmodium parasites provided many insights into the general biology of malaria infection. From the parasite perspective, it was shown that the complex Plasmodium life cycle is characterized by a high level of coordination in gene expression but at the same time parasites have a considerable capacity to alter their transcriptional profile as a response to external stimuli and/or adaptation to varying growth conditions in their host. In addition to transcriptional profiling, DNA microarrays were shown to be useful for quantitative analyses of Plasmodium genomic DNA including characterizations of sequence polymorphisms and copy number variants (CNV) as well as genomic loci associated with different chromatin factors (e.g., immunoprecipitated material (ChIP-on-chip)). Here, we present protocols for transcriptional profiling, comparative genomic hybridization (CGH), and ChIP-on-chip analyses that have been developed for the use of low-density long oligonucleotide DNA microarrays of Plasmodium species. Many of the presented procedures including RNA purification, DNA amplification, and chromatin immunoprecipitation are likely to be transferable to other genomic platforms such as other microarray technologies and new generation sequencing.

New Agilent platform DNA microarrays for transcriptome analysis of Plasmodium falciparum and Plasmodium berghei for the malaria research community

Background

DNA microarrays have been a valuable tool in malaria research for over a decade but remain in limited use in part due their relatively high cost, poor availability, and technical difficulty. With the aim of alleviating some of these factors next-generation DNA microarrays for genome-wide transcriptome analysis for both Plasmodium falciparum and Plasmodium berghei using the Agilent 8x15K platform were designed.

Methods

Probe design was adapted from previously published methods and based on the most current transcript predictions available at the time for P. falciparum or P. berghei. Array performance and transcriptome analysis was determined using dye-coupled, aminoallyl-labelled cDNA and streamlined methods for hybridization, washing, and array analysis were developed.

Results

The new array design marks a notable improvement in the number of transcripts covered and average number of probes per transcript. Array performance was excellent across a wide range of transcript abundance, with low inter-array and inter-probe variability for relative abundance measurements and it recapitulated previously observed transcriptional patterns. Additionally, improvements in sensitivity permitted a 20-fold reduction in necessary starting RNA amounts, further reducing experimental costs and widening the range of application.

Conclusions

DNA microarrays utilizing the Agilent 8x15K platform for genome-wide transcript analysis in P. falciparum and P. berghei mark an improvement in coverage and sensitivity, increased availability to the research community, and simplification of the experimental methods.

Analysis of subtelomeric virulence gene families in Plasmodium falciparum by comparative transcriptional profiling

The Plasmodium falciparum genome is equipped with several subtelomeric gene families that are implicated in parasite virulence and immune evasion. Members of these families are uniformly positioned within heterochromatic domains and are thus subject to variegated expression. The best-studied example is that of the var family encoding the major parasite virulence factor P. falciparum erythrocyte membrane protein 1 (PfEMP1). PfEMP1 undergoes antigenic variation through switches in mutually exclusive var gene transcription. var promoters function as crucial regulatory elements in the underlying epigenetic control strategy. Here, we analysed promoters of upsA, upsB and upsC var, rifA1-type rif, stevor, phist and pfmc-2tm genes and investigated their role in endogenous gene transcription by comparative genome-wide expression profiling of transgenic parasite lines. We find that the three major var promoter types are functionally equal and play an essential role in singular gene choice. Unlike var promoters, promoters of non-var families are not silenced by default, and transcription of non-var families is not subject to the same mode of mutually exclusive transcription as has been observed for var genes. Our findings identified a differential logic in the regulation of var and other subtelomeric virulence gene families, which will have important implications for our understanding and future analyses of phenotypic variation in malaria parasites.

Transcriptional variation in the malaria parasite Plasmodium falciparum

Malaria genetic variation has been extensively characterized, but the level of epigenetic plasticity remains largely unexplored. Here we provide a comprehensive characterization of transcriptional variation in the most lethal malaria parasite, Plasmodium falciparum, based on highly accurate transcriptional analysis of isogenic parasite lines grown under homogeneous conditions. This analysis revealed extensive transcriptional heterogeneity within genetically homogeneous clonal parasite populations. We show that clonally variant expression controlled at the epigenetic level is an intrinsic property of specific genes and gene families, the majority of which participate in host–parasite interactions. Intrinsic transcriptional variability is not restricted to genes involved in immune evasion, but also affects genes linked to lipid metabolism, protein folding, erythrocyte remodeling, or transcriptional regulation, among others, indicating that epigenetic variation results in both antigenic and functional variation. We observed a general association between heterochromatin marks and clonally variant expression, extending previous observations for specific genes to essentially all variantly expressed gene families. These results suggest that phenotypic variation of functionally unrelated P. falciparum gene families is mediated by a common mechanism based on reversible formation of H3K9me3-based heterochromatin. In changing environments, diversity confers fitness to a population. Our results support the idea that P. falciparum uses a bet-hedging strategy, as an alternative to directed transcriptional responses, to adapt to common fluctuations in its environment. Consistent with this idea, we found that transcriptionally different isogenic parasite lines markedly differed in their survival to heat-shock mimicking febrile episodes and adapted to periodic heat-shock with a pattern consistent with natural selection of pre-existing parasites.

Comparative gene expression profiling of P. falciparum malaria parasites exposed to three different histone deacetylase inhibitors

Histone deacetylase (HDAC) inhibitors are being intensively pursued as potential new drugs for a range of diseases, including malaria. HDAC inhibitors are also important tools for the study of epigenetic mechanisms, transcriptional control, and other important cellular processes. In this study the effects of three structurally related antimalarial HDAC inhibitors on P. falciparum malaria parasite gene expression were compared. The three hydroxamate-based compounds, trichostatin A (TSA), suberoylanilide hydroxamic acid (SAHA; Vorinostat®) and a 2-aminosuberic acid derivative (2-ASA-9), all caused profound transcriptional effects, with ∼2–21% of genes having >2-fold altered expression following 2 h exposure to the compounds. Only two genes, alpha tubulin II and a hydrolase, were up-regulated by all three compounds after 2 h exposure in all biological replicates examined. The transcriptional changes observed after 2 h exposure to HDAC inhibitors were found to be largely transitory, with only 1–5% of genes being regulated after removing the compounds and culturing for a further 2 h. Despite some structural similarity, the three inhibitors caused quite diverse transcriptional effects, possibly reflecting subtle differences in mode of action or cellular distribution. This dataset represents an important contribution to our understanding of how HDAC inhibitors act on malaria parasites and identifies alpha tubulin II as a potential transcriptional marker of HDAC inhibition in malaria parasites that may be able to be exploited for future development of HDAC inhibitors as new antimalarial agents.

Whole-genome analysis of Plasmodium spp. Utilizing a new agilent technologies DNA microarray platform

The application of DNA microarray technologies to malaria genomics has been widely used but has been limited by sample availability and technical variability. To address these issues, we present a microarray hybridization protocol that has been optimized for use with two new Agilent Technologies DNA microarrays for Plasmodium falciparum and P. berghei. Using the most recent genome sequences available for each species, we have designed ∼14,000 oligonucleotide probes representing ∼5,600 transcripts for each species. Included in each array design are numerous probes that allow for the identification of parasite developmental stages, common Plasmodium molecular markers used in genetic manipulation, and manufacturer probes that control for array consistency and quality. Overall, the Agilent Plasmodium spp. array designs and hybridization methodology provides a sensitive, easy-to-use, high-quality, cost-effective alternative to other currently available microarray platforms.

Quantitative proteomics reveals new insights into erythrocyte invasion by Plasmodium falciparum

Differential expression of ligands in the human malaria parasite Plasmodium falciparum enables it to recognize different receptors on the erythrocyte surface, thereby providing alternative invasion pathways. Switching of invasion from using sialated to nonsialated erythrocyte receptors has been linked to the transcriptional activation of a single parasite ligand. We have used quantitative proteomics to show that in addition to this single known change, there are a significant number of changes in the expression of merozoite proteins that are regulated independent of transcription during invasion pathway switching. These results demonstrate a so far unrecognized mechanism by which the malaria parasite is able to adapt to variations in the host cell environment by post-transcriptional regulation.

Artemisinin resistance in Plasmodium falciparum is associated with an altered temporal pattern of transcription

BACKGROUND

Artemisinin resistance in Plasmodium falciparum malaria has emerged in Western Cambodia. This is a major threat to global plans to control and eliminate malaria as the artemisinins are a key component of antimalarial treatment throughout the world. To identify key features associated with the delayed parasite clearance phenotype, we employed DNA microarrays to profile the physiological gene expression pattern of the resistant isolates.

RESULTS

In the ring and trophozoite stages, we observed reduced expression of many basic metabolic and cellular pathways which suggests a slower growth and maturation of these parasites during the first half of the asexual intraerythrocytic developmental cycle (IDC). In the schizont stage, there is an increased expression of essentially all functionalities associated with protein metabolism which indicates the prolonged and thus increased capacity of protein synthesis during the second half of the resistant parasite IDC. This modulation of the P. falciparum intraerythrocytic transcriptome may result from differential expression of regulatory proteins such as transcription factors or chromatin remodeling associated proteins. In addition, there is a unique and uniform copy number variation pattern in the Cambodian parasites which may represent an underlying genetic background that contributes to the resistance phenotype.

CONCLUSIONS

The decreased metabolic activities in the ring stages are consistent with previous suggestions of higher resilience of the early developmental stages to artemisinin. Moreover, the increased capacity of protein synthesis and protein turnover in the schizont stage may contribute to artemisinin resistance by counteracting the protein damage caused by the oxidative stress and/or protein alkylation effect of this drug. This study reports the first global transcriptional survey of artemisinin resistant parasites and provides insight to the complexities of the molecular basis of pathogens with drug resistance phenotypes in vivo.

Design of a variant surface antigen-supplemented microarray chip for whole transcriptome analysis of multiple Plasmodium falciparum cytoadherent strains, and identification of strain-transcendent rif and stevor genes

Background

The cytoadherence of Plasmodium falciparum is thought to be mediated by variant surface antigens (VSA), encoded by var, rif, stevor and pfmc-2tm genes. The last three families have rarely been studied in the context of cytoadherence. As most VSA genes are unique, the variability among sequences has impeded the functional study of VSA across different P. falciparum strains. However, many P. falciparum genomes have recently been sequenced, allowing the development of specific microarray probes for each VSA gene.

Methods

All VSA sequences from the HB3, Dd2 and IT/FCR3 genomes were extracted using HMMer software. Oligonucleotide probes were designed with OligoRankPick and added to the 3D7-based microarray chip. As a proof of concept, IT/R29 parasites were selected for and against rosette formation and the transcriptomes of isogenic rosetting and non-rosetting parasites were compared by microarray.

Results

From each parasite strain 50-56 var genes, 125-132 rif genes, 26-33 stevor genes and 3-8 pfmc-2tm genes were identified. Bioinformatic analysis of the new VSA sequences showed that 13 rif genes and five stevor genes were well-conserved across at least three strains (83-100% amino acid identity). The ability of the VSA-supplemented microarray chip to detect cytoadherence-related genes was assessed using P. falciparum clone IT/R29, in which rosetting is known to be mediated by PfEMP1 encoded by ITvar9. Whole transcriptome analysis showed that the most highly up-regulated gene in rosetting parasites was ITvar9 (19 to 429-fold up-regulated over six time points). Only one rif gene (IT4rifA_042) was up-regulated by more than four fold (five fold at 12 hours post-invasion), and no stevor or pfmc-2tm genes were up-regulated by more than two fold. 377 non-VSA genes were differentially expressed by three fold or more in rosetting parasites, although none was as markedly or consistently up-regulated as ITvar9.

Conclusions

Probes for the VSA of newly sequenced P. falciparum strains can be added to the 3D7-based microarray chip, allowing the analysis of the entire transcriptome of multiple strains. For the rosetting clone IT/R29, the striking transcriptional upregulation of ITvar9 was confirmed, and the data did not support the involvement of other VSA families in rosette formation.

Quantitative time-course profiling of parasite and host cell proteins in the human malaria parasite Plasmodium falciparum

Studies of the Plasmodium falciparum transcriptome have shown that the tightly controlled progression of the parasite through the intra-erythrocytic developmental cycle (IDC) is accompanied by a continuous gene expression cascade in which most expressed genes exhibit a single transcriptional peak. Because the biochemical and cellular functions of most genes are mediated by the encoded proteins, understanding the relationship between mRNA and protein levels is crucial for inferring biological activity from transcriptional gene expression data. Although studies on other organisms show that <50% of protein abundance variation may be attributable to corresponding mRNA levels, the situation in Plasmodium is further complicated by the dynamic nature of the cyclic gene expression cascade. In this study, we simultaneously determined mRNA and protein abundance profiles for P. falciparum parasites during the IDC at 2-hour resolution based on oligonucleotide microarrays and two-dimensional differential gel electrophoresis protein gels. We find that most proteins are represented by more than one isoform, presumably because of post-translational modifications. Like transcripts, most proteins exhibit cyclic abundance profiles with one peak during the IDC, whereas the presence of functionally related proteins is highly correlated. In contrast, the abundance of most parasite proteins peaks significantly later (median 11 h) than the corresponding transcripts and often decreases slowly in the second half of the IDC. Computational modeling indicates that the considerable and varied incongruence between transcript and protein abundance may largely be caused by the dynamics of translation and protein degradation. Furthermore, we present cyclic abundance profiles also for parasite-associated human proteins and confirm the presence of five human proteins with a potential role in antioxidant defense within the parasites. Together, our data provide fundamental insights into transcript-protein relationships in P. falciparum that are important for the correct interpretation of transcriptional data and that may facilitate the improvement and development of malaria diagnostics and drug therapy.

A major role for the Plasmodium falciparum ApiAP2 protein PfSIP2 in chromosome end biology

The heterochromatic environment and physical clustering of chromosome ends at the nuclear periphery provide a functional and structural framework for antigenic variation and evolution of subtelomeric virulence gene families in the malaria parasite Plasmodium falciparum. While recent studies assigned important roles for reversible histone modifications, silent information regulator 2 and heterochromatin protein 1 (PfHP1) in epigenetic control of variegated expression, factors involved in the recruitment and organization of subtelomeric heterochromatin remain unknown. Here, we describe the purification and characterization of PfSIP2, a member of the ApiAP2 family of putative transcription factors, as the unknown nuclear factor interacting specifically with cis-acting SPE2 motif arrays in subtelomeric domains. Interestingly, SPE2 is not bound by the full-length protein but rather by a 60kDa N-terminal domain, PfSIP2-N, which is released during schizogony. Our experimental re-definition of the SPE2/PfSIP2-N interaction highlights the strict requirement of both adjacent AP2 domains and a conserved bipartite SPE2 consensus motif for high-affinity binding. Genome-wide in silico mapping identified 777 putative binding sites, 94% of which cluster in heterochromatic domains upstream of subtelomeric var genes and in telomere-associated repeat elements. Immunofluorescence and chromatin immunoprecipitation (ChIP) assays revealed co-localization of PfSIP2-N with PfHP1 at chromosome ends. Genome-wide ChIP demonstrated the exclusive binding of PfSIP2-N to subtelomeric SPE2 landmarks in vivo but not to single chromosome-internal sites. Consistent with this specialized distribution pattern, PfSIP2-N over-expression has no effect on global gene transcription. Hence, contrary to the previously proposed role for this factor in gene activation, our results provide strong evidence for the first time for the involvement of an ApiAP2 factor in heterochromatin formation and genome integrity. These findings are highly relevant for our understanding of chromosome end biology and variegated expression in P. falciparum and other eukaryotes, and for the future analysis of the role of ApiAP2-DNA interactions in parasite biology.

Defining species specific genome differences in malaria parasites

Background

In recent years a number of genome sequences for different plasmodium species have become available. This has allowed the identification of numerous conserved genes across the different species and has significantly enhanced our understanding of parasite biology. In contrast little is known about species specific differences between the different genomes partly due to the lower sequence coverage and therefore relatively poor annotation of some of the draft genomes particularly the rodent malarias parasite species.

Results

To improve the current annotation and gene identification status of the draft genomes of P. berghei, P. chabaudi and P. yoelii, we performed genome-wide comparisons between these three species. Through analyses via comparative genome hybridizations using a newly designed pan-rodent array as well as in depth bioinformatics analysis, we were able to improve on the coverage of the draft rodent parasite genomes by detecting orthologous genes between these related rodent parasite species. More than 1,000 orthologs for P. yoelii were now newly associated with a P. falciparum gene. In addition to extending the current core gene set for all plasmodium species this analysis also for the first time identifies a relatively small number of genes that are unique to the primate malaria parasites while a larger gene set is uniquely conserved amongst the rodent malaria parasites.

Conclusions

These findings allow a more thorough investigation of the genes that are important for host specificity in malaria.

New insights into the blood-stage transcriptome of Plasmodium falciparum using RNA-Seq

Recent advances in high-throughput sequencing present a new opportunity to deeply probe an organism's transcriptome. In this study, we used Illumina-based massively parallel sequencing to gain new insight into the transcriptome (RNA-Seq) of the human malaria parasite, Plasmodium falciparum. Using data collected at seven time points during the intraerythrocytic developmental cycle, we (i) detect novel gene transcripts; (ii) correct hundreds of gene models; (iii) propose alternative splicing events; and (iv) predict 5' and 3' untranslated regions. Approximately 70% of the unique sequencing reads map to previously annotated protein-coding genes. The RNA-Seq results greatly improve existing annotation of the P. falciparum genome with over 10% of gene models modified. Our data confirm 75% of predicted splice sites and identify 202 new splice sites, including 84 previously uncharacterized alternative splicing events. We also discovered 107 novel transcripts and expression of 38 pseudogenes, with many demonstrating differential expression across the developmental time series. Our RNA-Seq results correlate well with DNA microarray analysis performed in parallel on the same samples, and provide improved resolution over the microarray-based method. These data reveal new features of the P. falciparum transcriptional landscape and significantly advance our understanding of the parasite's red blood cell-stage transcriptome.

Histone deacetylases play a major role in the transcriptional regulation of the Plasmodium falciparum life cycle

The apparent paucity of molecular factors of transcriptional control in the genomes of Plasmodium parasites raises many questions about the mechanisms of life cycle regulation in these malaria parasites. Epigenetic regulation has been suggested to play a major role in the stage specific gene expression during the Plasmodium life cycle. To address some of these questions, we analyzed global transcriptional responses of Plasmodium falciparum to a potent inhibitor of histone deacetylase activities (HDAC). The inhibitor apicidin induced profound transcriptional changes in multiple stages of the P. falciparum intraerythrocytic developmental cycle (IDC) that were characterized by rapid activation and repression of a large percentage of the genome. A major component of this response was induction of genes that are otherwise suppressed during that particular stage of the IDC or specific for the exo-erythrocytic stages. In the schizont stage, apicidin induced hyperacetylation of histone lysine residues H3K9, H4K8 and the tetra-acetyl H4 (H4Ac4) and demethylation of H3K4me3. Interestingly, we observed overlapping patterns of chromosomal distributions between H4K8Ac and H3K4me3 and between H3K9Ac and H4Ac4. There was a significant but partial association between the apicidin-induced gene expression and histone modifications, which included a number of stage specific transcription factors. Taken together, inhibition of HDAC activities leads to dramatic de-regulation of the IDC transcriptional cascade, which is a result of both disruption of histone modifications and up-regulation of stage specific transcription factors. These findings suggest an important role of histone modification and chromatin remodeling in transcriptional regulation of the Plasmodium life cycle. This also emphasizes the potential of P. falciparum HDACs as drug targets for malaria chemotherapy.

Plasmodium falciparum, a parasitic protozoan, causes the most lethal form of human malaria, killing more than 2 million people per year. It has a complex life cycle that involves distinct morphological stages accompanied by stage specific gene expression in both the mosquito and human hosts. The lack of a vaccine for malaria and widespread resistance highlights the urgency for new anti-malarial drugs that act on different parasite targets. We show that inhibition of histone deacetylase activities results in activation and repression of transcriptionally regulated genes in multiple stages of the P. falciparum asexual life cycle. We also show that inhibition disrupts the steady-state level of histone acetylation and methylation across the P. falciparum genome. Our data strongly implies that in P. falciparum, inhibition of histone deacetylase activity leads to a dramatic increase in global acetylation of histones and subsequently disruption of stage specific gene expression. This process then leads to a collapse of the transcriptional cascade of P. falciparum. Therefore, the essential role of histone deacetylases in Plasmodium parasites suggests their high potential as molecular targets for malaria intervention strategies.

Optimized detection of sequence variation in heterozygous genomes using DNA microarrays with isothermal-melting probes

The use of DNA microarrays to identify nucleotide variation is almost 20 years old. A variety of improvements in probe design and experimental conditions have brought this technology to the point that single-nucleotide differences can be efficiently detected in unmixed samples, although developing reliable methods for detection of mixed sequences (e.g., heterozygotes) remains challenging. Surprisingly, a comprehensive study of the probe design parameters and experimental conditions that optimize discrimination of single-nucleotide polymorphisms (SNPs) has yet to be reported, so the limits of this technology remain uncertain. By targeting 24,549 SNPs that differ between two Saccharomyces cerevisiae strains, we studied the effect of SNPs on hybridization efficiency to DNA microarray probes of different lengths under different hybridization conditions. We found that the critical parameter for optimization of sequence discrimination is the relationship between probe melting temperature (T(m)) and the temperature at which the hybridization reaction is performed. This relationship can be exploited through the design of microarrays containing probes of equal T(m) by varying the length of probes. We demonstrate using such a microarray that we detect >90% homozygous SNPs and >80% heterozygous SNPs using the SNPScanner algorithm. The optimized design and experimental parameters determined in this study should guide DNA microarray designs for applications that require sequence discrimination such as mutation detection, genotyping of unmixed and mixed samples, and allele-specific gene expression. Moreover, designing microarray probes with optimized sensitivity to mismatches should increase the accuracy of standard microarray applications such as copy-number variation detection and gene expression analysis.

Transcriptional profiling of growth perturbations of the human malaria parasite Plasmodium falciparum

Functions have yet to be defined for the majority of genes of Plasmodium falciparum, the agent responsible for the most serious form of human malaria. Here we report changes in P. falciparum gene expression induced by 20 compounds that inhibit growth of the schizont stage of the intraerythrocytic development cycle. In contrast with previous studies, which reported only minimal changes in response to chemically induced perturbations of P. falciparum growth, we find that approximately 59% of its coding genes display over three-fold changes in expression in response to at least one of the chemicals we tested. We use this compendium for guilt-by-association prediction of protein function using an interaction network constructed from gene co-expression, sequence homology, domain-domain and yeast two-hybrid data. The subcellular localizations of 31 of 42 proteins linked with merozoite invasion is consistent with their role in this process, a key target for malaria control. Our network may facilitate identification of novel antimalarial drugs and vaccines.

PrimerBank: a resource of human and mouse PCR primer pairs for gene expression detection and quantification

PrimerBank (http://pga.mgh.harvard.edu/primerbank/) is a public resource for the retrieval of human and mouse primer pairs for gene expression analysis by PCR and Quantitative PCR (QPCR). A total of 306,800 primers covering most known human and mouse genes can be accessed from the PrimerBank database, together with information on these primers such as T(m), location on the transcript and amplicon size. For each gene, at least one primer pair has been designed and in many cases alternative primer pairs exist. Primers have been designed to work under the same PCR conditions, thus facilitating high-throughput QPCR. There are several ways to search for primers for the gene(s) of interest, such as by: GenBank accession number, NCBI protein accession number, NCBI gene ID, PrimerBank ID, NCBI gene symbol or gene description (keyword). In all, 26,855 primer pairs covering most known mouse genes have been experimentally validated by QPCR, agarose gel analysis, sequencing and BLAST, and all validation data can be freely accessed from the PrimerBank web site.

Comparative transcriptional and genomic analysis of Plasmodium falciparum field isolates

Mechanisms for differential regulation of gene expression may underlie much of the phenotypic variation and adaptability of malaria parasites. Here we describe transcriptional variation among culture-adapted field isolates of Plasmodium falciparum, the species responsible for most malarial disease. It was found that genes coding for parasite protein export into the red cell cytosol and onto its surface, and genes coding for sexual stage proteins involved in parasite transmission are up-regulated in field isolates compared with long-term laboratory isolates. Much of this variability was associated with the loss of small or large chromosomal segments, or other forms of gene copy number variation that are prevalent in the P. falciparum genome (copy number variants, CNVs). Expression levels of genes inside these segments were correlated to that of genes outside and adjacent to the segment boundaries, and this association declined with distance from the CNV boundary. This observation could not be explained by copy number variation in these adjacent genes. This suggests a local-acting regulatory role for CNVs in transcription of neighboring genes and helps explain the chromosomal clustering that we observed here. Transcriptional co-regulation of physical clusters of adaptive genes may provide a way for the parasite to readily adapt to its highly heterogeneous and strongly selective environment.

Plasmodium falciparum heterochromatin protein 1 marks genomic loci linked to phenotypic variation of exported virulence factors

Epigenetic processes are the main conductors of phenotypic variation in eukaryotes. The malaria parasite Plasmodium falciparum employs antigenic variation of the major surface antigen PfEMP1, encoded by 60 var genes, to evade acquired immune responses. Antigenic variation of PfEMP1 occurs through in situ switches in mono-allelic var gene transcription, which is PfSIR2-dependent and associated with the presence of repressive H3K9me3 marks at silenced loci. Here, we show that P. falciparum heterochromatin protein 1 (PfHP1) binds specifically to H3K9me3 but not to other repressive histone methyl marks. Based on nuclear fractionation and detailed immuno-localization assays, PfHP1 constitutes a major component of heterochromatin in perinuclear chromosome end clusters. High-resolution genome-wide chromatin immuno-precipitation demonstrates the striking association of PfHP1 with virulence gene arrays in subtelomeric and chromosome-internal islands and a high correlation with previously mapped H3K9me3 marks. These include not only var genes, but also the majority of P. falciparum lineage-specific gene families coding for exported proteins involved in host-parasite interactions. In addition, we identified a number of PfHP1-bound genes that were not enriched in H3K9me3, many of which code for proteins expressed during invasion or at different life cycle stages. Interestingly, PfHP1 is absent from centromeric regions, implying important differences in centromere biology between P. falciparum and its human host. Over-expression of PfHP1 results in an enhancement of variegated expression and highlights the presence of well-defined heterochromatic boundaries. In summary, we identify PfHP1 as a major effector of virulence gene silencing and phenotypic variation. Our results are instrumental for our understanding of this widely used survival strategy in unicellular pathogens.

Clonally variant gene families in Plasmodium falciparum share a common activation factor

The genome of the malaria parasite Plasmodium falciparum contains several multicopy gene families, including var, rifin, stevor and Pfmc-2TM. These gene families undergo expression switching and appear to play a role in antigenic variation. It has recently been shown that forcing parasites to express high copy numbers of transcriptionally active, episomal var promoters led to gradual downregulation and eventual silencing of the entire var gene family, suggesting that a limiting titratable factor plays a role in var gene activation. Through similar experiments using rifin, stevor or Pfmc-2TM episomal promoters we show that promoter titration can be used as a general method to downregulate multicopy gene families in P. falciparum. Additionally, we show that promoter titration with var, rifin, stevor or Pfmc-2TM episomal promoters results in downregulation of expression not only of the family to which the episomal promoter belongs, but also members of the other gene families, suggesting that the var-specific titratable factor previously described is shared by all four families. Further, transcriptionally active promoters from different families colocalize within the same subnuclear expression site, indicating that the role that nuclear architecture plays in var gene regulation also likely applies to the other multicopy gene families of P. falciparum.

Microarrays and high-throughput transcriptomic analysis in species with incomplete availability of genomic sequences

Microarrays produce a measurement of gene expression based on the relative measures of dye intensities that correspond to the amount of target RNA. This technology is fast developing and its application is expanding from Homo sapiens to a wide number of species, where enough information on sequences and annotations exist. Anyway, the number of species for which a dedicated platform exists is not high. The use of heterologous array hybridization, screening for gene expression in one species using an array developed for another one, is still quite frequent, even though cross-species microarray hybridization has raised many arguments. Some methods which are high throughput and do not rely on knowledge of the DNA/RNA sequence exist, namely serial analysis of gene expression (SAGE), Massively Parallel Signature Sequencing (MPSS) and deep sequencing of full transcriptome. Although very powerful, particularly the latter, they are still quite costly and cumbersome methods. In some species where genome sequences are largely unknown, several anonymous sequences are deposited in gene banks as a result of Expressed Sequence Tags (ESTs) sequencing projects. The ESTs databases represent a valuable knowledge that can be exploited with some bioinformatic effort to build species-specific microarrays. We present here a method of high-density in situ synthesized microarrays starting from available EST sequences in, Ovis aries. Our data indicate that the method is very efficient and can be easily extended to other species of which genetic sequences are present in public databases, but neglected so far with advanced devices like microarrays. As a perspective, the approach can be applied also to species of which no sequences are available to date, thanks to high-throughput deep sequencing methods.

The Plasmodium falciparum STEVOR multigene family mediates antigenic variation of the infected erythrocyte

Modifications of the Plasmodium falciparum-infected red blood cell (iRBC) surface have been linked to parasite-associated pathology. Such modifications enable the parasite to establish long-lasting chronic infection by evading antibody mediate immune recognition and splenic clearance. With the exception of the well-demonstrated roles of var-encoded PfEMP1 in virulence and immune evasion, the biological significance of other variant surface antigens (rif and stevor) is largely unknown. While PfEMP1 and RIFIN have been located on the iRBC surface, recent studies have located STEVOR at the iRBC membrane where it may be exposed on the erythrocyte surface. To investigate the role of STEVOR in more detail, we have developed antibodies against two putative STEVOR proteins and used a combination of indirect immunofluorescence assays (IFA), live IFA, flow cytometry, as well as agglutination assays, which enable us to demonstrate that STEVOR is clonally variant at the surface of schizont stage parasites. Crucially, expression of different STEVOR on the surface of the iRBC changes the antigenic property of the parasite. Taken together, our data for the first time demonstrate that STEVOR plays a role in creating antigenic diversity of schizont stage parasites, thereby adding additional complexity to the immunogenic properties of the iRBC. Furthermore, it clearly demonstrates that to obtain a complete understanding of how parasite-induced pathology is linked to variation on the surface of the iRBC, focusing the interactions of multiple multigene families needs to be considered.

Design of 240,000 orthogonal 25mer DNA barcode probes

DNA barcodes linked to genetic features greatly facilitate screening these features in pooled formats using microarray hybridization, and new tools are needed to design large sets of barcodes to allow construction of large barcoded mammalian libraries such as shRNA libraries. Here we report a framework for designing large sets of orthogonal barcode probes. We demonstrate the utility of this framework by designing 240,000 barcode probes and testing their performance by hybridization. From the test hybridizations, we also discovered new probe design rules that significantly reduce cross-hybridization after their introduction into the framework of the algorithm. These rules should improve the performance of DNA microarray probe designs for many applications.

Model-based probe set optimization for high-performance microarrays

A major challenge in microarray design is the selection of highly specific oligonucleotide probes for all targeted genes of interest, while maintaining thermodynamic uniformity at the hybridization temperature. We introduce a novel microarray design framework (Thermodynamic Model-based Oligo Design Optimizer, TherMODO) that for the first time incorporates a number of advanced modelling features: (i) A model of position-dependent labelling effects that is quantitatively derived from experiment. (ii) Multi-state thermodynamic hybridization models of probe binding behaviour, including potential cross-hybridization reactions. (iii) A fast calibrated sequence-similarity-based heuristic for cross-hybridization prediction supporting large-scale designs. (iv) A novel compound score formulation for the integrated assessment of multiple probe design objectives. In contrast to a greedy search for probes meeting parameter thresholds, this approach permits an optimization at the probe set level and facilitates the selection of highly specific probe candidates while maintaining probe set uniformity. (v) Lastly, a flexible target grouping structure allows easy adaptation of the pipeline to a variety of microarray application scenarios. The algorithm and features are discussed and demonstrated on actual design runs. Source code is available on request.

Plasmodium vivax and the importance of the subtelomeric multigene vir superfamily

Plasmodium vivax is responsible for more than 100 million clinical cases yearly. Unlike P. falciparum, in which infected red blood cells cytoadhere via variant proteins, avoiding passage through the spleen, P.-vivax-infected reticulocytes seem not to cytoadhere. However, a variant subtelomeric multigene vir family has been identified in P. vivax. Thus, questions remain about how P. vivax circulates through the spleen and the role of Vir proteins. In this review, the importance of the vir multigene superfamily is reviewed in the light of the completion of the entire genome sequence of P. vivax and from data gathered from experimental infections in reticulocyte-prone non-lethal malaria parasites and natural P. vivax infections.

The transcriptome of Plasmodium vivax reveals divergence and diversity of transcriptional regulation in malaria parasites

Plasmodium vivax causes over 100 million clinical infections each year. Primarily because of the lack of a suitable culture system, our understanding of the biology of this parasite lags significantly behind that of the more deadly species P. falciparum. Here, we present the complete transcriptional profile throughout the 48-h intraerythrocytic cycle of three distinct P. vivax isolates. This approach identifies strain specific patterns of expression for subsets of genes predicted to encode proteins associated with virulence and host pathogen interactions. Comparison to P. falciparum revealed significant differences in the expression of genes involved in crucial cellular functions that underpin the biological differences between the two parasite species. These data provide insights into the biology of P. vivax and constitute an important resource for the development of therapeutic approaches.

Antibodies targeting the PfRH1 binding domain inhibit invasion of Plasmodium falciparum merozoites

Invasion by the malaria merozoite depends on recognition of specific erythrocyte surface receptors by parasite ligands. Plasmodium falciparum uses multiple ligands, including at least two gene families, reticulocyte binding protein homologues (RBLs) and erythrocyte binding proteins/ligands (EBLs). The combination of different RBLs and EBLs expressed in a merozoite defines the invasion pathway utilized and could also play a role in parasite virulence. The binding regions of EBLs lie in a conserved cysteine-rich domain while the binding domain of RBL is still not well characterized. Here, we identify the erythrocyte binding region of the P. falciparum reticulocyte binding protein homologue 1 (PfRH1) and show that antibodies raised against the functional binding region efficiently inhibit invasion. In addition, we directly demonstrate that changes in the expression of RBLs can constitute an immune evasion mechanism of the malaria merozoite.

Plasmodium falciparum STEVOR proteins are highly expressed in patient isolates and located in the surface membranes of infected red blood cells and the apical tips of merozoites

The human parasite Plasmodium falciparum has the potential to express a vast repertoire of variant proteins on the surface of the infected red blood cell (iRBC). Variation in the expression pattern of these proteins is linked to antigenic variation and thereby evasion of host antibody-mediated immunity. The genes in the stevor multigene family code for small variant antigens that are expressed in blood-stage parasites where they can be detected in membranous structures called Maurer's clefts (MC). Some studies have indicated that STEVOR protein may also be trafficked to the iRBC membrane. To address the location of STEVOR protein in more detail, we have analyzed expression in several cultured parasite lines and in parasites obtained directly from patients. We detected STEVOR expression in a higher proportion of parasites recently isolated from patients than in cultured parasite lines and show that STEVOR is trafficked in schizont-stage parasites from the MC to the RBC cytosol and the iRBC membrane. Furthermore, STEVOR protein is also detected at the apical end of merozoites. Importantly, we show that culture-adapted parasites do not require STEVOR for survival. These findings provide new insights into the role of the stevor multigene family during both the schizont and merozoite stages of the parasite and highlight the importance of studying freshly isolated parasites, rather than parasite lines maintained in culture, when investigating potential mediators of host-parasite interactions.