A Method for Amplicon Deep Sequencing of Drug Resistance Genes in Plasmodium falciparum Clinical Isolates from India

Lineage-informative microhaplotypes for recurrence classification and spatio-temporal surveillance of Plasmodium vivax malaria parasites

Challenges in classifying recurrent Plasmodium vivax infections constrain surveillance of antimalarial efficacy and transmission. Recurrent infections may arise from activation of dormant liver stages (relapse), blood-stage treatment failure (recrudescence) or reinfection. Molecular inference of familial relatedness (identity-by-descent or IBD) can help resolve the probable origin of recurrences. As whole genome sequencing of P. vivax remains challenging, targeted genotyping methods are needed for scalability. We describe a P. vivax marker discovery framework to identify and select panels of microhaplotypes (multi-allelic markers within small, amplifiable segments of the genome) that can accurately capture IBD. We evaluate panels of 50-250 microhaplotypes discovered in a global set of 615 P. vivax genomes. A candidate global 100-microhaplotype panel exhibits high marker diversity in the Asia-Pacific, Latin America and horn of Africa (median HE = 0.70-0.81) and identifies 89% of the polyclonal infections detected with genome-wide datasets. Data simulations reveal lower error in estimating pairwise IBD using microhaplotypes relative to traditional biallelic SNP barcodes. The candidate global panel also exhibits high accuracy in predicting geographic origin and captures local infection outbreak and bottlenecking events. Our framework is open-source enabling customised microhaplotype discovery and selection, with potential for porting to other species or data resources.

PHARE: a bioinformatics pipeline for compositional profiling of multiclonal Plasmodium falciparum infections from long-read Nanopore sequencing data

BACKGROUND

The emergence of drug-resistant clones of Plasmodium falciparum is a major public health concern, and the ability to detect and track the spread of these clones is crucial for effective malaria control and treatment. However, in endemic settings, malaria infected people often carry multiple P. falciparum clones simultaneously making it likely to miss drug-resistant clones using traditional molecular typing methods.

OBJECTIVES

Our goal was to develop a bioinformatics pipeline for compositional profiling in multiclonal P. falciparum samples, sequenced using the Oxford Nanopore Technologies MinION platform.

METHODS

We developed the 'Finding P. falciparum haplotypes with resistance mutations in polyclonal infections' (PHARE) pipeline using existing bioinformatics tools and custom scripts written in python. PHARE was validated on three control datasets containing P. falciparum DNA of four laboratory strains at varying mixing ratios. Additionally, the pipeline was tested on clinical samples from children admitted to a paediatric hospital in the Central African Republic.

RESULTS

The PHARE pipeline achieved high recall and accuracy rates in all control datasets. The pipeline can be used on any gene and was tested with amplicons of the P. falciparum drug resistance marker genes pfdhps, pfdhfr and pfK13.

CONCLUSIONS

The PHARE pipeline helps to provide a more complete picture of drug resistance in the circulating P. falciparum population and can help to guide treatment recommendations. PHARE is freely available under the GNU Lesser General Public License v.3.0 on GitHub: https://github.com/Fippu/PHARE.

Molecular tools are crucial for malaria elimination

The eradication of Plasmodium parasites, responsible for malaria, is a daunting global public health task. It requires a comprehensive approach that addresses symptomatic, asymptomatic, and submicroscopic cases. Overcoming this challenge relies on harnessing the power of molecular diagnostic tools, as traditional methods like microscopy and rapid diagnostic tests fall short in detecting low parasitaemia, contributing to the persistence of malaria transmission. By precisely identifying patients of all types and effectively characterizing malaria parasites, molecular tools may emerge as indispensable allies in the pursuit of malaria elimination. Furthermore, molecular tools can also provide valuable insights into parasite diversity, drug resistance patterns, and transmission dynamics, aiding in the implementation of targeted interventions and surveillance strategies. In this review, we explore the significance of molecular tools in the pursuit of malaria elimination, shedding light on their key contributions and potential impact on public health.

Drug resistance and vaccine target surveillance of Plasmodium falciparum using nanopore sequencing in Ghana

Malaria results in over 600,000 deaths annually, with the highest burden of deaths in young children living in sub-Saharan Africa. Molecular surveillance can provide important information for malaria control policies, including detection of antimalarial drug resistance. However, genome sequencing capacity in malaria-endemic countries is limited. We designed and implemented an end-to-end workflow to detect Plasmodium falciparum antimalarial resistance markers and diversity in the vaccine target circumsporozoite protein (csp) using nanopore sequencing in Ghana. We analysed 196 clinical samples and showed that our method is rapid, robust, accurate and straightforward to implement. Importantly, our method could be applied to dried blood spot samples, which are readily collected in endemic settings. We report that P. falciparum parasites in Ghana are mostly susceptible to chloroquine, with persistent sulfadoxine-pyrimethamine resistance and no evidence of artemisinin resistance. Multiple single nucleotide polymorphisms were identified in csp, but their significance is uncertain. Our study demonstrates the feasibility of nanopore sequencing for malaria genomic surveillance in endemic countries.

Metagenomic next generation sequencing for studying antibiotic resistance genes in the environment

Bacterial antimicrobial resistance (AMR) is a persisting and growing threat to human health. Characterization of antibiotic resistance genes (ARGs) in the environment is important to understand and control ARG-associated microbial risks. Numerous challenges exist in monitoring ARGs in the environment, due to the extraordinary diversity of ARGs, low abundance of ARGs with respect to the complex environmental microbiomes, difficulties in linking ARGs with bacterial hosts by molecular methods, difficulties in achieving quantification and high throughput simultaneously, difficulties in assessing mobility potential of ARGs, and difficulties in determining the specific AMR determinant genes. Advances in the next generation sequencing (NGS) technologies and related computational and bioinformatic tools are facilitating rapid identification and characterization ARGs in genomes and metagenomes from environmental samples. This chapter discusses NGS-based strategies, including amplicon-based sequencing, whole genome sequencing, bacterial population-targeted metagenome sequencing, metagenomic NGS, quantitative metagenomic sequencing, and functional/phenotypic metagenomic sequencing. Current bioinformatic tools for analyzing sequencing data for studying environmental ARGs are also discussed.

Antimalarial Drug Resistance Profiling of Plasmodium falciparum Infections in India Using Next-Generation Sequencing

Background

Tracking the emergence and spread of antimalarial drug resistance has become critical to sustaining progress towards the control and eventual elimination of malaria in South Asia, especially India.

Methods

An amplicon sequencing protocol was used for high-throughput molecular surveillance of antimalarial drug resistance in a total of 158 isolates at three sites in India: Chennai, Nadiad and Rourkela. Five genes of the Plasmodium falciparum implicated in antimalarial resistance were investigated here; Pfcrt for chloroquine resistance, Pfdhfr for pyrimethamine resistance, Pfdhps for sulfadoxine resistance, Pfk13 for artemisinin resistance and Pfmdr1 for resistance to multiple antimalarials.

Results

Mutations in the propeller domain of PfK13 were observed in two samples only, however these mutations are not validated for artemisinin resistance. A high proportion of parasites from the P. falciparum dominant site Rourkela showed wild-type Pfcrt and Pfdhfr haplotypes, while mutant Pfcrt and Pfdhfr haplotypes were fixed at the P. vivax dominant sites Chennai and Nadiad. The wild-type PfDHPS haplotype was predominant across all study sites. Finally, we observed the largest proportion of suspected multi-clonal infections at Rourkela, which has the highest transmission of P. falciparum among our study sites.

Conclusion

This is the first simultaneous high-throughput next generation sequencing of five complete P. falciparum genes from infected patients in India.

Insight into molecular diagnosis for antimalarial drug resistance of Plasmodium falciparum parasites: A review

Malaria is an infectious disease transmitted by the female Anopheles mosquito and poses a severe threat to human health. At present, antimalarial drugs are the primary treatment for malaria. The widespread use of artemisinin-based combination therapies (ACTs) has dramatically reduced the number of malaria-related deaths; however, the emergence of resistance has the potential to reverse this progress. Accurate and timely diagnosis of drug-resistant strains of Plasmodium parasites via detecting molecular markers (such as Pfnhe1, Pfmrp, Pfcrt, Pfmdr1, Pfdhps, Pfdhfr, and Pfk13) is essential for malaria control and elimination. Here, we review the current techniques which commonly used for molecular diagnosis of antimalarial resistance in P. falciparum and discuss their sensitivities and specificities for different drug resistance-associated molecular markers, with the aim of providing insights into possible directions for future precise point-of-care testing (POCT) of antimalarial drug resistance of malaria parasites.

Molecular assays for determining sulphadoxine-pyrimethamine drug resistance in India: a systematic review

A plethora of studies analyse the molecular markers of drug resistance and hence help in guiding the evidence-based malaria treatment policies in India. For reporting mutations, a number of techniques including DNA sequencing, restriction-fragment length polymorphism and mutation-specific polymerase chain reaction have been employed across numerous studies, including variations in the methodology used. However, there is no sufficient data from India comparing these methods as well as report the prevalence of polymorphisms in SP drug resistance molecular markers independently using such methods. Therefore, all data from Indian studies available for molecular marker studies of Plasmodium falciparum drug resistance to sulphadoxine-pyrimethamine was gathered, and a systematic review was performed. This systematic review identifies the molecular methods in use in India and compares each method for detecting sulphadoxine-pyrimethamine drug resistance marker. To delay the spread of drug-resistant parasite strains, a simplified and standardized molecular method is much needed which can be obtained by analysing the performance of each method in use and answering the necessity of newer methodological approaches.

Potential Opportunities and Challenges of Deploying Next Generation Sequencing and CRISPR-Cas Systems to Support Diagnostics and Surveillance Towards Malaria Control and Elimination in Africa

Recent developments in molecular biology and genomics have revolutionized biology and medicine mainly in the developed world. The application of next generation sequencing (NGS) and CRISPR-Cas tools is now poised to support endemic countries in the detection, monitoring and control of endemic diseases and future epidemics, as well as with emerging and re-emerging pathogens. Most low and middle income countries (LMICs) with the highest burden of infectious diseases still largely lack the capacity to generate and perform bioinformatic analysis of genomic data. These countries have also not deployed tools based on CRISPR-Cas technologies. For LMICs including Tanzania, it is critical to focus not only on the process of generation and analysis of data generated using such tools, but also on the utilization of the findings for policy and decision making. Here we discuss the promise and challenges of NGS and CRISPR-Cas in the context of malaria as Africa moves towards malaria elimination. These innovative tools are urgently needed to strengthen the current diagnostic and surveillance systems. We discuss ongoing efforts to deploy these tools for malaria detection and molecular surveillance highlighting potential opportunities presented by these innovative technologies as well as challenges in adopting them. Their deployment will also offer an opportunity to broadly build in-country capacity in pathogen genomics and bioinformatics, and to effectively engage with multiple stakeholders as well as policy makers, overcoming current workforce and infrastructure challenges. Overall, these ongoing initiatives will build the malaria molecular surveillance capacity of African researchers and their institutions, and allow them to generate genomics data and perform bioinformatics analysis in-country in order to provide critical information that will be used for real-time policy and decision-making to support malaria elimination on the continent.

Comparative analysis of targeted next-generation sequencing for Plasmodium falciparum drug resistance markers

Well-defined molecular resistance markers are available for a range of antimalarial drugs, and molecular surveillance is increasingly important for monitoring antimalarial drug resistance. Different genotyping platforms are available, but these have not been compared in detail. We compared Targeted Amplicon Deep sequencing (TADs) using Ion Torrent PGM with Illumina MiSeq for the typing of antimalarial drug resistance genes. We developed and validated protocols to type the molecular resistance markers pfcrt, pfdhfr, pfdhps, pfmdr1, pfkelch, and pfcytochrome b, in Plasmodium falciparum for the Ion Torrent PGM and Illumina MiSeq sequencing platforms. With P. falciparum 3D7 and K1 as reference strains, whole blood samples (N = 20) and blood spots from Rapid Diagnostic Test (RDT) samples (N = 5) from patients with uncomplicated falciparum malaria from Ubon Ratchathani were assessed on both platforms and compared for coverage (average reads per amplicon), sequencing accuracy, variant accuracy, false positive rate, false negative rate, and alternative allele detection, with conventional Sanger sequencing as the reference method for SNP calling. Both whole blood and RDT samples could be successfully sequenced using the Ion Torrent PGM and Illumina MiSeq platforms. Coverage of reads per amplicon was higher with Illumina MiSeq (28,886 reads) than with Ion Torrent PGM (1754 reads). In laboratory generated artificial mixed infections, the two platforms could detect the minor allele down to 1% density at 500X coverage. SNPs calls from both platforms were in complete agreement with conventional Sanger sequencing. The methods can be multiplexed with up to 96 samples per run, which reduces cost by 86% compared to conventional Sanger sequencing. Both platforms, using the developed TAD protocols, provide an accurate method for molecular surveillance of drug resistance markers in P. falciparum, but Illumina MiSeq provides higher coverage than Ion Torrent PGM.

Decrypting the complexity of the human malaria parasite biology through systems biology approaches

The human malaria parasite, Plasmodium falciparum, is a unicellular protozoan responsible for over half a million deaths annually. With a complex life cycle alternating between human and invertebrate hosts, this apicomplexan is notoriously adept at evading host immune responses and developing resistance to all clinically administered treatments. Advances in omics-based technologies, increased sensitivity of sequencing platforms and enhanced CRISPR based gene editing tools, have given researchers access to more in-depth and untapped information about this enigmatic micro-organism, a feat thought to be infeasible in the past decade. Here we discuss some of the most important scientific achievements made over the past few years with a focus on novel technologies and platforms that set the stage for subsequent discoveries. We also describe some of the systems-based methods applied to uncover gaps of knowledge left through single-omics applications with the hope that we will soon be able to overcome the spread of this life-threatening disease.

Advances and opportunities in malaria population genomics

Almost 20 years have passed since the first reference genome assemblies were published for Plasmodium falciparum, the deadliest malaria parasite, and Anopheles gambiae, the most important mosquito vector of malaria in sub-Saharan Africa. Reference genomes now exist for all human malaria parasites and nearly half of the ~40 important vectors around the world. As a foundation for genetic diversity studies, these reference genomes have helped advance our understanding of basic disease biology and drug and insecticide resistance, and have informed vaccine development efforts. Population genomic data are increasingly being used to guide our understanding of malaria epidemiology, for example by assessing connectivity between populations and the efficacy of parasite and vector interventions. The potential value of these applications to malaria control strategies, together with the increasing diversity of genomic data types and contexts in which data are being generated, raise both opportunities and challenges in the field. This Review discusses advances in malaria genomics and explores how population genomic data could be harnessed to further support global disease control efforts.

Should deep-sequenced amplicons become the new gold-standard for analysing malaria drug clinical trials?

Regulatory clinical trials are required to ensure the continued supply and deployment of effective antimalarial drugs. Patient follow-up in such trials typically lasts several weeks, as the drugs have long half-lives and new infections often occur during this period. “Molecular correction” is therefore used to distinguish drug failures from new infections. The current WHO-recommended method for molecular correction uses length-polymorphic alleles at highly diverse loci but is inherently poor at detecting low-density clones in polyclonal infections. This likely leads to substantial underestimates of failure rates, delaying the replacement of failing drugs with potentially lethal consequences. Deep-sequenced amplicons (AmpSeq) substantially increase the detectability of low-density clones and may offer a new “gold standard” for molecular correction. Pharmacological simulation of clinical trials was used to evaluate the suitability of AmpSeq for molecular correction. We investigated the impact of factors such as the number of amplicon loci analyzed, the informatics criteria used to distinguish genotyping “noise” from real low-density signals, the local epidemiology of malaria transmission, and the potential impact of genetic signals from gametocytes. AmpSeq greatly improved molecular correction and provided accurate drug failure rate estimates. The use of 3 to 5 amplicons was sufficient, and simple, nonstatistical criteria could be used to classify recurrent infections as drug failures or new infections. These results suggest AmpSeq is strongly placed to become the new standard for molecular correction in regulatory trials, with potential extension into routine surveillance once the requisite technical support becomes established.

The molecular basis of antimalarial drug resistance in Plasmodium vivax

Plasmodium vivax is the most geographically widespread cause of human malaria and is responsible for the majority of cases outside of the African continent. While great progress has been made towards eliminating human malaria, drug resistant parasite strains pose a threat towards continued progress. Resistance has arisen to multiple antimalarials in P. vivax, including to chloroquine, which is currently the first line therapy for P. vivax in most regions. Despite its importance, an understanding of the molecular mechanisms of drug resistance in this species remains elusive, in large part due to the complex biology of P. vivax and the lack of in vitro culture. In this review, we will cover the extent and challenges of measuring clinical and in vitro drug resistance in P. vivax. We will consider the roles of candidate drug resistance genes. We will highlight the development of molecular approaches for studying P. vivax biology that provide the opportunity to validate the role of putative drug resistance mutations as well as identify novel mechanisms of drug resistance in this understudied parasite. Validated molecular determinants and markers of drug resistance are essential for the rapid and cost-effective monitoring of drug resistance in P. vivax, and will be useful for optimizing drug regimens and for informing drug policy in control and elimination settings.

•

Drug resistance is emerging in Plasmodium vivax, an important cause of malaria.

•

The complex biology of P. vivax and the limited range of research tools make it difficult to identify drug resistance.

•

The molecular mechanisms of drug resistance in P. vivax remain elusive.

•

This review highlights the extent of drug resistance, the putative mechanisms of resistance and new technologies for the study of P. vivax drug resistance.

Antimalarial Drug Resistance and Implications for the WHO Global Technical Strategy

PURPOSE OF REVIEW

Five years have passed since the World Health Organization released its Global Technical Strategy for Malaria (GTS). In that time, progress against malaria has plateaued. This review focuses on the implications of antimalarial drug resistance for the GTS and how interim progress in parasite genomics and antimalarial pharmacology offer a bulwark against it.

RECENT FINDINGS

For the first time, drug resistance-conferring genes have been identified and validated before their global expansion in malaria parasite populations. More efficient methods for their detection and elaboration have been developed, although low-density infections and polyclonality remain a nuisance to be solved. Clinical trials of alternative regimens for multidrug-resistant malaria have delivered promising results. New agents continue down the development pipeline, while a nascent infrastructure in sub-Saharan Africa for conducting phase I trials and trials of transmission-blocking agents has come to fruition after years of preparation.

SUMMARY

These and other developments can help inform the GTS as the world looks ahead to the next two decades of its implementation. To remain ahead of the threat that drug resistance poses, wider application of genomic-based surveillance and optimization of existing and forthcoming antimalarial drugs are essential.

The effectiveness of malaria camps as part of the Durgama Anchalare Malaria Nirakaran (DAMaN) program in Odisha, India: study protocol for a cluster-assigned quasi-experimental study

The Indian state of Odisha has a longstanding battle with forest malaria. Many remote and rural villages have poor access to health care, a problem that is exacerbated during the rainy season when malaria transmission is at its peak. Approximately 62% of the rural population consists of tribal groups who are among the communities most negatively impacted by malaria. To address the persistently high rates of malaria in these remote regions, the Odisha State Malaria Control Program introduced 'malaria camps' in 2017 where teams of health workers visit villages to educate the population, enhance vector control methods, and perform village-wide screening and treatment. Malaria rates declined statewide, particularly in forested areas, following the introduction of the malaria camps, but the impact of the intervention is yet to be externally evaluated. This study protocol describes a cluster-assigned quasi-experimental stepped-wedge study with a pretest-posttest control group design that evaluates if malaria camps reduce the prevalence of malaria, compared to control villages which receive the usual malaria control interventions (e.g. IRS, ITNs), as detected by PCR.

Influx of diverse, drug resistant and transmissible Plasmodium falciparum into a malaria-free setting in Qatar

BACKGROUND

Successful control programs have impeded local malaria transmission in almost all Gulf Cooperation Council (GCC) countries: Qatar, Bahrain, Kuwait, Oman, the United Arab Emirates (UAE) and Saudi Arabia. Nevertheless, a prodigious influx of imported malaria via migrant workers sustains the threat of local transmission. Here we examine the origin of imported malaria in Qatar, assess genetic diversity and the prevalence of drug resistance genes in imported Plasmodium falciparum, and finally, address the potential for the reintroduction of local transmission.

METHODS

This study examined imported malaria cases reported in Qatar, between 2013 and 2016. We focused on P. falciparum infections and estimated both total parasite and gametocyte density, using qPCR and qRT-PCR, respectively. We also examined ten neutral microsatellites and four genes associated with drug resistance, Pfmrp1, Pfcrt, Pfmdr1, and Pfkelch13, to assess the genetic diversity of imported P. falciparum strains, and the potential for propagating drug resistance genotypes respectively.

RESULTS

The majority of imported malaria cases were P. vivax, while P. falciparum and mixed species infections (P. falciparum / P. vivax) were less frequent. The primary origin of P. vivax infection was the Indian subcontinent, while P. falciparum was mostly presented by African expatriates. Imported P. falciparum strains were highly diverse, carrying multiple genotypes, and infections also presented with early- and late-stage gametocytes. We observed a high prevalence of mutations implicated in drug resistance among these strains, including novel SNPs in Pfkelch13.

CONCLUSIONS

The influx of genetically diverse P. falciparum, with multiple drug resistance markers and a high capacity for gametocyte production, represents a threat for the reestablishment of drug-resistant malaria into GCC countries. This scenario highlights the impact of mass international migration on the reintroduction of malaria to areas with absent or limited local transmission.

Malaria Molecular Epidemiology: An Evolutionary Genetics Perspective

Malaria is a vector-borne disease that involves multiple parasite species in a variety of ecological settings. However, the parasite species causing the disease, the prevalence of subclinical infections, the emergence of drug resistance, the scale-up of interventions, and the ecological factors affecting malaria transmission, among others, are aspects that vary across areas where malaria is endemic. Such complexities have propelled the study of parasite genetic diversity patterns in the context of epidemiologic investigations. Importantly, molecular studies indicate that the time and spatial distribution of malaria cases reflect epidemiologic processes that cannot be fully understood without characterizing the evolutionary forces shaping parasite population genetic patterns. Although broad in scope, this review in the Microbiology Spectrum Curated Collection: Advances in Molecular Epidemiology highlights the need for understanding population genetic concepts when interpreting parasite molecular data. First, we discuss malaria complexity in terms of the parasite species involved. Second, we describe how molecular data are changing our understanding of malaria incidence and infectiousness. Third, we compare different approaches to generate parasite genetic information in the context of epidemiologically relevant questions related to malaria control. Finally, we describe a few Plasmodium genomic studies as evidence of how these approaches will provide new insights into the malaria disease dynamics. *This article is part of a curated collection.

Recent Advances in System Based Study for Anti-Malarial Drug Development Process

Background:

Presently, malaria is one of the most prevalent and deadly infectious disease across Africa, Asia, and America that has now started to spread in Europe. Despite large research being carried out in the field, still, there is a lack of efficient anti-malarial therapeutics. In this paper, we highlight the increasing efforts that are urgently needed towards the development and discovery of potential antimalarial drugs, which must be safe and affordable. The new drugs thus mentioned are also able to counter the spread of malaria parasites that have been resistant to the existing agents.

Objective:

The main objective of the review is to highlight the recent development in the use of system biologybased approaches towards the design and discovery of novel anti-malarial inhibitors.

Method:

A huge literature survey was performed to gain advance knowledge about the global persistence of malaria, its available treatment and shortcomings of the available inhibitors. Literature search and depth analysis were also done to gain insight into the use of system biology in drug discovery and how this approach could be utilized towards the development of the novel anti-malarial drug.

Results:

The system-based analysis has made easy to understand large scale sequencing data, find candidate genes expression during malaria disease progression further design of drug molecules those are complementary of the target proteins in term of shape and configuration.

Conclusion:

The review article focused on the recent computational advances in new generation sequencing, molecular modeling, and docking related to malaria disease and utilization of the modern system and network biology approach to antimalarial potential drug discovery and development.

Mosquitoes as a feasible sentinel group for anti-malarial resistance surveillance by Next Generation Sequencing of Plasmodium falciparum

Background

Plasmodium falciparum drug resistance surveillance is key to successful disease control and eradication. Contemporary methods that only allow determination of prevalence of resistance are expensive, time consuming and require ethical considerations. A newer method involving Next Generation Sequencing (NGS) permits obtaining frequency of resistance while allowing to detect minority variants in mixed infections. Here, NGS was tested for P. falciparum resistance marker detection in mosquito samples as a feasible and suitable alternative for molecular resistance surveillance. Anopheles funestus were collected in southern Mozambique using CDC light traps and manual collections. DNA was extracted from either whole mosquito, head-thorax and abdomen separately or pools of five mosquitoes. These samples were screened for P. falciparum and if positive for k13, pfcrt, pfmdr1, pfdhps and pfdhfr mutations related to anti-malarial drug resistance with Sanger sequencing and NGS.

Results

Among the 846 samples screened for P. falciparum, 122 were positive by 18S ssrDNA qPCR with an infection rate of 23.6%. No mutations were observed for k13 and pfcrt72-76 and almost zero for pfmdr86, but quintuple pfdhfr/pfdhps mutations were near fixation and about half of the isolates contained the pfmdr184F polymorphism. Similar allele frequencies of resistance markers were estimated with NGS in comparison with the prevalence of markers obtained with the gold standard Sanger sequencing.

Conclusions

Pooled deep sequencing of P. falciparum isolates extracted from mosquitoes is a promising, efficient and cost-effective method to quantify allele frequencies at population level which allows to detect known and unknown markers of resistance in single and mixed infections in a timelier manner. Using mosquitoes as sentinel group and focusing on allele frequency opposed to prevalence, permits active surveillance across a more homogeneous geographical range.

Review of applications of high-throughput sequencing in personalized medicine: barriers and facilitators of future progress in research and clinical application

There has been an exponential growth in the performance and output of sequencing technologies (omics data) with full genome sequencing now producing gigabases of reads on a daily basis. These data may hold the promise of personalized medicine, leading to routinely available sequencing tests that can guide patient treatment decisions. In the era of high-throughput sequencing (HTS), computational considerations, data governance and clinical translation are the greatest rate-limiting steps. To ensure that the analysis, management and interpretation of such extensive omics data is exploited to its full potential, key factors, including sample sourcing, technology selection and computational expertise and resources, need to be considered, leading to an integrated set of high-performance tools and systems. This article provides an up-to-date overview of the evolution of HTS and the accompanying tools, infrastructure and data management approaches that are emerging in this space, which, if used within in a multidisciplinary context, may ultimately facilitate the development of personalized medicine.

Deployment and utilization of next-generation sequencing of Plasmodium falciparum to guide anti-malarial drug policy decisions in sub-Saharan Africa: opportunities and challenges

Parasite resistance against anti-malarial drugs is a major threat to the ongoing malaria control and elimination strategies. This is especially true since resistance to the currently recommended artemisinins and partner drugs has been confirmed in South East Asia (SEA) and new anti-malarial compounds are not expected to be available in the near future. Spread from SEA or independent emergence of artemisinin resistance in sub-Saharan Africa (SSA) could reverse the achievements in malaria control that have been attained in the past two decades and derail the ongoing elimination strategies. The current surveillance of clinical efficacy and resistance to anti-malarial drugs is based on efficacy trials to assess the clinical performance of anti-malarials, in vivo/ex vivo assessment of parasite susceptibility to anti-malarials and prevalence of known molecular markers of drug resistance. Whereas clinical efficacy trials are restricted by cost and the complex logistics of patient follow-up, molecular detection of genetic mutations associated with resistance or reduced susceptibility to anti-malarials is by contrast a simple and powerful tool for early detection and monitoring of the prevalence of resistant parasites at population level. This provides needed information before clinical failure emerges, allowing policy makers to anticipate problems and respond. The various methods previously used in detection of molecular markers of drug resistance share some limitations: low-throughput, and high costs per sample and demanding infrastructure. However, recent technological advances including next-generation sequencing (NGS) methodologies promise greatly increased throughput and reduced costs, essentially providing unprecedented potential to address different research and operational questions of relevance for drug policy. This review assesses the potential role of NGS to provide comprehensive information that could guide drug policies in malaria endemic countries and looks at the foreseeable challenges facing the establishment of NGS approaches for routine surveillance of parasite resistance to anti-malarials in SSA.

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.

Molecular assays for antimalarial drug resistance surveillance: A target product profile

Antimalarial drug resistance is a major constraint for malaria control and elimination efforts. Artemisinin-based combination therapy is now the mainstay for malaria treatment. However, delayed parasite clearance following treatment with artemisinin derivatives has now spread in the Greater Mekong Sub region and may emerge or spread to other malaria endemic regions. This spread is of great concern for malaria control programmes, as no alternatives to artemisinin-based combination therapies are expected to be available in the near future. There is a need to strengthen surveillance systems for early detection and response to the antimalarial drug resistance threat. Current surveillance is mainly done through therapeutic efficacy studies; however these studies are complex and both time- and resource-intensive. For multiple common antimalarials, parasite drug resistance has been correlated with specific genetic mutations, and the molecular markers associated with antimalarial drug resistance offer a simple and powerful tool to monitor the emergence and spread of resistant parasites. Different techniques to analyse molecular markers associated with antimalarial drug resistance are available, each with advantages and disadvantages. However, procedures are not adequately harmonized to facilitate comparisons between sites. Here we describe the target product profiles for tests to analyse molecular markers associated with antimalarial drug resistance, discuss how use of current techniques can be standardised, and identify the requirements for an ideal product that would allow malaria endemic countries to provide useful spatial and temporal information on the spread of resistance.

Microbial sequence typing in the genomic era

Next-generation sequencing (NGS), also known as high-throughput sequencing, is changing the field of microbial genomics research. NGS allows for a more comprehensive analysis of the diversity, structure and composition of microbial genes and genomes compared to the traditional automated Sanger capillary sequencing at a lower cost. NGS strategies have expanded the versatility of standard and widely used typing approaches based on nucleotide variation in several hundred DNA sequences and a few gene fragments (MLST, MLVA, rMLST and cgMLST). NGS can now accommodate variation in thousands or millions of sequences from selected amplicons to full genomes (WGS, NGMLST and HiMLST). To extract signals from high-dimensional NGS data and make valid statistical inferences, novel analytic and statistical techniques are needed. In this review, we describe standard and new approaches for microbial sequence typing at gene and genome levels and guidelines for subsequent analysis, including methods and computational frameworks. We also present several applications of these approaches to some disciplines, namely genotyping, phylogenetics and molecular epidemiology.

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.

Next-Generation Sequencing and Bioinformatics Protocol for Malaria Drug Resistance Marker Surveillance

The recent advances in next-generation sequencing technologies provide a new and effective way of tracking malaria drug-resistant parasites. To take advantage of this technology, an end-to-end Illumina targeted amplicon deep sequencing (TADS) and bioinformatics pipeline for molecular surveillance of drug resistance in P. falciparum, called malaria resistance surveillance (MaRS), was developed. TADS relies on PCR enriching genomic regions, specifically target genes of interest, prior to deep sequencing. MaRS enables researchers to simultaneously collect data on allele frequencies of multiple full-length P. falciparum drug resistance genes (crt, mdr1, k13, dhfr, dhps, and the cytochrome b gene), as well as the mitochondrial genome. Information is captured at the individual patient level for both known and potential new single nucleotide polymorphisms associated with drug resistance. The MaRS pipeline was validated using 245 imported malaria cases that were reported to the Centers for Disease Control and Prevention (CDC). The chloroquine resistance crt CVIET genotype (mutations underlined) was observed in 42% of samples, the highly pyrimethamine-resistant dhpsIRN triple mutant in 92% of samples, and the sulfadoxine resistance dhps mutation SGEAA in 26% of samples. The mdr1 NFSND genotype was found in 40% of samples. With the exception of two cases imported from Cambodia, no artemisinin resistance k13 alleles were identified, and 99% of patients carried parasites susceptible to atovaquone-proguanil. Our goal is to implement MaRS at the CDC for routine surveillance of imported malaria cases in the United States and to aid in the adoption of this system at participating state public health laboratories, as well as by global partners.

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.

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

Background

Amplicon deep sequencing permits sensitive detection of minority clones and improves discriminatory power for genotyping multi-clone Plasmodium falciparum infections. New amplicon sequencing and data analysis protocols are needed for genotyping in epidemiological studies and drug efficacy trials of P. falciparum.

Methods

Targeted sequencing of molecular marker csp and novel marker cpmp was conducted in duplicate on mixtures of parasite culture strains and 37 field samples. A protocol allowing to multiplex up to 384 samples in a single sequencing run was applied. Software “HaplotypR” was developed for data analysis.

Results

Cpmp was highly diverse (H_e = 0.96) in contrast to csp (H_e = 0.57). Minority clones were robustly detected if their frequency was >1%. False haplotype calls owing to sequencing errors were observed below that threshold.

Conclusions

To reliably detect haplotypes at very low frequencies, experiments are best performed in duplicate and should aim for coverage of >10′000 reads/amplicon. When compared to length polymorphic marker msp2, highly multiplexed amplicon sequencing displayed greater sensitivity in detecting minority clones.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4260-y) contains supplementary material, which is available to authorized users.

malERA: An updated research agenda for diagnostics, drugs, vaccines, and vector control in malaria elimination and eradication

Since the turn of the century, a remarkable expansion has been achieved in the range and effectiveness of products and strategies available to prevent, treat, and control malaria, including advances in diagnostics, drugs, vaccines, and vector control. These advances have once again put malaria elimination on the agenda. However, it is clear that even with the means available today, malaria control and elimination pose a formidable challenge in many settings. Thus, currently available resources must be used more effectively, and new products and approaches likely to achieve these goals must be developed. This paper considers tools (both those available and others that may be required) to achieve and maintain malaria elimination. New diagnostics are needed to direct treatment and detect transmission potential; new drugs and vaccines to overcome existing resistance and protect against clinical and severe disease, as well as block transmission and prevent relapses; and new vector control measures to overcome insecticide resistance and more powerfully interrupt transmission. It is also essential that strategies for combining new and existing approaches are developed for different settings to maximise their longevity and effectiveness in areas with continuing transmission and receptivity. For areas where local elimination has been recently achieved, understanding which measures are needed to maintain elimination is necessary to prevent rebound and the reestablishment of transmission. This becomes increasingly important as more countries move towards elimination.

Malaria

Malaria is caused in humans by five species of single-celled eukaryotic Plasmodium parasites (mainly Plasmodium falciparum and Plasmodium vivax) that are transmitted by the bite of Anopheles spp. mosquitoes. Malaria remains one of the most serious infectious diseases; it threatens nearly half of the world's population and led to hundreds of thousands of deaths in 2015, predominantly among children in Africa. Malaria is managed through a combination of vector control approaches (such as insecticide spraying and the use of insecticide-treated bed nets) and drugs for both treatment and prevention. The widespread use of artemisinin-based combination therapies has contributed to substantial declines in the number of malaria-related deaths; however, the emergence of drug resistance threatens to reverse this progress. Advances in our understanding of the underlying molecular basis of pathogenesis have fuelled the development of new diagnostics, drugs and insecticides. Several new combination therapies are in clinical development that have efficacy against drug-resistant parasites and the potential to be used in single-dose regimens to improve compliance. This ambitious programme to eliminate malaria also includes new approaches that could yield malaria vaccines or novel vector control strategies. However, despite these achievements, a well-coordinated global effort on multiple fronts is needed if malaria elimination is to be achieved.

High throughput resistance profiling of Plasmodium falciparum infections based on custom dual indexing and Illumina next generation sequencing-technology

Genetic polymorphisms in P. falciparum can be used to indicate the parasite's susceptibility to antimalarial drugs as well as its geographical origin. Both of these factors are key to monitoring development and spread of antimalarial drug resistance. In this study, we combine multiplex PCR, custom designed dual indexing and Miseq sequencing for high throughput SNP-profiling of 457 malaria infections from Guinea-Bissau, at the cost of 10 USD per sample. By amplifying and sequencing 15 genetic fragments, we cover 20 resistance-conferring SNPs occurring in pfcrt, pfmdr1, pfdhfr, pfdhps, as well as the entire length of pfK13, and the mitochondrial barcode for parasite origin. SNPs of interest were sequenced with an average depth of 2,043 reads, and bases were called for the various SNP-positions with a p-value below 0.05, for 89.8-100% of samples. The SNP data indicates that artemisinin resistance-conferring SNPs in pfK13 are absent from the studied area of Guinea-Bissau, while the pfmdr1 86 N allele is found at a high prevalence. The mitochondrial barcodes are unanimous and accommodate a West African origin of the parasites. With this method, very reliable high throughput surveillance of antimalarial drug resistance becomes more affordable than ever before.

Molecular Epidemiology of Plasmodium falciparum kelch13 Mutations in Senegal Determined by Using Targeted Amplicon Deep Sequencing

The emergence of Plasmodium falciparum resistance to artemisinin in Southeast Asia threatens malaria control and elimination activities worldwide. Multiple polymorphisms in the P. falciparum kelch gene found in chromosome 13 (Pfk13) have been associated with artemisinin resistance. Surveillance of potential drug resistance loci within a population that may emerge under increasing drug pressure is an important public health activity. In this context, P. falciparum infections from an observational surveillance study in Senegal were genotyped using targeted amplicon deep sequencing (TADS) for Pfk13 polymorphisms. The results were compared to previously reported Pfk13 polymorphisms from around the world. A total of 22 Pfk13 propeller domain polymorphisms were identified in this study, of which 12 have previously not been reported. Interestingly, of the 10 polymorphisms identified in the present study that were also previously reported, all had a different amino acid substitution at these codon positions. Most of the polymorphisms were present at low frequencies and were confined to single isolates, suggesting they are likely transient polymorphisms that are part of naturally evolving parasite populations. The results of this study underscore the need to identify potential drug resistance loci existing within a population, which may emerge under increasing drug pressure.

Overlap Extension Barcoding for the Next Generation Sequencing and Genotyping of Plasmodium falciparum in Individual Patients in Western Kenya

Large-scale molecular epidemiologic studies of Plasmodium falciparum parasites have provided insights into parasite biology and transmission, can identify the spread of drug resistance, and are useful in assessing vaccine targets. The polyclonal nature infections in high transmission settings is problematic for traditional genotyping approaches. Next-generation sequencing (NGS) approaches to parasite genotyping allow sensitive detection of minority variants, disaggregation of complex parasite mixtures, and scalable processing of large samples sets. Therefore, we designed, validated, and applied to field parasites an approach that leverages sequencing of individually barcoded samples in a multiplex manner. We utilize variant barcodes, invariant linker sequences and modular template-specific primers to allow for the simultaneous generation of high-dimensional sequencing data of multiple gene targets. This modularity permits a cost-effective and reproducible way to query many genes at once. In mixtures of reference parasite genomes, we quantitatively detected unique haplotypes comprising as little as 2% of a polyclonal infection. We applied this genotyping approach to field-collected parasites collected in Western Kenya in order to simultaneously obtain parasites genotypes at three unlinked loci. In summary, we present a rapid, scalable, and flexible method for genotyping individual parasites that enables molecular epidemiologic studies of parasite evolution, population structure and transmission.

Dissecting malaria biology and epidemiology using population genetics and genomics

Molecular approaches have an increasingly recognized utility in surveillance of malaria parasite populations, not only in defining prevalence and incidence with higher sensitivity than traditional methods, but also in monitoring local and regional parasite transmission patterns. In this review, we provide an overview of population genetic and genomic studies of human-infecting Plasmodium species, highlighting recent advances in the field. In accordance with the renewed impetus for malaria eradication, many studies are now using genetic and genomic epidemiology to support local evidence-based intervention strategies. Microsatellite genotyping remains a popular approach for both Plasmodium falciparum and Plasmodium vivax. However, with the increasing availability of whole genome sequencing data enabling effective single nucleotide polymorphism-based panels tailored to a given study question and setting, this approach is gaining popularity. The availability of new reference genomes for Plasmodium malariae and Plasmodium ovale should see a surge in similar molecular studies on these currently neglected species. Genomic studies are revealing new insights into important adaptive mechanisms of the parasite including antimalarial drug resistance. The advent of new methodologies such as selective whole genome amplification for dealing with extensive human DNA in low density field isolates should see genome-wide approaches becoming routine for parasite surveillance once the economic costs outweigh the current cost benefits of targeted approaches.