Genomics and integrated systems biology in Plasmodium falciparum: a path to malaria control and eradication

Genetic polymorphism and evidence of signatures of selection in the Plasmodium falciparum circumsporozoite protein gene in Tanzanian regions with different malaria endemicity

BACKGROUND

In 2021 and 2023, the World Health Organization approved RTS,S/AS01 and R21/Matrix M malaria vaccines, respectively, for routine immunization of children in African countries with moderate to high transmission. These vaccines are made of Plasmodium falciparum circumsporozoite protein (PfCSP), but polymorphisms in the gene raise concerns regarding strain-specific responses and the long-term efficacy of these vaccines. This study assessed the Pfcsp genetic diversity, population structure and signatures of selection among parasites from areas of different malaria transmission intensities in Mainland Tanzania, to generate baseline data before the introduction of the malaria vaccines in the country.

METHODS

The analysis involved 589 whole genome sequences generated by and as part of the MalariaGEN Community Project. The samples were collected between 2013 and January 2015 from five regions of Mainland Tanzania: Morogoro and Tanga (Muheza) (moderate transmission areas), and Kagera (Muleba), Lindi (Nachingwea), and Kigoma (Ujiji) (high transmission areas). Wright's inbreeding coefficient (Fws), Wright's fixation index (FST), principal component analysis, nucleotide diversity, and Tajima's D were used to assess within-host parasite diversity, population structure and natural selection.

RESULTS

Based on Fws (< 0.95), there was high polyclonality (ranging from 69.23% in Nachingwea to 56.9% in Muheza). No population structure was detected in the Pfcsp gene in the five regions (mean FST = 0.0068). The average nucleotide diversity (π), nucleotide differentiation (K) and haplotype diversity (Hd) in the five regions were 4.19, 0.973 and 0.0035, respectively. The C-terminal region of Pfcsp showed high nucleotide diversity at Th2R and Th3R regions. Positive values for the Tajima's D were observed in the Th2R and Th3R regions consistent with balancing selection. The Pfcsp C-terminal sequences revealed 50 different haplotypes (H_1 to H_50), with only 2% of sequences matching the 3D7 strain haplotype (H_50). Conversely, with the NF54 strain, the Pfcsp C-terminal sequences revealed 49 different haplotypes (H_1 to H_49), with only 0.4% of the sequences matching the NF54 strain (Hap_49).

CONCLUSIONS

The findings demonstrate high diversity of the Pfcsp gene with limited population differentiation. The Pfcsp gene showed positive Tajima's D values, consistent with balancing selection for variants within Th2R and Th3R regions. The study observed differences between the intended haplotypes incorporated into the design of RTS,S and R21 vaccines and those present in natural parasite populations. Therefore, additional research is warranted, incorporating other regions and more recent data to comprehensively assess trends in genetic diversity within this important gene. Such insights will inform the choice of alleles to be included in the future vaccines.

Long-Read Sequencing and De Novo Genome Assembly Pipeline of Two Plasmodium falciparum Clones (Pf3D7, PfW2) Using Only the PromethION Sequencer from Oxford Nanopore Technologies without Whole-Genome Amplification

Antimalarial drug resistance has become a real public health problem despite WHO measures. New sequencing technologies make it possible to investigate genomic variations associated with resistant phenotypes at the genome-wide scale. Based on the use of hemisynthetic nanopores, the PromethION technology from Oxford Nanopore Technologies can produce long-read sequences, in contrast to previous short-read technologies used as the gold standard to sequence Plasmodium. Two clones of P. falciparum (Pf3D7 and PfW2) were sequenced in long-read using the PromethION sequencer from Oxford Nanopore Technologies without genomic amplification. This made it possible to create a processing analysis pipeline for human Plasmodium with ONT Fastq only. De novo assembly revealed N50 lengths of 18,488 kb and 17,502 kb for the Pf3D7 and PfW2, respectively. The genome size was estimated at 23,235,407 base pairs for the Pf3D7 clone and 21,712,038 base pairs for the PfW2 clone. The average genome coverage depth was estimated at 787X and 653X for the Pf3D7 and PfW2 clones, respectively. This study proposes an assembly processing pipeline for the human Plasmodium genome using software adapted to large ONT data and the high AT percentage of Plasmodium. This search provides all the parameters which were optimized for use with the software selected in the pipeline.

Genetic polymorphism and evidence of signatures of selection in the Plasmodium falciparum circumsporozoite protein gene in Tanzanian regions with different malaria endemicity

Background

In 2021 and 2023, the World Health Organization approved RTS, S/AS01 and R21/Matrix M malaria vaccines, respectively, for routine immunization of children in African countries with moderate to high transmission. These vaccines are made of Plasmodium falciparum circumsporozoite protein (Pfcsp) but polymorphisms in this gene raises concerns regarding strain-specific responses and the long-term efficacy of these vaccines. This study assessed the Pfcsp genetic diversity, population structure and signatures of selection among parasites from areas of different malaria transmission in mainland Tanzania, to generate baseline data before the introduction of the malaria vaccines in the country.

Methods

The analysis involved 589 whole genome sequences generated by and as part of the MalariaGEN Community Project. The samples were collected between 2013 and January 2015 from five regions of mainland Tanzania: Morogoro and Tanga (Muheza) (moderate transmission areas), and Kagera (Muleba), Lindi (Nachingwea), and Kigoma (Ujiji) (high transmission areas). Wright's inbreeding coefficient (Fws), Wright's fixation index (FST), principal component analysis, nucleotide diversity, and Tajima's D were used to assess within-host parasite diversity, population structure and natural selection.

Results

Based on Fws (< 0.95), there was high polyclonality (ranged from 69.23% in Nachingwea to 56.9% in Muheza). No population structure was detected in the Pfcsp gene in the five regions (mean FST= 0.0068). The average nucleotide diversity (π), nucleotide differentiation (K) and haplotype diversity (Hd) in the five regions were 4.19, 0.973 and 0.0035, respectively. The C-terminal region of Pfcsp showed high nucleotide diversity at Th2R and Th3R regions. Positive values for the Tajima's D were observed in the Th2R and Th3R regions consistent with balancing selection. The Pfcsp C-terminal sequences had 50 different haplotypes (H_1 to H_50) and only 2% of sequences matched the 3D7 strain haplotype (H_50).

Conclusions

The findings demonstrate high diversity of the Pfcsp gene with limited population differentiation. The Pfcsp gene showed positive Tajima's D values for parasite populations, consistent with balancing selection for variants within Th2R and Th3R regions. This data is consistent with other studies conducted across Africa and worldwide, which demonstrate low 3D7 haplotypes and little population structure. Therefore, additional research is warranted, incorporating other regions and more recent data to comprehensively assess trends in genetic diversity within this important gene. Such insights will inform the choice of alleles to be included in the future vaccines.

The potential of genomics for infectious disease forecasting

Genomic technologies have led to tremendous gains in understanding how pathogens function, evolve and interact. Pathogen diversity is now measurable at high precision and resolution, in part because over the past decade, sequencing technologies have increased in speed and capacity, at decreased cost. Alongside this, the use of models that can forecast emergence and size of infectious disease outbreaks has risen, highlighted by the coronavirus disease 2019 pandemic but also due to modelling advances that allow for rapid estimates in emerging outbreaks to inform monitoring, coordination and resource deployment. However, genomics studies have remained largely retrospective. While they contain high-resolution views of pathogen diversification and evolution in the context of selection, they are often not aligned with designing interventions. This is a missed opportunity because pathogen diversification is at the core of the most pressing infectious public health challenges, and interventions need to take the mechanisms of virulence and understanding of pathogen diversification into account. In this Perspective, we assess these converging fields, discuss current challenges facing both surveillance specialists and modellers who want to harness genomic data, and propose next steps for integrating longitudinally sampled genomic data with statistical learning and interpretable modelling to make reliable predictions into the future.

Temporal dynamics of Plasmodium falciparum population in Metehara, east-central Ethiopia

Background

Plasmodium falciparum is the most serious, genetically most complex and fastest-evolving malaria parasite. Information on genetic diversity of this parasite would guide policy decision and malaria elimination endeavors. This study explored the temporal dynamics of P. falciparum population in two time points in Metehara, east-central Ethiopia.

Methods

The participants were quantitative real-time polymerase chain reaction-confirmed patients who were recruited for uncomplicated falciparum malaria therapeutic efficacy test in 2015 and 2019. Dry blood spot samples were analysed by the nested PCR to genotype P. falciparum merozoite surface protein (msp1, msp2) and glutamate-rich protein (glurp) genes.

Results

While msp1, msp2 and glurp genotypes were successfully detected in 26(89.7%), 24(82.8%) and 14(48.3%) of 2015 samples (n = 29); the respective figures for 2019 (n = 41) were 31(68.3%), 39(95.1%), 25(61.0%). In 2015, the frequencies of K1, MAD20 and RO33 allelic families of msp1, and FC27 and IC/3D7 of msp2 were 19(73.1%), 8(30.6%), 14(53.8%), 21(87.5%), 12(50.5%); and in 2019 it was 15(48.4%), 19(61.3%), 15(48.4%), 30(76.9%), 27(69.2%) respectively. MAD20 has shown dominance over both K1 and RO33 in 2019 compared to the proportion in 2015. Similarly, although FC27 remained dominant, there was shifting trend in the frequency of IC/3D7 from 50.5% in 2015 to 69.2% in 2019. The multiplicity of infection (MOI) and expected heterozygosity index (He) in 2015 and 2019 were respectively [1.43 ± 0.84] and [1.15 ± 0.91], 0.3 and 0.03 for msp1. However, there was no significant association between MOI and age or parasitaemia in both time points.

Conclusion

The lower genetic diversity in P. falciparum population in the two time points and overall declining trend as demonstrated by the lower MOI and He may suggest better progress in malaria control in Metehara. But, the driving force and selective advantage of switching to MAD20 dominance over the other two msp1 allelic families, and the dynamics within msp2 alleles needs further investigation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-022-04277-5.

In-Silico Functional Annotation of Plasmodium falciparum Hypothetical Proteins to Identify Novel Drug Targets

Plasmodium falciparum is one of the plasmodium species responsible for the majority of life-threatening malaria cases. The current antimalarial therapies are becoming less effective due to growing drug resistance, leading to the urgent requirement for alternative and more effective antimalarial drugs or vaccines. To facilitate the novel drug discovery or vaccine development efforts, recent advances in sequencing technologies provide valuable information about the whole genome of the parasite, yet a lot more needs to be deciphered due to its incomplete proteome annotation. Surprisingly, out of the 5,389 proteins currently annotated in the Plasmodium falciparum 3D7 strain, 1,626 proteins (∼30% data) are annotated as hypothetical proteins. In parasite genomic studies, the challenge to annotate hypothetical proteins is often ignored, which may obscure the crucial information related to the pathogenicity of the parasite. In this study, we attempt to characterize hypothetical proteins of the parasite to identify novel drug targets using a computational pipeline. The study reveals that out of the overall pool of the hypothetical proteins, 266 proteins have conserved functional signatures. Furthermore, the pathway analysis of these proteins revealed that 23 proteins have an essential role in various biochemical, signalling and metabolic pathways. Additionally, all the proteins (266) were subjected to computational structure analysis. We could successfully model 11 proteins. We validated and checked the structural stability of the models by performing molecular dynamics simulation. Interestingly, eight proteins show stable conformations, and seven proteins are specific for Plasmodium falciparum, based on homology analysis. Lastly, mapping the seven shortlisted hypothetical proteins on the Plasmodium falciparum protein-protein interaction network revealed 3,299 nodes and 2,750,692 edges. Our study revealed interesting functional details of seven hypothetical proteins of the parasite, which help learn more about the less-studied molecules and their interactions, providing valuable clues to unravel the role of these proteins via future experimental validation.

Whole Genome Sequencing Contributions and Challenges in Disease Reduction Focused on Malaria

Malaria elimination remains an important goal that requires the adoption of sophisticated science and management strategies in the era of the COVID-19 pandemic. The advent of next generation sequencing (NGS) is making whole genome sequencing (WGS) a standard today in the field of life sciences, as PCR genotyping and targeted sequencing provide insufficient information compared to the whole genome. Thus, adapting WGS approaches to malaria parasites is pertinent to studying the epidemiology of the disease, as different regions are at different phases in their malaria elimination agenda. Therefore, this review highlights the applications of WGS in disease management, challenges of WGS in controlling malaria parasites, and in furtherance, provides the roles of WGS in pursuit of malaria reduction and elimination. WGS has invaluable impacts in malaria research and has helped countries to reach elimination phase rapidly by providing required information needed to thwart transmission, pathology, and drug resistance. However, to eliminate malaria in sub-Saharan Africa (SSA), with high malaria transmission, we recommend that WGS machines should be readily available and affordable in the region.

Characterizing the genomic variation and population dynamics of Plasmodium falciparum malaria parasites in and around Lake Victoria, Kenya

Characterising the genomic variation and population dynamics of Plasmodium falciparum parasites in high transmission regions of Sub-Saharan Africa is crucial to the long-term efficacy of regional malaria elimination campaigns and eradication. Whole-genome sequencing (WGS) technologies can contribute towards understanding the epidemiology and structural variation landscape of P. falciparum populations, including those within the Lake Victoria basin, a region of intense transmission. Here we provide a baseline assessment of the genomic diversity of P. falciparum isolates in the Lake region of Kenya, which has sparse genetic data. Lake region isolates are placed within the context of African-wide populations using Illumina WGS data and population genomic analyses. Our analysis revealed that P. falciparum isolates from Lake Victoria form a cluster within the East African parasite population. These isolates also appear to have distinct ancestral origins, containing genome-wide signatures from both Central and East African lineages. Known drug resistance biomarkers were observed at similar frequencies to those of East African parasite populations, including the S160N/T mutation in the pfap2mu gene, which has been associated with delayed clearance by artemisinin-based combination therapy. Overall, our work provides a first assessment of P. falciparum genetic diversity within the Lake Victoria basin, a region targeting malaria elimination.

Detecting sequence variants in clinically important protozoan parasites

Second and third generation sequencing methods are crucial for population genetic studies, and variant detection is a popular approach for exploiting this sequence data. While mini- and microsatellites are historically useful markers for studying important Protozoa such as Toxoplasma and Plasmodium spp., detecting non-repetitive variants such as those found in genes can be fundamental to investigating a pathogen's biology. These variants, namely single nucleotide polymorphisms and insertions and deletions, can help elucidate the genetic basis of an organism's pathogenicity, identify selective pressures, and resolve phylogenetic relationships. They also have the added benefit of possessing a comparatively low mutation rate, which contributes to their stability. However, there is a plethora of variant analysis tools with nuanced pipelines and conflicting recommendations for best practise, which can be confounding. This lack of standardisation means that variant analysis requires careful parameter optimisation, an understanding of its limitations, and the availability of high quality data. This review explores the value of variant detection when applied to non-model organisms such as clinically important protozoan pathogens. The limitations of current methods are discussed, including special considerations that require the end-users' attention to ensure that the results generated are reproducible, and the biological conclusions drawn are valid.

High Plasmodium falciparum genetic diversity and temporal stability despite control efforts in high transmission settings along the international border between Zambia and the Democratic Republic of the Congo

Background

While the utility of parasite genotyping for malaria elimination has been extensively documented in low to moderate transmission settings, it has been less well-characterized in holoendemic regions. High malaria burden settings have received renewed attention acknowledging their critical role in malaria elimination. Defining the role for parasite genomics in driving these high burden settings towards elimination will enhance future control programme planning.

Methods

Amplicon deep sequencing was used to characterize parasite population genetic diversity at polymorphic Plasmodium falciparum loci, Pfama1 and Pfcsp, at two timepoints in June–July 2016 and January–March 2017 in a high transmission region along the international border between Luapula Province, Zambia and Haut-Katanga Province, the Democratic Republic of the Congo (DRC).

Results

High genetic diversity was observed across both seasons and in both countries. No evidence of population structure was observed between parasite populations on either side of the border, suggesting that this region may be one contiguous transmission zone. Despite a decline in parasite prevalence at the sampling locations in Haut-Katanga Province, no genetic signatures of a population bottleneck were detected, suggesting that larger declines in transmission may be required to reduce parasite genetic diversity. Analysing rare variants may be a suitable alternative approach for detecting epidemiologically important genetic signatures in highly diverse populations; however, the challenge is distinguishing true signals from potential artifacts introduced by small sample sizes.

Conclusions

Continuing to explore and document the utility of various parasite genotyping approaches for understanding malaria transmission in holoendemic settings will be valuable to future control and elimination programmes, empowering evidence-based selection of tools and methods to address pertinent questions, thus enabling more efficient resource allocation.

Plasmodium Genomics and Genetics: New Insights into Malaria Pathogenesis, Drug Resistance, Epidemiology, and Evolution

Protozoan Plasmodium parasites are the causative agents of malaria, a deadly disease that continues to afflict hundreds of millions of people every year. Infections with malaria parasites can be asymptomatic, with mild or severe symptoms, or fatal, depending on many factors such as parasite virulence and host immune status. Malaria can be treated with various drugs, with artemisinin-based combination therapies (ACTs) being the first-line choice. Recent advances in genetics and genomics of malaria parasites have contributed greatly to our understanding of parasite population dynamics, transmission, drug responses, and pathogenesis. However, knowledge gaps in parasite biology and host-parasite interactions still remain. Parasites resistant to multiple antimalarial drugs have emerged, while advanced clinical trials have shown partial efficacy for one available vaccine. Here we discuss genetic and genomic studies of Plasmodium biology, host-parasite interactions, population structures, mosquito infectivity, antigenic variation, and targets for treatment and immunization. Knowledge from these studies will advance our understanding of malaria pathogenesis, epidemiology, and evolution and will support work to discover and develop new medicines and vaccines.

Performance evaluation of laboratory professionals on malaria microscopy at health facilities in Bahir Dar city administration, Northwest Ethiopia

BACKGROUND

Microscopic diagnosis of Giemsa stained thick and thin blood films by skilled microscopists has remained the gold standard laboratory method for the diagnosis of malaria. However, there is a scarcity of qualified laboratory professionals for correctly diagnosing malaria using microscopy. The aim of this study was to evaluate the performance of laboratory professionals on malaria microscopy at health facilities in Bahir Dar city administration, Northwest Ethiopia.

METHODS

A cross-sectional study was conducted from January to March 2017 in Bahir Dar City. A total of 87 laboratory professionals participated in the selected health facilities, with a response rate of (100%). Standardized pre-validated slide panels and questionnaire were distributed to laboratory professionals by the principal investigator. The panel slides were comprised of 5 positives and 3 negative blood films. The laboratory professionals were requested to report the parasite density using semi-quantitative (+) and per micro-liter methods. Their performances of slide readings were compared with the experts' readings. Agreement in detecting malaria parasites between laboratory professionals and expert was estimated using the kappa score.

RESULTS

The overall sensitivity and specificity of the laboratory professionals in detecting malaria parasites were 83% and 97%, respectively. Similarly, positive predictive values of 98.1% and negative predictive values of 77.7% were reported. The percent agreement between laboratory professionals and expert microscopist in the detection of malaria parasites was 88.5% with a Kappa index of 0.78. Percent agreement in species identification and reporting of Pf/Pv mixed infections were 27.2% and 22.4%, respectively. About 62.2% of the laboratory professionals reported parasite density using semi-quantitative method. While none of them reported per micro-liter method.

CONCLUSIONS

The current study showed that laboratory professionals had good performance in parasite detection. However, poor performance was seen in both species identification and reporting of Pf/Pv mixed infections.

Developmental Sensitivity in Schistosoma mansoni to Puromycin To Establish Drug Selection of Transgenic Schistosomes

Schistosomiasis is considered the most important disease caused by helminth parasites, in terms of morbidity and mortality. Tools to facilitate gain- and loss-of-function approaches can be expected to precipitate the discovery of novel interventions, and drug selection of transgenic schistosomes would facilitate the establishment of stable lines of engineered parasites. Sensitivity of developmental stages of schistosomes to the aminonucleoside antibiotic puromycin was investigated. For the schistosomulum and sporocyst stages, viability was quantified by fluorescence microscopy following dual staining with fluorescein diacetate and propidium iodine. By 6 days in culture, the 50% lethal concentration (LC₅₀) for schistosomula was 19 μg/ml whereas the sporocysts were 45-fold more resilient. Puromycin potently inhibited the development of in vitro-laid eggs (LC₅₀, 68 ng/ml) but was less effective against liver eggs (LC₅₀, 387 μg/ml). Toxicity for adult stages was evaluated using the xCELLigence-based, real-time motility assay (xWORM), which revealed LC₅₀s after 48 h of 4.9 and 17.3 μg/ml for male and female schistosomes, respectively. Also, schistosomula transduced with pseudotyped retrovirus encoding the puromycin resistance marker were partially rescued when cultured in the presence of the antibiotic. Together, these findings will facilitate selection on puromycin of transgenic schistosomes and the enrichment of cultures of transgenic eggs and sporocysts to facilitate the establishment of schistosome transgenic lines. Streamlining schistosome transgenesis with drug selection will open new avenues to understand parasite biology and hopefully lead to new interventions for this neglected tropical disease.

Quantitative chromatin proteomics reveals a dynamic histone post-translational modification landscape that defines asexual and sexual Plasmodium falciparum parasites

Gene expression in Plasmodia integrates post-transcriptional regulation with epigenetic marking of active genomic regions through histone post-translational modifications (PTMs). To generate insights into the importance of histone PTMs to the entire asexual and sexual developmental cycles of the parasite, we used complementary and comparative quantitative chromatin proteomics to identify and functionally characterise histone PTMs in 8 distinct life cycle stages of P. falciparum parasites. ~500 individual histone PTMs were identified of which 106 could be stringently validated. 46 individual histone PTMs and 30 co-existing PTMs were fully quantified with high confidence. Importantly, 15 of these histone PTMs are novel for Plasmodia (e.g. H3K122ac, H3K27me3, H3K56me3). The comparative nature of the data revealed a highly dynamic histone PTM landscape during life cycle development, with a set of histone PTMs (H3K4ac, H3K9me1 and H3K36me2) displaying a unique and conserved abundance profile exclusively during gametocytogenesis (P < 0.001). Euchromatic histone PTMs are abundant during schizogony and late gametocytes; heterochromatic PTMs mark early gametocytes. Collectively, this data provides the most accurate, complete and comparative chromatin proteomic analyses of the entire life cycle development of malaria parasites. A substantial association between histone PTMs and stage-specific transition provides insights into the intricacies characterising Plasmodial developmental biology.

Cis regulatory motifs and antisense transcriptional control in the apicomplexan Theileria parva

Background

Theileria parva is an intracellular parasite that causes a lymphoproliferative disease in cattle. It does so by inducing cancer-like phenotypes in the host cells it infects, although the molecular and regulatory mechanisms involved remain poorly understood. RNAseq data, and the resulting updated genome annotation now available for this parasite, offer an unprecedented opportunity to characterize the genomic features associated with gene regulation in this species. Our previous analyses revealed a T. parva genome even more gene-dense than previously thought, with many adjacent loci overlapping each other, not only at the level of untranslated sequences (UTRs) but even in coding sequences.

Results

Despite this compactness, Theileria intergenic regions show a pattern of size distribution indicative of monocistronic gene transcription. Three previously described motifs are conserved among Theileria species and highly prevalent in promoter regions near or at the transcription start sites. We found novel motifs at many transcription termination sites, as well as upstream of parasite genes thought to be critical for host transformation. Adjacent genes that could be regulated by antisense transcription from an overlapping transcriptional unit are syntenic between T. parva and P. falciparum at a frequency higher than expected by chance, suggesting the presence of common, and evolutionary old, regulatory mechanisms in the phylum Apicomplexa.

Conclusions

We propose a model of transcription with conserved sense and antisense transcription from a few taxonomically ubiquitous and several species-specific promoter motifs. Interestingly, the gene networks regulated by conserved promoters are themselves, in most cases, not conserved between species or genera.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2444-5) contains supplementary material, which is available to authorized users.

In Silico Designing and Analysis of Inhibitors against Target Protein Identified through Host-Pathogen Protein Interactions in Malaria

Malaria, a life-threatening blood disease, has been a major concern in the field of healthcare. One of the severe forms of malaria is caused by the parasite Plasmodium falciparum which is initiated through protein interactions of pathogen with the host proteins. It is essential to analyse the protein-protein interactions among the host and pathogen for better understanding of the process and characterizing specific molecular mechanisms involved in pathogen persistence and survival. In this study, a complete protein-protein interaction network of human host and Plasmodium falciparum has been generated by integration of the experimental data and computationally predicting interactions using the interolog method. The interacting proteins were filtered according to their biological significance and functional roles. α-tubulin was identified as a potential protein target and inhibitors were designed against it by modification of amiprophos methyl. Docking and binding affinity analysis showed two modified inhibitors exhibiting better docking scores of −10.5 kcal/mol and −10.43 kcal/mol and an improved binding affinity of −83.80 kJ/mol and −98.16 kJ/mol with the target. These inhibitors can further be tested and validated in vivo for their properties as an antimalarial drug.

The mechanics of malaria parasite invasion of the human erythrocyte - towards a reassessment of the host cell contribution

Despite decades of research, we still know little about the mechanics of Plasmodium host cell invasion. Fundamentally, while the essential or non‐essential nature of different parasite proteins is becoming clearer, their actual function and how each comes together to govern invasion are poorly understood. Furthermore, in recent years an emerging world view is shifting focus away from the parasite actin–myosin motor being the sole force responsible for entry to an appreciation of host cell dynamics and forces and their contribution to the process. In this review, we discuss merozoite invasion of the erythrocyte, focusing on the complex set of pre‐invasion events and how these might prime the red cell to facilitate invasion. While traditionally parasite interactions at this stage have been viewed simplistically as mediating adhesion only, recent work makes it apparent that by interacting with a number of host receptors and signalling pathways, combined with secretion of parasite‐derived lipid material, that the merozoite may initiate cytoskeletal re‐arrangements and biophysical changes in the erythrocyte that greatly reduce energy barriers for entry. Seen in this light Plasmodium invasion may well turn out to be a balance between host and parasite forces, much like that of other pathogen infection mechanisms.

Analysis of nucleosome positioning landscapes enables gene discovery in the human malaria parasite Plasmodium falciparum

Background

Plasmodium falciparum, the deadliest malaria-causing parasite, has an extremely AT-rich (80.7 %) genome. Because of high AT-content, sequence-based annotation of genes and functional elements remains challenging. In order to better understand the regulatory network controlling gene expression in the parasite, a more complete genome annotation as well as analysis tools adapted for AT-rich genomes are needed. Recent studies on genome-wide nucleosome positioning in eukaryotes have shown that nucleosome landscapes exhibit regular characteristic patterns at the 5’- and 3’-end of protein and non-protein coding genes. In addition, nucleosome depleted regions can be found near transcription start sites. These unique nucleosome landscape patterns may be exploited for the identification of novel genes. In this paper, we propose a computational approach to discover novel putative genes based exclusively on nucleosome positioning data in the AT-rich genome of P. falciparum.

Results

Using binary classifiers trained on nucleosome landscapes at the gene boundaries from two independent nucleosome positioning data sets, we were able to detect a total of 231 regions containing putative genes in the genome of Plasmodium falciparum, of which 67 highly confident genes were found in both data sets. Eighty-eight of these 231 newly predicted genes exhibited transcription signal in RNA-Seq data, indicative of active transcription. In addition, 20 out of 21 selected gene candidates were further validated by RT-PCR, and 28 out of the 231 genes showed significant matches using BLASTN against an expressed sequence tag (EST) database. Furthermore, 108 (47 %) out of the 231 putative novel genes overlapped with previously identified but unannotated long non-coding RNAs. Collectively, these results provide experimental validation for 163 predicted genes (70.6 %). Finally, 73 out of 231 genes were found to be potentially translated based on their signal in polysome-associated RNA-Seq representing transcripts that are actively being translated.

Conclusion

Our results clearly indicate that nucleosome positioning data contains sufficient information for novel gene discovery. As distinct nucleosome landscapes around genes are found in many other eukaryotic organisms, this methodology could be used to characterize the transcriptome of any organism, especially when coupled with other DNA-based gene finding and experimental methods (e.g., RNA-Seq).

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2214-9) contains supplementary material, which is available to authorized users.

Metabolomics in the fight against malaria

Metabolomics uses high-resolution mass spectrometry to provide a chemical fingerprint of thousands of metabolites present in cells, tissues or body fluids. Such metabolic phenotyping has been successfully used to study various biologic processes and disease states. High-resolution metabolomics can shed new light on the intricacies of host-parasite interactions in each stage of the Plasmodium life cycle and the downstream ramifications on the host's metabolism, pathogenesis and disease. Such data can become integrated with other large datasets generated using top-down systems biology approaches and be utilised by computational biologists to develop and enhance models of malaria pathogenesis relevant for identifying new drug targets or intervention strategies. Here, we focus on the promise of metabolomics to complement systems biology approaches in the quest for novel interventions in the fight against malaria. We introduce the Malaria Host-Pathogen Interaction Center (MaHPIC), a new systems biology research coalition. A primary goal of the MaHPIC is to generate systems biology datasets relating to human and non-human primate (NHP) malaria parasites and their hosts making these openly available from an online relational database. Metabolomic data from NHP infections and clinical malaria infections from around the world will comprise a unique global resource.

Heavy path mining of protein-protein associations in the malaria parasite

Annotating and understanding the function of proteins and other elements in a genome can be difficult in the absence of a well-studied and evolutionarily close relative. The causative agent of malaria, one of the oldest and most deadly global infectious diseases, is a good example of this problem. The burden of malaria is huge and there is a pressing need for new, more effective antimalarial strategies. However, techniques such as homology-dependent annotation transfer are severely impaired in this parasite because there are no well-understood close relatives. To circumvent this approach we developed a network-based method that uses a heavy path network-mining algorithm. We uncovered the protein-protein associations that are implicated in important cellular processes including genome integrity, DNA repair, transcriptional regulation, invasion, and pathogenesis, thus demonstrating the utility of this method. The URL of the source code for super-sequence mining method is http://www.cs.utsa.edu/~korkmaz/research/heavy-path-mining/.

Plasmodium knowlesi gene expression differs in ex vivo compared to in vitro blood-stage cultures

BACKGROUND

Plasmodium knowlesi is one of five Plasmodium species known to cause malaria in humans and can result in severe illness and death. While a zoonosis in humans, this simian malaria parasite species infects macaque monkeys and serves as an experimental model for in vivo, ex vivo and in vitro studies. It has underpinned malaria discoveries relating to host-pathogen interactions, the immune response and immune evasion strategies. This study investigated differences in P. knowlesi gene expression in samples from ex vivo and in vitro cultures.

METHODS

Gene expression profiles were generated using microarrays to compare the stage-specific transcripts detected for a clone of P. knowlesi propagated in the blood of a rhesus macaque host and then grown in an ex-vivo culture, and the same clone adapted to long-term in vitro culture. Parasite samples covering one blood-stage cycle were analysed at four-hour intervals. cDNA was generated and hybridized to an oligoarray representing the P. knowlesi genome. Two replicate experiments were developed from in vitro cultures. Expression values were filtered, normalized, and analysed using R and Perl language and applied to a sine wave model to determine changes in equilibrium and amplitude. Differentially expressed genes from ex vivo and in vitro time points were detected using limma R/Bioconductor and gene set enrichment analysis (GSEA).

RESULTS

Major differences were noted between the ex vivo and in vitro time courses in overall gene expression and the length of the cycle (25.5 hours ex vivo; 33.5 hours in vitro). GSEA of genes up-regulated ex vivo showed an enrichment of various genes including SICAvar, ribosomal- associated and histone acetylation pathway genes. In contrast, certain genes involved in metabolism and cell growth, such as porphobilinogen deaminase and tyrosine phosphatase, and one SICAvar gene, were significantly up-regulated in vitro.

CONCLUSIONS

This study demonstrates how gene expression in P. knowlesi blood-stage parasites can differ dramatically depending on whether the parasites are grown in vivo, with only one cycle of development ex vivo, or as an adapted isolate in long-term in vitro culture. These data bring emphasis to the importance of studying the parasite, its biology and disease manifestations in the context of the host.

Aspartic protease inhibitors containing tertiary alcohol transition-state mimics

Aspartic proteases (APs) are a class of enzymes engaged in the proteolytic digestion of peptide substrates. APs play important roles in physiological and infectious pathways, making them plausible drug targets. For instance in the treatment of HIV infections, access to an efficient combination of protease and reverse transcriptase inhibitors have changed a terminal illness to a chronic but manageable disease. However, the benefits have been limited due to the emergence of drug resistant viral strains, poor pharmacokinetic properties of peptidomimetic inhibitors and adverse effects associated with the treatment. In the 1980s, D. Rich and co-workers proposed a novel strategy for the development of AP inhibitors by replacing the secondary hydroxyl group with a tertiary alcohol as part of the transition state (TS) mimicking moiety. This strategy has been extensively explored over the last decade with a common belief that masking of the polar group, e.g. by intramolecular hydrogen bonding, has the potential to enhance transcellular transport. This is the first review presenting the advances of AP inhibitors comprising a tertiary hydroxyl group. The inhibitors have been classified into different tert-hydroxy TS mimics and their design strategies, synthesis, biological activities, structure-activity-relationships and X-ray structures are discussed.

Multiple dimensions of epigenetic gene regulation in the malaria parasite Plasmodium falciparum: gene regulation via histone modifications, nucleosome positioning and nuclear architecture in P. falciparum

Plasmodium falciparum is the most deadly human malarial parasite, responsible for an estimated 207 million cases of disease and 627,000 deaths in 2012. Recent studies reveal that the parasite actively regulates a large fraction of its genes throughout its replicative cycle inside human red blood cells and that epigenetics plays an important role in this precise gene regulation. Here, we discuss recent advances in our understanding of three aspects of epigenetic regulation in P. falciparum: changes in histone modifications, nucleosome occupancy and the three-dimensional genome structure. We compare these three aspects of the P. falciparum epigenome to those of other eukaryotes, and show that large-scale compartmentalization is particularly important in determining histone decomposition and gene regulation in P. falciparum. We conclude by presenting a gene regulation model for P. falciparum that combines the described epigenetic factors, and by discussing the implications of this model for the future of malaria research.

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.

Transcript maturation in apicomplexan parasites

The complex life cycles of apicomplexan parasites are associated with dynamic changes of protein repertoire. In Toxoplasma gondii, global analysis of gene expression demonstrates that dynamic changes in mRNA levels unfold in a serial cascade during asexual replication and up to 50% of encoded genes are unequally expressed in development. Recent studies indicate transcription and mRNA processing have important roles in fulfilling the 'just-in-time' delivery of proteins to parasite growth and development. The prominence of post-transcriptional mechanisms in the Apicomplexa was demonstrated by mechanistic studies of the critical RNA-binding proteins and regulatory kinases. However, it is still early in our understanding of how transcription and post-transcriptional mechanisms are balanced to produce adequate numbers of specialized forms that is required to complete the parasite life cycle.

A systems biology approach to the effect of aging, immunosenescence and vaccine response

Aging can lead to immunosenescence, which dramatically impairs the hosts' ability to develop protective immune responses to vaccine antigens. Reasons for this are not well understood. This topic's importance is reflected in the increases in morbidity and mortality due to infectious diseases among elderly persons, a population growing in size globally, and the significantly lower adaptive immune responses generated to vaccines in this population. Here, we endeavor to summarize the existing data on the genetic and immunologic correlates of immunosenescence with respect to vaccine response. We cover how the application of systems biology can advance our understanding of vaccine immunosenescence, with a view toward how such information could lead to strategies to overcome the lower immunogenicity of vaccines in the elderly.

Analysis of antibodies to newly described Plasmodium falciparum merozoite antigens supports MSPDBL2 as a predicted target of naturally acquired immunity

Prospective studies continue to identify malaria parasite genes with particular patterns of polymorphism which indicate they may be under immune selection, and the encoded proteins require investigation. Sixteen new recombinant protein reagents were designed to characterize three such polymorphic proteins expressed in Plasmodium falciparum schizonts and merozoites: MSPDBL1 (also termed MSP3.4) and MSPDBL2 (MSP3.8), which possess Duffy binding-like (DBL) domains, and SURFIN4.2, encoded by a member of the surface-associated interspersed (surf) multigene family. After testing the antigenicities of these reagents by murine immunization and parasite immunofluorescence, we analyzed naturally acquired antibody responses to the antigens in two cohorts in coastal Kenya in which the parasite was endemic (Chonyi [n = 497] and Ngerenya [n = 461]). As expected, the prevalence and levels of serum antibodies increased with age. We then investigated correlations with subsequent risk of clinical malaria among children <11 years of age during 6 months follow-up surveillance. Antibodies to the polymorphic central region of MSPDBL2 were associated with reduced risk of malaria in both cohorts, with statistical significance remaining for the 3D7 allelic type after adjustment for individuals' ages in years and antibody reactivity to whole-schizont extract (Chonyi, risk ratio, 0.51, and 95% confidence interval [CI], 0.28 to 0.93; Ngerenya, risk ratio, 0.38, and 95% CI, 0.18 to 0.82). For the MSPDBL1 Palo Alto allelic-type antigen, there was a protective association in one cohort (Ngerenya, risk ratio, 0.53, and 95% CI, 0.32 to 0.89), whereas the other antigens showed no protective associations after adjustment. These findings support the prediction that antibodies to the polymorphic region of MSPDBL2 contribute to protective immunity.

Plasmodium vivax: modern strategies to study a persistent parasite's life cycle

Plasmodium vivax has unique attributes to support its survival in varying ecologies and climates. These include hypnozoite forms in the liver, an invasion preference for reticulocytes, caveola-vesicle complex structures in the infected erythrocyte membrane and rapidly forming and circulating gametocytes. These characteristics make this species very different from P. falciparum. Plasmodium cynomolgi and other related simian species have identical biology and can serve as informative models of P. vivax infections. Plasmodium vivax and its model parasites can be grown in non-human primates (NHP), and in short-term ex vivo cultures. For P. vivax, in the absence of in vitro culture systems, these models remain highly relevant side by side with human clinical studies. While post-genomic technologies allow for greater exploration of P. vivax-infected blood samples from humans, these come with restrictions. Two advantages of NHP models are that infections can be experimentally tailored to address hypotheses, including genetic manipulation. Also, systems biology approaches can capitalise on computational biology combined with set experimental infection periods and protocols, which may include multiple sampling times, different types of samples, and the broad use of "omics" technologies. Opportunities for research on vivax malaria are increasing with the use of existing and new methodological strategies in combination with modern technologies.

Improving N-terminal protein annotation of Plasmodium species based on signal peptide prediction of orthologous proteins

BACKGROUND

Signal peptide is one of the most important motifs involved in protein trafficking and it ultimately influences protein function. Considering the expected functional conservation among orthologs it was hypothesized that divergence in signal peptides within orthologous groups is mainly due to N-terminal protein sequence misannotation. Thus, discrepancies in signal peptide prediction of orthologous proteins were used to identify misannotated proteins in five Plasmodium species.

METHODS

Signal peptide (SignalP) and orthology (OrthoMCL) were combined in an innovative strategy to identify orthologous groups showing discrepancies in signal peptide prediction among their protein members (Mixed groups). In a comparative analysis, multiple alignments for each of these groups and gene models were visually inspected in search of misannotated proteins and, whenever possible, alternative gene models were proposed. Thresholds for signal peptide prediction parameters were also modified to reduce their impact as a possible source of discrepancy among orthologs. Validation of new gene models was based on RT-PCR (few examples) or on experimental evidence already published (ApiLoc).

RESULTS

The rate of misannotated proteins was significantly higher in Mixed groups than in Positive or Negative groups, corroborating the proposed hypothesis. A total of 478 proteins were reannotated and change of signal peptide prediction from negative to positive was the most common. Reannotations triggered the conversion of almost 50% of all Mixed groups, which were further reduced by optimization of signal peptide prediction parameters.

CONCLUSIONS

The methodological novelty proposed here combining orthology and signal peptide prediction proved to be an effective strategy for the identification of proteins showing wrongly N-terminal annotated sequences, and it might have an important impact in the available data for genome-wide searching of potential vaccine and drug targets and proteins involved in host/parasite interactions, as demonstrated for five Plasmodium species.

Systems immunology of human malaria

Plasmodium falciparum malaria remains a global public health threat. Optimism that a highly effective malaria vaccine can be developed stems in part from the observation that humans can acquire immunity to malaria through experimental and natural P. falciparum infection. Recent advances in systems immunology could accelerate efforts to unravel the mechanisms of acquired immunity to malaria. Here, we review the tools of systems immunology, their current limitations in the context of human malaria research, and the human 'models' of malaria immunity to which these tools can be applied.