Var Genes Research Articles

Comparison of molecular surveillance methods to assess changes in the population genetics of Plasmodium falciparum in high transmission.

A major motivation for developing molecular methods for malaria surveillance is to measure the impact of control interventions on the population genetics of Plasmodium falciparum as a potential marker of progress towards elimination. Here we assess three established methods (i) single nucleotide polymorphism (SNP) barcoding (panel of 24-biallelic loci), (ii) microsatellite genotyping (panel of 12-multiallelic loci), and (iii) varcoding (fingerprinting var gene diversity, akin to microhaplotyping) to identify changes in parasite population genetics in response to a short-term indoor residual spraying (IRS) intervention. Typical of high seasonal transmission in Africa, multiclonal infections were found in 82.3% (median 3; range 1-18) and 57.8% (median 2; range 1-12) of asymptomatic individuals pre- and post-IRS, respectively, in Bongo District, Ghana. Since directly phasing multilocus haplotypes for population genetic analysis is not possible for biallelic SNPs and microsatellites, we chose ~200 low-complexity infections biased to single and double clone infections for analysis. Each genotyping method presented a different pattern of change in diversity and population structure as a consequence of variability in usable data and the relative polymorphism of the molecular markers (i.e., SNPs < microsatellites < var). Varcoding and microsatellite genotyping showed the overall failure of the IRS intervention to significantly change the population structure from pre-IRS characteristics (i.e., many diverse genomes of low genetic similarity). The 24-SNP barcode provided limited information for analysis, largely due to the biallelic nature of SNPs leading to a high proportion of double-allele calls and a view of more isolate relatedness compared to microsatellites and varcoding. Relative performance, suitability, and cost-effectiveness of the methods relevant to sample size and local malaria elimination in high-transmission endemic areas are discussed.

Molecular epidemiology of continued Plasmodium falciparum disease transmission after an outbreak in Ecuador.

To better understand the factors underlying the continued incidence of clinical episodes of falciparum malaria in E-2025 countries targeting elimination, we characterized the molecular epidemiology of Plasmodium falciparum disease transmission after a clonal outbreak in Ecuador. Here we study disease transmission by documenting the diversity and population structure of the major variant surface antigen of the blood stages of P. falciparum encoded by the var multigene family. We used a high-resolution genotyping method, "varcoding", involving targeted amplicon sequencing to fingerprint the DBLα encoding region of var genes to describe both antigenic var diversity and var repertoire similarity or relatedness in parasite isolates from clinical cases. We identified nine genetic varcodes in 58 P. falciparum isolates causing clinical disease in 2013-2015. Network analyses revealed that four of the varcodes were highly related to the outbreak varcode, with identification of possible diversification of the outbreak parasites by recombination as seen in three of those varcodes. The majority of clinical cases in Ecuador were associated with parasites with highly related or recombinant varcodes to the outbreak clone and due to local transmission rather than recent importation of parasites from other endemic countries. Sharing of types in Ecuadorian varcodes to those sampled in South American varcodes reflects historical parasite importation of some varcodes, especially from Colombia and Peru. Our findings highlight the translational application of varcoding for outbreak surveillance in epidemic/unstable malaria transmission, such as in E-2025 countries, and point to the need for surveillance of local reservoirs of infection in Ecuador to achieve the malaria elimination goal by 2025.

Genetic diversity of Plasmodium falciparum erythrocyte membrane protein 1 in field isolates from central Myanmar.

Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1), encoded by the polymorphic var multigene family, is a highly polymorphic antigen that plays a crucial role in the pathology of malaria. The contribution of the genetic diversity of var toward the immune escape of P. falciparum has not yet been fully elucidated. This study aimed to characterize the diversity of var repertoires by screening P. falciparum Duffy-binding-like α domain (PfDBLα) among field isolates from central Myanmar. Genetic analysis revealed that the D-H segments of var in Myanmar populations have an extensive polymorphic repertoire, with high numbers of unique sequence types in each individual. However, var genes from the global population, including Myanmar, shared close genetic lineages regardless of their geographic origins, indicating that they have not undergone rapid evolutionary changes.

Unravelling var complexity: Relationship between DBLα types and var genes in Plasmodium falciparum.

The enormous diversity and complexity of var genes that diversify rapidly by recombination has led to the exclusion of assembly of these genes from major genome initiatives (e.g., Pf6). A scalable solution in epidemiological surveillance of var genes is to use a small 'tag' region encoding the immunogenic DBLα domain as a marker to estimate var diversity. As var genes diversify by recombination, it is not clear the extent to which the same tag can appear in multiple var genes. This relationship between marker and gene has not been investigated in natural populations. Analyses of in vitro recombination within and between var genes have suggested that this relationship would not be exclusive. Using a dataset of publicly-available assembled var sequences, we test this hypothesis by studying DBLα-var relationships for four study sites in four countries: Pursat (Cambodia) and Mae Sot (Thailand), representing low malaria transmission, and Navrongo (Ghana) and Chikwawa (Malawi), representing high malaria transmission. In all study sites, DBLα-var relationships were shown to be predominantly 1-to-1, followed by a second largest proportion of 1-to-2 DBLα-var relationships. This finding indicates that DBLα tags can be used to estimate not just DBLα diversity but var gene diversity when applied in a local endemic area. Epidemiological applications of this result are discussed.

Patterns of Heterochromatin Transitions Linked to Changes in the Expression of Plasmodium falciparum Clonally Variant Genes.

The survival of malaria parasites in the changing human blood environment largely depends on their ability to alter gene expression by epigenetic mechanisms. The active state of Plasmodium falciparum clonally variant genes (CVGs) is associated with euchromatin characterized by the histone mark H3K9ac, whereas the silenced state is characterized by H3K9me3-based heterochromatin. Expression switches are linked to euchromatin-heterochromatin transitions, but these transitions have not been characterized for the majority of CVGs. To define the heterochromatin distribution patterns associated with the alternative transcriptional states of CVGs, we compared H3K9me3 occupancy at a genome-wide level among several parasite subclones of the same genetic background that differed in the transcriptional state of many CVGs. We found that de novo heterochromatin formation or the complete disruption of a heterochromatin domain is a relatively rare event, and for the majority of CVGs, expression switches can be explained by the expansion or retraction of heterochromatin domains. We identified different modalities of heterochromatin changes linked to transcriptional differences, but despite this complexity, heterochromatin distribution patterns generally enable the prediction of the transcriptional state of specific CVGs. We also found that in some subclones, several var genes were simultaneously in an active state. Furthermore, the heterochromatin levels in the putative regulatory region of the gdv1 antisense noncoding RNA, a regulator of sexual commitment, varied between parasite lines with different sexual conversion rates. IMPORTANCE The malaria parasite P. falciparum is responsible for more than half a million deaths every year. P. falciparum clonally variant genes (CVGs) mediate fundamental host-parasite interactions and play a key role in parasite adaptation to fluctuations in the conditions of the human host. The expression of CVGs is regulated at the epigenetic level by changes in the distribution of a type of chromatin called heterochromatin. Here, we describe at a genome-wide level the changes in the heterochromatin distribution associated with the different transcriptional states of CVGs. Our results also reveal a likely role for heterochromatin at a particular locus in determining the parasite investment in transmission to mosquitoes. Additionally, this data set will enable the prediction of the transcriptional state of CVGs from epigenomic data, which is important for the study of parasite adaptation to the conditions of the host in natural malaria infections.

A coordinated transcriptional switching network mediates antigenic variation of human malaria parasites.

Malaria parasites avoid immune clearance through their ability to systematically alter antigens exposed on the surface of infected red blood cells. This is accomplished by tightly regulated transcriptional control of individual members of a large, multicopy gene family called var and is the key to both the virulence and chronic nature of malaria infections. Expression of var genes is mutually exclusive and controlled epigenetically, however how large populations of parasites coordinate var gene switching to avoid premature exposure of the antigenic repertoire is unknown. Here, we provide evidence for a transcriptional network anchored by a universally conserved gene called var2csa that coordinates the switching process. We describe a structured switching bias that shifts overtime and could shape the pattern of var expression over the course of a lengthy infection. Our results provide an explanation for a previously mysterious aspect of malaria infections and shed light on how parasites possessing a relatively small repertoire of variant antigen-encoding genes can coordinate switching events to limit antigen exposure, thereby maintaining chronic infections.

Evaluation of an ultrasensitive HRP2–based rapid diagnostic test for detection of asymptomatic Plasmodium falciparum parasitaemia among children in western Kenya

BackgroundAccurate detection of asymptomatic malaria parasitaemia in children living in high transmission areas is important for malaria control and reduction programmes that employ screen-and-treat surveillance strategies. Relative to microscopy and conventional rapid diagnostic tests (RDTs), ultrasensitive RDTs (us-RDTs) have demonstrated reduced limits of detection with increased sensitivity to detect parasitaemia in symptomatic individuals. In this study, the performance of the NxTek™ Eliminate Malaria P.f test was compared with traditional microscopy and quantitative polymerase chain reaction (qPCR) testing methods of detection for P. falciparum parasitaemia among asymptomatic children aged 7–14 years living in an area of high malaria transmission intensity in western Kenya.MethodsIn October 2020, 240 healthy children without any reported malaria symptoms were screened for the presence of P. falciparum parasitaemia; 120 children were randomly selected to participate in a follow-up visit at 6–10 weeks. Malaria parasitaemia was assessed by blood-smear microscopy, us-RDT, and qPCR of a conserved var gene sequence from genomic DNA extracted from dried blood spots. Sensitivity, specificity, and predictive values were calculated for field diagnostic methods using qPCR as the gold standard. Comparison of detectable parasite density distributions and area under the curve were also calculated to determine the effectiveness of the us-RDT in detecting asymptomatic infections with low parasite densities.ResultsThe us-RDT detected significantly more asymptomatic P. falciparum infections than microscopy (42.5% vs. 32.2%, P = 0.002). The positive predictive value was higher for microscopy (92.2%) than for us-RDT (82.4%). However, false negative rates were high for microscopy and us-RDT, with negative predictive values of 53.7% and 54.6%, respectively. While us-RDT detected significantly more infections than microscopy overall, the density distribution of detectable infections did not differ (P = 0.21), and qPCR detected significantly more low-density infections than both field methods (P < 0.001, for both comparisons).ConclusionsUs-RDT is more sensitive than microscopy for detecting asymptomatic malaria parasitaemia in children. Though the detectable parasite density distributions by us-RDT in our specific study did not significantly differ from microscopy, the additional sensitivity of the us-RDT resulted in more identified asymptomatic infections in this important group of the population and makes the use of the us-RDT advisable compared to other currently available malaria field detection methods.

Evaluation of two commercial kits and two laboratory-developed qPCR assays compared to LAMP for molecular diagnosis of malaria

BackgroundMalaria is an infectious disease considered as one of the biggest causes of mortality in endemic areas. This life-threatening disease needs to be quickly diagnosed and treated. The standard diagnostic tools recommended by the World Health Organization are thick blood smears microscopy and immuno-chromatographic rapid diagnostic tests. However, these methods lack sensitivity especially in cases of low parasitaemia and non-falciparum infections. Therefore, the need for more accurate and reliable diagnostic tools, such as real-time polymerase chain reaction based methods which have proven greater sensitivity particularly in the screening of malaria, is prominent. This study was conducted at the French National Malaria Reference Centre to assess sensitivity and specificity of two commercial malaria qPCR kits and two in-house developed qPCRs compared to LAMP.Methods183 blood samples received for expertise at the FNMRC were included in this study and were subjected to four different qPCR methods: the Biosynex Ampliquick® Malaria test, the BioEvolution Plasmodium Typage test, the in-house HRM and the in-house TaqMan qPCRs. The specificity and sensitivity of each method and their confidence intervals were determined with the LAMP-based assay Alethia® Malaria as the reference for malaria diagnosis. The accuracy of species diagnosis of the Ampliquick® Malaria test and the two in-house qPCRs was also evaluated using the BioEvolution Plasmodium Typage test as the reference method for species identification.ResultsThe main results showed that when compared to LAMP, a test with excellent diagnostic performances, the two in-house developed qPCRs were the most sensitive (sensitivity at 100% for the in-house TaqMan qPCR and 98.1% for the in-house HRM qPCR), followed by the two commercial kits: the Biosynex Ampliquick® Malaria test (sensitivity at 97.2%) and the BioEvolution Plasmodium Typage (sensitivity at 95.4%). Additionally, with the in-house qPCRs we were able to confirm a Plasmodium falciparum infection in microscopically negative samples that were not detected by commercial qPCR kits. This demonstrates that the var genes of P. falciparum used in these in-house qPCRs are more reliable targets than the 18S sRNA commonly used in most of the developed qPCR methods for malaria diagnosis.ConclusionOverall, these results accentuate the role molecular methods could play in the screening of malaria. This may represent a helpful tool for other laboratories looking to implement molecular diagnosis methods in their routine analysis, which could be essential for the detection and treatment of malaria carriers and even for the eradication of this disease.

Differential Trafficking and Expression of PIR Proteins in Acute and Chronic Plasmodium Infections.

Plasmodium multigene families are thought to play important roles in the pathogenesis of malaria. Plasmodium interspersed repeat (pir) genes comprise the largest multigene family in many Plasmodium species. However, their expression pattern and localisation remain to be elucidated. Understanding protein subcellular localisation is fundamental to reveal the functional importance and cell-cell interactions of the PIR proteins. Here, we use the rodent malaria parasite, Plasmodium chabaudi chabaudi, as a model to investigate the localisation pattern of this gene family. We found that most PIR proteins are co-expressed in clusters during acute and chronic infection; members of the S7 clade are predominantly expressed during the acute-phase, whereas members of the L1 clade dominate the chronic-phase of infection. Using peptide antisera specific for S7 or L1 PIRS, we show that these PIRs have different localisations within the infected red blood cells. S7 PIRs are exported into the infected red blood cell cytoplasm where they are co-localised with parasite-induced host cell modifications termed Maurer’s clefts, whereas L1 PIRs are localised on or close to the parasitophorous vacuolar membrane. This localisation pattern changes following mosquito transmission and during progression from acute- to chronic-phase of infection. The presence of PIRs in Maurer’s clefts, as seen for Plasmodium falciparum RIFIN and STEVOR proteins, might suggest trafficking of the PIRs on the surface of the infected erythrocytes. However, neither S7 nor L1 PIR proteins detected by the peptide antisera are localised on the surface of infected red blood cells, suggesting that they are unlikely to be targets of surface variant-specific antibodies or to be directly involved in adhesion of infected red blood cells to host cells, as described for Plasmodium falciparum VAR proteins. The differences in subcellular localisation of the two major clades of Plasmodium chabaudi PIRs across the blood cycle, and the apparent lack of expression on the red cell surface strongly suggest that the function(s) of this gene family may differ from those of other multigene families of Plasmodium, such as the var genes of Plasmodium falciparum.

Ectopic Expression of Plasmodium vivax vir Genes in P. falciparum Affects Cytoadhesion via Increased Expression of Specific var Genes.

Plasmodium falciparum-infected erythrocytes (PfIEs) adhere to endothelial cell receptors (ECRs) of blood vessels mainly via PfEMP1 proteins to escape elimination via the spleen. Evidence suggests that P. vivax-infected reticulocytes (PvIRs) also bind to ECRs, presumably enabled by VIR proteins, as shown by inhibition experiments and studies with transgenic P. falciparum expressing vir genes. To test this hypothesis, our study investigated the involvement of VIR proteins in cytoadhesion using vir gene-expressing P. falciparum transfectants. Those VIR proteins with a putative transmembrane domain were present in Maurer’s clefts, and some were also present in the erythrocyte membrane. The VIR protein without a transmembrane domain (PVX_050690) was not exported. Five of the transgenic P. falciparum cell lines, including the one expressing PVX_050690, showed binding to CD36. We observed highly increased expression of specific var genes encoding PfEMP1s in all CD36-binding transfectants. These results suggest that ectopic vir expression regulates var expression through a yet unknown mechanism. In conclusion, the observed cytoadhesion of P. falciparum expressing vir genes depended on PfEMP1s, making this experimental unsuitable for characterizing VIR proteins.

Plasmodium falciparum PfRUVBL proteins bind at the TARE region and var gene promoter located in the subtelomeric region.

In order to survive and establish infection, the Plasmodium parasite employs various strategies to evade the host immune response. The var genes family, a repertoire of 60 genes, expresses parasite-specific protein PfEMP1, a variable surface antigen, on the membrane of infected erythrocytes, and by continuously switching the variants of PfEMP1, help the parasite to avoid detection and destruction by the host immune system during the intra-erythrocytic developmental cycle. Although chromatin modifications are recognised to be a prominent phenomenon in regulation of mono-allelic expression of these var genes, the precise histone codes and molecular players and mechanisms guiding these modifications have yet to be unravelled in depth. In this study, we have functionally characterised RUVBL proteins of Plasmodium falciparum and shown that PfMYST (an essential lysine acetyl transferase) and PfRUVBL protein complex occupy the TARE region and var gene promoter in the ring stage of the parasite. Further, we have demonstrated that the PfMYST/PfRUVBL complex interacts with core histones, H3 and H4. Overall the findings of this study add further information by identifying the potential role of epigenetic regulators, PfMYST and PfRUVBL, in the regulation of monoallelic expression of var genes in the malaria parasite.

Vivax Malaria and the Potential Role of the Subtelomeric Multigene vir Superfamily.

Vivax malaria, caused by Plasmodium vivax, remains a public health concern in Central and Southeast Asia and South America, with more than two billion people at risk of infection. Compared to Plasmodium falciparum, P. vivax is considered a benign infection. However, in recent decades, incidences of severe vivax malaria have been confirmed. The P. falciparum erythrocyte membrane protein 1 family encoded by var genes is known as a mediator of severe falciparum malaria by cytoadherence property. Correspondingly, the vir multigene superfamily has been identified as the largest multigene family in P. vivax and is implicated in cytoadherence to endothelial cells and immune response activation. In this review, the functions of vir genes are reviewed in the context of their potential roles in severe vivax malaria.

De Novo Assembly of Plasmodium knowlesi Genomes From Clinical Samples Explains the Counterintuitive Intrachromosomal Organization of Variant SICAvar and kir Multiple Gene Family Members.

Plasmodium knowlesi, a malaria parasite of Old World macaque monkeys, is used extensively to model Plasmodium biology. Recently, P. knowlesi was found in the human population of Southeast Asia, particularly Malaysia. P. knowlesi causes uncomplicated to severe and fatal malaria in the human host with features in common with the more prevalent and virulent malaria caused by Plasmodium falciparum. As such, P. knowlesi presents a unique opportunity to develop experimental translational model systems for malaria pathophysiology informed by clinical data from same-species human infections. Experimental lines of P. knowlesi represent well-characterized genetically stable parasites, and to maximize their utility as a backdrop for understanding malaria pathophysiology, genetically diverse contemporary clinical isolates, essentially wild-type, require comparable characterization. The Oxford Nanopore PCR-free long-read sequencing platform was used to sequence and de novo assemble P. knowlesi genomes from frozen clinical samples. The sequencing platform and assembly pipelines were designed to facilitate capturing data and describing, for the first time, P. knowlesi schizont-infected cell agglutination (SICA) var and Knowlesi-Interspersed Repeats (kir) multiple gene families in parasites acquired from nature. The SICAvar gene family members code for antigenically variant proteins analogous to the virulence-associated P. falciparum erythrocyte membrane protein (PfEMP1) multiple var gene family. Evidence presented here suggests that the SICAvar family members have arisen through a process of gene duplication, selection pressure, and variation. Highly evolving genes including PfEMP1family members tend to be restricted to relatively unstable sub-telomeric regions that drive change with core genes protected in genetically stable intrachromosomal locations. The comparable SICAvar and kir gene family members are counter-intuitively located across chromosomes. Here, we demonstrate that, in contrast to conserved core genes, SICAvar and kir genes occupy otherwise gene-sparse chromosomal locations that accommodate rapid evolution and change. The novel methods presented here offer the malaria research community not only new tools to generate comprehensive genome sequence data from small clinical samples but also new insight into the complexity of clinically important real-world parasites.

Transcriptome Analysis of Plasmodium falciparum Isolates From Benin Reveals Specific Gene Expression Associated With Cerebral Malaria.

Cerebral malaria (CM) is the severest form of Plasmodium falciparum infection. Children under 5 years old are those most vulnerable to CM, and they consequently have the highest risk of malaria-related death. Parasite-associated factors leading to CM are not yet fully elucidated. We therefore sought to characterize the gene expression profile associated with CM, using RNA sequencing data from 15 CM and 15 uncomplicated malaria isolates from Benin. Cerebral malaria parasites displayed reduced circulation times, possibly related to higher cytoadherence capacity. Consistent with the latter, we detected increased var genes abundance in CM isolates. Differential expression analyses showed that distinct transcriptome profiles are signatures of malaria severity. Genes involved in adhesion, excluding variant surface antigens, were dysregulated, supporting the idea of increased cytoadhesion capacity of CM parasites. Finally, we found dysregulated expression of genes in the entry into host pathway that may reflect greater erythrocyte invasion capacity of CM parasites.

Age-specific patterns of DBLα var diversity can explain why residents of high malaria transmission areas remain susceptible to Plasmodium falciparum blood stage infection throughout life

Immunity to Plasmodium falciparum is non-sterilising, thus individuals residing in malaria-endemic areas are at risk of infection throughout their lifetime. Here we seek to find a genomic epidemiological explanation for why residents of all ages harbour blood stage infections despite lifelong exposure to P. falciparum in areas of high transmission. We do this by exploring, for the first known time, the age-specific patterns of diversity of variant antigen encoding (var) genes in the reservoir of infection. Microscopic and submicroscopic P. falciparum infections were analysed at the end of the wet and dry seasons in 2012–2013 for a cohort of 1541 residents aged from 1 to 91 years in an area characterised by high seasonal malaria transmission in Ghana. By sequencing the near ubiquitous Duffy-binding-like alpha domain (DBLα) that encodes immunogenic domains, we defined var gene diversity in an estimated 1096 genomes detected in sequential wet and dry season sampling of this cohort. Unprecedented var (DBLα) diversity was observed in all ages with 42,399 unique var types detected. There was a high degree of maintenance of types between seasons (>40% seen more than once), with many of the same types, especially upsA, appearing multiple times in isolates from different individuals. Children and adolescents were found to be significant reservoirs of var DBLα diversity compared with adults. Var repertoires within individuals were highly variable, with children having more related var repertoires compared to adolescents and adults. Individuals of all ages harboured multiple genomes with var repertoires unrelated to those infecting other hosts. High turnover of parasites with diverse isolate var repertoires was also observed in all ages. These age-specific patterns are best explained by variant-specific immune selection. The observed level of var diversity for the population was then used to simulate the development of variant-specific immunity to the diverse var types under conservative assumptions. Simulations showed that the extent of observed var diversity with limited repertoire relatedness was sufficient to explain why adolescents and adults in this community remain susceptible to blood stage infection, even with multiple genomes.

Plasmodium falciparum Parasite Lines Expressing DC8 and Group A PfEMP1 Bind to Brain, Intestinal, and Kidney Endothelial Cells.

Cytoadhesion of Plasmodium falciparum-infected red blood cells is a virulence determinant associated with microvascular obstruction and organ complications. The gastrointestinal tract is a major site of sequestration in fatal cerebral malaria cases and kidney complications are common in severe malaria, but parasite interactions with these microvascular sites are poorly characterized. To study parasite tropism for different microvascular sites, we investigated binding of parasite lines to primary human microvascular endothelial cells from intestine (HIMEC) and peritubular kidney (HKMEC) sites. Of the three major host receptors for P. falciparum, CD36 had low or negligible expression; endothelial protein C receptor (EPCR) had the broadest constitutive expression; and intercellular adhesion molecule 1 (ICAM-1) was weakly expressed on resting cells and was strongly upregulated by TNF-α on primary endothelial cells from the brain, intestine, and peritubular kidney sites. By studying parasite lines expressing var genes linked to severe malaria, we provide evidence that both the DC8 and Group A EPCR-binding subsets of the P. falciparum erythrocyte membrane protein 1 (PfEMP1) family encodes binding affinity for brain, intestinal, and peritubular kidney endothelial cells, and that DC8 parasite adhesion was partially dependent on EPCR. Collectively, these findings raise the possibility of a brain-gut-kidney binding axis contributing to multi-organ complications in severe malaria.

Varia: a tool for prediction, analysis and visualisation of variable genes

BackgroundParasites use polymorphic gene families to evade the immune system or interact with the host. Assessing the diversity and expression of such gene families in pathogens can inform on the repertoire or host interaction phenotypes of clinical relevance. However, obtaining the sequences and quantifying their expression is a challenge. In Plasmodium falciparum, the highly polymorphic var genes encode the major virulence protein, PfEMP1, which bind a range of human receptors through varying combinations of DBL and CIDR domains. Here we present a tool, Varia, to predict near full-length gene sequences and domain compositions of query genes from database genes sharing short sequence tags. Varia generates output through two complementary pipelines. Varia_VIP returns all putative gene sequences and domain compositions of the query gene from any partial sequence provided, thereby enabling experimental validation of specific genes of interest and detailed assessment of their putative domain structure. Varia_GEM accommodates rapid profiling of var gene expression in complex patient samples from DBLα expression sequence tags (EST), by computing a sample overall transcript profile stratified by PfEMP1 domain types.ResultsVaria_VIP was tested querying sequence tags from all DBL domain types using different search criteria. On average 92% of query tags had one or more 99% identical database hits, resulting in the full-length query gene sequence being identified (> 99% identical DNA > 80% of query gene) among the five most prominent database hits, for ~ 33% of the query genes. Optimized Varia_GEM settings allowed correct prediction of > 90% of domains placed among the four most N-terminal domains, including the DBLα domain, and > 70% of C-terminal domains. With this accuracy, N-terminal domains could be predicted for > 80% of queries, whereas prediction rates of C-terminal domains dropped with the distance from the DBLα from 70 to 40%.ConclusionPrediction of var sequence and domain composition is possible from short sequence tags. Varia can be used to guide experimental validation of PfEMP1 sequences of interest and conduct high-throughput analysis of var type expression in patient samples.

An accurate method for identifying recent recombinants from unaligned sequences.

MotivationRecombination is a fundamental process in molecular evolution, and the identification of recombinant sequences is thus of major interest. However, current methods for detecting recombinants are primarily designed for aligned sequences. Thus, they struggle with analyses of highly diverse genes, such as the var genes of the malaria parasite Plasmodium falciparum, which are known to diversify primarily through recombination.ResultsWe introduce an algorithm to detect recent recombinant sequences from a dataset without a full multiple alignment. Our algorithm can handle thousands of gene-length sequences without the need for a reference panel. We demonstrate the accuracy of our algorithm through extensive numerical simulations; in particular, it maintains its effectiveness in the presence of insertions and deletions. We apply our algorithm to a dataset of 17 335 DBLα types in var genes from Ghana, observing that sequences belonging to the same ups group or domain subclass recombine amongst themselves more frequently, and that non-recombinant DBLα types are more conserved than recombinant ones.Availability and implementationSource code is freely available at https://github.com/qianfeng2/detREC_program.Supplementary information Supplementary data are available at Bioinformatics online.

Alteration of the expression of sirtuins and var genes by heat shock in the malaria parasite Plasmodium falciparum

BackgroundIn Plasmodium falciparum the monoallelic expression of var virulence genes is regulated through epigenetic mechanisms. A study in the Gambia showed that an increase in var genes commonly expressed in patients with severe malaria is associated with fever and high blood lactate. A strong association was demonstrated between the upregulation of PfSir2A and group B var genes. A subsequent study in Kenya extended this association to show a link between elevated expression of PfSir2A and overall var transcript levels. We investigate here the link between heat shock and/or lactate levels on sirtuin and var gene expression levels in vitro. MethodsIn vitro experiments were conducted using laboratory and recently-laboratory-adapted Kenyan isolates of P. falciparum. To investigate a potential cause-and-effect relationship between host stress factors and parasite gene expression, qPCR was used to measure the expression of sirtuins and var genes after highly synchronous cultured parasites had been exposed to 2 h or 6 h of heat shock at 40 °C or elevated lactate. ResultsHeat shock was shown to increase the expression ofPfSir2B in the trophozoites, whereas exposure to lactate was not. After the ring stages were exposed to heat shock and lactate, there was no alteration in the expression of sirtuins and severe-disease-associated upsA and upsB var genes. The association between high blood lactate and sirtuin/var gene expression that was previously observed in vivo appears to be coincidental rather than causative. ConclusionsThis study demonstrates that heat stress in a laboratory and recently-laboratory-adapted isolates of P. falciparum results in a small increase in PfSir2B transcripts in the trophozoite stages only. This finding adds to our understanding of how patient factors can influence the outcome of Plasmodium falciparum infections.

Resetting var Gene Transcription in Plasmodium falciparum.

One of the key mechanisms contributing to the virulence of Plasmodium falciparum is its ability to undergo antigenic switching among antigenically distinct variants of the PfEMP1 adhesive proteins, encoded by the var gene family. To avoid premature exposure of its antigenic repertoire, the parasite transcribes its var genes in a mutually exclusive manner, and switch expression at a very slow rate. This process is epigenetically regulated and it relies on "epigenetic memory," which imprints the single active var gene to remain active for multiple replication cycles. Erasing this epigenetic memory in parasites grown in culture resembles parasites, which egress from the liver. It could therefore be of interest for investigating var switching patterns at the onset of malaria infections. In addition, this procedure could be used for creating heterogeneity of var expression among parasite populations. The methodology described here for resetting of var gene expression is based on promoter titration, also known as molecular sponging.