RSpred, a set of Hidden Markov Models to detect and classify the RIFIN and STEVOR proteins of Plasmodium falciparum

<p>Abstract</p> <p>Background</p> <p>Many parasites use multicopy protein families to avoid their host's immune system through a strategy called antigenic variation. RIFIN and STEVOR proteins are variable surface antigens uniquely found in the malaria parasites <it>Plasmodium falciparum </it>and <it>P. reichenowi</it>. Although these two protein families are different, they have more similarity to each other than to any other proteins described to date. As a result, they have been grouped together in one Pfam domain. However, a recent study has described the sub-division of the RIFIN protein family into several functionally distinct groups. These sub-groups require phylogenetic analysis to sort out, which is not practical for large-scale projects, such as the sequencing of patient isolates and meta-genomic analysis.</p> <p>Results</p> <p>We have manually curated the <it>rif </it>and <it>stevor </it>gene repertoires of two <it>Plasmodium falciparum </it>genomes, isolates DD2 and HB3. We have identified 25% of mis-annotated and ~30 missing <it>rif </it>and <it>stevor </it>genes. Using these data sets, as well as sequences from the well curated reference genome (isolate 3D7) and field isolate data from Uniprot, we have developed a tool named RSpred. The tool, based on a set of hidden Markov models and an evaluation program, automatically identifies STEVOR and RIFIN sequences as well as the sub-groups: A-RIFIN, B-RIFIN, B1-RIFIN and B2-RIFIN. In addition to these groups, we distinguish a small subset of STEVOR proteins that we named STEVOR-like, as they either differ remarkably from typical STEVOR proteins or are too fragmented to reach a high enough score. When compared to Pfam and TIGRFAMs, RSpred proves to be a more robust and more sensitive method. We have applied RSpred to the proteomes of several <it>P. falciparum </it>strains, <it>P. reichenowi, P. vivax</it>, <it>P. knowlesi </it>and the rodent malaria species. All groups were found in the <it>P. falciparum </it>strains, and also in the <it>P. reichenowi </it>parasite, whereas none were predicted in the other species.</p> <p>Conclusions</p> <p>We have generated a tool for the sorting of RIFIN and STEVOR proteins, large antigenic variant protein groups, into homogeneous sub-families. Assigning functions to such protein families requires their subdivision into meaningful groups such as we have shown for the RIFIN protein family. RSpred removes the need for complicated and time consuming phylogenetic analysis methods. It will benefit both research groups sequencing whole genomes as well as others working with field isolates. RSpred is freely accessible via <url>http://www.ifm.liu.se/bioinfo/</url>.</p

Antigenic Variation in Plasmodium falciparum Malaria Involves a Highly Structured Switching Pattern

Many pathogenic bacteria, fungi, and protozoa achieve chronic infection through an immune evasion strategy known as antigenic variation. In the human malaria parasite Plasmodium falciparum, this involves transcriptional switching among members of the var gene family, causing parasites with different antigenic and phenotypic characteristics to appear at different times within a population. Here we use a genome-wide approach to explore this process in vitro within a set of cloned parasite populations. Our analyses reveal a non-random, highly structured switch pathway where an initially dominant transcript switches via a set of switch-intermediates either to a new dominant transcript, or back to the original. We show that this specific pathway can arise through an evolutionary conflict in which the pathogen has to optimise between safeguarding its limited antigenic repertoire and remaining capable of establishing infections in non-naïve individuals. Our results thus demonstrate a crucial role for structured switching during the early phases of infections and provide a unifying theory of antigenic variation in P. falciparum malaria as a balanced process of parasite-intrinsic switching and immune-mediated selection

Promoter regions of Plasmodium vivax are poorly or not recognized by Plasmodium falciparum

BACKGROUND: Heterologous promoter analysis in Plasmodium has revealed the existence of conserved cis regulatory elements as promoters from different species can drive expression of reporter genes in heterologous transfection assays. Here, the functional characterization of different Plasmodium vivax promoters in Plasmodium falciparum using luciferase as the reporter gene is presented. METHODS: Luciferase reporter plasmids harboring the upstream regions of the msp1, dhfr, and vir3 genes as well as the full-length intergenic regions of the vir23/24 and ef-1α genes of P. vivax were constructed and transiently transfected in P. falciparum. RESULTS: Only the constructs with the full-length intergenic regions of the vir23/24 and ef-1α genes were recognized by the P. falciparum transcription machinery albeit to values approximately two orders of magnitude lower than those reported by luc plasmids harbouring promoter regions from P. falciparum and Plasmodium berghei. A bioinformatics approach allowed the identification of a motif (GCATAT) in the ef-1α intergenic region that is conserved in five Plasmodium species but is degenerate (GCANAN) in P. vivax. Mutations of this motif in the P. berghei ef-1α promoter region decreased reporter expression indicating it is active in gene expression in Plasmodium. CONCLUSION: Together, this data indicates that promoter regions of P. vivax are poorly or not recognized by the P. falciparum transcription machinery suggesting the existence of P. vivax-specific transcription regulatory elements

The effects of symmetry on the dynamics of antigenic variation

In the studies of dynamics of pathogens and their interactions with a host immune system, an important role is played by the structure of antigenic variants associated with a pathogen. Using the example of a model of antigenic variation in malaria, we show how many of the observed dynamical regimes can be explained in terms of the symmetry of interactions between different antigenic variants. The results of this analysis are quite generic, and have wider implications for understanding the dynamics of immune escape of other parasites, as well as for the dynamics of multi-strain diseases.Comment: 21 pages, 4 figures; J. Math. Biol. (2012), Online Firs

A Statistically Rigorous Method for Determining Antigenic Switching Networks

Many vector-borne pathogens rely on antigenic variation to prolong infections and increase their likelihood of onward transmission. This immune evasion strategy often involves mutually exclusive switching between members of gene families that encode functionally similar but antigenically different variants during the course of a single infection. Studies of different pathogens have suggested that switching between variant genes is non-random and that genes have intrinsic probabilities of being activated or silenced. These factors could create a hierarchy of gene expression with important implications for both infection dynamics and the acquisition of protective immunity. Inferring complete switching networks from gene transcription data is problematic, however, because of the high dimensionality of the system and uncertainty in the data. Here we present a statistically rigorous method for analysing temporal gene transcription data to reconstruct an underlying switching network. Using artificially generated transcription profiles together with in vitro var gene transcript data from two Plasmodium falciparum laboratory strains, we show that instead of relying on data from long-term parasite cultures, accuracy can be greatly improved by using transcription time courses of several parasite populations from the same isolate, each starting with different variant distributions. The method further provides explicit indications about the reliability of the resulting networks and can thus be used to test competing hypotheses with regards to the underlying switching pathways. Our results demonstrate that antigenic switch pathways can be determined reliably from short gene transcription profiles assessing multiple time points, even when subject to moderate levels of experimental error. This should yield important new information about switching patterns in antigenically variable organisms and might help to shed light on the molecular basis of antigenic variation

An Upstream Open Reading Frame Controls Translation of var2csa, a Gene Implicated in Placental Malaria

Malaria, caused by the parasite Plasmodium falciparum, is responsible for substantial morbidity, mortality and economic losses in tropical regions of the world. Pregnant women are exceptionally vulnerable to severe consequences of the infection, due to the specific adhesion of parasite-infected erythrocytes in the placenta. This adhesion is mediated by a unique variant of PfEMP1, a parasite encoded, hyper-variable antigen placed on the surface of infected cells. This variant, called VAR2CSA, binds to chondroitin sulfate A on syncytiotrophoblasts in the intervillous space of placentas. VAR2CSA appears to only be expressed in the presence of a placenta, suggesting that its expression is actively repressed in men, children or non-pregnant women; however, the mechanism of repression is not understood. Using cultured parasite lines and reporter gene constructs, we show that the gene encoding VAR2CSA contains a small upstream open reading frame that acts to repress translation of the resulting mRNA, revealing a novel form of gene regulation in malaria parasites. The mechanism underlying this translational repression is reversible, allowing high levels of protein translation upon selection, thus potentially enabling parasites to upregulate expression of this variant antigen in the presence of the appropriate host tissue

Invasion of Wolbachia into Anopheles and Other Insect Germlines in an Ex vivo Organ Culture System

The common bacterial endosymbiont Wolbachia manipulates its host's reproduction to promote its own maternal transmission, and can interfere with pathogen development in many insects making it an attractive agent for the control of arthropod-borne disease. However, many important species, including Anopheles mosquitoes, are uninfected. Wolbachia can be artificially transferred between insects in the laboratory but this can be a laborious and sometimes fruitless process. We used a simple ex vivo culturing technique to assess the suitability of Wolbachia-host germline associations. Wolbachia infects the dissected germline tissue of multiple insect species when the host tissue and bacteria are cultured together. Ovary and testis infection occurs in a density-dependent manner. Wolbachia strains are more capable of invading the germline of their native or closely related rather than divergent hosts. The ability of Wolbachia to associate with the germline of novel hosts is crucial for the development of stably-transinfected insect lines. Rapid assessment of the suitability of a strain-host combination prior to transinfection may dictate use of a particular Wolbachia strain. Furthermore, the cultured germline tissues of two major Anopheline vectors of Plasmodium parasites are susceptible to Wolbachia infection. This finding further enhances the prospect of using Wolbachia for the biological control of malaria

Differential, Positional-Dependent Transcriptional Response of Antigenic Variation (var) Genes to Biological Stress in Plasmodium falciparum

1% of the genes of the human malaria causing agent Plasmodium falciparum belong to the heterogeneous var gene family which encodes P. falciparum erythrocyte membrane protein 1 (PFEMP1). This protein mediates part of the pathogenesis of the disease by causing adherence of infected erythrocytes (IE) to the host endothelium. At any given time, only one copy of the family is expressed on the IE surface. The cues which regulate the allelic exclusion of these genes are not known. We show the existence of a differential expression pattern of these genes upon exposure to biological stress in relation to their positional placement on the chromosome – expression of centrally located var genes is induced while sub-telomeric copies of the family are repressed - this phenomenon orchestrated by the histone deacetylase pfsir2. Moreover, stress was found to cause a switch in the pattern of the expressed var genes thus acting as a regulatory cue. By using pharmacological compounds which putatively affect pfsir2 activity, distinct changes of var gene expression patterns were achieved which may have therapeutic ramifications. As disease severity is partly associated with expression of particular var gene subtypes, manipulation of the IE environment may serve as a mechanism to direct transcription towards less virulent genes