Search

Abstract Malaria genomic surveillance often estimates parasite genetic relatedness using metrics such as Identity-By-Decent (IBD), yet strong positive selection stemming from antimalarial drug resistance or other interventions may bias IBD-based estimates. In this study, we use simulations, a true IBD inference algorithm, and empirical data sets from different malaria transmission settings to investigate the extent of this bias and explore potential correction strategies. We analyze whole genome sequence data generated from 640 new and 3089 publicly available Plasmodium falciparum clinical isolates. We demonstrate that positive selection distorts IBD distributions, leading to underestimated effective population size and blurred population structure. Additionally, we discover that the removal of IBD peak regions partially restores the accuracy of IBD-based inferences, with this effect contingent on the population’s background genetic relatedness and extent of inbreeding. Consequently, we advocate for selection correction for parasite populations undergoing strong, recent positive selection, particularly in high malaria transmission settings.

Abstract Background Early phase malaria vaccine field trials typically measure malaria infection by PCR or thick blood smear microscopy performed on serially sampled blood. Vaccine efficacy (VE) is the proportion reduction in an endpoint due to vaccination and is often calculated as VEHR = 1–hazard ratio or VERR = 1–risk ratio. Genotyping information can distinguish different clones and distinguish multiple infections over time, potentially increasing statistical power. This paper investigates two alternative VE endpoints incorporating genotyping information: VEmolFOI, the vaccine-induced proportion reduction in incidence of new clones acquired over time, and VEC, the vaccine-induced proportion reduction in mean number of infecting clones per exposure. Methods Power of VEmolFOI and VEC was compared to that of VEHR and VERR by simulations and analytic derivations, and the four VE methods were applied to three data sets: a Phase 3 trial of RTS,S malaria vaccine in 6912 African infants, a Phase 2 trial of PfSPZ Vaccine in 80 Burkina Faso adults, and a trial comparing Plasmodium vivax incidence in 466 Papua New Guinean children after receiving chloroquine + artemether lumefantrine with or without primaquine (as these VE methods can also quantify effects of other prevention measures). By destroying hibernating liver-stage P. vivax, primaquine reduces subsequent reactivations after treatment completion. Results In the trial of RTS,S vaccine, a significantly reduced number of clones at first infection was observed, but this was not the case in trials of PfSPZ Vaccine or primaquine, although the PfSPZ trial lacked power to show a reduction. Resampling smaller data sets from the large RTS,S trial to simulate phase 2 trials showed modest power gains from VEC compared to VEHR for data like those from RTS,S, but VEC is less powerful than VEHR for trials in which the number of clones at first infection is not reduced. VEmolFOI was most powerful in model-based simulations, but only the primaquine trial collected enough serial samples to precisely estimate VEmolFOI. The primaquine VEmolFOI estimate decreased after most control arm liver-stage infections reactivated (which mathematically resembles a waning vaccine), preventing VEmolFOI from improving power. Conclusions The power gain from the genotyping methods depends on the context. Because input parameters for early phase power calculations are often uncertain, these estimators are not recommended as primary endpoints for small trials unless supported by targeted data analysis. Trial registrations: NCT00866619, NCT02663700, NCT02143934.

Introduction Malaria, a devastating febrile illness caused by protozoan parasites, sickened 247,000,000 people in 2021 and killed 619,000, mostly children and pregnant women in sub-Saharan Africa. A highly effective vaccine is urgently needed, especially for Plasmodium falciparum (Pf), the deadliest human malaria parasite. Areas covered Sporozoites (SPZ), the parasite stage transmitted by Anopheles mosquitoes to humans, are the only vaccine immunogen achieving > 90% efficacy against Pf infection. This review describes > 30 clinical trials of PfSPZ vaccines in the U.S.A., Europe, Africa, and Asia, based on first-hand knowledge of the trials and PubMed searches of ‘sporozoites,’ ‘malaria,’ and ‘vaccines.’ Expert opinion First generation (radiation-attenuated) PfSPZ vaccines are safe, well tolerated, 80-100% efficacious against homologous controlled human malaria infection (CHMI) and provide 18-19 months protection without boosting in Africa. Second generation chemo-attenuated PfSPZ are more potent, 100% efficacious against stringent heterologous (variant strain) CHMI, but require a co-administered drug, raising safety concerns. Third generation, late liver stage-arresting, replication competent (LARC), genetically-attenuated PfSPZ are expected to be both safe and highly efficacious. Overall, PfSPZ vaccines meet safety, tolerability, and efficacy requirements for protecting pregnant women and travelers, with licensure for these populations possible within five years. Protecting children and mass vaccination programs to block transmission and eliminate malaria are long-term objectives.

Abstract Over the past two decades, a considerable expansion of malaria interventions has occurred at the national level in Angola, together with cross-border initiatives and regional efforts in southern Africa. Currently, Angola aims to consolidate malaria control and to accelerate the transition from control to pre-elimination, along with other country members of the Elimination 8 initiative. However, the tremendous heterogeneity in malaria prevalence among Angolan provinces, as well as internal population movements and migration across borders, represent major challenges for the Angolan National Malaria Control Programme. This review aims to contribute to the understanding of factors underlying the complex malaria situation in Angola and to encourage future research studies on transmission dynamics and population structure of Plasmodium falciparum, important areas to complement host epidemiological information and to help reenergize the goal of malaria elimination in the country.

Abstract Two malaria parasite species, Plasmodium falciparum (Pf) and P. vivax (Pv) are responsible for most of the disease burden caused by malaria. Vaccine development against this disease has focused mainly on Pf. Whole-sporozoite (WSp) vaccination, targeting pre-erythrocytic (PE) parasite stages, is a promising strategy for immunization against malaria and several PfWSp-based vaccine candidates are currently undergoing clinical evaluation. In contrast, no WSp candidates have been developed for Pv, mainly due to constraints in the production of Pv sporozoites in the laboratory. Recently, we developed a novel approach for WSp vaccination against Pf based on the use of transgenic rodent P. berghei (Pb) sporozoites expressing immunogens of this human-infective parasite. We showed that this platform can be used to deliver PE Pf antigens, eliciting both targeted humoral responses and cross-species cellular immune responses against Pf. Here we explored this WSp platform for the delivery of Pv antigens. As the Pv circumsporozoite protein (CSP) is a leading vaccine candidate antigen, we generated a transgenic Pb parasite, PbviVac, that, in addition to its endogenous PbCSP, expresses PvCSP under the control of a strictly PE promoter. Immunofluorescence microscopy analyses confirmed that both the PbCSP and the PvCSP antigens are expressed in PbviVac sporozoites and liver stages and that PbviVac sporozoite infectivity of hepatic cells is similar to that of its wild-type Pb counterpart. Immunization of mice with PbviVac sporozoites elicits the production of anti-PvCSP antibodies that efficiently recognize and bind to Pv sporozoites. Our results warrant further development and evaluation of PbviVac as a surrogate for WSp vaccination against Pv malaria.

Abstract Background The ability of malaria rapid diagnostic tests (RDTs) to effectively detect active infections is being compromised by the presence of malaria strains with genomic deletions at the hrp2 and hrp3 loci, encoding the antigens most commonly targeted in diagnostics for Plasmodium falciparum detection. The presence of such deletions can be determined in publically available P. falciparum whole genome sequencing (WGS) datasets. A computational approach was developed and validated, termed Gene Coverage Count and Classification (GC3), to analyse genome-wide sequence coverage data and provide informative outputs to assess presence and coverage profile of a target locus in WGS data. GC3 was applied to detect deletions at hrp2 and hrp3 (hrp2/3) and flanking genes in different geographic regions and across time points. Methods GC3 uses Python and R scripts to extract locus read coverage metrics from mapped WGS data according to user-defined parameters and generates relevant tables and figures. GC3 was tested using WGS data for laboratory reference strains with known hrp2/3 genotypes, and its results compared to those of a hrp2/3-specific qPCR assay. Samples with at least 25% of coding region positions with zero coverage were classified as having a deletion. Publicly available sequence data was analysed and compared with published deletion frequency estimates. Results GC3 results matched the expected coverage of known laboratory reference strains. Agreement between GC3 and a hrp2/3-specific qPCR assay reported for 19/19 (100%) hrp2 deletions and 18/19 (94.7%) hrp3 deletions. Among Cambodian (n = 127) and Brazilian (n = 20) WGS datasets, which had not been previously analysed for hrp2/3 deletions, GC3 identified hrp2 deletions in three and four samples, and hrp3 deletions in 10 and 15 samples, respectively. Plots of hrp2/3 coding regions, grouped by year of sample collection, showed a decrease in median standardized coverage among Malawian samples (n = 150) suggesting the importance of a careful, properly controlled follow up to determine if an increase in frequency of deletions has occurred between 2007–2008 and 2014–2015. Among Malian (n = 90) samples, median standardized coverage was lower in 2002 than 2010, indicating widespread deletions present at the gene locus in 2002. Conclusions The GC3 tool accurately classified hrp2/3 deletions and provided informative tables and figures to analyse targeted gene coverage. GC3 is an appropriate tool when performing preliminary and exploratory assessment of locus coverage data.

IntroductionHost gene and protein expression impact susceptibility to clinical malaria, but the balance of immune cell populations, cytokines and genes that contributes to protection, remains incompletely understood. Little is known about the determinants of host susceptibility to clinical malaria at a time when acquired immunity is developing.MethodsWe analyzed peripheral blood mononuclear cells (PBMCs) collected from children who differed in susceptibility to clinical malaria, all from a small town in Mali. PBMCs were collected from children aged 4-6 years at the start, peak and end of the malaria season. We characterized the immune cell composition and cytokine secretion for a subset of 20 children per timepoint (10 children with no symptomatic malaria age-matched to 10 children with >2 symptomatic malarial illnesses), and gene expression patterns for six children (three per cohort) per timepoint. ResultsWe observed differences between the two groups of children in the expression of genes related to cell death and inflammation; in particular, inflammatory genes such as CXCL10 and STAT1 and apoptotic genes such as XAF1 were upregulated in susceptible children before the transmission season began. We also noted higher frequency of HLA-DR+ CD4 T cells in protected children during the peak of the malaria season and comparable levels cytokine secretion after stimulation with malaria schizonts across all three time points. ConclusionThis study highlights the importance of baseline immune signatures in determining disease outcome. Our data suggests that differences in apoptotic and inflammatory gene expression patterns can serve as predictive markers of susceptibility to clinical malaria.

Here the authors show that controlled human malaria infection with a Brazilian parasite highly divergent from vaccine and West African field strains can provide estimates of vaccine efficacy in Mali, and could replace field testing, streamlining vaccine development.

Abstract Background RIFINs and STEVORs are variant surface antigens expressed by P. falciparum that play roles in severe malaria pathogenesis and immune evasion. These two highly diverse multigene families feature multiple paralogs, making their classification challenging using traditional bioinformatic methods. Results STRIDE (STevor and RIfin iDEntifier) is an HMM-based, command-line program that automates the identification and classification of RIFIN and STEVOR protein sequences in the malaria parasite Plasmodium falciparum. STRIDE is more sensitive in detecting RIFINs and STEVORs than available PFAM and TIGRFAM tools and reports RIFIN subtypes and the number of sequences with a FHEYDER amino acid motif, which has been associated with severe malaria pathogenesis. Conclusions STRIDE will be beneficial to malaria research groups analyzing genome sequences and transcripts of clinical field isolates, providing insight into parasite biology and virulence.

OBJECTIVES/GOALS: Recent in vitro evidence suggests that diverse parasite protein families called RIFINs and STEVORs are displayed on the surface of infected red blood cells and may have a role in severe malaria, but they remain sparsely studied in natural infections. We measured the RNA expression of these antigens in Malian children with severe or mild malaria illness. METHODS/STUDY POPULATION: We collected blood samples from Malian children aged six months to five years, including 14 with cerebral malaria, 10 with severe malarial anemia, and demographic-matched controls with mild, uncomplicated malaria. We extracted total RNA from each patient and used a custom capture array to selectively enrich Plasmodium falciparum parasite RNA. We then performed Illumina next-generation RNA sequencing and reconstructed parasite transcriptomes using reference-free de novo assembly. We identified RNA encoding RIFINs and STEVORs using an in-house classifier, then measured the diversity and abundance of gene expression for each infection. Expression diversity was defined as the number of unique variants transcribed. Expression abundance was calculated as transcripts per million (TPM). RESULTS/ANTICIPATED RESULTS: Cerebral malaria cases, but not severe malarial anemia cases, had higher diversity and abundance of RIFIN expression compared to mild infections. Type A RIFINs predominated over Type B RIFINs, and the same two RIFINs were predominantly expressed in all disease phenotypes. We anticipate that predominantly expressed RIFINs share high sequence homology with variants previously shown to bind blood antigens or immune inhibitory receptors. STEVOR expression was also higher in cerebral malaria compared to mild malaria, but STEVOR transcripts were sparse overall. DISCUSSION/SIGNIFICANCE: Elevated RIFIN expression in cerebral malaria over mild malaria supports a role for these antigens in pathogenesis. Severe malarial anemia may progress through a different pathogenic mechanism. Predominantly expressed RIFIN variants may be promising targets for vaccines and therapeutics to protect children against cerebral malaria.

Molina-Cruz et al. analyzed 4,971 sequences of the Plasmodium falciparum Pfs47 gene and identified a combination of two SNP haplotypes highly correlated with geographic origin and transmission potential by mosquito vectors. A simple High-Resolution Melting RT-PCR assay was developed to genotype these SNPs and establish the likely continent of origin and transmission risk of P. falciparum isolates.

Resistance to antimalarial drugs, and in particular to the artemisinin derivatives and their partner drugs, threatens recent progress toward regional malaria elimination and eventual global malaria eradication. Population-level studies utilizing whole-genome sequencing approaches have facilitated the identification of regions of the parasite genome associated with both clinical and in vitro drug-resistance phenotypes. However, the biological relevance of genes identified in these analyses and the establishment of a causal relationship between genotype and phenotype requires functional characterization. Here we examined data from population genomic and transcriptomic studies in the context of data generated from recent functional studies, using a new population genetic approach designed to identify potential favored mutations within the region of a selective sweep (iSAFE). We identified several genes functioning in pathways now known to be associated with artemisinin resistance that were supported in early population genomic studies, as well as potential new drug targets/pathways for further validation and consideration for treatment of artemisinin-resistant Plasmodium falciparum. In addition, we establish the utility of iSAFE in identifying positively-selected mutations in population genomic studies, potentially accelerating the time to functional validation of candidate genes.

Abstract Library preparation for high-throughput sequencing applications is a critical step in producing representative, unbiased sequencing data. The iGenomX Riptide High Throughput Rapid Library Prep Kit purports to provide high-quality sequencing data with lower costs compared to other Illumina library kits. To test these claims, we compared sequence data quality of Riptide libraries to libraries constructed with KAPA Hyper and NEBNext Ultra. Across several single-source genome samples, mapping performance and de novo assembly of Riptide libraries were similar to conventional libraries prepared with the same DNA. Poor performance of some libraries resulted in low sequencing depth. In particular, degraded DNA samples may be challenging to sequence with Riptide. There was little cross-well plate contamination with the overwhelming majority of reads belong to the proper source genomes. The sequencing of metagenome samples using different Riptide primer sets resulted in variable taxonomic assignment of reads. Increased adoption of the Riptide kit will decrease library preparation costs. However, this method might not be suitable for degraded DNA.

Abstract Background The apicomplexan parasite Theileria parva causes a livestock disease called East coast fever (ECF), with millions of animals at risk in sub-Saharan East and Southern Africa, the geographic distribution of T. parva. Over a million bovines die each year of ECF, with a tremendous economic burden to pastoralists in endemic countries. Comprehensive, accurate parasite genome annotation can facilitate the discovery of novel chemotherapeutic targets for disease treatment, as well as elucidate the biology of the parasite. However, genome annotation remains a significant challenge because of limitations in the quality and quantity of the data being used to inform the location and function of protein-coding genes and, when RNA data are used, the underlying biological complexity of the processes involved in gene expression. Here, we apply our recently published RNAseq dataset derived from the schizont life-cycle stage of T. parva to update structural and functional gene annotations across the entire nuclear genome. Results The re-annotation effort lead to evidence-supported updates in over half of all protein-coding sequence (CDS) predictions, including exon changes, gene merges and gene splitting, an increase in average CDS length of approximately 50 base pairs, and the identification of 128 new genes. Among the new genes identified were those involved in N-glycosylation, a process previously thought not to exist in this organism and a potentially new chemotherapeutic target pathway for treating ECF. Alternatively-spliced genes were identified, and antisense and multi-gene family transcription were extensively characterized. Conclusions The process of re-annotation led to novel insights into the organization and expression profiles of protein-coding sequences in this parasite, and uncovered a minimal N-glycosylation pathway that changes our current understanding of the evolution of this post-translational modification in apicomplexan parasites.

Abstract Background Owing to the large amount of host DNA in clinical samples, generation of high-quality Plasmodium falciparum whole genome sequencing (WGS) data requires enrichment for parasite DNA. Enrichment is often achieved by leukocyte depletion of infected blood prior to storage. However, leukocyte depletion is difficult in low-resource settings and limits analysis to prospectively-collected samples. As a result, approaches such as selective whole genome amplification (sWGA) are being used to enrich for parasite DNA. However, sWGA has had limited success in generating reliable sequencing data from low parasitaemia samples. In this study, enzymatic digestion with MspJI prior to sWGA and whole genome sequencing was evaluated to determine whether this approach improved genome coverage compared to sWGA alone. The potential of sWGA to cause amplification bias in polyclonal infections was also examined. Methods DNA extracted from laboratory-created dried blood spots was treated with a modification-dependent restriction endonuclease, MspJI, and filtered via vacuum filtration. Samples were then selectively amplified using a previously reported sWGA protocol and subjected to WGS. Genome coverage statistics were compared between the optimized sWGA approach and the previously reported sWGA approach performed in parallel. Differential amplification by sWGA was assessed by comparing WGS data generated from lab-created mixtures of parasite isolates, from the same geographical region, generated with or without sWGA. Results MspJI digestion did not enrich for parasite DNA. Samples that underwent vacuum filtration (without MspJI digestion) prior to sWGA had the highest parasite DNA concentration and displayed greater genome coverage compared to MspJI + sWGA and sWGA alone, particularly for low parasitaemia samples. The optimized sWGA (filtration + sWGA) approach was successfully used to generate WGS data from 218 non-leukocyte depleted field samples from Malawi. Sequences from lab-created mixtures of parasites did not show evidence of differential amplification of parasite strains compared to directly sequenced samples. Conclusion This optimized sWGA approach is a reliable method to obtain WGS data from non-leukocyte depleted, low parasitaemia samples. The absence of amplification bias in data generated from mixtures of isolates from the same geographic region suggests that this approach can be appropriately used for molecular epidemiological studies.

Abstract Background Plasmodium falciparum (Pf) whole-organism sporozoite vaccines have been shown to provide significant protection against controlled human malaria infection (CHMI) in clinical trials. Initial CHMI studies showed significantly higher durable protection against homologous than heterologous strains, suggesting the presence of strain-specific vaccine-induced protection. However, interpretation of these results and understanding of their relevance to vaccine efficacy have been hampered by the lack of knowledge on genetic differences between vaccine and CHMI strains, and how these strains are related to parasites in malaria endemic regions. Methods Whole genome sequencing using long-read (Pacific Biosciences) and short-read (Illumina) sequencing platforms was conducted to generate de novo genome assemblies for the vaccine strain, NF54, and for strains used in heterologous CHMI (7G8 from Brazil, NF166.C8 from Guinea, and NF135.C10 from Cambodia). The assemblies were used to characterize sequences in each strain relative to the reference 3D7 (a clone of NF54) genome. Strains were compared to each other and to a collection of clinical isolates (sequenced as part of this study or from public repositories) from South America, sub-Saharan Africa, and Southeast Asia. Results While few variants were detected between 3D7 and NF54, we identified tens of thousands of variants between NF54 and the three heterologous strains. These variants include SNPs, indels, and small structural variants that fall in regulatory and immunologically important regions, including transcription factors (such as PfAP2-L and PfAP2-G) and pre-erythrocytic antigens that may be key for sporozoite vaccine-induced protection. Additionally, these variants directly contributed to diversity in immunologically important regions of the genomes as detected through in silico CD8+ T cell epitope predictions. Of all heterologous strains, NF135.C10 had the highest number of unique predicted epitope sequences when compared to NF54. Comparison to global clinical isolates revealed that these four strains are representative of their geographic origin despite long-term culture adaptation; of note, NF135.C10 is from an admixed population, and not part of recently formed subpopulations resistant to artemisinin-based therapies present in the Greater Mekong Sub-region. Conclusions These results will assist in the interpretation of vaccine efficacy of whole-organism vaccines against homologous and heterologous CHMI.

ABSTRACT var genes encode Plasmodium falciparum erythrocyte membrane protein-1 (PfEMP1) antigens. These highly diverse antigens are displayed on the surface of infected erythrocytes and play a critical role in immune evasion and sequestration of infected erythrocytes. Studies of var expression using non-leukocyte-depleted blood are challenging because of the predominance of host genetic material and lack of conserved var segments. Our goal was to enrich for parasite RNA, allowing de novo assembly of var genes and detection of expressed novel variants. We used two overall approaches: (i) enriching for total mRNA in the sequencing library preparations and (ii) enriching for parasite RNA with a custom capture array based on Roche’s SeqCap EZ enrichment system. The capture array was designed with probes based on the whole 3D7 reference genome and an additional >4,000 full-length var gene sequences from other P. falciparum strains. We tested each method on RNA samples from Malian children with severe or uncomplicated malaria infections. All reads mapping to the human genome were removed, the remaining reads were assembled de novo into transcripts, and from these, var-like transcripts were identified and annotated. The capture array produced the longest maximum length and largest numbers of var gene transcripts in each sample, particularly in samples with low parasitemia. Identifying the most-expressed var gene sequences in whole-blood clinical samples without the need for extensive processing or generating sample-specific reference genome data is critical for understanding the role of PfEMP1s in malaria pathogenesis. IMPORTANCE Malaria parasites display antigens on the surface of infected red blood cells in the human host that facilitate attachment to blood vessels, contributing to the severity of infection. These antigens are highly variable, allowing the parasite to evade the immune system. Identifying these expressed antigens is critical to understanding the development of severe malarial disease. However, clinical samples contain limited amounts of parasite genetic material, a challenge for sequencing efforts further compounded by the extreme diversity of the parasite surface antigens. We present a method that enriches for these antigen sequences in clinical samples using a custom capture array, requiring minimal processing in the field. While our results are focused on the malaria parasite Plasmodium falciparum, this approach has broad applicability to other highly diverse antigens from other parasites and pathogens such as those that cause giardiasis and leishmaniasis.

Understanding genomic variation in Plasmodium falciparum parasites and inferring migration patterns can guide malaria elimination strategies. Using genome-wide data for 1722 parasites collected from 54 districts, the authors use identity-by-descent approaches to estimate regional parasite migration and spread of artemisinin drug resistance.

Failure to account for genetic diversity of antigens during vaccine design may lead to vaccine escape. To evaluate the vaccine escape potential of antigens used in vaccines currently in development or clinical testing, we surveyed the genetic diversity, measured population differentiation, and performed in silico prediction and analysis of T-cell epitopes of ten such Plasmodium falciparum pre-erythrocytic-stage antigens using whole-genome sequence data from 1010 field isolates. Of these, 699 were collected in Africa (Burkina Faso, Cameroon, Guinea, Kenya, Malawi, Mali, and Tanzania), 69 in South America (Brazil, Colombia, French Guiana, and Peru), 59 in Oceania (Papua New Guinea), and 183 in Asia (Cambodia, Myanmar, and Thailand). Antigens surveyed include cell-traversal protein for ookinetes and sporozoites, circumsporozoite protein, liver-stage antigens 1 and 3, sporozoite surface proteins P36 and P52, sporozoite asparagine-rich protein-1, sporozoite microneme protein essential for cell traversal-2, and upregulated-in-infectious-sporozoite 3 and 4 proteins. The analyses showed that a limited number of these protein variants, when combined, would be representative of worldwide parasite populations. Moreover, predicted T-cell epitopes were identified that could be further explored for immunogenicity and protective efficacy. Findings can inform the rational design of a multivalent malaria vaccine.

Babesia bovis is a hemoprotozoan parasite of cattle that has a complex life cycle within vertebrate and invertebrate hosts. In the mammalian host, B. bovis undergoes asexual reproduction while in the tick midgut, gametes are induced, fuse, and form zygotes. The zygote infects tick gut epithelial cells and transform into kinetes that are released into the hemolymph and invade other tick tissues such as the ovaries, resulting in transovarial transmission to tick offspring. To compare gene regulation between different B. bovis life stages, we collected parasites infecting bovine erythrocytes and tick hemolymph. Total RNA samples were isolated, and multiplexed libraries sequenced using paired-end 100 cycle reads of a HiSeq 2500. The data was normalized using the TMM method and analysed for significant differential expression using the generalized linear model likelihood ratio test (GLM LRT) in edgeR. To validate our datasets, ten genes were selected using NormFinder. Genes that had no significant fold change between the blood and tick stages in the RNA-Seq datasets were tested by quantitative PCR to determine their suitability as “housekeeping” genes. The normalized RNA-Seq data revealed genes upregulated during infection of the mammalian host or tick vector and six upregulated genes were validated by quantitative PCR. These datasets can help identify useful targets for controlling bovine babesiosis.

Theileria parva is a protozoan parasite that causes East Coast fever (ECF), an economically important disease of cattle in Africa. It is transmitted mainly by the tick Rhipicephalus appendiculatus. Research efforts to develop a subunit vaccine based on parasite neutralizing antibodies and cytotoxic T-lymphocytes have met with limited success. The molecular mechanisms underlying T. parva life cycle stages in the tick vector and bovine host are poorly understood, thus limiting progress toward an effective and efficient control of ECF. Transcriptomics has been used to identify candidate vaccine antigens or markers associated with virulence and disease pathology. Therefore, characterization of gene expression throughout the parasite's life cycle should shed light on host–pathogen interactions in ECF and identify genes underlying differences in parasite stages as well as potential, novel therapeutic targets. Recently, the first gene expression profiling of T. parva was conducted for the sporoblast, sporozoite, and schizont stages. The sporozoite is infective to cattle, whereas the schizont is the major pathogenic form of the parasite. The schizont can differentiate into piroplasm, which is infective to the tick vector. The present study was designed to extend the T. parva gene expression profiling to the piroplasm stage with reference to the schizont. Pairwise comparison revealed that 3,279 of a possible 4,084 protein coding genes were differentially expressed, with 1,623 (49%) genes upregulated and 1,656 (51%) downregulated in the piroplasm relative to the schizont. In addition, over 200 genes were stage-specific. In general, there were more molecular functions, biological processes, subcellular localizations, and pathways significantly enriched in the piroplasm than in the schizont. Using known antigens as benchmarks, we identified several new potential vaccine antigens, including TP04_0076 and TP04_0640, which were highly immunogenic in naturally T. parva-infected cattle. All the candidate vaccine antigens identified have yet to be investigated for their capacity to induce protective immune response against ECF.

Abstract Background VAR2CSA, a member of the Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1) family, mediates the binding of P. falciparum-infected erythrocytes to chondroitin sulfate A, a surface-associated molecule expressed in placental cells, and plays a central role in the pathogenesis of placental malaria. VAR2CSA is a target of naturally acquired immunity and, as such, is a leading vaccine candidate against placental malaria. This protein is very polymorphic and technically challenging to sequence. Published var2csa sequences, mostly limited to specific domains, have been generated through the sequencing of cloned PCR amplicons using capillary electrophoresis, a method that is both time consuming and costly, and that performs poorly when applied to clinical samples that are commonly polyclonal. A next-generation sequencing platform, Pacific Biosciences (PacBio), offers an alternative approach to overcome these issues. Methods PCR primers were designed that target a 5 kb segment in the 5′ end of var2csa and the resulting amplicons were sequenced using PacBio sequencing. The primers were optimized using two laboratory strains and were validated on DNA from 43 clinical samples, extracted from dried blood spots on filter paper or from cryopreserved P. falciparum-infected erythrocytes. Sequence reads were assembled using the SMRT-analysis ConsensusTools module. Results Here, a PacBio sequencing-based approach for recovering a segment encoding the majority of VAR2CSA’s extracellular region is described; this segment includes the totality of the first four domains in the 5′ end of var2csa (~5 kb), from clinical malaria samples. The feasibility of the method is demonstrated, showing a high success rate from cryopreserved samples and more limited success from dried blood spots stored at room temperature, and characterized the genetic variation of the var2csa locus. Conclusions This method will facilitate a detailed analysis of var2csa genetic variation and can be adapted to sequence other hypervariable P. falciparum genes.

Abstract Background Encoded by the var gene family, highly variable Plasmodium falciparum erythrocyte membrane protein-1 (PfEMP1) proteins mediate tissue-specific cytoadherence of infected erythrocytes, resulting in immune evasion and severe malaria disease. Sequencing and assembling the 40–60 var gene complement for individual infections has been notoriously difficult, impeding molecular epidemiological studies and the assessment of particular var elements as subunit vaccine candidates. Methods We developed and validated a novel algorithm, Exon-Targeted Hybrid Assembly (ETHA), to perform targeted assembly of var gene sequences, based on a combination of Pacific Biosciences and Illumina data. Results Using ETHA, we characterized the repertoire of var genes in 12 samples from uncomplicated malaria infections in children from a single Malian village and showed them to be as genetically diverse as vars from isolates from around the globe. The gene var2csa, a member of the var family associated with placental malaria pathogenesis, was present in each genome, as were vars previously associated with severe malaria. Conclusion ETHA, a tool to discover novel var sequences from clinical samples, will aid the understanding of malaria pathogenesis and inform the design of malaria vaccines based on PfEMP1. ETHA is available at: https://sourceforge.net/projects/etha/ .