Development of amplicon deep sequencing markers and data analysis pipeline for genotyping multi-clonal malaria infections

Sensitive, highly multiplexed sequencing of microhaplotypes from the Plasmodium falciparum heterozygome

Background

Targeted next-generation sequencing offers the potential for consistent, deep coverage of information-rich genomic regions to characterize polyclonal Plasmodium falciparum infections. However, methods to identify and sequence these genomic regions are currently limited.

Methods

A bioinformatic pipeline and multiplex methods were developed to identify and simultaneously sequence 100 targets and applied to dried blood spot (DBS) controls and field isolates from Mozambique. For comparison, whole-genome sequencing data were generated for the same controls.

Results

Using publicly available genomes, 4465 high-diversity genomic regions suited for targeted sequencing were identified, representing the P. falciparum heterozygome. For this study, 93 microhaplotypes with high diversity (median expected heterozygosity = 0.7) were selected along with 7 drug resistance loci. The sequencing method achieved very high coverage (median 99%), specificity (99.8%), and sensitivity (90% for haplotypes with 5% within sample frequency in dried blood spots with 100 parasites/µL). In silico analyses revealed that microhaplotypes provided much higher resolution to discriminate related from unrelated polyclonal infections than biallelic single-nucleotide polymorphism barcodes.

Conclusions

The bioinformatic and laboratory methods outlined here provide a flexible tool for efficient, low-cost, high-throughput interrogation of the P. falciparum genome, and can be tailored to simultaneously address multiple questions of interest in various epidemiological settings.

Microgeographic Epidemiology of Malaria Parasites in an Irrigated Area of Western Kenya by Deep Amplicon Sequencing

To improve food security, investments in irrigated agriculture are anticipated to increase throughout Africa. However, the extent to which environmental changes from water resource development will impact malaria epidemiology remains unclear. This study was designed to compare the sensitivity of molecular markers used in deep amplicon sequencing for evaluating malaria transmission intensities and to assess malaria transmission intensity at various proximities to an irrigation scheme. Compared to ama1, csp, and msp1 amplicons, cpmp required the smallest sample size to detect differences in infection complexity between transmission risk zones. Transmission intensity was highest within 5 km of the irrigation scheme by polymerase chain reaction positivity rate, infection complexity, and linkage disequilibrium. The irrigated area provided a source of parasite infections for the surrounding 2- to 10-km area. This study highlights the suitability of the cpmp amplicon as a measure for transmission intensities and the impact of irrigation on microgeographic epidemiology of malaria parasites.

Implementing parasite genotyping into national surveillance frameworks: feedback from control programmes and researchers in the Asia-Pacific region

The Asia-Pacific region faces formidable challenges in achieving malaria elimination by the proposed target in 2030. Molecular surveillance of Plasmodium parasites can provide important information on malaria transmission and adaptation, which can inform national malaria control programmes (NMCPs) in decision-making processes. In November 2019 a parasite genotyping workshop was held in Jakarta, Indonesia, to review molecular approaches for parasite surveillance and explore ways in which these tools can be integrated into public health systems and inform policy. The meeting was attended by 70 participants from 8 malaria-endemic countries and partners of the Asia Pacific Malaria Elimination Network. The participants acknowledged the utility of multiple use cases for parasite genotyping including: quantifying the prevalence of drug resistant parasites, predicting risks of treatment failure, identifying major routes and reservoirs of infection, monitoring imported malaria and its contribution to local transmission, characterizing the origins and dynamics of malaria outbreaks, and estimating the frequency of Plasmodium vivax relapses. However, the priority of each use case varies with different endemic settings. Although a one-size-fits-all approach to molecular surveillance is unlikely to be applicable across the Asia-Pacific region, consensus on the spectrum of added-value activities will help support data sharing across national boundaries. Knowledge exchange is needed to establish local expertise in different laboratory-based methodologies and bioinformatics processes. Collaborative research involving local and international teams will help maximize the impact of analytical outputs on the operational needs of NMCPs. Research is also needed to explore the cost-effectiveness of genetic epidemiology for different use cases to help to leverage funding for wide-scale implementation. Engagement between NMCPs and local researchers will be critical throughout this process.

The use of spatial and genetic tools to assess Plasmodium falciparum transmission in Lusaka, Zambia between 2011 and 2015

BACKGROUND

Zambia has set itself the ambitious target of eliminating malaria by 2021. To continue tracking transmission to zero, new interventions, tools and approaches are required.

METHODS

Urban reactive case detection (RCD) was performed in Lusaka city from 2011 to 2015 to better understand the location and drivers of malaria transmission. Briefly, index cases were followed to their home and all consenting individuals living in the index house and nine proximal houses were tested with a malaria rapid diagnostic test and treated if positive. A brief survey was performed and for certain responses, a dried blood spot sample collected for genetic analysis. Aggregate health facility data, individual RCD response data and genetic results were analysed spatially and against environmental correlates.

RESULTS

Total number of malaria cases remained relatively constant, while the average age of incident cases and the proportion of incident cases reporting recent travel both increased. The estimated R₀ in Lusaka was < 1 throughout the study period. RCD responses performed within 250 m of uninhabited/vacant land were associated with a higher probability of identifying additional infections.

CONCLUSIONS

Evidence suggests that the majority of malaria infections are imported from outside Lusaka. However there remains some level of local transmission occurring on the periphery of urban settlements, namely in the wet season. Unfortunately, due to the higher-than-expected complexity of infections and the small number of samples tested, genetic analysis was unable to identify any meaningful trends in the data.

PCR correction strategies for malaria drug trials: updates and clarifications

Malaria drug trials conducted in endemic areas face a major challenge in their analysis because it is difficult to establish whether parasitaemia in blood samples collected after treatment indicate drug failure or a new infection acquired after treatment. It is therefore vital to reliably distinguish drug failures from new infections in order to obtain accurate estimates of drug failure rates. This distinction can be achieved for Plasmodium falciparum by comparing parasite genotypes obtained at the time of treatment (the baseline) and on the day of recurring parasitaemia. Such PCR correction is required to obtain accurate failure rates, even for new effective drugs. Despite the routine use of PCR correction in surveillance of drug resistance and in clinical drug trials, limitations inherent to the molecular genotyping methods have led some researchers to question the validity of current PCR correction strategies. Here we describe and discuss recent developments in these genotyping approaches, with a particular focus on method validation and limitations of the genotyping strategies. Our aim is to update scientists from public and private bodies who are working on the development, deployment, and surveillance of new malaria drugs. We aim to promote discussion around these issues and argue for the adoption of improved standardised PCR correction methodologies.

Strengthening Surveillance Systems for Malaria Elimination by Integrating Molecular and Genomic Data

Unprecedented efforts in malaria control over the last 15 years have led to a substantial decrease in both morbidity and mortality in most endemic settings. However, these progresses have stalled over recent years, and resurgence may cause dramatic impact on both morbidity and mortality. Nevertheless, elimination efforts are currently going on with the objective of reducing malaria morbidity and mortality by 90% and malaria elimination in at least 35 countries by 2030. Strengthening surveillance systems is of paramount importance to reach those targets, and the integration of molecular and genomic techniques into routine surveillance could substantially improve the quality and robustness of data. Techniques such as polymerase chain reaction (PCR) and quantitative PCR (qPCR) are increasingly available in malaria endemic countries, whereas others such as sequencing are already available in a few laboratories. However, sequencing, especially next-generation sequencing (NGS), requires sophisticated infrastructure with adequate computing power and highly trained personnel for data analysis that require substantial investment. Different techniques will be required for different applications, and cost-effective planning must ensure the appropriate use of available resources. The development of national and sub-regional reference laboratories could help in minimizing the resources required in terms of equipment and trained staff. Concerted efforts from different stakeholders at national, sub-regional, and global level are needed to develop the required framework to establish and maintain these reference laboratories.

Amplicon deep sequencing improves Plasmodium falciparum genotyping in clinical trials of antimalarial drugs

Clinical trials monitoring malaria drug resistance require genotyping of recurrent Plasmodium falciparum parasites to distinguish between treatment failure and new infection occurring during the trial follow up period. Because trial participants usually harbour multi-clonal P. falciparum infections, deep amplicon sequencing (AmpSeq) was employed to improve sensitivity and reliability of minority clone detection. Paired samples from 32 drug trial participants were Illumina deep-sequenced for five molecular markers. Reads were analysed by custom-made software HaplotypR and trial outcomes compared to results from the previous standard genotyping method based on length-polymorphic markers. Diversity of AmpSeq markers in pre-treatment samples was comparable or higher than length-polymorphic markers. AmpSeq was highly reproducible with consistent quantification of co-infecting parasite clones within a host. Outcomes of the three best-performing markers, cpmp, cpp and ama1-D3, agreed in 26/32 (81%) of patients. Discordance between the three markers performed per sample was much lower by AmpSeq (six patients) compared to length-polymorphic markers (eleven patients). Using AmpSeq for discrimination of recrudescence and new infection in antimalarial drug trials provides highly reproducible and robust characterization of clone dynamics during trial follow-up. AmpSeq overcomes limitations inherent to length-polymorphic markers. Regulatory clinical trials of antimalarial drugs will greatly benefit from this unbiased typing method.

Asymptomatic malaria in the clinical and public health context

Introduction: Historically, the global community has focused on the control of symptomatic malaria. However, interest in asymptomatic malaria has been growing, particularly in the context of malaria elimination.Areas covered: We undertook a comprehensive PubMed literature review on asymptomatic malaria as it relates to detection and elimination with emphasis between 2014 and 2019. Diagnostic tools with a low limit of detection (LOD) have allowed us to develop a more detailed understanding of asymptomatic malaria and its impact. These highly sensitive diagnostics have demonstrated that the prevalence of asymptomatic malaria is greater than previously thought. In addition, it is now possible to detect the malaria reservoir in the community, something that was previously not feasible. Asymptomatic malaria has previously not been treated, but research has begun to examine whether treating individuals with asymptomatic malaria may lead to health benefits. Finally, we have begun to understand the importance of asymptomatic malaria in ongoing transmission.Expert opinion: Therefore, with malaria elimination back on the agenda, asymptomatic malaria can no longer be ignored, especially in light of new ultra-sensitive diagnostic tools.

Applying next-generation sequencing to track falciparum malaria in sub-Saharan Africa

Next-generation sequencing (NGS) technologies are increasingly being used to address a diverse range of biological and epidemiological questions. The current understanding of malaria transmission dynamics and parasite movement mainly relies on the analyses of epidemiologic data, e.g. case counts and self-reported travel history data. However, travel history data are often not routinely collected or are incomplete, lacking the necessary level of accuracy. Although genetic data from routinely collected field samples provides an unprecedented opportunity to track the spread of malaria parasites, it remains an underutilized resource for surveillance due to lack of local awareness and capacity, limited access to sensitive laboratory methods and associated computational tools and difficulty in interpreting genetic epidemiology data. In this review, the potential roles of NGS in better understanding of transmission patterns, accurately tracking parasite movement and addressing the emerging challenges of imported malaria in low transmission settings of sub-Saharan Africa are discussed. Furthermore, this review highlights the insights gained from malaria genomic research and challenges associated with integrating malaria genomics into existing surveillance tools to inform control and elimination strategies.

Detection of low-density Plasmodium falciparum infections using amplicon deep sequencing

BACKGROUND

Deep sequencing of targeted genomic regions is becoming a common tool for understanding the dynamics and complexity of Plasmodium infections, but its lower limit of detection is currently unknown. Here, a new amplicon analysis tool, the Parallel Amplicon Sequencing Error Correction (PASEC) pipeline, is used to evaluate the performance of amplicon sequencing on low-density Plasmodium DNA samples. Illumina-based sequencing of two Plasmodium falciparum genomic regions (CSP and SERA2) was performed on two types of samples: in vitro DNA mixtures mimicking low-density infections (1-200 genomes/μl) and extracted blood spots from a combination of symptomatic and asymptomatic individuals (44-653,080 parasites/μl). Three additional analysis tools-DADA2, HaplotypR, and SeekDeep-were applied to both datasets and the precision and sensitivity of each tool were evaluated.

RESULTS

Amplicon sequencing can contend with low-density samples, showing reasonable detection accuracy down to a concentration of 5 Plasmodium genomes/μl. Due to increased stochasticity and background noise, however, all four tools showed reduced sensitivity and precision on samples with very low parasitaemia (< 5 copies/μl) or low read count (< 100 reads per amplicon). PASEC could distinguish major from minor haplotypes with an accuracy of 90% in samples with at least 30 Plasmodium genomes/μl, but only 61% at low Plasmodium concentrations (< 5 genomes/μl) and 46% at very low read counts (< 25 reads per amplicon). The four tools were additionally used on a panel of extracted parasite-positive blood spots from natural malaria infections. While all four identified concordant patterns of complexity of infection (COI) across four sub-Saharan African countries, COI values obtained for individual samples differed in some cases.

CONCLUSIONS

Amplicon deep sequencing can be used to determine the complexity and diversity of low-density Plasmodium infections. Despite differences in their approach, four state-of-the-art tools resolved known haplotype mixtures with similar sensitivity and precision. Researchers can therefore choose from multiple robust approaches for analysing amplicon data, however, error filtration approaches should not be uniformly applied across samples of varying parasitaemia. Samples with very low parasitaemia and very low read count have higher false positive rates and call for read count thresholds that are higher than current default recommendations.

Genetic diversity and drug resistance surveillance of Plasmodium falciparum for malaria elimination: is there an ideal tool for resource-limited sub-Saharan Africa?

The intensification of malaria control interventions has resulted in its global decline, but it remains a significant public health burden especially in sub-Saharan Africa (sSA). Knowledge on the parasite diversity, its transmission dynamics, mechanisms of adaptation to environmental and interventional pressures could help refine or develop new control and elimination strategies. Critical to this is the accurate assessment of the parasite’s genetic diversity and monitoring of genetic markers of anti-malarial resistance across all susceptible populations. Such wide molecular surveillance will require selected tools and approaches from a variety of ever evolving advancements in technology and the changing epidemiology of malaria. The choice of an effective approach for specific endemic settings remains challenging, particularly for countries in sSA with limited access to advanced technologies. This article examines the current strategies and tools for Plasmodium falciparum genetic diversity typing and resistance monitoring and proposes how the different tools could be employed in resource-poor settings. Advanced approaches enabling targeted deep sequencing is valued as a sensitive method for assessing drug resistance and parasite diversity but remains out of the reach of most laboratories in sSA due to the high cost of development and maintenance. It is, however, feasible to equip a limited number of laboratories as Centres of Excellence in Africa (CEA), which will receive and process samples from a network of peripheral laboratories in the continent. Cheaper, sensitive and portable real-time PCR methods can be used in peripheral laboratories to pre-screen and select samples for targeted deep sequence or genome wide analyses at these CEAs.

Using parasite genetic and human mobility data to infer local and cross-border malaria connectivity in Southern Africa

Local and cross-border importation remain major challenges to malaria elimination and are difficult to measure using traditional surveillance data. To address this challenge, we systematically collected parasite genetic data and travel history from thousands of malaria cases across northeastern Namibia and estimated human mobility from mobile phone data. We observed strong fine-scale spatial structure in local parasite populations, providing positive evidence that the majority of cases were due to local transmission. This result was largely consistent with estimates from mobile phone and travel history data. However, genetic data identified more detailed and extensive evidence of parasite connectivity over hundreds of kilometers than the other data, within Namibia and across the Angolan and Zambian borders. Our results provide a framework for incorporating genetic data into malaria surveillance and provide evidence that both strengthening of local interventions and regional coordination are likely necessary to eliminate malaria in this region of Southern Africa.

Longitudinal tracking and quantification of individual Plasmodium falciparum clones in complex infections

Longitudinal tracking of individual Plasmodium falciparum strains in multi-clonal infections is essential for investigating infection dynamics of malaria. The traditional genotyping techniques did not permit tracking changes in individual clone density during persistent natural infections. Amplicon deep sequencing (Amp-Seq) offers a tool to address this knowledge gap. The sensitivity of Amp-Seq for relative quantification of clones was investigated using three molecular markers, ama1-D2, ama1-D3, and cpmp. Amp-Seq and length-polymorphism based genotyping were compared for their performance in following minority clones in longitudinal samples from Papua New Guinea. Amp-Seq markers were superior to length-polymorphic marker msp2 in detecting minority clones (sensitivity Amp-Seq: 95%, msp2: 85%). Multiplicity of infection (MOI) by Amp-Seq was 2.32 versus 1.73 for msp2. The higher sensitivity had no effect on estimates of force of infection because missed minority clones were detected in preceding or succeeding bleeds. Individual clone densities were tracked longitudinally by Amp-Seq despite MOI > 1, thus providing an additional parameter for investigating malaria infection dynamics. Amp-Seq based genotyping of longitudinal samples improves detection of minority clones and estimates of MOI. Amp-Seq permits tracking of clone density over time to study clone competition or the dynamics of specific, i.e. resistance-associated genotypes.

Multiplicity of Infection Is a Poor Predictor of Village-Level Plasmodium vivax and P. falciparum Population Prevalence in the Southwest Pacific

Across 8101 individuals in 46 villages, the proportion of Plasmodium spp. multiple clone infections (0%–53.8%) did not reflect prevalence by quantitative polymerase chain reaction (qPCR; 1.9%–38.4%), except for P. vivax in Solomon Islands (P < .001). Thus this parameter was not informative to identify transmission foci. In contrast, prevalence by microscopy and qPCR correlated well.

Multiplicity and molecular epidemiology of Plasmodium vivax and Plasmodium falciparum infections in East Africa

Background

Parasite genetic diversity and multiplicity of infection (MOI) affect clinical outcomes, response to drug treatment and naturally-acquired or vaccine-induced immunity. Traditional methods often underestimate the frequency and diversity of multiclonal infections due to technical sensitivity and specificity. Next-generation sequencing techniques provide a novel opportunity to study complexity of parasite populations and molecular epidemiology.

Methods

Symptomatic and asymptomatic Plasmodium vivax samples were collected from health centres/hospitals and schools, respectively, from 2011 to 2015 in Ethiopia. Similarly, both symptomatic and asymptomatic Plasmodium falciparum samples were collected, respectively, from hospitals and schools in 2005 and 2015 in Kenya. Finger-pricked blood samples were collected and dried on filter paper. Long amplicon (> 400 bp) deep sequencing of merozoite surface protein 1 (msp1) gene was conducted to determine multiplicity and molecular epidemiology of P. vivax and P. falciparum infections. The results were compared with those based on short amplicon (117 bp) deep sequencing.

Results

A total of 139 P. vivax and 222 P. falciparum samples were pyro-sequenced for pvmsp1 and pfmsp1, yielding a total of 21 P. vivax and 99 P. falciparum predominant haplotypes. The average MOI for P. vivax and P. falciparum were 2.16 and 2.68, respectively, which were significantly higher than that of microsatellite markers and short amplicon (117 bp) deep sequencing. Multiclonal infections were detected in 62.2% of the samples for P. vivax and 74.8% of the samples for P. falciparum. Four out of the five subjects with recurrent P. vivax malaria were found to be a relapse 44–65 days after clearance of parasites. No difference was observed in MOI among P. vivax patients of different symptoms, ages and genders. Similar patterns were also observed in P. falciparum except for one study site in Kenyan lowland areas with significantly higher MOI.

Conclusions

The study used a novel method to evaluate Plasmodium MOI and molecular epidemiological patterns by long amplicon ultra-deep sequencing. The complexity of infections were similar among age groups, symptoms, genders, transmission settings (spatial heterogeneity), as well as over years (pre- vs. post-scale-up interventions). This study demonstrated that long amplicon deep sequencing is a useful tool to investigate multiplicity and molecular epidemiology of Plasmodium parasite infections.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2337-y) contains supplementary material, which is available to authorized users.

Molecular approaches to determine the multiplicity of Plasmodium infections

Multiplicity of infection (MOI), also termed complexity of infection (COI), is defined as the number of genetically distinct parasite strains co-infecting a single host, which is an important indicator of malaria epidemiology. PCR-based genotyping often underestimates MOI. Next generation sequencing technologies provide much more accurate and genome-wide characterization of polyclonal infections. However, complete haplotype characterization of multiclonal infections remains a challenge due to PCR artifacts and sequencing errors, and requires efficient computational tools. In this review, the advantages and limitations of current molecular approaches to determine multiplicity of malaria parasite infection are discussed.

Evolutionary ARMS Race: Antimalarial Resistance Molecular Surveillance

Molecular surveillance of antimalarial drug resistance markers has become an important part of resistance detection and containment. In the current climate of multidrug resistance, including resistance to the global front-line drug artemisinin, there is a consensus to upscale molecular surveillance. The most salient limitation to current surveillance efforts is that skill and infrastructure requirements preclude many regions. This includes sub-Saharan Africa, where Plasmodium falciparum is responsible for most of the global malaria disease burden. New molecular and data technologies have emerged with an emphasis on accessibility. These may allow surveillance to be conducted in broad settings where it is most needed, including at the primary healthcare level in endemic countries, and extending to the village health worker.

Malaria Epidemiology at the Clone Level

Genotyping to distinguish between parasite clones is nowadays a standard in many molecular epidemiological studies of malaria. It has become crucial in drug trials and to follow individual clones in epidemiological studies, and to understand how drug resistance emerges and spreads. Here, we review the applications of the increasingly available genotyping tools and whole-genome sequencing data, and argue for a better integration of population genetics findings into malaria-control strategies.