Indels, structural variation, and recombination drive genomic diversity in Plasmodium falciparum

Co-transmission of Related Malaria Parasite Lineages Shapes Within-Host Parasite Diversity

In high-transmission regions, we expect parasite lineages within complex malaria infections to be unrelated due to parasite inoculations from different mosquitoes. This project was designed to test this prediction. We generated 485 single-cell genome sequences from fifteen P. falciparum malaria patients from Chikhwawa, Malawi-an area of intense transmission. Patients harbored up to seventeen unique parasite lineages. Surprisingly, parasite lineages within infections tend to be closely related, suggesting that superinfection by repeated mosquito bites is rarer than co-transmission of parasites from a single mosquito. Both closely and distantly related parasites comprise an infection, suggesting sequential transmission of complex infections between multiple hosts. We identified tetrads and reconstructed parental haplotypes, which revealed the inbred ancestry of infections and non-Mendelian inheritance. Our analysis suggests strong barriers to secondary infection and outbreeding amongst malaria parasites from a high transmission setting, providing unexpected insights into the biology and transmission of malaria.

Polymorphism and natural selection in the merozoite surface protein 3F2 (PVX_97710) locus of Plasmodium vivax among field isolates

Plasmodium vivax, the chronic relapsing human malaria parasite with the most widespread distribution, possesses proteins associated with the merozoite surface that could be targets for host immune responses and potential vaccine candidates. Of these, the merozoite surface protein 3 of P. vivax (PvMSP3) is an attractive vaccine target as well as a genetic marker for epidemiological surveillance. PvMSP3 comprises a group of protein members encoded by a multigene family. Although some protein members, i.e. PvMSP3α and PvMSP3β, have been targets for molecular and immunological investigations, the most abundantly expressed protein member during late asexual erythrocytic stages, PvMSP3F2 (PVX_97710), remains unexplored. To address domain organization and evolution of this locus, the complete coding sequences of 31 P. vivax isolates from diverse malaria endemic areas of Thailand were analyzed and compared with 10 previously reported sequences. Results revealed that all PvMSP3F2 sequences differed but could be divided into 5 repeat-containing domains flanked by 6 non-repeat domains. Repeat domains II and IV at the 5' portion and domain X at the 3' portion exhibited extensive sequence and length variation whereas repeat domains VI and VIII located at the central region were relatively conserved. Despite a repertoire of PvMSP3F2 variants, predicted coiled-coil tertiary structure and predicted B-cell epitopes seem to be maintained. Evidence of intragenic recombination has been detected among field isolates in Thailand that could enhance sequence diversity at this locus. Non-repeat domains I and IX located at the 5' end and at the 3' portion, respectively, seem to have evolved under purifying selection. Evidence of positive selection was found in non-repeat domains III, V and VII where a number of predicted HLA class I epitopes were identified. Amino acid substitutions in these predicted epitopes could alter predicted peptide binding affinity or abolish peptide epitope property, suggesting that polymorphism in these epitopes conferred host immune evasion. Further studies on PvMSP3F2 are warranted, particularly on interaction with host immune system and the potential role of this PvMSP3 protein member as a vaccine target.

Genetic mapping of fitness determinants across the malaria parasite Plasmodium falciparum life cycle

Determining the genetic basis of fitness is central to understanding evolution and transmission of microbial pathogens. In human malaria parasites (Plasmodium falciparum), most experimental work on fitness has focused on asexual blood stage parasites, because this stage can be easily cultured, although the transmission of malaria requires both female Anopheles mosquitoes and vertebrate hosts. We explore a powerful approach to identify the genetic determinants of parasite fitness across both invertebrate and vertebrate life-cycle stages of P. falciparum. This combines experimental genetic crosses using humanized mice, with selective whole genome amplification and pooled sequencing to determine genome-wide allele frequencies and identify genomic regions under selection across multiple lifecycle stages. We applied this approach to genetic crosses between artemisinin resistant (ART-R, kelch13-C580Y) and ART-sensitive (ART-S, kelch13-WT) parasites, recently isolated from Southeast Asian patients. Two striking results emerge: we observed (i) a strong genome-wide skew (>80%) towards alleles from the ART-R parent in the mosquito stage, that dropped to ~50% in the blood stage as selfed ART-R parasites were selected against; and (ii) repeatable allele specific skews in blood stage parasites with particularly strong selection (selection coefficient (s) ≤ 0.18/asexual cycle) against alleles from the ART-R parent at loci on chromosome 12 containing MRP2 and chromosome 14 containing ARPS10. This approach robustly identifies selected loci and has strong potential for identifying parasite genes that interact with the mosquito vector or compensatory loci involved in drug resistance.

Plasmodium vivax chloroquine resistance links to pvcrt transcription in a genetic cross

Mainstay treatment for Plasmodium vivax malaria has long relied on chloroquine (CQ) against blood-stage parasites plus primaquine against dormant liver-stage forms (hypnozoites), however drug resistance confronts this regimen and threatens malaria control programs. Understanding the basis of P. vivax chloroquine resistance (CQR) will inform drug discovery and malaria control. Here we investigate the genetics of P. vivax CQR by a cross of parasites differing in drug response. Gametocytogenesis, mosquito infection, and progeny production are performed with mixed parasite populations in nonhuman primates, as methods for P. vivax cloning and in vitro cultivation remain unavailable. Linkage mapping of progeny surviving >15 mg/kg CQ identifies a 76 kb region in chromosome 1 including pvcrt, an ortholog of the Plasmodium falciparum CQR transporter gene. Transcriptional analysis supports upregulated pvcrt expression as a mechanism of CQR.

Immunogenic Evaluation of Ribosomal P-Protein Antigen P0, P1, and P2 and Pentameric Protein Complex P0-(P1-P2)2 of Plasmodium falciparum in a Mouse Model

Malaria remains one the most infectious and destructive protozoan diseases worldwide. Plasmodium falciparum, a protozoan parasite with a complex life cycle and high genetic variability responsible for the difficulties in vaccine development, is implicated in most malaria-related deaths. In the course of study, we prepared a set of antigens based on P-proteins from P. falciparum and determined their immunogenicity in an in vivo assay on a mouse model. The pentameric complex P0-(P1-P2)₂ was prepared along with individual P1, P2, and P0 antigens. We determined the level of cellular- and humoral-type immunological response followed by development of specific immunological memory. We have shown that the number of Tc cells increased significantly after the first immunization with P2 and after the second immunization with P1 and P0-(P1-P2)₂, which highly correlated with the number of Th1 cells. P0 appeared as a poor inducer of cellular response. After the third boost with P1, P2, or P0-(P1-P2)₂, the initially high cellular response dropped to the control level accompanied by elevation of the number of activated Treg cells and a high level of suppressive TGF-β. Subsequently, the humoral response against the examined antigens was activated. Although the titers of specific IgG were increasing during the course of immunization for all antigens used, P2 and P0-(P1-P2)₂ were found to be significantly stronger than P1 and P0. A positive correlation between the Th2 cell abundance and the level of IL-10 was observed exclusively after immunization with P0-(P1-P2)₂. An in vitro exposure of spleen lymphocytes from the immunized mice especially to the P1, P2, and P0-(P1-P2)₂ protein caused 2-3-fold higher cell proliferation than that in the case of lymphocytes from the nonimmunized animals, suggesting development of immune memory. Our results demonstrate for the first time that the native-like P-protein pentameric complex represents much stronger immune potential than individual P-antigens.

Estimating Relatedness Between Malaria Parasites

Understanding the relatedness of individuals within or between populations is a common goal in biology. Increasingly, relatedness features in genetic epidemiology studies of pathogens. These studies are relatively new compared to those in humans and other organisms, but are important for designing interventions and understanding pathogen transmission. Only recently have researchers begun to routinely apply relatedness to apicomplexan eukaryotic malaria parasites, and to date have used a range of different approaches on an ad hoc basis. Therefore, it remains unclear how to compare different studies and which measures to use. Here, we systematically compare measures based on identity-by-state (IBS) and identity-by-descent (IBD) using a globally diverse data set of malaria parasites, Plasmodium falciparum and P. vivax, and provide marker requirements for estimates based on IBD. We formally show that the informativeness of polyallelic markers for relatedness inference is maximized when alleles are equifrequent. Estimates based on IBS are sensitive to allele frequencies, which vary across populations and by experimental design. For portability across studies, we thus recommend estimates based on IBD. To generate estimates with errors below an arbitrary threshold of 0.1, we recommend ∼100 polyallelic or 200 biallelic markers. Marker requirements are immediately applicable to haploid malaria parasites and other haploid eukaryotes. C.I.s facilitate comparison when different marker sets are used. This is the first attempt to provide rigorous analysis of the reliability of, and requirements for, relatedness inference in malaria genetic epidemiology. We hope it will provide a basis for statistically informed prospective study design and surveillance strategies.

Plasmodium Genomics and Genetics: New Insights into Malaria Pathogenesis, Drug Resistance, Epidemiology, and Evolution

Protozoan Plasmodium parasites are the causative agents of malaria, a deadly disease that continues to afflict hundreds of millions of people every year. Infections with malaria parasites can be asymptomatic, with mild or severe symptoms, or fatal, depending on many factors such as parasite virulence and host immune status. Malaria can be treated with various drugs, with artemisinin-based combination therapies (ACTs) being the first-line choice. Recent advances in genetics and genomics of malaria parasites have contributed greatly to our understanding of parasite population dynamics, transmission, drug responses, and pathogenesis. However, knowledge gaps in parasite biology and host-parasite interactions still remain. Parasites resistant to multiple antimalarial drugs have emerged, while advanced clinical trials have shown partial efficacy for one available vaccine. Here we discuss genetic and genomic studies of Plasmodium biology, host-parasite interactions, population structures, mosquito infectivity, antigenic variation, and targets for treatment and immunization. Knowledge from these studies will advance our understanding of malaria pathogenesis, epidemiology, and evolution and will support work to discover and develop new medicines and vaccines.

Whole-genome sequencing reveals high complexity of copy number variation at insecticide resistance loci in malaria mosquitoes

Polymorphisms in genetic copy number can influence gene expression, coding sequence, and zygosity, making them powerful actors in the evolutionary process. Copy number variants (CNVs) are however understudied, being more difficult to detect than single-nucleotide polymorphisms. We take advantage of the intense selective pressures on the major malaria vector Anopheles gambiae, caused by the widespread use of insecticides for malaria control, to investigate the role of CNVs in the evolution of insecticide resistance. Using the whole-genome sequencing data from 1142 samples in the An. gambiae 1000 genomes project, we identified 250 gene-containing CNVs, encompassing a total of 267 genes of which 28 were in gene families linked to metabolic insecticide resistance, representing significant enrichment of these families. The five major gene clusters for metabolic resistance all contained CNVs, with 44 different CNVs being found across these clusters and multiple CNVs frequently covering the same genes. These 44 CNVs are widespread (45% of individuals carry at least one of them) and have been spreading through positive selection, indicated by their high local frequencies and extended haplotype homozygosity. Our results demonstrate the importance of CNVs in the response to selection, highlighting the urgent need to identify the contribution of each CNV to insecticide resistance and to track their spread as the use of insecticides in malaria endemic countries intensifies and as the operational deployment of next-generation bed nets targeting metabolic resistance gathers pace. Our detailed descriptions of CNVs found across the species range provide the tools to do so.

The origins and relatedness structure of mixed infections vary with local prevalence of P. falciparum malaria

Individual malaria infections can carry multiple strains of Plasmodium falciparum with varying levels of relatedness. Yet, how local epidemiology affects the properties of such mixed infections remains unclear. Here, we develop an enhanced method for strain deconvolution from genome sequencing data, which estimates the number of strains, their proportions, identity-by-descent (IBD) profiles and individual haplotypes. Applying it to the Pf3k data set, we find that the rate of mixed infection varies from 29% to 63% across countries and that 51% of mixed infections involve more than two strains. Furthermore, we estimate that 47% of symptomatic dual infections contain sibling strains likely to have been co-transmitted from a single mosquito, and find evidence of mixed infections propagated over successive infection cycles. Finally, leveraging data from the Malaria Atlas Project, we find that prevalence correlates within Africa, but not Asia, with both the rate of mixed infection and the level of IBD.

De Novo Mutations Resolve Disease Transmission Pathways in Clonal Malaria

Detecting de novo mutations in viral and bacterial pathogens enables researchers to reconstruct detailed networks of disease transmission and is a key technique in genomic epidemiology. However, these techniques have not yet been applied to the malaria parasite, Plasmodium falciparum, in which a larger genome, slower generation times, and a complex life cycle make them difficult to implement. Here, we demonstrate the viability of de novo mutation studies in P. falciparum for the first time. Using a combination of sequencing, library preparation, and genotyping methods that have been optimized for accuracy in low-complexity genomic regions, we have detected de novo mutations that distinguish nominally identical parasites from clonal lineages. Despite its slower evolutionary rate compared with bacterial or viral species, de novo mutation can be detected in P. falciparum across timescales of just 1–2 years and evolutionary rates in low-complexity regions of the genome can be up to twice that detected in the rest of the genome. The increased mutation rate allows the identification of separate clade expansions that cannot be found using previous genomic epidemiology approaches and could be a crucial tool for mapping residual transmission patterns in disease elimination campaigns and reintroduction scenarios.

An analysis of large structural variation in global Plasmodium falciparum isolates identifies a novel duplication of the chloroquine resistance associated gene

The evolution of genetic mechanisms for host immune evasion and anti-malarial resistance has enabled the Plasmodium falciparum malaria parasite to inflict high morbidity and mortality on human populations. Most studies of P. falciparum genetic diversity have focused on single-nucleotide polymorphisms (SNPs), assisting the identification of drug resistance-associated loci such as the chloroquine related crt and sulfadoxine-pyrimethamine related dhfr. Whilst larger structural variants are known to impact adaptation, for example, mdr1 duplications with anti-malarial resistance, no large-scale, genome-wide study on clinical isolates has been undertaken using whole genome sequencing data. By applying a structural variant detection pipeline across whole genome sequence data from 2,855 clinical isolates in 21 malaria-endemic countries, we identified >70,000 specific deletions and >600 duplications. Most structural variants are rare (48.5% of deletions and 94.7% of duplications are found in single isolates) with 2.4% of deletions and 0.2% of duplications found in >5% of global isolates. A subset of variants was present at high frequency in drug-resistance related genes including mdr1, the gch1 promoter region, and a putative novel duplication of crt. Regional-specific variants were identified, and a companion visualisation tool has been developed to assist web-based investigation of these polymorphisms by the wider scientific community.

Whole-Genome Resequencing Identifies KIT New Alleles That Affect Coat Color Phenotypes in Pigs

The Duroc × (Landrace × Large White) hybrid pig (DLY) is the most popular commercial pig used in the Chinese pig industry. DLY pigs are usually white but sometimes show colored phenotypes. Colored DLY pigs are not favored by slaughterhouses and retailers, thus causing certain economic losses to farmers in China. In this study, we first conducted a genome-wide association study and RNA sequencing to demonstrate that KIT variants are responsible for diversifying coat color phenotypes segregating in a DLY population. We then defined the precise sizes and locations of four duplications (DUP1-4), four candidate causative mutations at the KIT locus, in the pig reference genome using the whole-genome sequence data of representative colored individuals. The sequence data also enabled us to identify a list of new KIT alleles. By investigating the association between these new alleles and coat color phenotypes, we provide further evidence that DUP2 is another causative mutation for the solid white coat color in pigs. DUP1 (the KIT gene duplication), DUP2 and the splice mutation are all required for the manifestation of a solid white coat color. DUP4 had a more significant effect on the formation of the belt phenotype compared with DUP3. Given the necessity of DUP2 for the solid white coat color, we detected I^N/I^N homozygotes lacking DUP2 in Large White and Landrace pigs and found that French Landrace pigs had the highest frequency (8.98%) of I^N/I^N individuals. This study not only advances our understanding of the molecular mechanism of the color phenotype in pigs, but also establishes a simple and accurate method for the screening of KIT I^N/I^N homozygotes in Large White and Landrace that would cause colored DLY pigs.

A forward genetic screen reveals a primary role for Plasmodium falciparum Reticulocyte Binding Protein Homologue 2a and 2b in determining alternative erythrocyte invasion pathways

Invasion of human erythrocytes is essential for Plasmodium falciparum parasite survival and pathogenesis, and is also a complex phenotype. While some later steps in invasion appear to be invariant and essential, the earlier steps of recognition are controlled by a series of redundant, and only partially understood, receptor-ligand interactions. Reverse genetic analysis of laboratory adapted strains has identified multiple genes that when deleted can alter invasion, but how the relative contributions of each gene translate to the phenotypes of clinical isolates is far from clear. We used a forward genetic approach to identify genes responsible for variable erythrocyte invasion by phenotyping the parents and progeny of previously generated experimental genetic crosses. Linkage analysis using whole genome sequencing data revealed a single major locus was responsible for the majority of phenotypic variation in two invasion pathways. This locus contained the PfRh2a and PfRh2b genes, members of one of the major invasion ligand gene families, but not widely thought to play such a prominent role in specifying invasion phenotypes. Variation in invasion pathways was linked to significant differences in PfRh2a and PfRh2b expression between parasite lines, and their role in specifying alternative invasion was confirmed by CRISPR-Cas9-mediated genome editing. Expansion of the analysis to a large set of clinical P. falciparum isolates revealed common deletions, suggesting that variation at this locus is a major cause of invasion phenotypic variation in the endemic setting. This work has implications for blood-stage vaccine development and will help inform the design and location of future large-scale studies of invasion in clinical isolates.

Plasmodium parasites cause more than 200 million cases of malaria each year. All the symptoms of malaria are caused after Plasmodium parasites invade human red blood cells. Once inside, they grow, multiply and break open the red blood cells to release new parasites. This cycle is repeated every 48 hours, rapidly amplifying the number of parasites and causing severe anemia and other complications. Plasmodium falciparum, the parasite species responsible for almost all malaria deaths, can use multiple different pathways to invade human red blood cells, but the relative importance of each is not well understood. We tested the invasion pathways used by a collection of closely related parasites and compared their genome sequences to identify the genes responsible. This analysis revealed that expression differences in two neighboring genes of the Reticulocyte Binding Homologue family are responsible for most of the variation in two invasion pathways. P. falciparum may use variation in these genes to avoid the immune system or adapt to specific blood groups, which has important implications for vaccine development against malaria.

Highlights on the Application of Genomics and Bioinformatics in the Fight Against Infectious Diseases: Challenges and Opportunities in Africa

Genomics and bioinformatics are increasingly contributing to our understanding of infectious diseases caused by bacterial pathogens such as Mycobacterium tuberculosis and parasites such as Plasmodium falciparum. This ranges from investigations of disease outbreaks and pathogenesis, host and pathogen genomic variation, and host immune evasion mechanisms to identification of potential diagnostic markers and vaccine targets. High throughput genomics data generated from pathogens and animal models can be combined with host genomics and patients’ health records to give advice on treatment options as well as potential drug and vaccine interactions. However, despite accounting for the highest burden of infectious diseases, Africa has the lowest research output on infectious disease genomics. Here we review the contributions of genomics and bioinformatics to the management of infectious diseases of serious public health concern in Africa including tuberculosis (TB), dengue fever, malaria and filariasis. Furthermore, we discuss how genomics and bioinformatics can be applied to identify drug and vaccine targets. We conclude by identifying challenges to genomics research in Africa and highlighting how these can be overcome where possible.

Multi-population genomic analysis of malaria parasites indicates local selection and differentiation at the gdv1 locus regulating sexual development

Parasites infect hosts in widely varying environments, encountering diverse challenges for adaptation. To identify malaria parasite genes under locally divergent selection across a large endemic region with a wide spectrum of transmission intensity, genome sequences were obtained from 284 clinical Plasmodium falciparum infections from four newly sampled locations in Senegal, The Gambia, Mali and Guinea. Combining these with previous data from seven other sites in West Africa enabled a multi-population analysis to identify discrete loci under varying local selection. A genome-wide scan showed the most exceptional geographical divergence to be at the early gametocyte gene locus gdv1 which is essential for parasite sexual development and transmission. We identified a major structural dimorphism with alternative 1.5 kb and 1.0 kb sequence deletions at different positions of the 3'-intergenic region, in tight linkage disequilibrium with the most highly differentiated single nucleotide polymorphism, one of the alleles being very frequent in Senegal and The Gambia but rare in the other locations. Long non-coding RNA transcripts were previously shown to include the entire antisense of the gdv1 coding sequence and the portion of the intergenic region with allelic deletions, suggesting adaptive regulation of parasite sexual development and transmission in response to local conditions.

Whole genome sequencing and microsatellite analysis of the Plasmodium falciparum E5 NF54 strain show that the var, rifin and stevor gene families follow Mendelian inheritance

Background

Plasmodium falciparum exhibits a high degree of inter-isolate genetic diversity in its variant surface antigen (VSA) families: P. falciparum erythrocyte membrane protein 1, repetitive interspersed family (RIFIN) and subtelomeric variable open reading frame (STEVOR). The role of recombination for the generation of this diversity is a subject of ongoing research. Here the genome of E5, a sibling of the 3D7 genome strain is presented. Short and long read whole genome sequencing (WGS) techniques (Ilumina, Pacific Bioscience) and a set of 84 microsatellites (MS) were employed to characterize the 3D7 and non-3D7 parts of the E5 genome. This is the first time that VSA genes in sibling parasites were analysed with long read sequencing technology.

Results

Of the 5733 E5 genes only 278 genes, mostly var and rifin/stevor genes, had no orthologues in the 3D7 genome. WGS and MS analysis revealed that chromosomal crossovers occurred at a rate of 0–3 per chromosome. var, stevor and rifin genes were inherited within the respective non-3D7 or 3D7 chromosomal context. 54 of the 84 MS PCR fragments correctly identified the respective MS as 3D7- or non-3D7 and this correlated with var and rifin/stevor gene inheritance in the adjacent chromosomal regions. E5 had 61 var and 189 rifin/stevor genes. One large non-chromosomal recombination event resulted in a new var gene on chromosome 14. The remainder of the E5 3D7-type subtelomeric and central regions were identical to 3D7.

Conclusions

The data show that the rifin/stevor and var gene families represent the most diverse compartments of the P. falciparum genome but that the majority of var genes are inherited without alterations within their respective parental chromosomal context. Furthermore, MS genotyping with 54 MS can successfully distinguish between two sibling progeny of a natural P. falciparum cross and thus can be used to investigate identity by descent in field isolates.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2503-2) contains supplementary material, which is available to authorized users.

Global genetic diversity of var2csa in Plasmodium falciparum with implications for malaria in pregnancy and vaccine development

Malaria infection during pregnancy, caused by the sequestering of Plasmodium falciparum parasites in the placenta, leads to high infant mortality and maternal morbidity. The parasite-placenta adherence mechanism is mediated by the VAR2CSA protein, a target for natural occurring immunity. Currently, vaccine development is based on its ID1-DBL2Xb domain however little is known about the global genetic diversity of the encoding var2csa gene, which could influence vaccine efficacy. In a comprehensive analysis of the var2csa gene in >2,000 P. falciparum field isolates across 23 countries, we found that var2csa is duplicated in high prevalence (>25%), African and Oceanian populations harbour a much higher diversity than other regions, and that insertions/deletions are abundant leading to an underestimation of the diversity of the locus. Further, ID1-DBL2Xb haplotypes associated with adverse birth outcomes are present globally, and African-specific haplotypes exist, which should be incorporated into vaccine design.

High-Resolution Single-Cell Sequencing of Malaria Parasites

Single-cell genomics is a powerful tool for determining the genetic architecture of complex communities of unicellular organisms. In areas of high transmission, malaria patients are often challenged by the activities of multiple Plasmodium falciparum lineages, which can potentiate pathology, spread drug resistance loci, and also complicate most genetic analysis. Single-cell sequencing of P. falciparum would be key to understanding infection complexity, though efforts are hampered by the extreme nucleotide composition of its genome (∼80% AT-rich). To counter the low coverage achieved in previous studies, we targeted DNA-rich late-stage parasites by Fluorescence-Activated Cell Sorting and whole genome sequencing. Our method routinely generates accurate, near-complete capture of the 23 Mb P. falciparum genome (mean breadth of coverage 90.7%) at high efficiency. Data from 48 single-cell genomes derived from a polyclonal infection sampled in Chikhwawa, Malawi allowed for unambiguous determination of haplotype diversity and recent meiotic events, information that will aid public health efforts.

Fitness Costs and the Rapid Spread of kelch13-C580Y Substitutions Conferring Artemisinin Resistance

Fitness costs are key determinants of whether drug resistance alleles establish and how fast they spread within populations. More than 125 different kelch13 alleles, each containing a different amino acid substitution, have arisen in Southeast Asian malaria parasite (Plasmodium falciparum) populations under artemisinin selection over the past 15 years in a dramatic example of a soft selective event.

Fitness costs are key determinants of whether drug resistance alleles establish and how fast they spread within populations. More than 125 different kelch13 alleles, each containing a different amino acid substitution, have arisen in Southeast Asian malaria parasite (Plasmodium falciparum) populations under artemisinin selection over the past 15 years in a dramatic example of a soft selective event. However, just one of these alleles (C580Y) is now outcompeting other alleles in multiple different countries and is spreading toward fixation. Here we examine the fitness consequences of C580Y, relative to another less successful kelch13 mutation (R561H), to try to explain the distinctive dynamics of C580Y. We hypothesized that C580Y will show lower fitness costs than other kelch13 substitutions in the absence of artemisinin treatment. We used CRISPR/Cas9 methods to introduce single mutations (C580Y or R561H) or synonymous control edits into a wild-type parasite isolated on the Thailand-Myanmar border, conducted replicated head-to-head competition assays, and determined the outcome of competition using deep sequencing of kelch13 amplicons. Contrary to our predictions, these experiments reveal that C580Y carries higher fitness costs (s [selection coefficient] = 0.15 ± 0.008 [1 standard error {SE}]) than R561H (s = 0.084 ± 0.005). Furthermore, R561H outcompetes C580Y in direct competition (s = 0.065 ± 0.004). We conclude that fitness costs of C580Y in isolation are unlikely to explain the rapid spread of this substitution.

Panoptes: web-based exploration of large scale genome variation data

Motivation

The size and complexity of modern large-scale genome variation studies demand novel approaches for exploring and sharing the data. In order to unlock the potential of these data for a broad audience of scientists with various areas of expertise, a unified exploration framework is required that is accessible, coherent and user-friendly.

Results

Panoptes is an open-source software framework for collaborative visual exploration of large-scale genome variation data and associated metadata in a web browser. It relies on technology choices that allow it to operate in near real-time on very large datasets. It can be used to browse rich, hybrid content in a coherent way, and offers interactive visual analytics approaches to assist the exploration. We illustrate its application using genome variation data of Anopheles gambiae, Plasmodium falciparum and Plasmodium vivax.

Availability and implementation

Freely available at https://github.com/cggh/panoptes, under the GNU Affero General Public License.

Contact

paul.vauterin@gmail.com.

Commit, hide and escape: the story of Plasmodium gametocytes

Malaria is the major cause of mortality and morbidity in tropical countries. The causative agent, Plasmodium sp., has a complex life cycle and is armed with various mechanisms which ensure its continuous transmission. Gametocytes represent the sexual stage of the parasite and are indispensable for the transmission of the parasite from the human host to the mosquito. Despite its vital role in the parasite's success, it is the least understood stage in the parasite's life cycle. The presence of gametocytes in asymptomatic populations and induction of gametocytogenesis by most antimalarial drugs warrants further investigation into its biology. With a renewed focus on malaria elimination and advent of modern technology available to biologists today, the field of gametocyte biology has developed swiftly, providing crucial insights into the molecular mechanisms driving sexual commitment. This review will summarise key current findings in the field of gametocyte biology and address the associated challenges faced in malaria detection, control and elimination.

Cross-Contamination Explains "Inter and Intraspecific Horizontal Genetic Transfers" between Asexual Bdelloid Rotifers

A few metazoan lineages are thought to have persisted for millions of years without sexual reproduction. If so, they would offer important clues to the evolutionary paradox of sex itself [1, 2]. Most "ancient asexuals" are subject to ongoing doubt because extant populations continue to invest in males [3-9]. However, males are famously unknown in bdelloid rotifers, a class of microscopic invertebrates comprising hundreds of species [10-12]. Bdelloid genomes have acquired an unusually high proportion of genes from non-metazoans via horizontal transfer [13-17]. This well-substantiated finding has invited speculation [13] that homologous horizontal transfer between bdelloid individuals also may occur, perhaps even "replacing" sex [14]. In 2016, Current Biology published an article claiming to supply evidence for this idea. Debortoli et al. [18] sampled rotifers from natural populations and sequenced one mitochondrial and four nuclear loci. Species assignments were incongruent among loci for several samples, which was interpreted as evidence of "interspecific horizontal genetic transfers." Here, we use sequencing chromatograms supplied by the authors to demonstrate that samples treated as individuals actually contained two or more highly divergent mitochondrial and ribosomal sequences, revealing cross-contamination with DNA from multiple animals of different species. Other chromatograms indicate contamination with DNA from conspecific animals, explaining genetic and genomic evidence for "intraspecific horizontal exchanges" reported in the same study. Given the clear evidence of contamination, the data and findings of Debortoli et al. [18] provide no reliable support for their conclusions that DNA is transferred horizontally between or within bdelloid species.

Identity-by-descent analyses for measuring population dynamics and selection in recombining pathogens

Identification of genomic regions that are identical by descent (IBD) has proven useful for human genetic studies where analyses have led to the discovery of familial relatedness and fine-mapping of disease critical regions. Unfortunately however, IBD analyses have been underutilized in analysis of other organisms, including human pathogens. This is in part due to the lack of statistical methodologies for non-diploid genomes in addition to the added complexity of multiclonal infections. As such, we have developed an IBD methodology, called isoRelate, for analysis of haploid recombining microorganisms in the presence of multiclonal infections. Using the inferred IBD status at genomic locations, we have also developed a novel statistic for identifying loci under positive selection and propose relatedness networks as a means of exploring shared haplotypes within populations. We evaluate the performance of our methodologies for detecting IBD and selection, including comparisons with existing tools, then perform an exploratory analysis of whole genome sequencing data from a global Plasmodium falciparum dataset of more than 2500 genomes. This analysis identifies Southeast Asia as having many highly related isolates, possibly as a result of both reduced transmission from intensified control efforts and population bottlenecks following the emergence of antimalarial drug resistance. Many signals of selection are also identified, most of which overlap genes that are known to be associated with drug resistance, in addition to two novel signals observed in multiple countries that have yet to be explored in detail. Additionally, we investigate relatedness networks over the selected loci and determine that one of these sweeps has spread between continents while the other has arisen independently in different countries. IBD analysis of microorganisms using isoRelate can be used for exploring population structure, positive selection and haplotype distributions, and will be a valuable tool for monitoring disease control and elimination efforts of many diseases.

There are growing concerns over the emergence of antimicrobial drug resistance, which threatens the efficacy of treatments for infectious diseases such as malaria. As such, it is important to understand the dynamics of resistance by investigating population structure, natural selection and disease transmission in microorganisms. The study of disease dynamics has been hampered by the lack of suitable statistical models for analysis of isolates containing multiple infections. We introduce a statistical model that uses population genomic data to identify genomic regions (loci) that are inherited from a common ancestor, in the presence of multiple infections. We demonstrate its potential for biological discovery using a global Plasmodium falciparum dataset. We identify low genetic diversity in isolates from Southeast Asia, possibly from clonal expansion following intensified control efforts after the emergence of artemisinin resistance. We also identify loci under positive selection, most of which contain genes that have been associated with antimalarial drug resistance. We discover two loci under strong selection in multiple countries throughout Southeast Asia and Africa where the selection pressure is currently unknown. We find that the selection pressure at one of these loci has originated from gene flow, while the other loci has originated from multiple independent events.

hmmIBD: software to infer pairwise identity by descent between haploid genotypes

Background

A number of recent malaria studies have used identity by descent (IBD) to study epidemiological processes relevant to malaria control. In this paper, a software package, hmmIBD, is introduced for estimating pairwise IBD between haploid genomes, such as those of the malaria parasite, sampled from one or two populations. Source code is freely available.

Methods

The performance of hmmIBD was verified using simulated data and benchmarked against an existing method for detecting IBD within populations. Code for all tests is freely available. The utility of hmmIBD for detecting IBD across populations was demonstrated using Plasmodium falciparum data from Cambodia and Ghana.

Results

Alongside an existing method, hmmIBD was highly accurate, sensitive and specific. It is fast, requiring only 70 s on average to analyse 50 whole genome sequences on a laptop computer, and scales linearly in the number of pairwise comparisons. Treatment of different populations under hmmIBD improves detection of IBD across populations.

Conclusion

Fast and accurate software for detecting IBD in malaria parasite genetic data sampled from one or two populations is presented. The latter will likely be a useful feature for malaria elimination efforts, since it could facilitate identification of imported malaria cases. Software is robust to possible misspecification of the genotyping error and the recombination rate. However, exclusion of data in regions whose rates vary greatly from their genome-wide average is recommended.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2349-7) contains supplementary material, which is available to authorized users.

Long read assemblies of geographically dispersed Plasmodium falciparum isolates reveal highly structured subtelomeres

Background: Although thousands of clinical isolates of Plasmodium falciparum are being sequenced and analysed by short read technology, the data do not resolve the highly variable subtelomeric regions of the genomes that contain polymorphic gene families involved in immune evasion and pathogenesis. There is also no current standard definition of the boundaries of these variable subtelomeric regions.

Methods: Using long-read sequence data (Pacific Biosciences SMRT technology), we assembled and annotated the genomes of 15 P. falciparum isolates, ten of which are newly cultured clinical isolates. We performed comparative analysis of the entire genome with particular emphasis on the subtelomeric regions and the internal var genes clusters.

  Results: The nearly complete sequence of these 15 isolates has enabled us to define a highly conserved core genome, to delineate the boundaries of the subtelomeric regions, and to compare these across isolates. We found highly structured variable regions in the genome. Some exported gene families purportedly involved in release of merozoites show copy number variation. As an example of ongoing genome evolution, we found a novel CLAG gene in six isolates.  We also found a novel gene that was relatively enriched in the South East Asian isolates compared to those from Africa.

Conclusions: These 15 manually curated new reference genome sequences with their nearly complete subtelomeric regions and fully assembled genes are an important new resource for the malaria research community. We report the overall conserved structure and pattern of important gene families and the more clearly defined subtelomeric regions.

Global analysis of Plasmodium falciparum histidine-rich protein-2 (pfhrp2) and pfhrp3 gene deletions using whole-genome sequencing data and meta-analysis

Many rapid diagnostic tests (RDT) used on suspected malaria cases are based on the detection of the protein encoded by the Plasmodium falciparum histidine-rich protein-2 (pfhrp2) gene, which shares a high sequence homology with pfhrp3 in the 3D7 reference genome. Parasite isolates showing pfhrp2 and pfhrp3 gene deletions have been emerging over the years, but a comprehensive genetic analysis of these variants is still lacking. With this purpose, genomic data from experimental P. falciparum genetic crosses between different laboratory lines (3D7, HB3, DD2, 7G8 and GB4) were first analysed (n = 98). The frequency of pfhrp2 deletions was consistent with a Mendelian prediction in HB3 × DD2 (56.7%; 95%CI = (39.5%-72.9%)). Moreover, the pfhrp2 and pfhrp3 deletions segregated independently of each other in the same genetic cross. Analysis of 3D7 × HB3 and 7G8 × GB4 estimated the probability of spontaneously generating a pfhrp2 deletion during sexual recombination to be up to 6.2%. Next, whole genome sequence data from 1970 P. falciparum isolates collected globally were analysed. Nine samples displayed depth of coverage consistent with pfhrp2 deletions (0.5%), but the corresponding split-read analysis could not confirm deletions in seven of these samples. Twenty-eight isolates had evidence of pfhrp3 deletions (1.4%), which are widespread in Southeast Asia. Finally, a meta-analysis of published data revealed a positive mean association between the frequencies of pfhrp2 and pfhrp3 deletions in Africa and South America. This result suggested a shared selective pressure acting on these genetic variants. In conclusion, evidence of genetic selection on both pfhrp2 and pfhrp3 deletions was presented, but experimental crosses do not provide evidence of a fitness cost of these variants. Further work is urgently needed to accurately determine the prevalence and the degree of association between these genetic variants, and the respective impact on diagnostic accuracy of many in-use RDT.

Putting evolution in elimination: Winning our ongoing battle with evolving malaria mosquitoes and parasites

Since 2000, the world has made significant progress in reducing malaria morbidity and mortality, and several countries in Africa, South America and South-East Asia are working hard to eliminate the disease. These elimination efforts continue to rely heavily on antimalarial drugs and insecticide-based interventions, which remain the cornerstones of malaria treatment and prevention. However, resistance has emerged against nearly every antimalarial drug and insecticide that is available. In this review we discuss the evolutionary consequences of the way we currently implement antimalarial interventions, which is leading to resistance and may ultimately lead to control failure, but also how evolutionary principles can be applied to extend the lifespan of current and novel interventions. A greater understanding of the general evolutionary principles that are at the core of emerging resistance is urgently needed if we are to develop improved resistance management strategies with the ultimate goal to achieve a malaria-free world.

Repetitive sequences in malaria parasite proteins

Five species of parasite cause malaria in humans with the most severe disease caused by Plasmodium falciparum. Many of the proteins encoded in the P. falciparum genome are unusually enriched in repetitive low-complexity sequences containing a limited repertoire of amino acids. These repetitive sequences expand and contract dynamically and are among the most rapidly changing sequences in the genome. The simplest repetitive sequences consist of single amino acid repeats such as poly-asparagine tracts that are found in approximately 25% of P. falciparum proteins. More complex repeats of two or more amino acids are also common in diverse parasite protein families. There is no universal explanation for the occurrence of repetitive sequences and it is possible that many confer no function to the encoded protein and no selective advantage or disadvantage to the parasite. However, there are increasing numbers of examples where repetitive sequences are important for parasite protein function. We discuss the diverse roles of low-complexity repetitive sequences throughout the parasite life cycle, from mediating protein-protein interactions to enabling the parasite to evade the host immune system.

Origins of the current outbreak of multidrug-resistant malaria in southeast Asia: a retrospective genetic study

BACKGROUND

Antimalarial resistance is rapidly spreading across parts of southeast Asia where dihydroartemisinin-piperaquine is used as first-line treatment for Plasmodium falciparum malaria. The first published reports about resistance to antimalarial drugs came from western Cambodia in 2013. Here, we analyse genetic changes in the P falciparum population of western Cambodia in the 6 years before those reports.

METHODS

We analysed genome sequence data on 1492 P falciparum samples from 11 locations across southeast Asia, including 464 samples collected in western Cambodia between 2007 and 2013. Different epidemiological origins of resistance were identified by haplotypic analysis of the kelch13 artemisinin resistance locus and the plasmepsin 2-3 piperaquine resistance locus.

FINDINGS

We identified more than 30 independent origins of artemisinin resistance, of which the KEL1 lineage accounted for 140 (91%) of 154 parasites resistant to dihydroartemisinin-piperaquine. In 2008, KEL1 combined with PLA1, the major lineage associated with piperaquine resistance. By 2013, the KEL1/PLA1 co-lineage had reached a frequency of 63% (24/38) in western Cambodia and had spread to northern Cambodia.

INTERPRETATION

The KEL1/PLA1 co-lineage emerged in the same year that dihydroartemisinin-piperaquine became the first-line antimalarial drug in western Cambodia and spread rapidly thereafter, displacing other artemisinin-resistant parasite lineages. These findings have important implications for management of the global health risk associated with the current outbreak of multidrug-resistant malaria in southeast Asia.

FUNDING

Wellcome Trust, Bill & Melinda Gates Foundation, Medical Research Council, UK Department for International Development, and the Intramural Research Program of the National Institute of Allergy and Infectious Diseases.

A set of microsatellite markers to differentiate Plasmodium falciparum progeny of four genetic crosses

Background

Four Plasmodium falciparum genetic crosses (HB3×3D7, HB3×Dd2, 7G8×GB4, and 803×GB4) have produced sets of recombinant progeny that are widely used for malaria research, including investigations of anti-malarial drug resistance. It is critical to maintain the progeny free from cross-contamination. Microsatellite polymorphisms can be used to validate parasite identity.

Results

A set of 12 markers was developed that differentiates the parents of the four P. falciparum crosses. This typing set identified distinguishing patterns of inheritance (fingerprints) in segregant collections of 15 progeny clones from HB3×3D7, 32 from HB3×Dd2, 33 from 7G8×GB4, and 81 from 803×GB4. Stronger amplification was observed with shorter relative to longer alleles of individual microsatellites. In experiments with mixed parental DNAs, electropherograms showed that signals of cross-contamination can be missed when minor peaks less than 1/4 or 1/3 the height of the major peak are disregarded by threshold settings commonly used for population studies.

Conclusions

Microsatellite typing is an effective method to check the identity of P. falciparum lines and detect parasite cross-contamination in cultures; however, care must be taken not to ignore minor peaks that can be overlooked. The 12 microsatellite markers presented here provide a rapid and efficient means to distinguish the segregants of laboratory crosses. Fingerprint patterns from these markers are useful to maintain the integrity of diverse parasite lines in and between research laboratories.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2210-z) contains supplementary material, which is available to authorized users.

Disrupting CD147-RAP2 interaction abrogates erythrocyte invasion by Plasmodium falciparum

Effective vaccines against malaria caused by Plasmodium falciparum are still lacking, and the molecular mechanism of the host-parasite interaction is not fully understood. Here we demonstrate that the interaction of RAP2, a parasite-secreted rhoptry protein that functions in the parasitophorous vacuole formation stage of the invasion, and CD147 on the host erythrocyte is essential for erythrocyte invasion by P falciparum and is independent from all previously identified interactions involved. Importantly, the blockade of the CD147-RAP2 interaction by HP6H8, a humanized CD147 antibody, completely abolished the parasite invasion with both cure and preventative functions in a humanized mouse model. Together with its long half-life on human red blood cells and its safety profile in cynomolgus monkeys, HP6H8 is the first antibody that offers an advantageous approach by targeting a more conserved late-stage parasite ligand for preventing as well as treating severe malaria.

Comparative analysis of low complexity regions in Plasmodia

Low complexity regions (LCRs) are a common feature shared by many genomes, but their evolutionary and functional significance remains mostly unknown. At the core of the uncertainty is a poor understanding of the mechanisms that regulate their retention in genomes, whether driven by natural selection or neutral evolution. Applying a comparative approach of LCRs to multiple strains and species is a powerful approach to identify patterns of conservation in these regions. Using this method, we investigate the evolutionary history of LCRs in the genus Plasmodium based on orthologous protein coding genes shared by 11 species and strains from primate and rodent-infecting pathogens. We find multiple lines of evidence in support of natural selection as a major evolutionary force shaping the composition and conservation of LCRs through time and signatures that their evolutionary paths are species specific. Our findings add a comparative analysis perspective to the debate on the evolution of LCRs and harness the power of sequence comparisons to identify potential functionally important LCR candidates.

Deconvolution of multiple infections in Plasmodium falciparum from high throughput sequencing data

Motivation

The presence of multiple infecting strains of the malarial parasite Plasmodium falciparum affects key phenotypic traits, including drug resistance and risk of severe disease. Advances in protocols and sequencing technology have made it possible to obtain high-coverage genome-wide sequencing data from blood samples and blood spots taken in the field. However, analyzing and interpreting such data is challenging because of the high rate of multiple infections present.

Results

We have developed a statistical method and implementation for deconvolving multiple genome sequences present in an individual with mixed infections. The software package DEploid uses haplotype structure within a reference panel of clonal isolates as a prior for haplotypes present in a given sample. It estimates the number of strains, their relative proportions and the haplotypes presented in a sample, allowing researchers to study multiple infection in malaria with an unprecedented level of detail.

Availability and implementation

The open source implementation DEploid is freely available at https://github.com/mcveanlab/DEploid under the conditions of the GPLv3 license. An R version is available at https://github.com/mcveanlab/DEploid-r.

Contact

joe.zhu@bdi.ox.ac.uk or gil.mcvean@bdi.ox.ac.uk.

Supplementary information

Supplementary data are available at Bioinformatics online.

Scale up of a Plasmodium falciparum elimination program and surveillance system in Kayin State, Myanmar

Background: Myanmar has one of the largest malaria burdens in the Greater Mekong Subregion (GMS). Throughout the GMS, Plasmodium falciparum parasites are increasingly resistant to artemisinin combination therapies. Given that there are no current alternative treatment therapies, one proposed solution to the threat of untreatable P. falciparum malaria is to eliminate the parasite from the region. Several small-scale elimination projects have been piloted in the GMS, including along the Myanmar-Thailand border. Following the success of the pilot elimination project along the Myanmar-Thailand border, there was a scale up to a broad area of Eastern Kayin State, Myanmar. Here we describe the establishment of the scale up elimination project in Easter Kayin State.

Methods: The scale up relied on geographic reconnaissance and a geographic information system, community engagement, generalized access to community-based early diagnosis and treatment, near real-time epidemiological surveillance, cross sectional malaria prevalence surveys and targeted mass drug administration in villages with high prevalence of P. falciparum malaria. Molecular markers of drug resistance were also monitored in individuals with symptomatic and asymptomatic infections.

Discussion: This protocol illustrates the establishment of an elimination project and operational research in a remote, rural area encompassing several armed groups, multiple political organizations and a near-absent health care infrastructure. The establishment of the project relied on a strong rapport with the target community, on-the-ground knowledge (through geographic surveys and community engagement), rapid decision making and an approach that was flexible enough to quickly adapt to a complex landscape. The elimination project is ongoing, now over three years in operation, and assessment of the impact of this operational research will follow. This project has relevance not only for other malaria elimination projects but also for operational research aimed at eliminating other diseases.

Scale up of a Plasmodium falciparum elimination program and surveillance system in Kayin State, Myanmar

Background: Myanmar has one of the largest malaria burdens in the Greater Mekong Subregion (GMS). Throughout the GMS, Plasmodium falciparum parasites are increasingly resistant to artemisinin combination therapies. Given that there are no current alternative treatment therapies, one proposed solution to the threat of untreatable P. falciparum malaria is to eliminate the parasite from the region. Several small-scale elimination projects have been piloted in the GMS, including along the Myanmar-Thailand border. Following the success of the pilot elimination project along the Myanmar-Thailand border, there was a scale up to a broad area of Eastern Kayin State, Myanmar. Here we describe the establishment of the scale up elimination project in Easter Kayin State. Methods: The scale up relied on geographic reconnaissance and a geographic information system, community engagement, generalized access to community-based early diagnosis and treatment, near real-time epidemiological surveillance, cross sectional malaria prevalence surveys and targeted mass drug administration in villages with high prevalence of P. falciparum malaria. Molecular markers of drug resistance were also monitored in individuals with symptomatic and asymptomatic infections. Discussion: This protocol illustrates the establishment of an elimination project and operational research in a remote, rural area encompassing several armed groups, multiple political organizations and a near-absent health care infrastructure. The establishment of the project relied on a strong rapport with the target community, on-the-ground knowledge (through geographic surveys and community engagement), rapid decision making and an approach that was flexible enough to quickly adapt to a complex landscape. The elimination project is ongoing, now over three years in operation, and assessment of the impact of this operational research will follow. This project has relevance not only for other malaria elimination projects but also for operational research aimed at eliminating other diseases.

Extreme mutation bias and high AT content in Plasmodium falciparum

For reasons that remain unknown, the Plasmodium falciparum genome has an exceptionally high AT content compared to other Plasmodium species and eukaryotes in general - nearly 80% in coding regions and approaching 90% in non-coding regions. Here, we examine how this phenomenon relates to genome-wide patterns of de novo mutation. Mutation accumulation experiments were performed by sequential cloning of six P. falciparum isolates growing in human erythrocytes in vitro for 4 years, with 279 clones sampled for whole genome sequencing at different time points. Genome sequence analysis of these samples revealed a significant excess of G:C to A:T transitions compared to other types of nucleotide substitution, which would naturally cause AT content to equilibrate close to the level seen across the P. falciparum reference genome (80.6% AT). These data also uncover an extremely high rate of small indel mutation relative to other species, primarily associated with repetitive AT-rich sequences, in addition to larger-scale structural rearrangements focused in antigen-coding var genes. In conclusion, high AT content in P. falciparum is driven by a systematic mutational bias and ultimately leads to an unusual level of microstructural plasticity, raising the question of whether this contributes to adaptive evolution.

Quantifying connectivity between local Plasmodium falciparum malaria parasite populations using identity by descent

With the rapidly increasing abundance and accessibility of genomic data, there is a growing interest in using population genetic approaches to characterize fine-scale dispersal of organisms, providing insight into biological processes across a broad range of fields including ecology, evolution and epidemiology. For sexually recombining haploid organisms such as the human malaria parasite P. falciparum, however, there have been no systematic assessments of the type of data and methods required to resolve fine scale connectivity. This analytical gap hinders the use of genomics for understanding local transmission patterns, a crucial goal for policy makers charged with eliminating this important human pathogen. Here we use data collected from four clinics with a catchment area spanning approximately 120 km of the Thai-Myanmar border to compare the ability of divergence (FST) and relatedness based on identity by descent (IBD) to resolve spatial connectivity between malaria parasites collected from proximal clinics. We found no relationship between inter-clinic distance and FST, likely due to sampling of highly related parasites within clinics, but a significant decline in IBD-based relatedness with increasing inter-clinic distance. This association was contingent upon the data set type and size. We estimated that approximately 147 single-infection whole genome sequenced parasite samples or 222 single-infection parasite samples genotyped at 93 single nucleotide polymorphisms (SNPs) were sufficient to recover a robust spatial trend estimate at this scale. In summary, surveillance efforts cannot rely on classical measures of genetic divergence to measure P. falciparum transmission on a local scale. Given adequate sampling, IBD-based relatedness provides a useful alternative, and robust trends can be obtained from parasite samples genotyped at approximately 100 SNPs.

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.

MinorityReport, software for generalized analysis of causal genetic variants

Background

The widespread availability of next generation genome sequencing technologies has enabled a wide range of variant detection applications, especially in cancer and inborn genetic disorders. For model systems and microorganisms, the same technology may be used to discover the causative mutations for any phenotype, including those generated in response to chemical perturbation. In the case of pathogenic organisms, these approaches have allowed the determination of drug targets by means of resistance selection followed by genome sequencing.

Results

MinorityReport is open source software written in python that facilitates the comparison of any two sets of genome alignments for the purpose of rapidly identifying the spectrum of nonsynonymous changes, insertions or deletions, and copy number variations in a presumed mutant relative to its parent. Specifically, MinorityReport relates mapped sequence reads in SAM format output from any alignment tool for both the mutant and parent genome, relative to a reference genome, and produces the set of variants that distinguishes the mutant from the parent, all presented in an intuitive, straightforward report format. MinorityReport features tunable parameters for evaluating evidence and a scoring system that prioritizes reported variants based on relative proportions of read counts supporting the variant in the mutant versus parent data sets. The utility of MinorityReport is demonstrated using previously published publicly available data sets to find the determinants of resistance for novel anti-malarial drugs.

Conclusions

MinorityReport is readily available (github: JeremyHorst/MinorityReport) to identify the genetic mechanisms of drug resistance in Plasmodium, genotype-phenotype relationships in human diads, or genomic variations between any two related organisms.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-017-1730-2) contains supplementary material, which is available to authorized users.

Genetic markers associated with dihydroartemisinin-piperaquine failure in Plasmodium falciparum malaria in Cambodia: a genotype-phenotype association study

BACKGROUND

As the prevalence of artemisinin-resistant Plasmodium falciparum malaria increases in the Greater Mekong subregion, emerging resistance to partner drugs in artemisinin combination therapies seriously threatens global efforts to treat and eliminate this disease. Molecular markers that predict failure of artemisinin combination therapy are urgently needed to monitor the spread of partner drug resistance, and to recommend alternative treatments in southeast Asia and beyond.

METHODS

We did a genome-wide association study of 297 P falciparum isolates from Cambodia to investigate the relationship of 11 630 exonic single-nucleotide polymorphisms (SNPs) and 43 copy number variations (CNVs) with in-vitro piperaquine 50% inhibitory concentrations (IC50s), and tested whether these genetic variants are markers of treatment failure with dihydroartemisinin-piperaquine. We then did a survival analysis of 133 patients to determine whether candidate molecular markers predicted parasite recrudescence following dihydroartemisinin-piperaquine treatment.

FINDINGS

Piperaquine IC50s increased significantly from 2011 to 2013 in three Cambodian provinces (2011 vs 2013 median IC50s: 20·0 nmol/L [IQR 13·7-29·0] vs 39·2 nmol/L [32·8-48·1] for Ratanakiri, 19·3 nmol/L [15·1-26·2] vs 66·2 nmol/L [49·9-83·0] for Preah Vihear, and 19·6 nmol/L [11·9-33·9] vs 81·1 nmol/L [61·3-113·1] for Pursat; all p≤10-3; Kruskal-Wallis test). Genome-wide analysis of SNPs identified a chromosome 13 region that associates with raised piperaquine IC50s. A non-synonymous SNP (encoding a Glu415Gly substitution) in this region, within a gene encoding an exonuclease, associates with parasite recrudescence following dihydroartemisinin-piperaquine treatment. Genome-wide analysis of CNVs revealed that a single copy of the mdr1 gene on chromosome 5 and a novel amplification of the plasmepsin 2 and plasmepsin 3 genes on chromosome 14 also associate with raised piperaquine IC50s. After adjusting for covariates, both exo-E415G and plasmepsin 2-3 markers significantly associate (p=3·0 × 10-8 and p=1·7 × 10-7, respectively) with decreased treatment efficacy (survival rates 0·38 [95% CI 0·25-0·51] and 0·41 [0·28-0·53], respectively).

INTERPRETATION

The exo-E415G SNP and plasmepsin 2-3 amplification are markers of piperaquine resistance and dihydroartemisinin-piperaquine failures in Cambodia, and can help monitor the spread of these phenotypes into other countries of the Greater Mekong subregion, and elucidate the mechanism of piperaquine resistance. Since plasmepsins are involved in the parasite's haemoglobin-to-haemozoin conversion pathway, targeted by related antimalarials, plasmepsin 2-3 amplification probably mediates piperaquine resistance.

FUNDING

Intramural Research Program of the US National Institute of Allergy and Infectious Diseases, National Institutes of Health, Wellcome Trust, Bill & Melinda Gates Foundation, Medical Research Council, and UK Department for International Development.

Culture adaptation of malaria parasites selects for convergent loss-of-function mutants

Cultured human pathogens may differ significantly from source populations. To investigate the genetic basis of laboratory adaptation in malaria parasites, clinical Plasmodium falciparum isolates were sampled from patients and cultured in vitro for up to three months. Genome sequence analysis was performed on multiple culture time point samples from six monoclonal isolates, and single nucleotide polymorphism (SNP) variants emerging over time were detected. Out of a total of five positively selected SNPs, four represented nonsense mutations resulting in stop codons, three of these in a single ApiAP2 transcription factor gene, and one in SRPK1. To survey further for nonsense mutants associated with culture, genome sequences of eleven long-term laboratory-adapted parasite strains were examined, revealing four independently acquired nonsense mutations in two other ApiAP2 genes, and five in Epac. No mutants of these genes exist in a large database of parasite sequences from uncultured clinical samples. This implicates putative master regulator genes in which multiple independent stop codon mutations have convergently led to culture adaptation, affecting most laboratory lines of P. falciparum. Understanding the adaptive processes should guide development of experimental models, which could include targeted gene disruption to adapt fastidious malaria parasite species to culture.