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.

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.

SummaryBackgroundWestern 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.MethodsWe 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.FindingsWhole-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%.InterpretationPiperaquine 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.FundingInstitut Pasteur in Cambodia, Institut Pasteur Paris, National Institutes of Health, WHO, Agence Nationale de la Recherche, Investissement d’Avenir programme, Laboratoire d’Excellence In...

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

BackgroundThe 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.MethodsThe 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).ResultsAll 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).ConclusionPSA 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.Electronic supplementary materialThe online version of this article (doi:10.1186/s12916-015-0539-5) contains supplementary material, which is available to authorized users.

Tolerance is not detected by current assays and represents a major threat to antimalarial drug policy.

We present here the first genome-wide characterization of linkage disequilibrium (LD) in the French Blonde d'Aquitaine (BLA) breed, a well-muscled breed renowned for producing high-yielding beef carcasses. To assess the pattern and extent of LD, we used a sample of 30 unrelated bulls and 36 923 single nucleotide polymorphisms (SNPs) covering all cattle autosomes. The squared correlation of the alleles at two loci (r(2) ) was used as a measure of LD. The analysis of adjacent marker pairs revealed that the level of LD decreases rapidly with physical distance between SNPs. Overall mean r(2) was 0.205 (±0.262). Strong LD (r(2) > 0.8) and useful LD (measured as r(2 ) > 0.2) were observed within genomic regions of up to 720 and 724 kb, respectively. We analysed the genetic structure of the BLA population and found stratification. The observed genetic sub-structuring is consistent with the known recent demographic history that occurred during BLA breed formation. Our results indicate that LD mapping of phenotypic traits in the BLA population is feasible; however, because of this sub-structuring, special care is needed to reduce the likelihood of false-positive associations between marker loci and traits of interest.

MotivationMost 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.ResultsWe 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 ImplementationThe R-package LEANR implementing LEAN is supplied as supplementary material and available on CRAN (https://cran.r-project.org).Supplementary information Supplementary data are available at Bioinformatics online.