Approche génomique et bioinformatique de l'émergence et de la diffusion des résistances chez Plasmodium au Cambodge

Malaria, a protozoan vector-borne disease, is mainly prevalent in tropical areas, where nearly 40% of the world population is residing and remains one of the most concerns for public health worldwide. In Cambodia, the five Plasmodium species known to cause malaria in humans are present. The main feature of this country is that it is recognized as the epicenter of the emergence of multi-resistant P. falciparum parasites (to chloroquine, sulfadoxine-pyrimethamine, mefloquine, and artemisinin), a very significant menace to public health in the Mekong region that could impact the worldwide strategy to fight malaria. The thesis presented here, entitled “Genomics and Bioinformatics in the emergence and spread of resistant Plasmodium in Cambodia” aimed to develop new molecular and biological tools for:1) improving our understanding of the collateral impact of the strategies implemented to fight against falciparum malaria on the other Plasmodium species; 2) defining the molecular epidemiology of antimalarial resistant parasites, especially resistance to quinine and artemisinin derivatives;3) studying and defining the structure of P. falciparum parasite populations circulating in Cambodia to estimate areas at risk of spread of artemisinin resistance, using genomic approaches and bioinformatics. This thesis was carried out in the Malaria Molecular Epidemiology Unit at Pasteur Institute in Cambodia (IPC) under the co-direction of Dr. Didier Ménard (Head of the Unit, IP) and Pr. Emmanuel Cornillot (Professor, University of Montpellier I). The first objective of this work was to study the impact of drug used to treat falciparum malaria on the dynamics of other Plasmodium species. In a first step, we evaluated the polymorphism in gene associated to pyrimethamine resistance (dhfr gene, dihydrofolate reductase) in Plasmodium malariae and in Plasmodium ovale (article 1 and manuscript in preparation 1) and the polymorphism in mdr-1 gene (multidrug resistance 1 gene) associated to mefloquine resistance in P. vivax (article 2). Secondly, in collaboration with Pasteur Institute in Madagascar, we investigated the association between the polymorphism in Plasmodium falciparum Na + / H + exchanger gene (Pfnhe-1) and quinine resistance defined either by clinical or in vitro phenotypes (articles 3 and 4). The second objective was focused on the development of novel biological and molecular tools to assess the resistance of P. falciparum to artemisinin derivatives. The three papers presented (articles 5, 6 and 7) describe an original approach combining genomics, biological, clinical and epidemiological studies, which lead to the discovery of a molecular marker (mutations Kelch 13 gene) associated to artemisinin resistance.The third and final objective was devoted to the development of the PCR-LDR-FMA technology applied to the detection of a panel of 24 SNPs to characterize a "barcode" of P. falciparum isolates. This technic coupled with bioinformatics and statistical analysis allowed us to study and define the structure of the parasite populations circulating in Cambodia for estimating areas at risk of spread of artemisinin resistance (manuscript in preparation 2). Through this work, we have tried to show the usefulness of available molecular biology methods coupled with genomic and bioinformatics approaches to improve our understanding of the dynamics of the malaria parasite populations. This work has been mainly focused on the emergence and spread of antimalarial resistant parasites, keeping in mind that the ultimate goal of this work was to improve strategies implemented to achieve the ambitious goal of malaria elimination.Le paludisme, maladie parasitaire et vectorielle, sévissant principalement dans les zones intertropicales où vit près de 40% de la population mondiale, reste un problème majeur en santé publique. Les cinq espèces de Plasmodium connues infectées le paludisme chez l'homme sont présentes au Cambodge, qui est reconnu comme l'épicentre de l'émergence de souches de P. falciparum multi-résistantes (chloroquine, sulfadoxine- pyriméthamine, méfloquine, artémisinine), pouvant entraver les progrès accomplis depuis plus d'une décennie. Le travail de thèse intitulé « Approche génomique et bio-informatique de l'émergence et de la diffusion des résistances chez Plasmodium au Cambodge » avait pour objectif de développer de nouveaux outils moléculaires et biologiques pour 1) une meilleure compréhension de l'impact des stratégies mises en place pour lutter contre le paludisme à P. falciparum sur les autres espèces de Plasmodium, 2) la mise en place d'outils biologique et moléculaire, permettant de mieux définir l'épidémiologie des parasites résistants, en particulier la résistance à la quinine et aux dérivés de l'artémisinine, 3) l'étude et la définition des sous-populations parasitaires circulant au Cambodge afin d'estimer les risques associés à la diffusion de la résistance à l'artémisinine. Cette thèse a été réalisée dans l'Unité d'Epidémiologie Moléculaire du Paludisme à l'Institut Pasteur du Cambodge (IPC) sous la codirection du Dr. Didier Ménard (Chef de laboratoire à l'IP) et du Pr. Emmanuel Cornillot (Professeur à l'Université Montpellier I). Le premier objectif visait à étudier l'impact de la pression médicamenteuse sur la dynamique d'évolution des populations parasitaires. Nous avons d'abord évalué le polymorphisme des gènes associés à la résistance à la pyriméthamine (gène dhfr, dihydrofolate reductase) chez Plasmodium malariae et Plasmodium ovale (article 1 et manuscrit en préparation1), et le polymorphisme du gène mdr-1 (multidrug resistance 1) associé à la résistance à la mefloquine chez P. vivax (article 2). De plus, en collaboration avec l'Institut Pasteur de Madagascar, nous avons étudié le lien pouvant exister entre le polymorphisme du gène candidat Plasmodium falciparum Na+/H+ exchanger (Pfnhe-1) et la résistance (clinique et in vitro) de P. falciparum à la quinine (articles 3 et 4).Le deuxième objectif s'est interessé au développement d'outils biologiques et moléculaires permettant d'évaluer la résistance des souches de P. falciparum aux dérivés de l'artémisinine. Les 3 articles présentés (articles 5, 6 et 7) decrivent la méthodologie d'approche originale utilisée associant la génomique, la biologie, la clinique et l'épidémiologie, qui a permis d'aboutir à la découverte d'un marqueur moléculaire (mutations au sein du gène Kelch 13) fiable pour identifier les souches résistantes aux dérivés de l'artémisinine.Le dernier objectif était consacré au développement de la technique PCR-LDR-FMA appliqué à la détection d'un panel de 24 SNPs permettant de caractériser par un « barcode » chaque isolat de P. falciparum. Cette technique couplée avec une analyse bio-informatique et statistique des données nous a permis d'étudier et de définir la structuration des populations parasitaires circulant au Cambodge afin d'estimer les zones à risque de diffusion de la résistance à l'artémisinine (manuscrit en préparation 2).A travers ce travail de thèse, nous nous sommes efforcés de montrer la puissance des techniques de biologie moléculaire disponibles couplées avec des approches génomique et bio-informatique pour améliorer notre compréhension de la dynamique d'évolution des populations parasitaires. Ce travail s'est essentiellement concentré sur les phénomènes liés à l'émergence et de la diffusion des parasites résistants aux antipaludiques, le but final de ce travail étant d'améliorer les stratégies de lutte mises en place pour atteindre l'ambitieux objectif d'élimination du paludisme

Large-scale malaria survey in Cambodia: Novel insights on species distribution and risk factors

BACKGROUND: In Cambodia, estimates of the malaria burden rely on a public health information system that does not record cases occurring among remote populations, neither malaria cases treated in the private sector nor asymptomatic carriers. A global estimate of the current malaria situation and associated risk factors is, therefore, still lacking. METHODS: A large cross-sectional survey was carried out in three areas of multidrug resistant malaria in Cambodia, enrolling 11,652 individuals. Fever and splenomegaly were recorded. Malaria prevalence, parasite densities and spatial distribution of infection were determined to identify parasitological profiles and the associated risk factors useful for improving malaria control programmes in the country. RESULTS: Malaria prevalence was 3.0%, 7.0% and 12.3% in Sampovloun, Koh Kong and Preah Vihear areas. Prevalences and Plasmodium species were heterogeneously distributed, with higher Plasmodium vivax rates in areas of low transmission. Malaria-attributable fevers accounted only for 10–33% of malaria cases, and 23–33% of parasite carriers were febrile. Multivariate multilevel regression analysis identified adults and males, mostly involved in forest activities, as high risk groups in Sampovloun, with additional risks for children in forest-fringe villages in the other areas along with an increased risk with distance from health facilities. CONCLUSION: These observations point to a more complex malaria situation than suspected from official reports. A large asymptomatic reservoir was observed. The rates of P. vivax infections were higher than recorded in several areas. In remote areas, malaria prevalence was high. This indicates that additional health facilities should be implemented in areas at higher risk, such as remote rural and forested parts of the country, which are not adequately served by health services. Precise malaria risk mapping all over the country is needed to assess the extensive geographical heterogeneity of malaria endemicity and risk populations, so that current malaria control measures can be reinforced accordingly

Performance of Ultrasensitive Rapid Diagnostic Tests for Detecting Asymptomatic Plasmodium falciparum.

Proposed interventions for eliminating drug-resistant Plasmodium falciparum malaria include the targeting of asymptomatic carriers through screening and treatment. We report on the diagnostic performance of the recently developed ultrasensitive rapid diagnostic test (uRDT) compared with screening with conventional RDTs (cRDT) and polymerase chain reaction (PCR) under field conditions in Cambodia in a total of 2,729 individuals. The P. falciparum positivity by quantitative PCR (qPCR) was 3.8% (26/678) in those screened during active case detection and 0.5% (10/2,051) in the cross-sectional survey. Compared with qPCR, the sensitivity of the uRDTs was 53.8% (95% CI: 33.4-73.4%) when used in active case detection and 60.0% (95% CI: 26.2-87.8%) in the cross-sectional survey. The uRDTs did not show a significant improvement in diagnostic performance over cRDTs when used for active case detection and for a malaria prevalence survey in the context of this low-transmission setting

Plasmodium falciparum parasite population structure and gene flow associated to anti-malarial drugs resistance in Cambodia

Background: Western Cambodia is recognized as the epicentre of emergence of Plasmodium falciparum multi-drug resistance. The emergence of artemisinin resistance has been observed in this area since 2008–2009 and molecular signatures associated to artemisinin resistance have been characterized in k13 gene. At present, one of the major threats faced, is the possible spread of Asian artemisinin resistant parasites over the world threatening millions of people and jeopardizing malaria elimination programme efforts. To anticipate the diffusion of artemisinin resistance, the identification of the P. falciparum population structure and the gene flow among the parasite population in Cambodia are essential. Methods: To this end, a mid-throughput PCR-LDR-FMA approach based on LUMINEX technology was developed to screen for genetic barcode in 533 blood samples collected in 2010–2011 from 16 health centres in malaria endemics areas in Cambodia. Results: Based on successful typing of 282 samples, subpopulations were characterized along the borders of the country. Each 11-loci barcode provides evidence supporting allele distribution gradient related to subpopulations and gene flow. The 11-loci barcode successfully identifies recently emerging parasite subpopulations in western Cambodia that are associated with the C580Y dominant allele for artemisinin resistance in k13 gene. A subpopulation was identified in northern Cambodia that was associated to artemisinin (R539T resistant allele of k13 gene) and mefloquine resistance. Conclusions: The gene flow between these subpopulations might have driven the spread of artemisinin resistance over Cambodia

Functional analysis of Plasmodium falciparum subpopulations associated with artemisinin resistance in Cambodia

Background: Plasmodium falciparum malaria is one of the most widespread parasitic infections in humans and remains a leading global health concern. Malaria elimination efforts are threatened by the emergence and spread of resistance to artemisinin-based combination therapy, the first-line treatment of malaria. Promising molecular markers and pathways associated with artemisinin drug resistance have been identified, but the underlying molecular mechanisms of resistance remains unknown. The genomic data from early period of emergence of artemisinin resistance (2008–2011) was evaluated, with aim to define k13 associated genetic background in Cambodia, the country identified as epicentre of anti-malarial drug resistance, through characterization of 167 parasite isolates using a panel of 21,257 SNPs. Results: Eight subpopulations were identified suggesting a process of acquisition of artemisinin resistance consistent with an emergence-selection-diffusion model, supported by the shifting balance theory. Identification of population specific mutations facilitated the characterization of a core set of 57 background genes associated with artemisinin resistance and associated pathways. The analysis indicates that the background of artemisinin resistance was not acquired after drug pressure, rather is the result of fixation followed by selection on the daughter subpopulations derived from the ancestral population. Conclusions: Functional analysis of artemisinin resistance subpopulations illustrates the strong interplay between ubiquitination and cell division or differentiation in artemisinin resistant parasites. The relationship of these pathways with the P. falciparum resistant subpopulation and presence of drug resistance markers in addition to k13, highlights the major role of admixed parasite population in the diffusion of artemisinin resistant background. The diffusion of resistant genes in the Cambodian admixed population after selection resulted from mating of gametocytes of sensitive and resistant parasite populations. (Résumé d'auteur

Reactive case-detection of malaria in Pailin Province, Western Cambodia: lessons from a year-long evaluation in a pre-elimination setting.

BACKGROUND: As momentum towards malaria elimination grows, strategies are being developed for scale-up in elimination settings. One prominent strategy, reactive case detection (RACD), involves screening and treating individuals living in close proximity to passively detected, or "index" cases. This study aims to use RACD to quantify Plasmodium parasitaemia in households of index cases, and identify risk factors for infection; these data could inform reactive screening approaches and identify target risk groups. METHODS: This study was conducted in the Western Cambodian province of Pailin between May 2013 and March 2014 among 440 households. Index participants/index cases (n = 270) and surrounding households (n = 110) were screened for Plasmodium infection with rapid diagnostic tests (RDT), microscopy and real-time polymerase chain reaction (PCR). Participants were interviewed to identify risk factors. A comparison group of 60 randomly-selected households was also screened, to compare infection levels of RACD and non-RACD households. In order to identify potential risk factors that would inform screening approaches and identify risk groups, multivariate logistic regression models were applied. RESULTS: Nine infections were identified in households of index cases (RACD approach) through RDT screening of 1898 individuals (seven Plasmodium vivax, two Plasmodium falciparum); seven were afebrile. Seventeen infections were identified through PCR screening of 1596 individuals (15 P. vivax, and 22 % P. falciparum/P. vivax mixed infections). In the control group, 25 P. falciparum infections were identified through PCR screening of 237 individuals, and no P. vivax was found. Plasmodium falciparum infection was associated with fever (p = 0.013), being a member of a control household (p ≤ 0.001), having a history of malaria infection (p = 0.041), and sleeping without a mosquito net (p = 0.011). Significant predictors of P. vivax infection, as diagnosed by PCR, were fever (p = 0.058, borderline significant) and history of malaria infection (p ≤ 0.001). CONCLUSION: This study found that RACD identified very few secondary infections when targeting index and neighbouring households for screening. The results suggest RACD is not appropriate, where exposure to malaria occurs away from the community, and there is a high level of treatment-seeking from the private sector. Piloting RACD in a range of transmission settings would help to identify the ideal environment for feasible and effective reactive screening methods

Plasmodium knowlesi Infection in Humans, Cambodia, 2007–2010

Two cases of Plasmodium knowlesi infection in humans were identified in Cambodia by 3 molecular detection assays and sequencing. This finding confirms the widespread distribution of P. knowlesi malaria in humans in Southeast Asia. Further wide-scale studies are required to assess the public health relevance of this zoonotic malaria parasite

African isolates show a high proportion of multiple copies of the Plasmodium falciparum plasmepsin-2 gene, a piperaquine resistance marker

Background: Today, the development of new and well-tolerated anti-malarial drugs is strongly justifed by the emer‑ gence of Plasmodium falciparum resistance. In 2014–2015, a phase 2b clinical study was conducted to evaluate the efcacy of a single oral dose of Artefenomel (OZ439)–piperaquine (PPQ) in Asian and African patients presenting with uncomplicated falciparum malaria. Methods: Blood samples collected before treatment ofered the opportunity to investigate the proportion of mul‑ tidrug resistant parasite genotypes, including P. falciparum kelch13 mutations and copy number variation of both P. falciparum plasmepsin 2 (Pfpm2) and P. falciparum multidrug resistance 1 (Pfmdr1) genes. Results: Validated kelch13 resistance mutations including C580Y, I543T, P553L and V568G were only detected in parasites from Vietnamese patients. In Africa, isolates with multiple copies of the Pfmdr1 gene were shown to be more frequent than previously reported (21.1%, range from 12.4% in Burkina Faso to 27.4% in Uganda). More strikingly, high proportions of isolates with multiple copies of the Pfpm2 gene, associated with piperaquine (PPQ) resistance, were frequently observed in the African sites, especially in Burkina Faso and Uganda (>30%). Conclusions: These fndings were considered to sharply contrast with the recent description of increased sensitivity to PPQ of Ugandan parasite isolates. This emphasizes the necessity to investigate in vitro susceptibility profles to PPQ of African isolates with multiple copies of the Pfpm2 gene and estimate the risk of development of PPQ resistance in Africa

Field-based evidence of fast and global increase of Plasmodium falciparum drug-resistance by DNA-microarrays and PCR/RFLP in Niger

<p>Abstract</p> <p>Background</p> <p>Over the last years, significant progress has been made in the comprehension of the molecular mechanism of malaria resistance to drugs. Together with <it>in vivo </it>tests, the molecular monitoring is now part of the survey strategy of the <it>Plasmodium </it>sensitivity. Currently, DNA-microarray analysis allows the simultaneous study of many single nucleotide polymorphisms (SNP) of <it>Plasmodium </it>isolates. In December 2005, the International Federation of the Red Cross distributed two million three hundred thousand long-lasting insecticide nets to pregnant women and mothers of under five years children in the whole Niger. Then, Niger adopted artemisinin-based combination therapy as first-line treatment.</p> <p>Methods</p> <p>Thirty four SNPs of <it>pfcrt, pfdhfr, pfdhps, pfmdr </it>and <it>pfATPase </it>were analysed by DNA-microarray and PCR/RFLP in two villages – Zindarou and Banizoumbou – with different durations of malaria transmission. The main objective of the study was to measure the dynamics <it>of Plasmodium falciparum </it>resistant strains and associated factors.</p> <p>Results</p> <p>This study shows a global and clear increase of the drug-resistance associated molecular markers frequencies during a relatively short-time period of four years. Markers associated with resistance to chloroquine and sulphonamids were more frequently found in the short transmission zone than in the long transmission one. The <it>pfcrt76T </it>mutation is significantly more present at Banizoumbou than Zindarou (38.3% vs 25.2%, p = 0.013).</p> <p>This work allowed the screening of several field strains for five SNPs of <it>PfATPase6 </it>gene. The <it>pfATPase6S769N</it>, candidate mutation of resistance to artemisinin was not found. However the <it>pfATPsaeA623E </it>mutation was found in 4.7% of samples.</p> <p>Conclusion</p> <p>A significant increase of several SNPs frequencies was highlighted over a four-year period. The polymorphism of five <it>PfATPase6 </it>gene SNPs was described. The global, large and fast increase of the molecular resistance is discussed in the context of current changes of health policy and malaria control in Niger.</p

Spatial clustering and risk factors of malaria infections in Ratanakiri Province, Cambodia

Background: Malaria incidence worldwide has steadily declined over the past decades. Consequently, increasingly more countries will proceed from control to elimination. The malaria distribution in low incidence settings appears patchy, and local transmission hotspots are a continuous source of infection. In this study, species-specific clusters and associated risk factors were identified based on malaria prevalence data collected in the north-east of Cambodia. In addition, Plasmodium falciparum genetic diversity, population structure and gene flows were studied.Method: In 2012, blood samples from 5793 randomly selected individuals living in 117 villages were collected from Ratanakiri province, Cambodia. Malariometric data of each participant were simultaneously accumulated using a standard questionnaire. A two-step PCR allowed for species-specific detection of malaria parasites, and SNPgenotyping of P. falciparum was performed. SaTScan was used to determine species-specific areas of elevated risk to infection, and univariate and multivariate risk analyses were carried out.Result: PCR diagnosis found 368 positive individuals (6.4%) for malaria parasites, of which 22% contained mixed species infections. The occurrence of these co-infections was more frequent than expected. Specific areas with elevated risk of infection were detected for all Plasmodium species. The clusters for Falciparum, Vivax and Ovale malaria appeared in the north of the province along the main river, while the cluster for Malariae malaria was situated elsewhere. The relative risk to be a malaria parasite carrier within clusters along the river was twice that outside the area. The main risk factor associated with three out of four malaria species was overnight stay in the plot hut, a human behaviour associated with indigenous farming. Haplotypes did not show clear geographical population structure, but pairwise Fst value comparison indicated higher parasite flow along the river.Discussion: Spatial aggregation of malaria parasite carriers, and the identification of malaria species-specific risk factors provide key insights in malaria epidemiology in low transmission settings, which can guide targeted supplementary interventions. Consequently, future malaria programmes in the province should implement additional specific policies targeting households staying overnight at their farms outside the village, in addition to migrants and forest workers