ClermonTyping: an easy-to-use and accurate in silico method for Escherichia genus strain phylotyping.

The genus Escherichia is composed of Escherichia albertii, E. fergusonii, five cryptic Escherichia clades and E. coli sensu stricto. Furthermore, the E. coli species can be divided into seven main phylogroups termed A, B1, B2, C, D, E and F. As specific lifestyles and/or hosts can be attributed to these species/phylogroups, their identification is meaningful for epidemiological studies. Classical phenotypic tests fail to identify non-sensu stricto E. coli as well as phylogroups. Clermont and colleagues have developed PCR assays that allow the identification of most of these species/phylogroups, the triplex/quadruplex PCR for E. coli phylogroup determination being the most popular. With the growing availability of whole genome sequences, we have developed the ClermonTyping method and its associated web-interface, the ClermonTyper, that allows a given strain sequence to be assigned to E. albertii, E. fergusonii, Escherichia clades I-V, E. coli sensu stricto as well as to the seven main E. coli phylogroups. The ClermonTyping is based on the concept of in vitro PCR assays and maintains the principles of ease of use and speed that prevailed during the development of the in vitro assays. This in silico approach shows 99.4 % concordance with the in vitro PCR assays and 98.8 % with the Mash genome-clustering tool. The very few discrepancies result from various errors occurring mainly from horizontal gene transfers or SNPs in the primers. We propose the ClermonTyper as a freely available resource to the scientific community at: http://clermontyping.iame-research.center/.

A clone of the emergent Streptococcus pyogenes emm89 clade responsible for a large outbreak in a post-surgery oncology unit in France.

An outbreak of nosocomial infections due to Streptococcus pyogenes (Group A Streptococcus; GAS) occurred in a post-surgery oncology unit and concerned more than 60 patients and lasted 20 months despite enhanced infection control and prophylaxis measures. All GAS strains were characterized (emm genotype, toxin gene profile and pulse-field gel electrophoresis subtype). Selected strains were sequenced and phylogenetic relationship established. Capacity to form biofilm and interaction with human pulmonary epithelial cells and macrophages were determined. Twenty-six GAS strains responsible for invasive infections (II) and 57 for non-II or colonization were isolated from patients (n = 66) or healthcare workers (n = 13). Seventy strains shared the same molecular markers and 69 the same PFGE pattern; 56 were sequenced. They all belonged to the emerging emm89 clade 3; all but 1 were clonal. Whole genome sequencing identified 43 genetic profiles with sporadic mutations in regulatory genes and acquired mutations in 2 structural genes. Except for two regulatory gene mutants, all strains tested had the same biofilm formation capacity and displayed similar adherence and invasion of pulmonary epithelial cells and phagocytosis and survival in human macrophages. This large outbreak of GAS infection in a post-surgery oncology unit, a setting that contains highly susceptible patients, arose from a strain of the emergent emm89 clade. No relationship between punctual or acquired mutations, invasive status, and strain phenotypic characteristics was found. Noteworthy, the phenotypic characteristics of this clone account for its emergence and its remarkable capacity to elicit outbreaks.

Network-based analysis of omics data: the LEAN method.

Motivation: Most computational approaches for the analysis of omics data in the context of interaction networks have very long running times, provide single or partial, often heuristic, solutions and/or contain user-tuneable parameters. Results: We introduce local enrichment analysis (LEAN) for the identification of dysregulated subnetworks from genome-wide omics datasets. By substituting the common subnetwork model with a simpler local subnetwork model, LEAN allows exact, parameter-free, efficient and exhaustive identification of local subnetworks that are statistically dysregulated, and directly implicates single genes for follow-up experiments.Evaluation on simulated and biological data suggests that LEAN generally detects dysregulated subnetworks better, and reflects biological similarity between experiments more clearly than standard approaches. A strong signal for the local subnetwork around Von Willebrand Factor (VWF), a gene which showed no change on the mRNA level, was identified by LEAN in transcriptome data in the context of the genetic disease Cerebral Cavernous Malformations (CCM). This signal was experimentally found to correspond to an unexpected strong cellular effect on the VWF protein. LEAN can be used to pinpoint statistically significant local subnetworks in any genome-scale dataset. Availability and Implementation: The R-package LEANR implementing LEAN is supplied as supplementary material and available on CRAN ( https://cran.r-project.org ). Contacts: benno@pasteur.fr or tournier-lasserve@univ-paris-diderot.fr. Supplementary information: Supplementary data are available at Bioinformatics online.

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".

A Worldwide Map of Plasmodium falciparum K13-Propeller Polymorphisms.

BACKGROUND: Recent gains in reducing the global burden of malaria are threatened by the emergence of Plasmodium falciparum resistance to artemisinins. The discovery that mutations in portions of a P. falciparum gene encoding kelch (K13)-propeller domains are the major determinant of resistance has provided opportunities for monitoring such resistance on a global scale. METHODS: We analyzed the K13-propeller sequence polymorphism in 14,037 samples collected in 59 countries in which malaria is endemic. Most of the samples (84.5%) were obtained from patients who were treated at sentinel sites used for nationwide surveillance of antimalarial resistance. We evaluated the emergence and dissemination of mutations by haplotyping neighboring loci. RESULTS: We identified 108 nonsynonymous K13 mutations, which showed marked geographic disparity in their frequency and distribution. In Asia, 36.5% of the K13 mutations were distributed within two areas--one in Cambodia, Vietnam, and Laos and the other in western Thailand, Myanmar, and China--with no overlap. In Africa, we observed a broad array of rare nonsynonymous mutations that were not associated with delayed parasite clearance. The gene-edited Dd2 transgenic line with the A578S mutation, which expresses the most frequently observed African allele, was found to be susceptible to artemisinin in vitro on a ring-stage survival assay. CONCLUSIONS: No evidence of artemisinin resistance was found outside Southeast Asia and China, where resistance-associated K13 mutations were confined. The common African A578S allele was not associated with clinical or in vitro resistance to artemisinin, and many African mutations appear to be neutral. (Funded by Institut Pasteur Paris and others.).

Plasmodium copy number variation scan: gene copy numbers evaluation in haploid genomes.

BACKGROUND: In eukaryotic genomes, deletion or amplification rates have been estimated to be a thousand more frequent than single nucleotide variation. In Plasmodium falciparum, relatively few transcription factors have been identified, and the regulation of transcription is seemingly largely influenced by gene amplification events. Thus copy number variation (CNV) is a major mechanism enabling parasite genomes to adapt to new environmental changes. METHODS: Currently, the detection of CNVs is based on quantitative PCR (qPCR), which is significantly limited by the relatively small number of genes that can be analysed at any one time. Technological advances that facilitate whole-genome sequencing, such as next generation sequencing (NGS) enable deeper analyses of the genomic variation to be performed. Because the characteristics of Plasmodium CNVs need special consideration in algorithms and strategies for which classical CNV detection programs are not suited a dedicated algorithm to detect CNVs across the entire exome of P. falciparum was developed. This algorithm is based on a custom read depth strategy through NGS data and called PlasmoCNVScan. RESULTS: The analysis of CNV identification on three genes known to have different levels of amplification and which are located either in the nuclear, apicoplast or mitochondrial genomes is presented. The results are correlated with the qPCR experiments, usually used for identification of locus specific amplification/deletion. CONCLUSIONS: This tool will facilitate the study of P. falciparum genomic adaptation in response to ecological changes: drug pressure, decreased transmission, reduction of the parasite population size (transition to pre-elimination endemic area).

Plasmodium falciparum dihydroartemisinin-piperaquine failures in Cambodia are associated with mutant K13 parasites presenting high survival rates in novel piperaquine in vitro assays: retrospective and prospective investigations.

BACKGROUND: The declining efficacy of dihydroartemisinin-piperaquine against Plasmodium falciparum in Cambodia, along with increasing numbers of recrudescent cases, suggests resistance to both artemisinin and piperaquine. Available in vitro piperaquine susceptibility assays do not correlate with treatment outcome. A novel assay using a pharmacologically relevant piperaquine dose/time exposure was designed and its relevance explored in retrospective and prospective studies. METHODS: The piperaquine survival assay (PSA) exposed parasites to 200 nM piperaquine for 48 hours and monitored survival 24 hours later. The retrospective study tested 32 culture-adapted, C580Y-K13 mutant parasites collected at enrolment from patients treated with a 3-day course of dihydroartemisinin-piperaquine and having presented or not with a recrudescence at day 42 (registered ACTRN12615000793516). The prospective study assessed ex vivo PSA survival rate alongside K13 polymorphism of isolates collected from patients enrolled in an open-label study with dihydroartemisinin-piperaquine for uncomplicated P. falciparum malaria in Cambodia (registered ACTRN12615000696594). RESULTS: All parasites from recrudescent cases had in vitro or ex vivo PSA survival rates ≥10%, a relevant cut-off value for piperaquine-resistance. Ex vivo PSA survival rates were higher for recrudescent than non-recrudescent cases (39.2% vs. 0.17%, P <1 × 10(-7)). Artemisinin-resistant K13 mutants with ex vivo PSA survival rates ≥10% were associated with 32-fold higher risk of recrudescence (95% CI, 4.5-224; P = 0.0005). CONCLUSION: PSA adequately captures the piperaquine resistance/recrudescence phenotype, a mainstay to identify molecular marker(s) and evaluate efficacy of alternative drugs. Combined ex vivo PSA and K13 genotyping provides a convenient monitor for both artemisinin and piperaquine resistance where dihydroartemisinin-piperaquine is used.

Induction of Multidrug Tolerance in Plasmodium falciparum by Extended Artemisinin Pressure.

Plasmodium falciparum resistance to artemisinin derivatives in Southeast Asia threatens global malaria control strategies. Whether delayed parasite clearance, which exposes larger parasite numbers to artemisinins for longer times, selects higher-grade resistance remains unexplored. We investigated whether long-lasting artemisinin pressure selects a novel multidrug-tolerance profile. Although 50% inhibitory concentrations for 10 antimalarial drugs tested were unchanged, drug-tolerant parasites showed higher recrudescence rates for endoperoxides, quinolones, and an antifolate, including partner drugs of recommended combination therapies, but remained susceptible to atovaquone. Moreover, the age range of intraerythrocytic stages able to resist artemisinin was extended to older ring forms and trophozoites. Multidrug tolerance results from drug-induced quiescence, which enables parasites to survive exposure to unrelated antimalarial drugs that inhibit a variety of metabolic pathways. This novel resistance pattern should be urgently monitored in the field because this pattern is not detected by current assays and represents a major threat to antimalarial drug policy.

A molecular marker of artemisinin-resistant Plasmodium falciparum malaria.

Plasmodium falciparum resistance to artemisinin derivatives in southeast Asia threatens malaria control and elimination activities worldwide. To monitor the spread of artemisinin resistance, a molecular marker is urgently needed. Here, using whole-genome sequencing of an artemisinin-resistant parasite line from Africa and clinical parasite isolates from Cambodia, we associate mutations in the PF3D7_1343700 kelch propeller domain ('K13-propeller') with artemisinin resistance in vitro and in vivo. Mutant K13-propeller alleles cluster in Cambodian provinces where resistance is prevalent, and the increasing frequency of a dominant mutant K13-propeller allele correlates with the recent spread of resistance in western Cambodia. Strong correlations between the presence of a mutant allele, in vitro parasite survival rates and in vivo parasite clearance rates indicate that K13-propeller mutations are important determinants of artemisinin resistance. K13-propeller polymorphism constitutes a useful molecular marker for large-scale surveillance efforts to contain artemisinin resistance in the Greater Mekong Subregion and prevent its global spread.