Whole-Genome Sequencing to Evaluate the Resistance Landscape Following Antimalarial Treatment Failure With Fosmidomycin-Clindamycin

Turning the needle into the haystack: Culture-independent amplification of complex microbial genomes directly from their native environment

High-throughput sequencing (HTS) has revolutionized microbiology, but many microbes exist at low abundance in their natural environment and/or are difficult, if not impossible, to culture in the laboratory. This makes it challenging to use HTS to study the genomes of many important microbes and pathogens. In this review, we discuss the development and application of selective whole genome amplification (SWGA) to allow whole or partial genomes to be sequenced for low abundance microbes directly from complex biological samples. We highlight ways in which genomic data generated by SWGA have been used to elucidate the population dynamics of important human pathogens and monitor development of antimicrobial resistance and the emergence of potential outbreaks. We also describe the limitations of this method and propose some potential innovations that could be used to improve the quality of SWGA and lower the barriers to using this method across a wider range of infectious pathogens.

Beyond the MEP Pathway: A novel kinase required for prenol utilization by malaria parasites

A proposed treatment for malaria is a combination of fosmidomycin and clindamycin. Both compounds inhibit the methylerythritol 4-phosphate (MEP) pathway, the parasitic source of farnesyl and geranylgeranyl pyrophosphate (FPP and GGPP, respectively). Both FPP and GGPP are crucial for the biosynthesis of several essential metabolites such as ubiquinone and dolichol, as well as for protein prenylation. Dietary prenols, such as farnesol (FOH) and geranylgeraniol (GGOH), can rescue parasites from MEP inhibitors, suggesting the existence of a missing pathway for prenol salvage via phosphorylation. In this study, we identified a gene in the genome of P. falciparum, encoding a transmembrane prenol kinase (PolK) involved in the salvage of FOH and GGOH. The enzyme was expressed in Saccharomyces cerevisiae, and its FOH/GGOH kinase activities were experimentally validated. Furthermore, conditional knockout parasites (Δ-PolK) were created to investigate the biological importance of the FOH/GGOH salvage pathway. Δ-PolK parasites were viable but displayed increased susceptibility to fosmidomycin. Their sensitivity to MEP inhibitors could not be rescued by adding prenols. Additionally, Δ-PolK parasites lost their capability to utilize prenols for protein prenylation. Experiments using culture medium supplemented with whole/delipidated human plasma in transgenic parasites revealed that human plasma has components that can diminish the effectiveness of fosmidomycin. Mass spectrometry tests indicated that both bovine supplements used in culture and human plasma contain GGOH. These findings suggest that the FOH/GGOH salvage pathway might offer an alternate source of isoprenoids for malaria parasites when de novo biosynthesis is inhibited. This study also identifies a novel kind of enzyme related to isoprenoid metabolism.

A fast machine-learning-guided primer design pipeline for selective whole genome amplification

Addressing many of the major outstanding questions in the fields of microbial evolution and pathogenesis will require analyses of populations of microbial genomes. Although population genomic studies provide the analytical resolution to investigate evolutionary and mechanistic processes at fine spatial and temporal scales-precisely the scales at which these processes occur-microbial population genomic research is currently hindered by the practicalities of obtaining sufficient quantities of the relatively pure microbial genomic DNA necessary for next-generation sequencing. Here we present swga2.0, an optimized and parallelized pipeline to design selective whole genome amplification (SWGA) primer sets. Unlike previous methods, swga2.0 incorporates active and machine learning methods to evaluate the amplification efficacy of individual primers and primer sets. Additionally, swga2.0 optimizes primer set search and evaluation strategies, including parallelization at each stage of the pipeline, to dramatically decrease program runtime. Here we describe the swga2.0 pipeline, including the empirical data used to identify primer and primer set characteristics, that improve amplification performance. Additionally, we evaluate the novel swga2.0 pipeline by designing primer sets that successfully amplify Prevotella melaninogenica, an important component of the lung microbiome in cystic fibrosis patients, from samples dominated by human DNA.

Selective whole-genome amplification reveals population genetics of Leishmania braziliensis directly from patient skin biopsies

In Brazil, Leishmania braziliensis is the main causative agent of the neglected tropical disease, cutaneous leishmaniasis (CL). CL presents on a spectrum of disease severity with a high rate of treatment failure. Yet the parasite factors that contribute to disease presentation and treatment outcome are not well understood, in part because successfully isolating and culturing parasites from patient lesions remains a major technical challenge. Here we describe the development of selective whole genome amplification (SWGA) for Leishmania and show that this method enables culture-independent analysis of parasite genomes obtained directly from primary patient skin samples, allowing us to circumvent artifacts associated with adaptation to culture. We show that SWGA can be applied to multiple Leishmania species residing in different host species, suggesting that this method is broadly useful in both experimental infection models and clinical studies. SWGA carried out directly on skin biopsies collected from patients in Corte de Pedra, Bahia, Brazil, showed extensive genomic diversity. Finally, as a proof-of-concept, we demonstrated that SWGA data can be integrated with published whole genome data from cultured parasite isolates to identify variants unique to specific geographic regions in Brazil where treatment failure rates are known to be high. SWGA provides a relatively simple method to generate Leishmania genomes directly from patient samples, unlocking the potential to link parasite genetics with host clinical phenotypes.

Isoprenoid alcohols utilization by malaria parasites

Plasmodium falciparum is the etiological agent of human malaria, one of the most widespread diseases in tropical and subtropical regions. Drug resistance is one of the biggest problems in controlling the disease, which leads to the need to discover new antimalarial compounds. One of the most promissory drugs purposed is fosmidomycin, an inhibitor of the biosynthesis of isoprene units by the methylerythritol 4-phosphate (MEP) pathway, which in some cases failed in clinical studies. Once formed, isoprene units are condensed to form longer structures such as farnesyl and geranylgeranyl pyrophosphate, which are necessary for Heme O and A formation, ubiquinone, and dolichyl phosphate biosynthesis as well as for protein isoprenylation. Even though the natural substrates of polyprenyl transferases and synthases are polyprenyl pyrophosphates, it was already demonstrated that isoprenoid alcohols (polyprenols) such as farnesol (FOH) and geranylgeraniol (GGOH) can rescue parasites from fosmidomycin. This study better investigated how this rescue phenomenon occurs by performing drug-rescue assays. Similarly, to FOH and GGOH, it was observed that phytol (POH), a 20-carbon plant isoprenoid, as well as unsaponifiable lipid extracts from foods rescue parasites from the antimalarial effect of fosmidomycin. Contrarily, neither dolichols nor nonaprenol rescue parasites from fosmidomycin. Considering this, here we characterized the transport of FOH, GGOH, and POH. Once incorporated, it was observed that these substances are phosphorylated, condensed into longer isoprenoid alcohols, and incorporated into proteins and dolichyl phosphates. Through proteomic and radiolabelling approaches, it was found that prenylated proteins are naturally attached to several isoprenoids, derived from GGOH, dolichol, and POH if exogenously added. Furthermore, the results suggest the presence of at least two promiscuous protein prenyltransferases in the parasite: one enzyme which can use FPP among other unidentified substrates and another enzyme that can use GGPP, phytyl pyrophosphate (PPP), and dolichols, among other substrates not identified here. Thus, further evidence was obtained for dolichols and other isoprenoid products attached to proteins. This study helps to better understand the apicoplast-targeting antimalarial mechanism of action and a novel post-translational modification of proteins in P. falciparum.

5WBF: a low-cost and straightforward whole blood filtration method suitable for whole-genome sequencing of Plasmodium falciparum clinical isolates

Background

Whole-genome sequencing (WGS) is becoming increasingly helpful to assist malaria control programmes. A major drawback of this approach is the large amount of human DNA compared to parasite DNA extracted from unprocessed whole blood. As red blood cells (RBCs) have a diameter of about 7–8 µm and exhibit some deformability, it was hypothesized that cheap and commercially available 5 µm filters might retain leukocytes but much less of Plasmodium falciparum-infected RBCs. This study aimed to test the hypothesis that such a filtration method, named 5WBF (for 5 µm Whole Blood Filtration), may provide highly enriched parasite material suitable for P. falciparum WGS.

Methods

Whole blood was collected from five patients experiencing a P. falciparum malaria episode (ring-stage parasitaemia range: 0.04–5.5%) and from mock samples obtained by mixing synchronized, ring-stage cultured P. falciparum 3D7 parasites with uninfected human whole blood (final parasitaemia range: 0.02–1.1%). These whole blood samples (50 to 400 µL) were diluted in RPMI 1640 medium or PBS 1× buffer and filtered with a syringe connected to a 5 µm commercial filter. DNA was extracted from 5WBF-treated and unfiltered counterpart blood samples using a commercial kit. The 5WBF method was evaluated on the ratios of parasite:human DNA assessed by qPCR and by sequencing depth and percentages of coverage from WGS data (Illumina NextSeq 500). As a comparison, the popular selective whole-genome amplification (sWGA) method, which does not rely on blood filtration, was applied to the unfiltered counterpart blood samples.

Results

After applying 5WBF, qPCR indicated an average of twofold loss in the amount of parasite template DNA (Pf ARN18S gene) and from 4096- to 65,536-fold loss of human template DNA (human β actin gene). WGS analyses revealed that > 95% of the  parasite nuclear and organellar genomes were all covered at ≥ 10× depth for all samples tested. In sWGA counterparts, the organellar genomes were poorly covered and from 47.7 to 82.1% of the nuclear genome was covered at ≥ 10× depth depending on parasitaemia. Sequence reads were homogeneously distributed across gene sequences for 5WBF-treated samples (n = 5460 genes; mean coverage: 91×; median coverage: 93×; 5th percentile: 70×; 95th percentile: 103×), allowing the identification of gene copy number variations such as for gch1. This later analysis was not possible for sWGA-treated samples, as a much more heterogeneous distribution of reads across gene sequences was observed (mean coverage: 80×; median coverage: 51×; 5th percentile: 7×; 95th percentile: 245×).

Conclusions

The novel 5WBF leucodepletion method is simple to implement and based on commercially available, standardized 5 µm filters which cost from 1.0 to 1.7€ per unit depending on suppliers. 5WBF permits extensive genome-wide analysis of P. falciparum ring-stage isolates from minute amounts of whole blood even with parasitaemias as low as 0.02%.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-022-04073-1.

Diversity and Within-Host Evolution of Leishmania donovani from Visceral Leishmaniasis Patients with and without HIV Coinfection in Northern Ethiopia

Visceral leishmaniasis (VL) is a fatal disease and a growing public health problem in East Africa, where Ethiopia has one of the highest VL burdens. The largest focus of VL in Ethiopia is driven by high prevalence in migrant agricultural workers and associated with a high rate of coinfection with HIV. This coinfection makes VL more difficult to treat successfully and is associated with a high rate of relapse, with VL/HIV patients frequently experiencing many relapses of VL before succumbing to this infection. We present genome-wide data on Leishmania donovani isolates from a longitudinal study of cohorts of VL and VL/HIV patients reporting to a single clinic in Ethiopia. Extensive clinical data allow us to investigate the influence of coinfection and relapse on the populations of parasites infecting these patients. We find that the same parasite population is responsible for both VL and VL/HIV infections and that, in most cases, disease relapse is caused by recrudescence of the population of parasites that caused primary VL. Complex, multiclonal infections are present in both primary and relapse cases, but the infrapopulation of parasites within a patient loses genetic diversity between primary disease presentation and subsequent relapses, presumably due to a population bottleneck induced by treatment. These data suggest that VL/HIV relapses are not caused by genetically distinct parasite infections or by reinfection. Treatment of VL does not lead to sterile cure, and in VL/HIV, the infecting parasites are able to reestablish after clinically successful treatment, leading to repeated relapse of VL. IMPORTANCE Visceral leishmaniasis (VL) is the second largest cause of deaths due to parasite infections and a growing problem in East Africa. In Ethiopia, it is particularly associated with migrant workers moving from regions of nonendemicity for seasonal agricultural work and is frequently found as a coinfection with HIV, which leads to frequent VL relapse following treatment. Insight into the process of relapse in these patients is thus key to controlling the VL epidemic in Ethiopia. We show that there is little genetic differentiation between the parasites infecting HIV-positive and HIV-negative VL patients. Moreover, we provide evidence that relapses are caused by the initially infecting parasite population and that treatment induces a loss of genetic diversity in this population. We propose that restoring functioning immunity and improving antiparasitic treatment may be key in breaking the cycle of relapsing VL in VL/HIV patients.

Optimization of whole-genome sequencing of Plasmodium falciparum from low-density dried blood spot samples

Background

Whole-genome sequencing (WGS) is becoming increasingly useful to study the biology, epidemiology, and ecology of malaria parasites. Despite ease of sampling, DNA extracted from dried blood spots (DBS) has a high ratio of human DNA compared to parasite DNA, which poses a challenge for downstream genetic analyses. The effects of multiple methods for DNA extraction, digestion of methylated DNA, and amplification were evaluated on the quality and fidelity of WGS data recovered from DBS.

Methods

Low parasite density mock DBS samples were created, extracted either with Tween-Chelex or QIAamp, treated with or without McrBC, and amplified with one of three different amplification techniques (two sWGA primer sets and one rWGA). Extraction conditions were evaluated on performance of sequencing depth, percentiles of coverage, and expected SNP concordance.

Results

At 100 parasites/μL, Chelex-Tween-McrBC samples had higher coverage (5 × depth = 93% genome) than QIAamp extracted samples (5 × depth = 76% genome). The two evaluated sWGA primer sets showed minor differences in overall genome coverage and SNP concordance, with a newly proposed combination of 20 primers showing a modest improvement in coverage over those previously published.

Conclusions

Overall, Tween-Chelex extracted samples that were treated with McrBC digestion and are amplified using 6A10AD sWGA conditions had minimal dropout rate, higher percentages of coverage at higher depth, and more accurate SNP concordance than QiaAMP extracted samples. These findings extend the results of previously reported methods, making whole genome sequencing accessible to a larger number of low density samples that are commonly encountered in cross-sectional surveys.

Proteome size reduction in Apicomplexans is linked with loss of DNA repair and host redundant pathways

Apicomplexans are alveolate parasites which include Plasmodium falciparum, the main cause of malaria, one of the world's biggest killers from infectious disease. Apicomplexans are characterized by a reduction in proteome size, which appears to result from metabolic and functional simplification, commensurate with their parasitic lifestyle. However, other factors may also help to explain gene loss such as population bottlenecks experienced during transmission, and the effect of reducing the overall genomic information content. The latter constitutes an 'informational constraint', which is proposed to exert a selective pressure to evolve and maintain genes involved in informational fidelity and error correction, proportional to the quantity of information in the genome (which approximates to proteome size). The dynamics of gene loss was examined in 41 Apicomplexan genomes using orthogroup analysis. We show that loss of genes involved in amino acid metabolism and steroid biosynthesis can be explained by metabolic redundancy with the host. We also show that there is a marked tendency to lose DNA repair genes as proteome size is reduced. This may be explained by a reduction in size of the informational constraint and can help to explain elevated mutation rates in pathogens with reduced genome size. Multiple Sequentially Markovian Coalescent (MSMC) analysis indicates a recent bottleneck, consistent with predictions generated using allele-based population genetics approaches, implying that relaxed selection pressure due to reduced population size might have contributed to gene loss. However, the non-randomness of pathways that are lost challenges this scenario. Lastly, we identify unique orthogroups in malaria-causing Plasmodium species that infect humans, with a high proportion of membrane associated proteins. Thus, orthogroup analysis appears useful for identifying novel candidate pathogenic factors in parasites, when there is a wide sample of genomes available.

Modulating Isoprenoid Biosynthesis Increases Lipooligosaccharides and Restores Acinetobacter baumannii Resistance to Host and Antibiotic Stress

Acinetobacter baumannii is a leading cause of ventilator-associated pneumonia and a critical threat due to multidrug resistance. The A. baumannii outer membrane is an asymmetric lipid bilayer composed of inner leaflet glycerophospholipids and outer leaflet lipooligosaccharides. Deleting mlaF of the maintenance of lipid asymmetry (Mla) system causes A. baumannii to become more susceptible to pulmonary surfactants and antibiotics and decreases bacterial survival in the lungs of mice. Spontaneous suppressor mutants isolated from infected mice contain an ISAba11 insertion upstream of the ispB initiation codon, an essential isoprenoid biosynthesis gene. The insertion restores antimicrobial resistance and virulence to ΔmlaF. The suppressor strain increases lipooligosaccharides, suggesting that the mechanism involves balancing the glycerophospholipids/lipooligosaccharides ratio on the bacterial surface. An identical insertion exists in an extensively drug-resistant A. baumannii isolate, demonstrating its clinical relevance. These data show that the stresses bacteria encounter during infection select for genomic rearrangements that increase resistance to antimicrobials.

Efficacy and Safety of Fosmidomycin-Piperaquine as Nonartemisinin-Based Combination Therapy for Uncomplicated Falciparum Malaria: A Single-Arm, Age De-escalation Proof-of-Concept Study in Gabon

Background

Fosmidomycin–piperaquine is being developed as nonartemisinin-based combination therapy to meet the challenge of emerging artemisinin resistance.

Methods

The study was a phase 2, single-arm, proof-of-concept study of the efficacy, tolerability, and safety of fosmidomycin–piperaquine for the treatment of uncomplicated Plasmodium falciparum monoinfection in Gabon. Adults and children of both sexes with initial parasite counts between 1000 and 150000/µL received oral treatment with fosmidomycin (twice daily doses of 30 mg/kg) and piperaquine (once daily dose of 16 mg/kg) for 3 days and followed-up for 63 days. The primary efficacy endpoint was the per-protocol polymerase chain reaction (PCR)–corrected day 28 adequate clinical and parasitological response (ACPR).

Results

One hundred patients were enrolled. The PCR-corrected day 28 ACPR rate was 83/83, or 100% (95% confidence interval, 96–100). Fourteen patients had asexual parasitaemia between day 28 and day 63; all were typed by PCR as new infections. Fosmidomycin–piperaquine therapy led to rapid parasite clearance (median, 36 hours; interquartile range [IQR], 6–60) and fever clearance time (median, 12 hours; IQR, 6–48). The electrocardiogram assessments showed 2 patients with prolonged QT interval >500 msec following study drug administration. The majority of adverse events affected the gastrointestinal and respiratory tracts and were transient and mild to moderate in severity.

Conclusions

This is the first report of the use of the combination fosmidomycin–piperaquine. The combination appeared to have high efficacy and be safe and well tolerated despite observed transient changes in electrocardiogram with prolongation of the QT interval.

Clinical Trials Registration. NCT02198807.

Genetic mapping of fitness determinants across the malaria parasite Plasmodium falciparum life cycle

Determining the genetic basis of fitness is central to understanding evolution and transmission of microbial pathogens. In human malaria parasites (Plasmodium falciparum), most experimental work on fitness has focused on asexual blood stage parasites, because this stage can be easily cultured, although the transmission of malaria requires both female Anopheles mosquitoes and vertebrate hosts. We explore a powerful approach to identify the genetic determinants of parasite fitness across both invertebrate and vertebrate life-cycle stages of P. falciparum. This combines experimental genetic crosses using humanized mice, with selective whole genome amplification and pooled sequencing to determine genome-wide allele frequencies and identify genomic regions under selection across multiple lifecycle stages. We applied this approach to genetic crosses between artemisinin resistant (ART-R, kelch13-C580Y) and ART-sensitive (ART-S, kelch13-WT) parasites, recently isolated from Southeast Asian patients. Two striking results emerge: we observed (i) a strong genome-wide skew (>80%) towards alleles from the ART-R parent in the mosquito stage, that dropped to ~50% in the blood stage as selfed ART-R parasites were selected against; and (ii) repeatable allele specific skews in blood stage parasites with particularly strong selection (selection coefficient (s) ≤ 0.18/asexual cycle) against alleles from the ART-R parent at loci on chromosome 12 containing MRP2 and chromosome 14 containing ARPS10. This approach robustly identifies selected loci and has strong potential for identifying parasite genes that interact with the mosquito vector or compensatory loci involved in drug resistance.

Evaluating 18s-rRNA LAMP and selective whole genome amplification (sWGA) assay in detecting asymptomatic Plasmodium falciparum infections in blood donors

Background

Undesirable consequences of donor Plasmodium falciparum parasitaemia on stored donor blood have been reported. Therefore, it is imperative that all prospective blood donors are screened for P. falciparum infections using sensitive techniques. In this study, the sensitivities of microscopy, rapid diagnostic test (RDT), loop-mediated isothermal amplification (LAMP) assay and selective whole genome amplification (sWGA) technique in detecting P. falciparum infections in blood donors was assessed.

Methods

Randomly selected blood donors from 5 districts in Greater Accra Region of Ghana were screened for asymptomatic P. falciparum infections. Each donor sample was screened with SD Bioline RDT kit for P. falciparum histidine rich protein 2 and Plasmodium lactate dehydrogenase antigens, sWGA and 18s-rRNA LAMP. Crude DNA LAMP (crDNA-LAMP) was compared to purified DNA LAMP (pDNA-LAMP).

Results

A total of 771 blood donors were screened. The respective overall prevalence of P. falciparum in Ghana by microscopy, RDT, crDNA-LAMP, pDNA-LAMP and sWGA was 7.4%, 11.8%, 16.9%, 17.5% and 18.0%. Using sWGA as the reference test, the sensitivities of microscopy, RDT, crDNA-LAMP and pDNA-LAMP were 41.0% (95% CI 32.7–49.7), 65.5% (95% CI 56.9–73.3), 82.6% (95% CI 75.8–88.3) and 95.7% (95% CI 90.1–98.4), respectively. There was near perfect agreement between LAMP and sWGA (sWGA vs. crDNA-LAMP, κ = 0.87; sWGA vs. pDNA-LAMP, κ = 0.96), while crDNA-LAMP and pDNA-LAMP agreed perfectly (κ = 0.91). Goodness of fit test indicated non-significant difference between the performance of LAMP and sWGA (crDNA-LAMP vs. sWGA: x² = 0.71, p = 0.399 and pDNA-LAMP vs. sWGA: x² = 0.14, p = 0.707). Finally, compared to sWGA, the performance of LAMP did not differ in detecting sub-microscopic parasitaemia (sWGA vs. crDNA-LAMP: x² = 1.12, p = 0.290 and sWGA vs. pDNA-LAMP: x² = 0.22, p = 0.638).

Conclusions

LAMP assay agreed near perfectly with sWGA with non-significant differences in their ability to detect asymptomatic P. falciparum parasitaemia in blood donors. Therefore, it is recommended that LAMP based assays are employed to detect P. falciparum infections in blood donors due to its high sensitivity, simplicity, cost-effectiveness and user-friendliness.

swga: a primer design toolkit for selective whole genome amplification

Motivation

Population genomic analyses are often hindered by difficulties in obtaining sufficient numbers of genomes for analysis by DNA sequencing. Selective whole-genome amplification (SWGA) provides an efficient approach to amplify microbial genomes from complex backgrounds for sequence acquisition. However, the process of designing sets of primers for this method has many degrees of freedom and would benefit from an automated process to evaluate the vast number of potential primer sets.

Results

Here, we present swga, a program that identifies primer sets for SWGA and evaluates them for efficiency and selectivity. We used swga to design and test primer sets for the selective amplification of Wolbachia pipientis genomic DNA from infected Drosophila melanogaster and Mycobacterium tuberculosis from human blood. We identify primer sets that successfully amplify each against their backgrounds and describe a general method for using swga for arbitrary targets. In addition, we describe characteristics of primer sets that correlate with successful amplification, and present guidelines for implementation of SWGA to detect new targets.

Availability and Implementation

Source code and documentation are freely available on https://www.github.com/eclarke/swga. The program is implemented in Python and C and licensed under the GNU Public License.

Supplementary information

Supplementary data are available at Bioinformatics online.

Direct whole-genome sequencing of Plasmodium falciparum specimens from dried erythrocyte spots

BACKGROUND

Plasmodium falciparum malaria remains a major health burden and genomic research represents one of the necessary approaches for continued progress towards malaria control and elimination. Sample acquisition for this purpose is troublesome, with the majority of malaria-infected individuals living in rural areas, away from main infrastructure and the electrical grid. The aim of this study was to describe a low-tech procedure to sample P. falciparum specimens for direct whole genome sequencing (WGS), without use of electricity and cold-chain.

METHODS

Venous blood samples were collected from malaria patients in Bandim, Guinea-Bissau and leukocyte-depleted using Plasmodipur filters, the enriched parasite sample was spotted on Whatman paper and dried. The samples were stored at ambient temperatures and subsequently used for DNA-extraction. Ratios of parasite:human content of the extracted DNA was assessed by qPCR, and five samples with varying parasitaemia, were sequenced. Sequencing data were used to analyse the sample content, as well as sample coverage and depth as compared to the 3d7 reference genome.

RESULTS

qPCR revealed that 73% of the 199 samples were applicable for WGS, as defined by a minimum ratio of parasite:human DNA of 2:1. WGS revealed an even distribution of sequence data across the 3d7 reference genome, regardless of parasitaemia. The acquired read depths varied from 16 to 99×, and coverage varied from 87.5 to 98.9% of the 3d7 reference genome. SNP-analysis of six genes, for which amplicon sequencing has been performed previously, confirmed the reliability of the WGS-data.

CONCLUSION

This study describes a simple filter paper based protocol for sampling P. falciparum from malaria patients for subsequent direct WGS, enabling acquisition of samples in remote settings with no access to electricity.

Small molecule inhibition of apicomplexan FtsH1 disrupts plastid biogenesis in human pathogens

The malaria parasite Plasmodium falciparum and related apicomplexan pathogens contain an essential plastid organelle, the apicoplast, which is a key anti-parasitic target. Derived from secondary endosymbiosis, the apicoplast depends on novel, but largely cryptic, mechanisms for protein/lipid import and organelle inheritance during parasite replication. These critical biogenesis pathways present untapped opportunities to discover new parasite-specific drug targets. We used an innovative screen to identify actinonin as having a novel mechanism-of-action inhibiting apicoplast biogenesis. Resistant mutation, chemical-genetic interaction, and biochemical inhibition demonstrate that the unexpected target of actinonin in P. falciparum and Toxoplasma gondii is FtsH1, a homolog of a bacterial membrane AAA+ metalloprotease. PfFtsH1 is the first novel factor required for apicoplast biogenesis identified in a phenotypic screen. Our findings demonstrate that FtsH1 is a novel and, importantly, druggable antimalarial target. Development of FtsH1 inhibitors will have significant advantages with improved drug kinetics and multistage efficacy against multiple human parasites.

Selective Whole-Genome Amplification Is a Robust Method That Enables Scalable Whole-Genome Sequencing of Plasmodium vivax from Unprocessed Clinical Samples

Whole-genome sequencing (WGS) of microbial pathogens from clinical samples is a highly sensitive tool used to gain a deeper understanding of the biology, epidemiology, and drug resistance mechanisms of many infections. However, WGS of organisms which exhibit low densities in their hosts is challenging due to high levels of host genomic DNA (gDNA), which leads to very low coverage of the microbial genome. WGS of Plasmodium vivax, the most widely distributed form of malaria, is especially difficult because of low parasite densities and the lack of an ex vivo culture system. Current techniques used to enrich P. vivax DNA from clinical samples require significant resources or are not consistently effective. Here, we demonstrate that selective whole-genome amplification (SWGA) can enrich P. vivax gDNA from unprocessed human blood samples and dried blood spots for high-quality WGS, allowing genetic characterization of isolates that would otherwise have been prohibitively expensive or impossible to sequence. We achieved an average genome coverage of 24×, with up to 95% of the P. vivax core genome covered by ≥5 reads. The single-nucleotide polymorphism (SNP) characteristics and drug resistance mutations seen were consistent with those of other P. vivax sequences from a similar region in Peru, demonstrating that SWGA produces high-quality sequences for downstream analysis. SWGA is a robust tool that will enable efficient, cost-effective WGS of P. vivax isolates from clinical samples that can be applied to other neglected microbial pathogens.

Malaria is a disease caused by Plasmodium parasites that caused 214 million symptomatic cases and 438,000 deaths in 2015. Plasmodium vivax is the most widely distributed species, causing the majority of malaria infections outside sub-Saharan Africa. Whole-genome sequencing (WGS) of Plasmodium parasites from clinical samples has revealed important insights into the epidemiology and mechanisms of drug resistance of malaria. However, WGS of P. vivax is challenging due to low parasite levels in humans and the lack of a routine system to culture the parasites. Selective whole-genome amplification (SWGA) preferentially amplifies the genomes of pathogens from mixtures of target and host gDNA. Here, we demonstrate that SWGA is a simple, robust method that can be used to enrich P. vivax genomic DNA (gDNA) from unprocessed human blood samples and dried blood spots for cost-effective, high-quality WGS.

Whole genome sequencing of Plasmodium falciparum from dried blood spots using selective whole genome amplification

Background

Translating genomic technologies into healthcare applications for the malaria parasite Plasmodium falciparum has been limited by the technical and logistical difficulties of obtaining high quality clinical samples from the field. Sampling by dried blood spot (DBS) finger-pricks can be performed safely and efficiently with minimal resource and storage requirements compared with venous blood (VB). Here, the use of selective whole genome amplification (sWGA) to sequence the P. falciparum genome from clinical DBS samples was evaluated, and the results compared with current methods that use leucodepleted VB.

Methods

Parasite DNA with high (>95%) human DNA contamination was selectively amplified by Phi29 polymerase using short oligonucleotide probes of 8–12 mers as primers. These primers were selected on the basis of their differential frequency of binding the desired (P. falciparum DNA) and contaminating (human) genomes.

Results

Using sWGA method, clinical samples from 156 malaria patients, including 120 paired samples for head-to-head comparison of DBS and leucodepleted VB were sequenced. Greater than 18-fold enrichment of P. falciparum DNA was achieved from DBS extracts. The parasitaemia threshold to achieve >5× coverage for 50% of the genome was 0.03% (40 parasites per 200 white blood cells). Over 99% SNP concordance between VB and DBS samples was achieved after excluding missing calls.

Conclusion

The sWGA methods described here provide a reliable and scalable way of generating P. falciparum genome sequence data from DBS samples. The current data indicate that it will be possible to get good quality sequence on most if not all drug resistance loci from the majority of symptomatic malaria patients. This technique overcomes a major limiting factor in P. falciparum genome sequencing from field samples, and paves the way for large-scale epidemiological applications.

Electronic supplementary material

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