Climate, demography, immunology, and virology combine to drive two decades of dengue virus dynamics in Cambodia.

The incidence of dengue virus disease has increased globally across the past half-century, with highest number of cases ever reported in 2019 and again in 2023. We analyzed climatological, epidemiological, and phylogenomic data to investigate drivers of two decades of dengue in Cambodia, an understudied endemic setting. Using epidemiological models fit to a 19-y dataset, we first demonstrate that climate-driven transmission alone is insufficient to explain three epidemics across the time series. We then use wavelet decomposition to highlight enhanced annual and multiannual synchronicity in dengue cycles between provinces in epidemic years, suggesting a role for climate in homogenizing dynamics across space and time. Assuming reported cases correspond to symptomatic secondary infections, we next use an age-structured catalytic model to estimate a declining force of infection for dengue through time, which elevates the mean age of reported cases in Cambodia. Reported cases in >70-y-old individuals in the 2019 epidemic are best explained when also allowing for waning multitypic immunity and repeat symptomatic infections in older patients. We support this work with phylogenetic analysis of 192 dengue virus (DENV) genomes that we sequenced between 2019 and 2022, which document emergence of DENV-2 Cosmopolitan Genotype-II into Cambodia. This lineage demonstrates phylogenetic homogeneity across wide geographic areas, consistent with invasion behavior and in contrast to high phylogenetic diversity exhibited by endemic DENV-1. Finally, we simulate an age-structured, mechanistic model of dengue dynamics to demonstrate how expansion of an antigenically distinct lineage that evades preexisting multitypic immunity effectively reproduces the older-age infections witnessed in our data.

Discovering disease-causing pathogens in resource-scarce Southeast Asia using a global metagenomic pathogen monitoring system.

SignificanceMetagenomic pathogen sequencing offers an unbiased approach to characterizing febrile illness. In resource-scarce settings with high biodiversity, it is critical to identify disease-causing pathogens in order to understand burden and to prioritize efforts for control. Here, metagenomic next-generation sequencing (mNGS) characterization of the pathogen landscape in Cambodia revealed diverse vector-borne and zoonotic pathogens irrespective of age and gender as risk factors. Identification of key pathogens led to changes in national program surveillance. This study is a "real world" example of the use of mNGS surveillance of febrile individuals, executed in-country, to identify outbreaks of vector-borne, zoonotic, and other emerging pathogens in a resource-scarce setting.

Immunity to Non-Dengue Flaviviruses Impacts Dengue Virus Immunoglobulin G Enzyme-Linked Immunosorbent Assay Specificity in Cambodia.

BACKGROUND: Seroprevalence studies are the standard for disease surveillance, and serology determined eligibility for the first dengue vaccine. Expanding flavivirus co-circulation and vaccination complicate testing. We evaluate the accuracy of a common dengue virus serological assay, examine immunity to non-dengue flaviviruses as a contributor to decreased performance, and assess whether alternative cut points may improve assay performance. METHODS: Children (n = 770) aged 2-9 years in Kampong Speu, Cambodia were enrolled in a prospective longitudinal study, and PanBio indirect dengue virus immunoglobulin G (IgG) enzyme-linked immunosorbent assay (ELISA) was performed. Plaque reduction neutralization tests (PRNTs) using dengue viruses were performed on a subset to assess the accuracy of the IgG ELISA, and PRNTs with Zika, Japanese encephalitis, and West Nile viruses evaluated immunity to non-dengue flaviviruses. Receiver operating curve analysis identified an alternative cut point to improve IgG ELISA accuracy. RESULTS: The dengue IgG ELISA had a lower specificity than previously reported (58% vs 93%-100%). Of those with false-positive IgG results, 46% had detectable neutralizing antibodies against other flaviviruses including 14% against West Nile virus. A higher IgG cut point improved the test accuracy in this population. CONCLUSIONS: Physicians and public health authorities should be alert for West Nile in Cambodia. Immunity to non-dengue flaviviruses can impact dengue surveillance. CLINICAL TRIALS REGISTRATION: NCT03534245.

Bacterial Genomics for National Antimicrobial Resistance Surveillance in Cambodia.

BACKGROUND: Antimicrobial resistance (AMR) surveillance in low- and middle-income countries (LMICs) often relies on poorly resourced laboratory processes. Centralized sequencing was combined with cloud-based, open-source bioinformatics solutions for national AMR surveillance in Cambodia. METHODS: Blood cultures growing gram-negative bacteria were collected at six Cambodian hospitals (January 2021 - October 2022). Isolates were obtained from pure plate growth and shotgun DNA sequencing performed in-country. Using public nucleotide and protein databases, reads were aligned for pathogen identification and AMR gene characterization. Multilocus sequence typing was performed on whole genome assemblies and haplotype clusters compared against published genomes. FINDINGS: Genes associated with acquired resistance to fluoroquinolones were identified in 59%, TMP/SMX in 45%, and aminoglycosides in 52% of 715 isolates. Extended-spectrum beta-lactamase encoding genes were identified in 34% isolates, most commonly blaCTX-M-15, blaCTX-M-27, and blaCTX-M-55 in E. coli sequence types 131 and 1193. Carbapenemase genes were identified in 12% isolates, most commonly blaOXA-23, blaNDM-1, blaOXA-58 and blaOXA-66 in Acinetobacter species. Phylogenetic analysis revealed clonal strains of A. baumannii, representing suspected nosocomial outbreaks, and genetic clusters of quinolone-resistant typhoidal Salmonella and ESBL E. coli cases suggesting community transmission. INTERPRETATION: With accessible sequencing platforms and bioinformatics solutions, bacterial genomics can supplement AMR surveillance in LMICs. FUNDING: Research was supported by the National Institute of Allergy and Infectious Diseases and the Bill and Melinda Gates Foundation [OPP1211806].

National dengue surveillance, Cambodia 2002-2020.

Global dengue incidence has increased dramatically over the past few decades from approximately 500 000 reported cases in 2000 to over 5 million in 2019. This trend has been attributed to population growth in endemic areas, rapid unplanned urbanization, increasing global connectivity, and climate change expanding the geographic range of the Aedes spp. mosquito, among other factors. Reporting dengue surveillance data is key to understanding the scale of the problem, identifying important changes in the landscape of disease, and developing policies for clinical management, vector control and vaccine rollout. However, surveillance practices are not standardized, and data may be difficult to interpret particularly in low- and middle-income countries with fragmented health-care systems. The latest national dengue surveillance data for Cambodia was published in 2010. Since its publication, the country experienced marked changes in health policies, population demographics, climate and urbanization. How these changes affected dengue control remains unknown. In this article, we summarize two decades of policy changes, published literature, country statistics, and dengue case data collected by the Cambodia National Dengue Control Programme to: (i) identify important changes in the disease landscape; and (ii) derive lessons to inform future surveillance and disease control strategies. We report that while dengue case morbidity and mortality rates in Cambodia fell between 2002 and 2020, dengue incidence doubled and age at infection increased. Future national surveillance, disease prevention and treatment, and vector control policies will have to account for these changes to optimize disease control.

Le taux d'incidence de la dengue dans le monde a considérablement augmenté au cours des dernières décennies, passant d'environ 500 000 cas notifiés en 2000 à plus de 5 millions en 2019. Cette tendance est attribuée à la croissance démographique dans les zones d'endémie, à l'urbanisation rapide non planifiée, au développement de la connectivité à l'échelle internationale, ainsi qu'au changement climatique, qui agrandit le territoire géographique du moustique Aedes spp., entre autres. La communication des données de surveillance de la dengue est essentielle pour comprendre l'étendue du problème, identifier les principales variations de contexte entourant la maladie et mettre au point des politiques pour la prise en charge clinique, la lutte contre les vecteurs et le déploiement des vaccins. Les pratiques en matière de surveillance ne sont toutefois pas standardisées et les données peuvent être difficiles à interpréter, surtout dans les pays à revenu faible et intermédiaire où les systèmes de soins de santé sont fragmentés. Les données de surveillance les plus récentes concernant la dengue au Cambodge ont été publiées en 2010. Depuis leur publication, le pays a subi de profondes mutations au niveau des politiques de santé, de l'évolution démographique, du climat et de l'urbanisation. L'impact de ces mutations sur la lutte contre la dengue reste à établir. Dans le présent article, nous résumons deux décennies d'amendements politiques, de documentation, de statistiques nationales et d'informations collectées sur les cas par le programme cambodgien de lutte contre la dengue afin de: (i) définir les changements importants survenus dans le contexte entourant la maladie; mais aussi (ii) tirer des leçons en vue d'élaborer, à l'avenir, des stratégies de surveillance et de lutte contre la maladie. Nous signalons qu'en dépit d'une baisse des taux de morbidité et de mortalité liés aux cas de dengue entre 2002 et 2020 au Cambodge, son incidence a doublé et l'âge des patients au moment de l'infection a augmenté. Les futures politiques nationales de surveillance, de prévention et de traitement de la dengue, mais aussi de lutte contre ses vecteurs, devront tenir compte de ces changements de façon à mieux maîtriser la maladie.

La incidencia del dengue a nivel mundial ha aumentado considerablemente en las últimas décadas, desde aproximadamente 500 000 casos notificados en el año 2000 a más de 5 millones en 2019. Esta tendencia se ha atribuido al crecimiento de la población en zonas endémicas, a una urbanización rápida y no planificada, al aumento de la conectividad a nivel mundial y al cambio climático, que está permitiendo una distribución geográfica más amplia del mosquito Aedes spp., entre otros factores. Para comprender la magnitud del problema resulta clave la notificación de datos sobre vigilancia del dengue, la identificación de cambios importantes dentro del escenario de la enfermedad, la creación de políticas enfocadas a la gestión clínica, así como el control de vectores y la implantación de la vacuna. Sin embargo, las prácticas sobre vigilancia no están estandarizadas y es posible que sea difícil interpretar los datos, especialmente en países con ingresos medios y bajos, que cuentan con sistemas fragmentados de atención sanitaria. Los datos nacionales más recientes sobre vigilancia del dengue en Camboya se publicaron en 2010. Desde su publicación, el país experimentó cambios significativos en las políticas sanitarias, la demografía de la población, el clima y la urbanización. Aún no se sabe cómo afectaron dichos cambios al control del dengue. En el presente artículo, resumimos dos décadas de cambios políticos, de bibliografía publicada, de datos estadísticos a nivel nacional y datos sobre casos de dengue recopilados por el programa nacional de control de dengue en Camboya, con el fin de: (i) identificar cambios importantes en el escenario de la enfermedad; y (ii) extraer conclusiones para orientar futuras estrategias sobre vigilancia y control de la enfermedad. Informamos de que, aunque las tasas de morbilidad y mortalidad de los casos de dengue en Camboya descendieron entre 2002 y 2020, la incidencia del dengue se duplicó y la edad de infección aumentó. Las futuras políticas nacionales sobre vigilancia, prevención y tratamiento de la enfermedad y control de vectores deberán tener en cuenta estos cambios para optimizar el control de la enfermedad.

Development of Inapparent Dengue Associated With Increased Antibody Levels to Aedes aegypti Salivary Proteins: A Longitudinal Dengue Cohort in Cambodia.

BACKGROUND: We established the first prospective cohort to understand how infection with dengue virus is influenced by vector-specific determinants such as humoral immunity to Aedes aegypti salivary proteins. METHODS: Children aged 2-9 years were enrolled in the PAGODAS (Pediatric Assessment Group of Dengue and Aedes Saliva) cohort with informed consent by their guardians. Children were followed semi-annually for antibodies to dengue and to proteins in Ae. aegypti salivary gland homogenate using enzyme-linked immunosorbent assays and dengue-specific neutralization titers. Children presented with fever at any time for dengue testing. RESULTS: From 13 July to 30 August 2018, we enrolled 771 children. At baseline, 22% (173/770) had evidence of neutralizing antibodies to 1 or more dengue serotypes. By April 2020, 51 children had symptomatic dengue while 148 dengue-naive children had inapparent dengue defined by neutralization assays. In a multivariate model, individuals with higher antibodies to Ae. aegypti salivary proteins were 1.5 times more likely to have dengue infection (hazard ratio [HR], 1.47 [95% confidence interval {CI}, 1.05-2.06]; P = .02), particularly individuals with inapparent dengue (HR, 1.64 [95% CI, 1.12-2.41]; P = .01). CONCLUSIONS: High levels of seropositivity to Ae. aegypti salivary proteins are associated with future development of dengue infection, primarily inapparent, in dengue-naive Cambodian children. CLINICAL TRIALS REGISTRATION: NCT03534245.

SARS-CoV-2 Cross-Reactivity in Prepandemic Serum from Rural Malaria-Infected Persons, Cambodia.

Inhabitants of the Greater Mekong Subregion in Cambodia are exposed to pathogens that might influence serologic cross-reactivity with severe acute respiratory syndrome coronavirus 2. A prepandemic serosurvey of 528 malaria-infected persons demonstrated higher-than-expected positivity of nonneutralizing IgG to spike and receptor-binding domain antigens. These findings could affect interpretation of large-scale serosurveys.

Single-Dose Tafenoquine to Prevent Relapse of Plasmodium vivax Malaria.

BACKGROUND: Treatment of Plasmodium vivax malaria requires the clearing of asexual parasites, but relapse can be prevented only if dormant hypnozoites are cleared from the liver (a treatment termed "radical cure"). Tafenoquine is a single-dose 8-aminoquinoline that has recently been registered for the radical cure of P. vivax. METHODS: This multicenter, double-blind, double-dummy, parallel group, randomized, placebo-controlled trial was conducted in Ethiopia, Peru, Brazil, Cambodia, Thailand, and the Philippines. We enrolled 522 patients with microscopically confirmed P. vivax infection (>100 to <100,000 parasites per microliter) and normal glucose-6-phosphate dehydrogenase (G6PD) activity (with normal activity defined as ≥70% of the median value determined at each trial site among 36 healthy male volunteers who were otherwise not involved in the trial). All patients received a 3-day course of chloroquine (total dose of 1500 mg). In addition, patients were assigned to receive a single 300-mg dose of tafenoquine on day 1 or 2 (260 patients), placebo (133 patients), or a 15-mg dose of primaquine once daily for 14 days (129 patients). The primary outcome was the Kaplan-Meier estimated percentage of patients who were free from recurrence at 6 months, defined as P. vivax clearance without recurrent parasitemia. RESULTS: In the intention-to-treat population, the percentage of patients who were free from recurrence at 6 months was 62.4% in the tafenoquine group (95% confidence interval [CI], 54.9 to 69.0), 27.7% in the placebo group (95% CI, 19.6 to 36.6), and 69.6% in the primaquine group (95% CI, 60.2 to 77.1). The hazard ratio for the risk of recurrence was 0.30 (95% CI, 0.22 to 0.40) with tafenoquine as compared with placebo (P<0.001) and 0.26 (95% CI, 0.18 to 0.39) with primaquine as compared with placebo (P<0.001). Tafenoquine was associated with asymptomatic declines in hemoglobin levels, which resolved without intervention. CONCLUSIONS: Single-dose tafenoquine resulted in a significantly lower risk of P. vivax recurrence than placebo in patients with phenotypically normal G6PD activity. (Funded by GlaxoSmithKline and Medicines for Malaria Venture; DETECTIVE ClinicalTrials.gov number, NCT01376167 .).

Distribution and Temporal Dynamics of Plasmodium falciparum Chloroquine Resistance Transporter Mutations Associated With Piperaquine Resistance in Northern Cambodia.

BACKGROUND: Newly emerged mutations within the Plasmodium falciparum chloroquine resistance transporter (PfCRT) can confer piperaquine resistance in the absence of amplified plasmepsin II (pfpm2). In this study, we estimated the prevalence of co-circulating piperaquine resistance mutations in P. falciparum isolates collected in northern Cambodia from 2009 to 2017. METHODS: The sequence of pfcrt was determined for 410 P. falciparum isolates using PacBio amplicon sequencing or whole genome sequencing. Quantitative polymerase chain reaction was used to estimate pfpm2 and pfmdr1 copy number. RESULTS: Newly emerged PfCRT mutations increased in prevalence after the change to dihydroartemisinin-piperaquine in 2010, with >98% of parasites harboring these mutations by 2017. After 2014, the prevalence of PfCRT F145I declined, being outcompeted by parasites with less resistant, but more fit PfCRT alleles. After the change to artesunate-mefloquine, the prevalence of parasites with amplified pfpm2 decreased, with nearly half of piperaquine-resistant PfCRT mutants having single-copy pfpm2. CONCLUSIONS: The large proportion of PfCRT mutants that lack pfpm2 amplification emphasizes the importance of including PfCRT mutations as part of molecular surveillance for piperaquine resistance in this region. Likewise, it is critical to monitor for amplified pfmdr1 in these PfCRT mutants, as increased mefloquine pressure could lead to mutants resistant to both drugs.

Selection of Cytochrome b Mutants Is Rare among Plasmodium falciparum Patients Failing Treatment with Atovaquone-Proguanil in Cambodia.

Atovaquone-proguanil remains effective against multidrug-resistant Plasmodium falciparum in Southeast Asia, but resistance is mediated by a single point mutation in cytochrome b (cytb) that can arise during treatment. Among 14 atovaquone-proguanil treatment failures in a clinical trial in Cambodia, only one recrudescence harbored the cytb mutation Y268C. Deep sequencing did not detect the mutation at baseline or in the first 3 days of treatment, suggesting that it arose de novo Further sequencing across cytb similarly found no low-frequency cytb mutations that were up-selected from baseline to recrudescence. Copy number amplification in dihydroorotate dehydrogenase (DHODH) and cytb as markers of atovaquone tolerance was also absent. Cytb mutation played a minor role in atovaquone-proguanil treatment failures in an active comparator clinical trial.

Update on Shigella and Nontyphoidal Salmonella Antimicrobial Drug Resistance: Implications on Empirical Treatment of Acute Infectious Diarrhea in Cambodia.

Information on causative diarrheal pathogens and their associated antimicrobial susceptibility remains limited for Cambodia. This study describes antimicrobial resistance patterns for Shigella and nontyphoidal Salmonella isolates collected in Cambodia over a 5-year period. Multidrug resistance was shown in 98% of Shigella isolates, with 70%, 11%, and 29% of isolates being resistant to fluoroquinolones, azithromycin, and cephalosporin, respectively. As many as 11% of Shigella isolates were resistant to nearly all oral and parenteral drugs typically used for shigellosis, demonstrating extreme drug resistance phenotypes. Although a vast majority of nontyphoidal Salmonella isolates remained susceptible to cephalosporins (99%) and macrolides (98%), decreased susceptibility to ciprofloxacin was found in 67% of isolates, which is notably higher than previous reports. In conclusion, increasing antimicrobial resistance of Shigella and nontyphoidal Salmonella is a major concern for selecting empirical treatment of acute infectious diarrhea in Cambodia. Treatment practices should be updated and follow local antimicrobial resistance data for the identified pathogens.

Efficient Transmission of Mixed Plasmodium falciparum/vivax Infections From Humans to Mosquitoes.

BACKGROUND: In Southeast Asia, people are often coinfected with different species of malaria (Plasmodium falciparum [Pf] and Plasmodium vivax [Pv]) as well as with multiple clones of the same species. Whether particular species or clones within mixed infections are more readily transmitted to mosquitoes remains unknown. METHODS: Laboratory-reared Anopheles dirus were fed on blood from 119 Pf-infected Cambodian adults, with 5950 dissected to evaluate for transmitted infection. Among 12 persons who infected mosquitoes, polymerase chain reaction and amplicon deep sequencing were used to track species and clone-specific transmission to mosquitoes. RESULTS: Seven of 12 persons that infected mosquitoes harbored mixed Pf/Pv infection. Among these 7 persons, all transmitted Pv with 2 transmitting both Pf and Pv, leading to Pf/Pv coinfection in 21% of infected mosquitoes. Up to 4 clones of each species were detected within persons. Shifts in clone frequency were detected during transmission. However, in general, all parasite clones in humans were transmitted to mosquitoes, with individual mosquitoes frequently carrying multiple transmitted clones. CONCLUSIONS: Malaria diversity in human hosts was maintained in the parasite populations recovered from mosquitoes fed on their blood. However, in persons with mixed Pf/Pv malaria, Pv appears to be transmitted more readily, in association with more prevalent patent gametocytemia.

Correction: Quantification of glucose-6-phosphate dehydrogenase activity by spectrophotometry: A systematic review and meta-analysis.

[This corrects the article DOI: 10.1371/journal.pmed.1003084.].

Quantification of glucose-6-phosphate dehydrogenase activity by spectrophotometry: A systematic review and meta-analysis.

BACKGROUND: The radical cure of Plasmodium vivax and P. ovale requires treatment with primaquine or tafenoquine to clear dormant liver stages. Either drug can induce haemolysis in individuals with glucose-6-phosphate dehydrogenase (G6PD) deficiency, necessitating screening. The reference diagnostic method for G6PD activity is ultraviolet (UV) spectrophotometry; however, a universal G6PD activity threshold above which these drugs can be safely administered is not yet defined. Our study aimed to quantify assay-based variation in G6PD spectrophotometry and to explore the diagnostic implications of applying a universal threshold. METHODS AND FINDINGS: Individual-level data were pooled from studies that used G6PD spectrophotometry. Studies were identified via PubMed search (25 April 2018) and unpublished contributions from contacted authors (PROSPERO: CRD42019121414). Studies were excluded if they assessed only individuals with known haematological conditions, were family studies, or had insufficient details. Studies of malaria patients were included but analysed separately. Included studies were assessed for risk of bias using an adapted form of the Quality Assessment of Diagnostic Accuracy Studies-2 (QUADAS-2) tool. Repeatability and intra- and interlaboratory variability in G6PD activity measurements were compared between studies and pooled across the dataset. A universal threshold for G6PD deficiency was derived, and its diagnostic performance was compared to site-specific thresholds. Study participants (n = 15,811) were aged between 0 and 86 years, and 44.4% (7,083) were women. Median (range) activity of G6PD normal (G6PDn) control samples was 10.0 U/g Hb (6.3-14.0) for the Trinity assay and 8.3 U/g Hb (6.8-15.6) for the Randox assay. G6PD activity distributions varied significantly between studies. For the 13 studies that used the Trinity assay, the adjusted male median (AMM; a standardised metric of 100% G6PD activity) varied from 5.7 to 12.6 U/g Hb (p < 0.001). Assay precision varied between laboratories, as assessed by variance in control measurements (from 0.1 to 1.5 U/g Hb; p < 0.001) and study-wise mean coefficient of variation (CV) of replicate measures (from 1.6% to 14.9%; p < 0.001). A universal threshold of 100% G6PD activity was defined as 9.4 U/g Hb, yielding diagnostic thresholds of 6.6 U/g Hb (70% activity) and 2.8 U/g Hb (30% activity). These thresholds diagnosed individuals with less than 30% G6PD activity with study-wise sensitivity from 89% (95% CI: 81%-94%) to 100% (95% CI: 96%-100%) and specificity from 96% (95% CI: 89%-99%) to 100% (100%-100%). However, when considering intermediate deficiency (<70% G6PD activity), sensitivity fell to a minimum of 64% (95% CI: 52%-75%) and specificity to 35% (95% CI: 24%-46%). Our ability to identify underlying factors associated with study-level heterogeneity was limited by the lack of availability of covariate data and diverse study contexts and methodologies. CONCLUSIONS: Our findings indicate that there is substantial variation in G6PD measurements by spectrophotometry between sites. This is likely due to variability in laboratory methods, with possible contribution of unmeasured population factors. While an assay-specific, universal quantitative threshold offers robust diagnosis at the 30% level, inter-study variability impedes performance of universal thresholds at the 70% level. Caution is advised in comparing findings based on absolute G6PD activity measurements across studies. Novel handheld quantitative G6PD diagnostics may allow greater standardisation in the future.

The risk of Plasmodium vivax parasitaemia after P. falciparum malaria: An individual patient data meta-analysis from the WorldWide Antimalarial Resistance Network.

BACKGROUND: There is a high risk of Plasmodium vivax parasitaemia following treatment of falciparum malaria. Our study aimed to quantify this risk and the associated determinants using an individual patient data meta-analysis in order to identify populations in which a policy of universal radical cure, combining artemisinin-based combination therapy (ACT) with a hypnozoitocidal antimalarial drug, would be beneficial. METHODS AND FINDINGS: A systematic review of Medline, Embase, Web of Science, and the Cochrane Database of Systematic Reviews identified efficacy studies of uncomplicated falciparum malaria treated with ACT that were undertaken in regions coendemic for P. vivax between 1 January 1960 and 5 January 2018. Data from eligible studies were pooled using standardised methodology. The risk of P. vivax parasitaemia at days 42 and 63 and associated risk factors were investigated by multivariable Cox regression analyses. Study quality was assessed using a tool developed by the Joanna Briggs Institute. The study was registered in the International Prospective Register of Systematic Reviews (PROSPERO: CRD42018097400). In total, 42 studies enrolling 15,341 patients were included in the analysis, including 30 randomised controlled trials and 12 cohort studies. Overall, 14,146 (92.2%) patients had P. falciparum monoinfection and 1,195 (7.8%) mixed infection with P. falciparum and P. vivax. The median age was 17.0 years (interquartile range [IQR] = 9.0-29.0 years; range = 0-80 years), with 1,584 (10.3%) patients younger than 5 years. 2,711 (17.7%) patients were treated with artemether-lumefantrine (AL, 13 studies), 651 (4.2%) with artesunate-amodiaquine (AA, 6 studies), 7,340 (47.8%) with artesunate-mefloquine (AM, 25 studies), and 4,639 (30.2%) with dihydroartemisinin-piperaquine (DP, 16 studies). 14,537 patients (94.8%) were enrolled from the Asia-Pacific region, 684 (4.5%) from the Americas, and 120 (0.8%) from Africa. At day 42, the cumulative risk of vivax parasitaemia following treatment of P. falciparum was 31.1% (95% CI 28.9-33.4) after AL, 14.1% (95% CI 10.8-18.3) after AA, 7.4% (95% CI 6.7-8.1) after AM, and 4.5% (95% CI 3.9-5.3) after DP. By day 63, the risks had risen to 39.9% (95% CI 36.6-43.3), 42.4% (95% CI 34.7-51.2), 22.8% (95% CI 21.2-24.4), and 12.8% (95% CI 11.4-14.5), respectively. In multivariable analyses, the highest rate of P. vivax parasitaemia over 42 days of follow-up was in patients residing in areas of short relapse periodicity (adjusted hazard ratio [AHR] = 6.2, 95% CI 2.0-19.5; p = 0.002); patients treated with AL (AHR = 6.2, 95% CI 4.6-8.5; p < 0.001), AA (AHR = 2.3, 95% CI 1.4-3.7; p = 0.001), or AM (AHR = 1.4, 95% CI 1.0-1.9; p = 0.028) compared with DP; and patients who did not clear their initial parasitaemia within 2 days (AHR = 1.8, 95% CI 1.4-2.3; p < 0.001). The analysis was limited by heterogeneity between study populations and lack of data from very low transmission settings. Study quality was high. CONCLUSIONS: In this meta-analysis, we found a high risk of P. vivax parasitaemia after treatment of P. falciparum malaria that varied significantly between studies. These P. vivax infections are likely attributable to relapses that could be prevented with radical cure including a hypnozoitocidal agent; however, the benefits of such a novel strategy will vary considerably between geographical areas.

Pharmacogenetic assessment of tafenoquine efficacy in patients with Plasmodium vivax malaria.

Plasmodium vivax has the largest geographic range of human malaria species and is challenging to manage and eradicate due to its ability to establish a dormant liver stage, the hypnozoite, which can reactivate leading to relapse. Until recently, the only treatment approved to kill hypnozoites was the 8-aminoquinoline, primaquine, requiring daily treatment for 14 days. Tafenoquine, an 8-aminoquinoline single-dose treatment with activity against P. vivax hypnozoites, has recently been approved by the US Food and Drug Administration and Australian Therapeutic Goods Administration for the radical cure of P. vivax malaria in patients 16 years and older. We conducted an exploratory pharmacogenetic analysis (GSK Study 208099) to assess the role of host genome-wide variation on tafenoquine efficacy in patients with P. vivax malaria using data from three GSK clinical trials, GATHER and DETECTIVE Part 1 and Part 2. Recurrence-free efficacy at 6 and 4 months and time to recurrence up to 6 months postdosing were analyzed in 438 P. vivax malaria patients treated with tafenoquine. Among the approximately 10.6 million host genetic variants analyzed, two signals reached genome-wide significance (P value ≤ 5 × 10). rs62103056, and variants in a chromosome 12 intergenic region, were associated with recurrence-free efficacy at 6 and 4 months, respectively. Neither of the signals has an obvious biological rationale and would need replication in an independent population. This is the first genome-wide association study to evaluate genetic influence on response to tafenoquine in P. vivax malaria.

Piperaquine resistant Cambodian Plasmodium falciparum clinical isolates: in vitro genotypic and phenotypic characterization.

BACKGROUND: High rates of dihydroartemisinin-piperaquine (DHA-PPQ) treatment failures have been documented for uncomplicated Plasmodium falciparum in Cambodia. The genetic markers plasmepsin 2 (pfpm2), exonuclease (pfexo) and chloroquine resistance transporter (pfcrt) genes are associated with PPQ resistance and are used for monitoring the prevalence of drug resistance and guiding malaria drug treatment policy. METHODS: To examine the relative contribution of each marker to PPQ resistance, in vitro culture and the PPQ survival assay were performed on seventeen P. falciparum isolates from northern Cambodia, and the presence of E415G-Exo and pfcrt mutations (T93S, H97Y, F145I, I218F, M343L, C350R, and G353V) as well as pfpm2 copy number polymorphisms were determined. Parasites were then cloned by limiting dilution and the cloned parasites were tested for drug susceptibility. Isobolographic analysis of several drug combinations for standard clones and newly cloned P. falciparum Cambodian isolates was also determined. RESULTS: The characterization of culture-adapted isolates revealed that the presence of novel pfcrt mutations (T93S, H97Y, F145I, and I218F) with E415G-Exo mutation can confer PPQ-resistance, in the absence of pfpm2 amplification. In vitro testing of PPQ resistant parasites demonstrated a bimodal dose-response, the existence of a swollen digestive vacuole phenotype, and an increased susceptibility to quinine, chloroquine, mefloquine and lumefantrine. To further characterize drug sensitivity, parental parasites were cloned in which a clonal line, 14-B5, was identified as sensitive to artemisinin and piperaquine, but resistant to chloroquine. Assessment of the clone against a panel of drug combinations revealed antagonistic activity for six different drug combinations. However, mefloquine-proguanil and atovaquone-proguanil combinations revealed synergistic antimalarial activity. CONCLUSIONS: Surveillance for PPQ resistance in regions relying on DHA-PPQ as the first-line treatment is dependent on the monitoring of molecular markers of drug resistance. P. falciparum harbouring novel pfcrt mutations with E415G-exo mutations displayed PPQ resistant phenotype. The presence of pfpm2 amplification was not required to render parasites PPQ resistant suggesting that the increase in pfpm2 copy number alone is not the sole modulator of PPQ resistance. Genetic background of circulating field isolates appear to play a role in drug susceptibility and biological responses induced by drug combinations. The use of latest field isolates may be necessary for assessment of relevant drug combinations against P. falciparum strains and when down-selecting novel drug candidates.

Detecting geospatial patterns of Plasmodium falciparum parasite migration in Cambodia using optimized estimated effective migration surfaces.

BACKGROUND: Understanding the genetic structure of natural populations provides insight into the demographic and adaptive processes that have affected those populations. Such information, particularly when integrated with geospatial data, can have translational applications for a variety of fields, including public health. Estimated effective migration surfaces (EEMS) is an approach that allows visualization of the spatial patterns in genomic data to understand population structure and migration. In this study, we developed a workflow to optimize the resolution of spatial grids used to generate EEMS migration maps and applied this optimized workflow to estimate migration of Plasmodium falciparum in Cambodia and bordering regions of Thailand and Vietnam. METHODS: The optimal density of EEMS grids was determined based on a new workflow created using density clustering to define genomic clusters and the spatial distance between genomic clusters. Topological skeletons were used to capture the spatial distribution for each genomic cluster and to determine the EEMS grid density; i.e., both genomic and spatial clustering were used to guide the optimization of EEMS grids. Model accuracy for migration estimates using the optimized workflow was tested and compared to grid resolutions selected without the optimized workflow. As a test case, the optimized workflow was applied to genomic data generated from P. falciparum sampled in Cambodia and bordering regions, and migration maps were compared to estimates of malaria endemicity, as well as geographic properties of the study area, as a means of validating observed migration patterns. RESULTS: Optimized grids displayed both high model accuracy and reduced computing time compared to grid densities selected in an unguided manner. In addition, EEMS migration maps generated for P. falciparum using the optimized grid corresponded to estimates of malaria endemicity and geographic properties of the study region that might be expected to impact malaria parasite migration, supporting the validity of the observed migration patterns. CONCLUSIONS: Optimized grids reduce spatial uncertainty in the EEMS contours that can result from user-defined parameters, such as the resolution of the spatial grid used in the model. This workflow will be useful to a broad range of EEMS users as it can be applied to analyses involving other organisms of interest and geographic areas.

Selective sweep suggests transcriptional regulation may underlie Plasmodium vivax resilience to malaria control measures in Cambodia.

Cambodia, in which both Plasmodium vivax and Plasmodium falciparum are endemic, has been the focus of numerous malaria-control interventions, resulting in a marked decline in overall malaria incidence. Despite this decline, the number of P vivax cases has actually increased. To understand better the factors underlying this resilience, we compared the genetic responses of the two species to recent selective pressures. We sequenced and studied the genomes of 70 P vivax and 80 P falciparum isolates collected between 2009 and 2013. We found that although P falciparum has undergone population fracturing, the coendemic P vivax population has grown undisrupted, resulting in a larger effective population size, no discernable population structure, and frequent multiclonal infections. Signatures of selection suggest recent, species-specific evolutionary differences. Particularly, in contrast to P falciparum, P vivax transcription factors, chromatin modifiers, and histone deacetylases have undergone strong directional selection, including a particularly strong selective sweep at an AP2 transcription factor. Together, our findings point to different population-level adaptive mechanisms used by P vivax and P falciparum parasites. Although population substructuring in P falciparum has resulted in clonal outgrowths of resistant parasites, P vivax may use a nuanced transcriptional regulatory approach to population maintenance, enabling it to preserve a larger, more diverse population better suited to facing selective threats. We conclude that transcriptional control may underlie P vivax's resilience to malaria control measures. Novel strategies to target such processes are likely required to eradicate P vivax and achieve malaria elimination.

Association of a Novel Mutation in the Plasmodium falciparum Chloroquine Resistance Transporter With Decreased Piperaquine Sensitivity.

Background: Amplified copy number in the plasmepsin II/III genes within Plasmodium falciparum has been associated with decreased sensitivity to piperaquine. To examine this association and test whether additional loci might also contribute, we performed a genome-wide association study of ex vivo P. falciparum susceptibility to piperaquine. Methods: Plasmodium falciparum DNA from 183 samples collected primarily from Cambodia was genotyped at 33716 genome-wide single nucleotide polymorphisms (SNPs). Linear mixed models and random forests were used to estimate associations between parasite genotypes and piperaquine susceptibility. Candidate polymorphisms were evaluated for their association with dihydroartemisinin-piperaquine treatment outcomes in an independent dataset. Results: Single nucleotide polymorphisms on multiple chromosomes were associated with piperaquine 90% inhibitory concentrations (IC90) in a genome-wide analysis. Fine-mapping of genomic regions implicated in genome-wide analyses identified multiple SNPs in linkage disequilibrium with each other that were significantly associated with piperaquine IC90, including a novel mutation within the gene encoding the P. falciparum chloroquine resistance transporter, PfCRT. This mutation (F145I) was associated with dihydroartemisinin-piperaquine treatment failure after adjusting for the presence of amplified plasmepsin II/III, which was also associated with decreased piperaquine sensitivity. Conclusions: Our data suggest that, in addition to plasmepsin II/III copy number, other loci, including pfcrt, may also be involved in piperaquine resistance.