Monocyte Single-Cell Multimodal Profiling in Cardiovascular Disease Risk States.

BACKGROUND: Monocytes are a critical innate immune system cell type that serves homeostatic and immunoregulatory functions. They have been identified historically by the cell surface expression of CD14 and CD16. However, recent single-cell studies have revealed that they are much more heterogeneous than previously realized. METHODS: We utilized cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) and single-cell RNA sequencing to describe the comprehensive transcriptional and phenotypic landscape of 437 126 monocytes. RESULTS: This high-dimensional multimodal approach identified vast phenotypic diversity and functionally distinct subsets, including IFN-responsive, MHCIIhi (major histocompatibility complex class II), monocyte-platelet aggregates, as well as nonclassical, and several subpopulations of classical monocytes. Using flow cytometry, we validated the existence of MHCII+CD275+ MHCIIhi, CD42b+ monocyte-platelet aggregates, CD16+CD99- nonclassical monocytes, and CD99+ classical monocytes. Each subpopulation exhibited unique characteristics, developmental trajectories, transcriptional regulation, and tissue distribution. In addition, alterations associated with cardiovascular disease risk factors, including race, smoking, and hyperlipidemia were identified. Moreover, the effect of hyperlipidemia was recapitulated in mouse models of elevated cholesterol. CONCLUSIONS: This integrative and cross-species comparative analysis provides a new perspective on the comparison of alterations in monocytes in pathological conditions and offers insights into monocyte-driven mechanisms in cardiovascular disease and the potential for monocyte subpopulation targeted therapies.

Atherosclerosis Is a Smooth Muscle Cell-Driven Tumor-Like Disease.

BACKGROUND: Atherosclerosis, a leading cause of cardiovascular disease, involves the pathological activation of various cell types, including immunocytes (eg, macrophages and T cells), smooth muscle cells (SMCs), and endothelial cells. Accumulating evidence suggests that transition of SMCs to other cell types, known as phenotypic switching, plays a central role in atherosclerosis development and complications. However, the characteristics of SMC-derived cells and the underlying mechanisms of SMC transition in disease pathogenesis remain poorly understood. Our objective is to characterize tumor cell-like behaviors of SMC-derived cells in atherosclerosis, with the ultimate goal of developing interventions targeting SMC transition for the prevention and treatment of atherosclerosis. METHODS: We used SMC lineage tracing mice and human tissues and applied a range of methods, including molecular, cellular, histological, computational, human genetics, and pharmacological approaches, to investigate the features of SMC-derived cells in atherosclerosis. RESULTS: SMC-derived cells in mouse and human atherosclerosis exhibit multiple tumor cell-like characteristics, including genomic instability, evasion of senescence, hyperproliferation, resistance to cell death, invasiveness, and activation of comprehensive cancer-associated gene regulatory networks. Specific expression of the oncogenic mutant KrasG12D in SMCs accelerates phenotypic switching and exacerbates atherosclerosis. Furthermore, we provide proof of concept that niraparib, an anticancer drug targeting DNA damage repair, attenuates atherosclerosis progression and induces regression of lesions in advanced disease in mouse models. CONCLUSIONS: Our findings demonstrate that atherosclerosis is an SMC-driven tumor-like disease, advancing our understanding of its pathogenesis and opening prospects for innovative precision molecular strategies aimed at preventing and treating atherosclerotic cardiovascular disease.

High-Dimensional Single-Cell Multimodal Landscape of Human Carotid Atherosclerosis.

BACKGROUND: Atherosclerotic plaques are complex tissues composed of a heterogeneous mixture of cells. However, our understanding of the comprehensive transcriptional and phenotypic landscape of the cells within these lesions is limited. METHODS: To characterize the landscape of human carotid atherosclerosis in greater detail, we combined cellular indexing of transcriptomes and epitopes by sequencing and single-cell RNA sequencing to classify all cell types within lesions (n=21; 13 symptomatic) to achieve a comprehensive multimodal understanding of the cellular identities of atherosclerosis and their association with clinical pathophysiology. RESULTS: We identified 25 cell populations, each with a unique multiomic signature, including macrophages, T cells, NK (natural killer) cells, mast cells, B cells, plasma cells, neutrophils, dendritic cells, endothelial cells, fibroblasts, and smooth muscle cells (SMCs). Among the macrophages, we identified 2 proinflammatory subsets enriched in IL-1B (interleukin-1B) or C1Q expression, 2 TREM2-positive foam cells (1 expressing inflammatory genes), and subpopulations with a proliferative gene signature and SMC-specific gene signature with fibrotic pathways upregulated. Further characterization revealed various subsets of SMCs and fibroblasts, including SMC-derived foam cells. These foamy SMCs were localized in the deep intima of coronary atherosclerotic lesions. Utilizing cellular indexing of transcriptomes and epitopes by sequencing data, we developed a flow cytometry panel, using cell surface proteins CD29, CD142, and CD90, to isolate SMC-derived cells from lesions. Lastly, we observed reduced proportions of efferocytotic macrophages, classically activated endothelial cells, and contractile and modulated SMC-derived cells, while inflammatory SMCs were enriched in plaques of clinically symptomatic versus asymptomatic patients. CONCLUSIONS: Our multimodal atlas of cell populations within atherosclerosis provides novel insights into the diversity, phenotype, location, isolation, and clinical relevance of the unique cellular composition of human carotid atherosclerosis. These findings facilitate both the mapping of cardiovascular disease susceptibility loci to specific cell types and the identification of novel molecular and cellular therapeutic targets for the treatment of the disease.

Comparative 3D genome organization in apicomplexan parasites.

The positioning of chromosomes in the nucleus of a eukaryotic cell is highly organized and has a complex and dynamic relationship with gene expression. In the human malaria parasite Plasmodium falciparum, the clustering of a family of virulence genes correlates with their coordinated silencing and has a strong influence on the overall organization of the genome. To identify conserved and species-specific principles of genome organization, we performed Hi-C experiments and generated 3D genome models for five Plasmodium species and two related apicomplexan parasites. Plasmodium species mainly showed clustering of centromeres, telomeres, and virulence genes. In P. falciparum, the heterochromatic virulence gene cluster had a strong repressive effect on the surrounding nuclear space, while this was less pronounced in Plasmodium vivax and Plasmodium berghei, and absent in Plasmodium yoelii In Plasmodium knowlesi, telomeres and virulence genes were more dispersed throughout the nucleus, but its 3D genome showed a strong correlation with gene expression. The Babesia microti genome showed a classical Rabl organization with colocalization of subtelomeric virulence genes, while the Toxoplasma gondii genome was dominated by clustering of the centromeres and lacked virulence gene clustering. Collectively, our results demonstrate that spatial genome organization in most Plasmodium species is constrained by the colocalization of virulence genes. P. falciparum and P. knowlesi, the only two Plasmodium species with gene families involved in antigenic variation, are unique in the effect of these genes on chromosome folding, indicating a potential link between genome organization and gene expression in more virulent pathogens.

In vitro resistance selections for Plasmodium falciparum dihydroorotate dehydrogenase inhibitors give mutants with multiple point mutations in the drug-binding site and altered growth.

Malaria is a preventable and treatable disease; yet half of the world's population lives at risk of infection, and an estimated 660,000 people die of malaria-related causes every year. Rising drug resistance threatens to make malaria untreatable, necessitating both the discovery of new antimalarial agents and the development of strategies to identify and suppress the emergence and spread of drug resistance. We focused on in-development dihydroorotate dehydrogenase (DHODH) inhibitors. Characterizing resistance pathways for antimalarial agents not yet in clinical use will increase our understanding of the potential for resistance. We identified resistance mechanisms of Plasmodium falciparum (Pf) DHODH inhibitors via in vitro resistance selections. We found 11 point mutations in the PfDHODH target. Target gene amplification and unknown mechanisms also contributed to resistance, albeit to a lesser extent. These mutant parasites were often hypersensitive to other PfDHODH inhibitors, which immediately suggested a novel combination therapy approach to preventing resistance. Indeed, a combination of wild-type and mutant-type selective inhibitors led to resistance far less often than either drug alone. The effects of point mutations in PfDHODH were corroborated with purified recombinant wild-type and mutant-type PfDHODH proteins, which showed the same trends in drug response as the cognate cell lines. Comparative growth assays demonstrated that two mutant parasites grew less robustly than their wild-type parent, and the purified protein of those mutants showed a decrease in catalytic efficiency, thereby suggesting a reason for the diminished growth rate. Co-crystallography of PfDHODH with three inhibitors suggested that hydrophobic interactions are important for drug binding and selectivity.

Identification of the drug/metabolite transporter 1 as a marker of quinine resistance in a NF54×Cam3.II P. falciparum genetic cross.

The genetic basis of Plasmodium falciparum resistance to quinine (QN), a drug used to treat severe malaria, has long been enigmatic. To gain further insight, we used FRG-NOD human liver-chimeric mice to conduct a P. falciparum genetic cross between QN-sensitive and QN-resistant parasites, which also differ in their susceptibility to chloroquine (CQ). By applying different selective conditions to progeny pools prior to cloning, we recovered 120 unique recombinant progeny. These progeny were subjected to drug profiling and QTL analyses with QN, CQ, and monodesethyl-CQ (md-CQ, the active metabolite of CQ), which revealed predominant peaks on chromosomes 7 and 12, consistent with a multifactorial mechanism of resistance. A shared chromosome 12 region mapped to resistance to all three antimalarials and was preferentially co-inherited with pfcrt. We identified an ATP-dependent zinc metalloprotease (FtsH1) as one of the top candidates and observed using CRISPR/Cas9 SNP-edited lines that ftsh1 is a potential mediator of QN resistance and a modulator of md-CQ resistance. As expected, CQ and md-CQ resistance mapped to a chromosome 7 region harboring pfcrt. However, for QN, high-grade resistance mapped to a chromosome 7 peak centered 295kb downstream of pfcrt. We identified the drug/metabolite transporter 1 (DMT1) as the top candidate due to its structural similarity to PfCRT and proximity to the peak. Deleting DMT1 in QN-resistant Cam3.II parasites significantly sensitized the parasite to QN but not to the other drugs tested, suggesting that DMT1 mediates QN response specifically. We localized DMT1 to structures associated with vesicular trafficking, as well as the parasitophorous vacuolar membrane, lipid bodies, and the digestive vacuole. We also observed that mutant DMT1 transports more QN than the wild-type isoform in vitro. Our study demonstrates that DMT1 is a novel marker of QN resistance and a new chromosome 12 locus associates with CQ and QN response, with ftsh1 is a potential candidate, suggesting these genes should be genotyped in surveillance and clinical settings.

Atherosclerosis is a smooth muscle cell-driven tumor-like disease.

Atherosclerosis, the leading cause of cardiovascular disease, is a chronic inflammatory disease involving pathological activation of multiple cell types, such as immunocytes (e.g., macrophage, T cells), smooth muscle cells (SMCs), and endothelial cells. Multiple lines of evidence have suggested that SMC "phenotypic switching" plays a central role in atherosclerosis development and complications. Yet, SMC roles and mechanisms underlying the disease pathogenesis are poorly understood. Here, employing SMC lineage tracing mice, comprehensive molecular, cellular, histological, and computational profiling, coupled to genetic and pharmacological studies, we reveal that atherosclerosis, in terms of SMC behaviors, share extensive commonalities with tumors. SMC-derived cells in the disease show multiple characteristics of tumor cell biology, including genomic instability, replicative immortality, malignant proliferation, resistance to cell death, invasiveness, and activation of comprehensive cancer-associated gene regulatory networks. SMC-specific expression of oncogenic KrasG12D accelerates SMC phenotypic switching and exacerbates atherosclerosis. Moreover, we present a proof of concept showing that niraparib, an anti-cancer drug targeting DNA damage repair, attenuates atherosclerosis progression and induces regression of lesions in advanced disease in mouse models. Our work provides systematic evidence that atherosclerosis is a tumor-like disease, deepening the understanding of its pathogenesis and opening prospects for novel precision molecular strategies to prevent and treat atherosclerotic cardiovascular disease.

High-Dimensional Single-Cell Multimodal Landscape of Human Carotid Atherosclerosis.

Background: Atherosclerotic plaques are complex tissues composed of a heterogeneous mixture of cells. However, we have limited understanding of the comprehensive transcriptional and phenotypical landscape of the cells within these lesions. Methods: To characterize the landscape of human carotid atherosclerosis in greater detail, we combined cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) and single-cell RNA sequencing (scRNA-seq) to classify all cell types within lesions (n=21; 13 symptomatic) to achieve a comprehensive multimodal understanding of the cellular identities of atherosclerosis and their association with clinical pathophysiology. Results: We identified 25 distinct cell populations each having a unique multi-omic signature, including macrophages, T cells, NK cells, mast cells, B cells, plasma cells, neutrophils, dendritic cells, endothelial cells, fibroblasts, and smooth muscle cells (SMCs). Within the macrophage populations, we identified 2 proinflammatory subsets that were enriched in IL1B or C1Q expression, 2 distinct TREM2 positive foam cell subsets, one of which also expressed inflammatory genes, as well as subpopulations displaying a proliferative gene expression signature and one expressing SMC-specific genes and upregulation of fibrotic pathways. An in-depth characterization uncovered several subsets of SMCs and fibroblasts, including a SMC-derived foam cell. We localized this foamy SMC to the deep intima of coronary atherosclerotic lesions. Using CITE-seq data, we also developed the first flow cytometry panel, using cell surface proteins CD29, CD142, and CD90, to isolate SMC-derived cells from lesions. Last, we found that the proportion of efferocytotic macrophages, classically activated endothelial cells, contractile and modulated SMC-derived cell types were reduced, and inflammatory SMCs were enriched in plaques of clinically symptomatic vs. asymptomatic patients. Conclusions: Our multimodal atlas of cell populations within atherosclerosis provides novel insights into the diversity, phenotype, location, isolation, and clinical relevance of the unique cellular composition of human carotid atherosclerosis. This facilitates both the mapping of cardiovascular disease susceptibility loci to specific cell types as well as the identification of novel molecular and cellular therapeutic targets for treatment of the disease.

Mapping the genomic landscape of multidrug resistance in Plasmodium falciparum and its impact on parasite fitness.

Drug-resistant Plasmodium falciparum parasites have swept across Southeast Asia and now threaten Africa. By implementing a P. falciparum genetic cross using humanized mice, we report the identification of key determinants of resistance to artemisinin (ART) and piperaquine (PPQ) in the dominant Asian KEL1/PLA1 lineage. We mapped k13 as the central mediator of ART resistance and identified secondary markers. Applying bulk segregant analysis, quantitative trait loci mapping and gene editing, our data reveal an epistatic interaction between mutant PfCRT and multicopy plasmepsins 2/3 in mediating high-grade PPQ resistance. Susceptibility and parasite fitness assays implicate PPQ as a driver of selection for KEL1/PLA1 parasites. Mutant PfCRT enhanced susceptibility to lumefantrine, the first-line partner drug in Africa, highlighting a potential benefit of opposing selective pressures with this drug and PPQ. We also identified that the ABCI3 transporter can operate in concert with PfCRT and plasmepsins 2/3 in mediating multigenic resistance to antimalarial agents.

Mapping the genomic landscape of multidrug resistance in Plasmodium falciparum and its impact on parasite fitness.

Drug-resistant Plasmodium falciparum parasites have swept across Southeast Asia and now threaten Africa. By implementing a P. falciparum genetic cross using humanized mice, we report the identification of key determinants of resistance to artemisinin (ART) and piperaquine (PPQ) in the dominant Asian KEL1/PLA1 lineage. We mapped k13 as the central mediator of ART resistance in vitro and identified secondary markers. Applying bulk segregant analysis, quantitative trait loci mapping using 34 recombinant haplotypes, and gene editing, our data reveal an epistatic interaction between mutant PfCRT and multicopy plasmepsins 2/3 in mediating high-grade PPQ resistance. Susceptibility and parasite fitness assays implicate PPQ as a driver of selection for KEL1/PLA1 parasites. Mutant PfCRT enhanced susceptibility to lumefantrine, the first-line partner drug in Africa, highlighting a potential benefit of opposing selective pressures with this drug and PPQ. We also identified that the ABCI3 transporter can operate in concert with PfCRT and plasmepsins 2/3 in mediating multigenic resistance to antimalarial agents.

Artemisinin-resistant K13 mutations rewire Plasmodium falciparum's intra-erythrocytic metabolic program to enhance survival.

The emergence and spread of artemisinin resistance, driven by mutations in Plasmodium falciparum K13, has compromised antimalarial efficacy and threatens the global malaria elimination campaign. By applying systems-based quantitative transcriptomics, proteomics, and metabolomics to a panel of isogenic K13 mutant or wild-type P. falciparum lines, we provide evidence that K13 mutations alter multiple aspects of the parasite's intra-erythrocytic developmental program. These changes impact cell-cycle periodicity, the unfolded protein response, protein degradation, vesicular trafficking, and mitochondrial metabolism. K13-mediated artemisinin resistance in the Cambodian Cam3.II line was reversed by atovaquone, a mitochondrial electron transport chain inhibitor. These results suggest that mitochondrial processes including damage sensing and anti-oxidant properties might augment the ability of mutant K13 to protect P. falciparum against artemisinin action by helping these parasites undergo temporary quiescence and accelerated growth recovery post drug elimination.

Plasmodium falciparum K13 mutations in Africa and Asia impact artemisinin resistance and parasite fitness.

The emergence of mutant K13-mediated artemisinin (ART) resistance in Plasmodium falciparum malaria parasites has led to widespread treatment failures across Southeast Asia. In Africa, K13-propeller genotyping confirms the emergence of the R561H mutation in Rwanda and highlights the continuing dominance of wild-type K13 elsewhere. Using gene editing, we show that R561H, along with C580Y and M579I, confer elevated in vitro ART resistance in some African strains, contrasting with minimal changes in ART susceptibility in others. C580Y and M579I cause substantial fitness costs, which may slow their dissemination in high-transmission settings, in contrast with R561H that in African 3D7 parasites is fitness neutral. In Cambodia, K13 genotyping highlights the increasing spatio-temporal dominance of C580Y. Editing multiple K13 mutations into a panel of Southeast Asian strains reveals that only the R561H variant yields ART resistance comparable to C580Y. In Asian Dd2 parasites C580Y shows no fitness cost, in contrast with most other K13 mutations tested, including R561H. Editing of point mutations in ferredoxin or mdr2, earlier associated with resistance, has no impact on ART susceptibility or parasite fitness. These data underline the complex interplay between K13 mutations, parasite survival, growth and genetic background in contributing to the spread of ART resistance.

Chemoprotective antimalarials identified through quantitative high-throughput screening of Plasmodium blood and liver stage parasites.

The spread of Plasmodium falciparum parasites resistant to most first-line antimalarials creates an imperative to enrich the drug discovery pipeline, preferably with curative compounds that can also act prophylactically. We report a phenotypic quantitative high-throughput screen (qHTS), based on concentration-response curves, which was designed to identify compounds active against Plasmodium liver and asexual blood stage parasites. Our qHTS screened over 450,000 compounds, tested across a range of 5 to 11 concentrations, for activity against Plasmodium falciparum asexual blood stages. Active compounds were then filtered for unique structures and drug-like properties and subsequently screened in a P. berghei liver stage assay to identify novel dual-active antiplasmodial chemotypes. Hits from thiadiazine and pyrimidine azepine chemotypes were subsequently prioritized for resistance selection studies, yielding distinct mutations in P. falciparum cytochrome b, a validated antimalarial drug target. The thiadiazine chemotype was subjected to an initial medicinal chemistry campaign, yielding a metabolically stable analog with sub-micromolar potency. Our qHTS methodology and resulting dataset provides a large-scale resource to investigate Plasmodium liver and asexual blood stage parasite biology and inform further research to develop novel chemotypes as causal prophylactic antimalarials.

Elucidating Mechanisms of Drug-Resistant Plasmodium falciparum.

Intensified treatment and control efforts since the early 2000s have dramatically reduced the burden of Plasmodium falciparum malaria. However, drug resistance threatens to derail this progress. In this review, we present four antimalarial resistance case studies that differ in timeline, technical approaches, mechanisms of action, and categories of resistance: chloroquine, sulfadoxine-pyrimethamine, artemisinin, and piperaquine. Lessons learned from prior losses of treatment efficacy, drug combinations, and control strategies will help advance mechanistic research into how P. falciparum parasites acquire resistance to current first-line artemisinin-based combination therapies. Understanding resistance in the clinic and laboratory is essential to prolong the effectiveness of current antimalarial drugs and to optimize the pipeline of future medicines.

Emerging Southeast Asian PfCRT mutations confer Plasmodium falciparum resistance to the first-line antimalarial piperaquine.

The widely used antimalarial combination therapy dihydroartemisinin + piperaquine (DHA + PPQ) has failed in Cambodia. Here, we perform a genomic analysis that reveals a rapid increase in the prevalence of novel mutations in the Plasmodium falciparum chloroquine resistance transporter PfCRT following DHA + PPQ implementation. These mutations occur in parasites harboring the K13 C580Y artemisinin resistance marker. By introducing PfCRT mutations into sensitive Dd2 parasites or removing them from resistant Cambodian isolates, we show that the H97Y, F145I, M343L, or G353V mutations each confer resistance to PPQ, albeit with fitness costs for all but M343L. These mutations sensitize Dd2 parasites to chloroquine, amodiaquine, and quinine. In Dd2 parasites, multicopy plasmepsin 2, a candidate molecular marker, is not necessary for PPQ resistance. Distended digestive vacuoles were observed in pfcrt-edited Dd2 parasites but not in Cambodian isolates. Our findings provide compelling evidence that emerging mutations in PfCRT can serve as a molecular marker and mediator of PPQ resistance.

Changes in genome organization of parasite-specific gene families during the Plasmodium transmission stages.

The development of malaria parasites throughout their various life cycle stages is coordinated by changes in gene expression. We previously showed that the three-dimensional organization of the Plasmodium falciparum genome is strongly associated with gene expression during its replication cycle inside red blood cells. Here, we analyze genome organization in the P. falciparum and P. vivax transmission stages. Major changes occur in the localization and interactions of genes involved in pathogenesis and immune evasion, host cell invasion, sexual differentiation, and master regulation of gene expression. Furthermore, we observe reorganization of subtelomeric heterochromatin around genes involved in host cell remodeling. Depletion of heterochromatin protein 1 (PfHP1) resulted in loss of interactions between virulence genes, confirming that PfHP1 is essential for maintenance of the repressive center. Our results suggest that the three-dimensional genome structure of human malaria parasites is strongly connected with transcriptional activity of specific gene families throughout the life cycle.

A surrogate marker of piperaquine-resistant Plasmodium falciparum malaria: a phenotype-genotype association study.

BACKGROUND: Western Cambodia is the epicentre of Plasmodium falciparum multidrug resistance and is facing high rates of dihydroartemisinin-piperaquine treatment failures. Genetic tools to detect the multidrug-resistant parasites are needed. Artemisinin resistance can be tracked using the K13 molecular marker, but no marker exists for piperaquine resistance. We aimed to identify genetic markers of piperaquine resistance and study their association with dihydroartemisinin-piperaquine treatment failures. METHODS: We obtained blood samples from Cambodian patients infected with P falciparum and treated with dihydroartemisinin-piperaquine. Patients were followed up for 42 days during the years 2009-15. We established in-vitro and ex-vivo susceptibility profiles for a subset using piperaquine survival assays. We determined whole-genome sequences by Illumina paired-reads sequencing, copy number variations by qPCR, RNA concentrations by qRT-PCR, and protein concentrations by immunoblotting. Fisher's exact and non-parametric Wilcoxon rank-sum tests were used to identify significant differences in single-nucleotide polymorphisms or copy number variants, respectively, for differential distribution between piperaquine-resistant and piperaquine-sensitive parasite lines. FINDINGS: Whole-genome exon sequence analysis of 31 culture-adapted parasite lines associated amplification of the plasmepsin 2-plasmepsin 3 gene cluster with in-vitro piperaquine resistance. Ex-vivo piperaquine survival assay profiles of 134 isolates correlated with plasmepsin 2 gene copy number. In 725 patients treated with dihydroartemisinin-piperaquine, multicopy plasmepsin 2 in the sample collected before treatment was associated with an adjusted hazard ratio (aHR) for treatment failure of 20·4 (95% CI 9·1-45·5, p<0·0001). Multicopy plasmepsin 2 predicted dihydroartemisinin-piperaquine failures with 0·94 (95% CI 0·88-0·98) sensitivity and 0·77 (0·74-0·81) specificity. Analysis of samples collected across the country from 2002 to 2015 showed that the geographical and temporal increase of the proportion of multicopy plasmepsin 2 parasites was highly correlated with increasing dihydroartemisinin-piperaquine treatment failure rates (r=0·89 [95% CI 0·77-0·95], p<0·0001, Spearman's coefficient of rank correlation). Dihydroartemisinin-piperaquine efficacy at day 42 fell below 90% when the proportion of multicopy plasmepsin 2 parasites exceeded 22%. INTERPRETATION: Piperaquine resistance in Cambodia is strongly associated with amplification of plasmepsin 2-3, encoding haemoglobin-digesting proteases, regardless of the location. Multicopy plasmepsin 2 constitutes a surrogate molecular marker to track piperaquine resistance. A molecular toolkit combining plasmepsin 2 with K13 and mdr1 monitoring should provide timely information for antimalarial treatment and containment policies. FUNDING: Institut Pasteur in Cambodia, Institut Pasteur Paris, National Institutes of Health, WHO, Agence Nationale de la Recherche, Investissement d'Avenir programme, Laboratoire d'Excellence Integrative "Biology of Emerging Infectious Diseases".