Rapid profiling of Plasmodium parasites from genome sequences to assist malaria control.

BACKGROUND: Malaria continues to be a major threat to global public health. Whole genome sequencing (WGS) of the underlying Plasmodium parasites has provided insights into the genomic epidemiology of malaria. Genome sequencing is rapidly gaining traction as a diagnostic and surveillance tool for clinical settings, where the profiling of co-infections, identification of imported malaria parasites, and detection of drug resistance are crucial for infection control and disease elimination. To support this informatically, we have developed the Malaria-Profiler tool, which rapidly (within minutes) predicts Plasmodium species, geographical source, and resistance to antimalarial drugs directly from WGS data. RESULTS: The online and command line versions of Malaria-Profiler detect ~ 250 markers from genome sequences covering Plasmodium speciation, likely geographical source, and resistance to chloroquine, sulfadoxine-pyrimethamine (SP), and other anti-malarial drugs for P. falciparum, but also providing mutations for orthologous resistance genes in other species. The predictive performance of the mutation library was assessed using 9321 clinical isolates with WGS and geographical data, with most being single-species infections (P. falciparum 7152/7462, P. vivax 1502/1661, P. knowlesi 143/151, P. malariae 18/18, P. ovale ssp. 5/5), but co-infections were identified (456/9321; 4.8%). The accuracy of the predicted geographical profiles was high to both continental (96.1%) and regional levels (94.6%). For P. falciparum, markers were identified for resistance to chloroquine (49.2%; regional range: 24.5% to 100%), sulfadoxine (83.3%; 35.4- 90.5%), pyrimethamine (85.4%; 80.0-100%) and combined SP (77.4%). Markers associated with the partial resistance of artemisinin were found in WGS from isolates sourced from Southeast Asia (30.6%). CONCLUSIONS: Malaria-Profiler is a user-friendly tool that can rapidly and accurately predict the geographical regional source and anti-malarial drug resistance profiles across large numbers of samples with WGS data. The software is flexible with modifiable bioinformatic pipelines. For example, it is possible to select the sequencing platform, display specific variants, and customise the format of outputs. With the increasing application of next-generation sequencing platforms on Plasmodium DNA, Malaria-Profiler has the potential to be integrated into point-of-care and surveillance settings, thereby assisting malaria control. Malaria-Profiler is available online (bioinformatics.lshtm.ac.uk/malaria-profiler) and as standalone software ( https://github.com/jodyphelan/malaria-profiler ).

Anopheles diversity, biting behaviour and transmission potential in forest and farm environments of Gia Lai province, Vietnam.

BACKGROUND: Despite recent reductions in Vietnam, malaria transmission persists in some areas in forests and farmlands where a high density of Anopheles mosquitoes relative to other environments occurs. To inform effective malaria control measures, it is important to understand vector bionomics and the malaria transmission role of Anopheles spp. in the highland regions of Vietnam. This study was conducted to quantify the abundance, composition and biting behaviour of the Anopheles mosquito population, and the proportion of Plasmodium spp. infected mosquitoes collected from forest and agricultural farm sites in Gia Lai province, Vietnam. METHODS: Forest and agricultural farm sites in Gia Lai province were selected for mosquito collections (total eight sites). Mosquito collection was performed by Human-baited Double Net Trap (HDNT), animal-baited traps (ABT) using cattle, and CDC light traps. Captured mosquitoes were identified morphologically, and salivary glands of Anopheles mosquitoes were examined for sporozoites using microscopy. Plasmodium infection was determined by Polymerase Chain Reaction (PCR), and identification of blood meal type was determined by PCR and diffuse serum agglutination assay. RESULTS: A total of 1815 Anopheles mosquitoes belonging to 19 species were collected by ABT (n = 1169), HDNT (n = 471) and CDC light trap (n = 175). Anopheles abundance and diversity varied by district and environment. Capture by HDNT of Anopheles of vectorial concern was observed between early evening and early morning. Plasmodium vivax infection was determined by PCR in two Anopheles dirus specimens captured by HDNT in forest sites. Blood from a range of hosts could, including human blood, could be detected in species considered primary and secondary vectors An. dirus, and Anopheles aconitus, and Anopheles maculatus, respectively. CONCLUSIONS: A low number of Anopheles spp. considered primary vectors of concern and very low numbers of Plasmodium spp. infected Anopheles mosquitoes were captured at the end of the rainy season in the Central Highlands of Vietnam. However, capture species of vectorial concern by HDNT throughout the early to late evening demonstrates that use of additional personal protective measures could supplement current preventative measures, such as bed nets to prevent exposure to vectors of concern in this region.

Surveillance of pfhrp2 and pfhrp3 gene deletions among symptomatic Plasmodium falciparum malaria patients in Central Vietnam.

BACKGROUND: Malaria rapid diagnostic tests (RDTs) remain the main point-of-care tests for diagnosis of symptomatic Plasmodium falciparum malaria in endemic areas. However, parasites with gene deletions in the most common RDT target, histidine rich protein 2 (pfhrp2/HRP2), can produce false-negative RDT results leading to inadequate case management. The objective of this study was to determine the prevalence of hrp2/3 deletions causing false-negative RDT results in Vietnam (Gia Lai and Dak Lak provinces). METHODS: Individuals presenting with malaria symptoms at health facilities were screened for P. falciparum infection using light microscopy and HRP2-RDT (SD Bioline Malaria Antigen Pf/Pv RDT, Abbott). Microscopically confirmed P. falciparum infections were analysed for parasite species by 18S rRNA qPCR, and pfhrp2 and pfhrp3 exon2 deletions were investigated by nested PCR. pfhrp2 amplicons were sequenced by the Sanger method and HRP2 plasma levels were determined by enzyme-linked immunosorbent assay (ELISA). RESULTS: The prevalence of false-negative RDT results among symptomatic cases was 5.6% (15/270). No pfhrp2 and pfhrp3 deletions were identified. False-negative RDT results were associated with lower parasite density (p = 0.005) and lower HRP2 plasma concentrations (p < 0.001), as compared to positive RDT. CONCLUSIONS: The absence of hrp2/3 deletions detected in this survey suggests that HRP2-based malaria RDTs remain effective for the diagnosis of symptomatic P. falciparum malaria in Central Vietnam.

High Proportion of Genome-Wide Homology and Increased Pretreatment pvcrt Levels in Plasmodium vivax Late Recurrences: a Chloroquine Therapeutic Efficacy Study.

Chloroquine (CQ) is the first-line treatment for Plasmodium vivax malaria in most countries where malaria is endemic. Monitoring P. vivax CQ resistance (CQR) is critical but remains challenged by the difficulty to distinguish real treatment failure from reinfection or liver relapse. The therapeutic efficacy of CQ against uncomplicated P. vivax malaria was evaluated in Gia Lai Province, Vietnam. Sixty-seven patients were enrolled and followed for 42 days using microscopy and quantitative PCR. Adequate clinical and parasitological response (ACPR) was 100% (66/66) on day 28 but 75.4% (49/65) on day 42. Eighteen recurrences (27.7%) were detected, with a median time to recurrence of 42 days (interquartile range [IQR], 35 to 42) and blood CQ concentration of <100 ng/ml. Primary infections leading to recurrence occurred in younger individuals (median age for ACPR = 25 years [IQR, 20 to 28]; recurrences = 18 [16 to 21]; P = 0.002) had a longer parasite clearance time (PCT for ACPR = 47.5 h [IQR, 36.2 to 59.8 h]; recurrences = 54.2 [48.4 to 62.0]; P = 0.035) and higher pvcrt gene expression (median relative expression ratio for ACPR = 0.09 [IQR, 0.05 to 0.22]; recurrences = 0.20 [0.15 to 0.56]; P = 0.002), but showed no differences in ex vivo CQ sensitivity. Parasite genotyping by microsatellites, single nucleotide polymorphism (SNP) barcoding, and whole-genome sequencing (WGS) identified a majority of homologous recurrences, with 80% (8/10) showing >98% identity by descent to paired day 0 samples. This study shows that CQ remained largely efficacious to treat P. vivax in Gia Lai; i.e., recurrences occurred late (>day 28) and in the presence of low blood CQ concentrations. However, the combination of both WGS and gene expression analysis (pvcrt) data with clinical data (PCT) allowed us to identify potential emergence of low-grade CQR, which should be closely monitored. (This study has been registered at ClinicalTrials.gov under identifier NCT02610686.).

Efficacy of dihydroartemisinin/piperaquine and artesunate monotherapy for the treatment of uncomplicated Plasmodium falciparum malaria in Central Vietnam.

BACKGROUND: Artemisinin-based combination therapies (ACTs) have significantly contributed to reduce Plasmodium falciparum malaria burden in Vietnam, but their efficacy is challenged by treatment failure of dihydroartemisinin/piperaquine ACT in Southern provinces. OBJECTIVES: To assess the efficacy of dihydroartemisinin/piperaquine for uncomplicated P. falciparum malaria in Gia Lai, Central Vietnam, and determine parasite resistance to artemisinin (ClinicalTrials.gov identifier NCT02604966). METHODS: Sixty patients received either dihydroartemisinin/piperaquine (4 mg/kg/day, 3 days; n = 33) or artesunate monotherapy (4 mg/kg/day, 3 days; n = 27) followed by dihydroartemisinin/piperaquine (AS + DHA/PPQ). Clinical phenotypes were determined during a 42 day follow-up and analysed together with ex vivo susceptibility to antimalarials and molecular markers of drug resistance. RESULTS: Day 3 positivity rate was significantly higher in the AS + DHA/PPQ arm compared with dihydroartemisinin/piperaquine (70.4% versus 39.4%, P = 0.016). Parasite clearance time was 95.2 h (AS + DHA/PPQ) versus 71.9 h (dihydroartemisinin/piperaquine, P = 0.063) and parasite clearance half-life was 7.4 h (AS + DHA/PPQ) versus 7.0 h (dihydroartemisinin/piperaquine, P = 0.140). Adequate clinical and parasitological response at Day 42 was 100% in both arms. By RT-qPCR, 36% (19/53) patients remained positive until Day 7. No recurrences were detected. kelch13 artemisinin resistance mutations were found in 87% (39/45) of isolates and 50% (20/40) were KEL1/C580Y. The piperaquine resistance marker plasmepsin-2 was duplicated in 10.4% (5/48). Isolates from Day 3-positive patients (n = 18) had higher ex vivo survival rates to artemisinin compounds (P < 0.048) and prevalence of kelch13 mutations (P = 0.005) than Day 3-negative patients (n = 5). The WHO definition of artemisinin resistance was fulfilled in 60% (24/40) of cases. CONCLUSIONS: Although dihydroartemisinin/piperaquine remained effective to treat P. falciparum, the high Day 3 positivity rate and prevalence of KEL1 strains calls for continuous monitoring of dihydroartemisinin/piperaquine efficacy in Central Vietnam.

Transcriptome Sequencing and Analysis of Changes Associated with Insecticide Resistance in the Dengue Mosquito (Aedes aegypti) in Vietnam.

The mosquito Aedes aegypti is a transmission vector for dangerous epidemic diseases in humans. Insecticides have been used as the most general vector control method in the world. However, Ae. aegypti have developed many resistant mechanisms such as reduced neuronal sensitivity to insecticides (target-site resistance), enhanced insecticide metabolism (metabolic resistance), altered transport, sequestration, and other mechanisms. It has become a major problem for vector control programs. Transcriptome sequencing and bioinformatic analysis were used to compare transcription levels between a susceptible strain (Bora7) and a resistant strain (KhanhHoa7) collected from the field. A total of 161 million Illumina reads, including 66,076,678 reads from the Bora7 strain and 69,606,654 reads from the KhanhHoa7 strain, were generated and assembled into 11,174 genes. A comparison of the KhanhHoa7 transcriptome to that of Bora7 showed 672 upregulated genes and 488 downregulated genes. We identified the highly upregulated genes: cytochrome P450 4C1, 4C3, 4C21, 4D1, 4D1 isoform X2, 4D2, 4D2 isoform X2, 4G15, 6A2, 6A8, 6D3, and 9E2; Glutathione S transferase (GST1), UGT1-3, 1-7, 2B15, and 2B37; binding cassette transporter (ABC) transporter F family member 4 and ABC transporter G family member 20. Interestingly, there was a significant increase in the expression of the genes such as CYP9E2 (8.3-fold), CYP6A8 (5.9-fold), CYP6D3 (5.4-fold), CYP4C21 (5.4-fold), CYP4G15 (5.2-fold), GST1 (3.5-fold), and ABC transporter 4 (2.1-fold). Our results suggested a potential relationship between the expression of the genes in metabolic processes and insecticide resistance in the studied strain. These results may contribute to the understanding of the mechanisms of insecticide resistance in Ae. aegypti.