#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

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


Authors: Xue Li aff001;  Sudhir Kumar aff002;  Marina McDew-White aff001;  Meseret Haile aff002;  Ian H. Cheeseman aff001;  Scott Emrich aff003;  Katie Button-Simons aff003;  Francois Nosten aff005;  Stefan H. I. Kappe aff002;  Michael T. Ferdig aff003;  Tim J. C. Anderson aff001;  Ashley M. Vaughan aff002
Authors place of work: Texas Biomedical Research Institute, San Antonio, Texas, United States of America aff001;  Center for Global Infectious Disease Research, Seattle Children’s Research Institute, Seattle, Washington, United States of America aff002;  Eck Institute for Global Health, Department of Biological Sciences, University of Notre Dame, Notre Dame, Indiana, United States of America aff003;  Electrical Engineering and Computer Science, University of Tennessee, Knoxville, Tennessee, United States of America aff004;  Shoklo Malaria Research Unit, Mahidol-Oxford Tropical Medicine Research Unit, Faculty of Tropical Medicine, Mahidol University, Mae Sot, Thailand aff005;  Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine Research building, University of Oxford Old Road campus, Oxford, United Kingdom aff006;  Department of Global Health, University of Washington, Seattle, Washington, United states of America aff007
Published in the journal: Genetic mapping of fitness determinants across the malaria parasite Plasmodium falciparum life cycle. PLoS Genet 15(10): e32767. doi:10.1371/journal.pgen.1008453
Category: Research Article
doi: https://doi.org/10.1371/journal.pgen.1008453

Summary

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.

Keywords:

Quantitative trait loci – Natural selection – Mammalian genomics – Blood – Parasitic diseases – Plasmodium – Malarial parasites – Parasitic life cycles

Introduction

Parasitic organisms frequently use multiple hosts and have several morphologically and transcriptionally distinctive life cycle stages. Within each host, parasites must circumvent immune defenses and navigate to new tissues. There are frequently extreme bottlenecks in parasite numbers during transmission [1], with rapid proliferative growth within hosts, and intense competition between co-infecting parasite genotypes. For example, the life cycle of malaria parasites involves successive infection of two hosts: female Anopheles mosquitoes, where gamete fusion, meiosis and recombination occurs, and humans in which parasites travel from the skin, develop in the liver and then proliferate asexually in the blood stream. Ideally, we would like to understand how natural selection operates across the complete life cycle and document the genes subject to selection pressures at each life cycle stage: during erythrocytic growth, gametocyte production, oocyst development in the mosquito midgut, migration of sporozoites to the salivary glands, transmission from the salivary glands, sporozoite survival in the skin, and establishment and parasite growth during liver stage development and exoerythrocytic merozoite release.

Selection can be directly measured by examining changes in allele frequency across these developmental stages. Shifts in allele frequencies in populations of thousands of progeny generated by experimental genetic crosses provide locus-specific readouts of competitive fitness. For example, deep sequencing of bulk populations containing thousands of recombinants identified genes selected under different regimens in yeast and Caenorhabditis elegans [25]. Bulk segregant analysis (BSA) has also been successfully applied to studies of several different parasitic organisms including coccidia (Eimeria tenella) and the human blood fluke Schistosoma mansoni [6, 7]. Our work was inspired by an exciting series of papers applying pooled sequencing approaches (termed linkage group selection in the malaria literature) for mapping genes of interest in rodent malaria parasites [813].

Most studies of Plasmodium falciparum to date focus only on the asexual erythrocytic stages [8, 1417], because they can be easily cultured in vitro in red blood cells, circumventing the need for humans or great apes, the natural hosts for this parasite. Two new research tools now allow us to examine selection across the complete life cycle of P. falciparum. First, we can maintain the complete life cycle of P. falciparum in a laboratory setting by using humanized mice [18] in place of splenectomized chimpanzees or human volunteers. These mice contain human hepatocytes and are therefore able to support liver stage development of P. falciparum. Hence, we can stage genetic crosses between different P. falciparum parasites, including parasites recently isolated from infected patients, and sample multiple parasite life stages for measurement of allele frequency changes throughout the life cycle. Second, selective whole genome amplification (sWGA) provides a simple and effective way to enrich Plasmodium DNA from contaminating host tissues. This is critical because Plasmodium DNA constitutes a very small fraction of DNA present in malaria-infected mosquitoes; likewise, Plasmodium DNA makes up a very small fraction of DNA extracted from malaria-infected livers (Table 1). sWGA uses short 8–12 mer oligonucleotide probes that preferentially bind to the target genome, rather than random hexamers used in normal whole genome amplification. This approach was pioneered by Leichty and Brisson [19], and protocols for sWGA have been successfully developed to amplify and sequence malaria parasite genomes from contaminating host tissues [2023].

Tab. 1. Sample collection and sequence statistics.
Sample collection and sequence statistics.

Artemisinin resistance is currently spreading across Southeast Asia [24]. SNPs in Kelch13 (PF3D7_1343700) locus on chromosome (chr) 13 underlie resistance and greater than 124 independent alleles have been recorded in dramatic example of a soft select sweep [25, 26]. One particular allele (Kelch13-C580Y) is currently replacing other resistant alleles and spreading toward fixation in independent transmission foci in western Cambodia/Laos/Vietnam and the Thailand-Myanmar border [2729]. Several studies have suggested that mutations within loci other than Kelch13 may provide a permissive background for evolution of artemisinin resistance or play a compensatory role [30, 31], but the role of such accessory loci is poorly understood.

In this study, we measured skews in allele frequencies across the genome in the progeny of a genetic cross between artemisinin resistant (ART-R, kelch13-C580Y) and ART sensitive (ART-S, kelch13-WT) parasites throughout the life cycle to identify genes that influence parasite fitness in parasite stages infecting both the mosquito and vertebrate host. We used ART-R and ART-S parental parasites in order to examine loci contributing to fitness and compensation for deleterious effects of ART-R alleles [16]. We used the humanized mouse model to allow parasite liver stage development of the genetic cross progeny, sWGA to enrich parasite DNA from host contaminations and pooled sequencing to determine temporal change in allele frequency and characterize genomic regions under selection. Our results demonstrate pervasive selection across the parasite genome over the course of a single parasite generation, selection against progeny produced from selfed matings, and strong locus-specific selection against parasite loci on chr 12 and 14.

Results

Identification of high-confidence SNPs between parents

P. falciparum NHP1337 and MKK2835 were cloned by limiting dilution and used as parents for genetic crosses. Both parasites are from the Thailand-Myanmar border. MKK2835 (ART-S) is a kelch13 wild-type ART-susceptible parasite collected from a patient visited the clinic in 2003 prior to the spread of ART resistance [32]. NHP1337 is a recent cloned ART-R parasite, that cleared slowly (Clearance half-life (T½P) = 7.84 h) from the blood of a patient treated with artemisinin combination therapy in 2011 and carries the C580Y kelch13 mutant. Parasites with C580Y mutation have been rapidly spreading in Southeast Asia and are replacing other ART-resistant kelch13 alleles [25, 33]. We detected 9,462 high confidence SNPs– 1 SNP per 2.43kb–between the two parental strains from 21 Mb core genome (defined in [34]) (S1 Fig). Comparison with single clone Southeast Asian parasites from Sanger pf3k project (ftp://ngs.sanger.ac.uk/production/pf3k/release_5/), shows that the parental parasites both fall into group designated as KH1 [35] (S7 Fig).

Genetic cross and generation of segregant pools

To generate segregant pools of progeny, we crossed NHP1337 and MKK2835 (Fig 1). We fed 500 mosquitoes with a ~50:50 gametocyte mixture of the two parental parasites. Recombinant progeny are generated after gametes fuse to form a diploid zygote that then rapidly transforms into a short-lived tetraploid ookinete which migrates to the basal lamina of the mosquito midgut and transforms into an oocyst. Mitotic division of the 4 meiotic products ultimately leads to the generation of on average over 3000 haploid sporozoites within each oocyst [36]. Oocyst prevalence was 80% with an average burden of three oocysts per mosquito midgut (range: 0–6), giving an estimate of 12 (3×4) recombinant genotypes per mosquito. We dissected a proportion of the infected mosquitoes to collect midguts (48 at each time point) for monitoring allele frequencies during oocyst development. Salivary gland sporozoites from 204 mosquitoes were pooled together and injected in to a single FRG huHep mouse.

Genetic mapping of parasite competition through <i>Plasmodium falciparum</i> life cycle.
Fig. 1. Genetic mapping of parasite competition through Plasmodium falciparum life cycle.
We generated genetic crosses using Anopheles stephensi mosquitoes and FRG huHep mice. We collected midgut and salivary glands from infected mosquitoes, infected mouse liver and emerging merozoites from in vivo blood, and recovered aliquots of in vitro cultured progeny parasites at intervals of 30 days (marked with arrows, parasite stage 1–6). Cross generation and sample collection were completed in two months (marked with green letters). For samples with host contamination or small amounts of DNA isolated (blue arrows, Table 1), selective whole genome amplification (sWGA) was performed before Illumina whole-genome sequencing (WGA). We used amplicon sequencing to trace biases in mtDNA transmission in those samples. For in vitro blood samples (pink arrow), we performed sequencing both before and after sWGA to evaluate the accuracy of allele frequency after sWGA.

We collected samples for allele frequency analysis from infected mosquito midguts, infected mosquito salivary glands, infected humanized mouse livers and infected blood (both mouse blood and injected human red blood cells) after the liver stage-to-blood stage transition. We then recovered aliquots of in vitro cultured progeny parasites at two-four day intervals over 30 days (Fig 1). We also set up cultures to enrich gametocytes from the in vitro cultures. These samples represent the important developmental stages across the parasite life cycle, including early oocyst, maturing oocyst, sporozoites, liver stage schizonts, transitioned blood stage parasites, fifteen asexual cycles in blood stage culture and reproductive gametocytes, required for transmission to the mosquito (Table 1). Our experiment examines the impact of selection across one complete P. falciparum life cycle.

We measured the total number of parasite genome copies and the amount of host DNA contamination for these segregant pools using qPCR. At the early midgut oocyst stage (4 days after mosquito infection), we isolated ~8,000 copies of the P. falciparum genome from 48 mosquito midguts. The parasite DNA represented approximately 0.01% of the total DNA within these isolated midguts. The percentage reached 1.80% after 10 days of mosquito infection, which indicating a 196-fold increase of parasite DNA in the six days following initial midgut isolation. The percentage of parasite DNA found in samples from mosquito salivary gland containing sporozoites, liver containing liver stage parasites and liver stage-to-blood stage transitioned in vivo blood samples were 3%, 3% and 30%, respectively (Table 1).

sWGA-WGS, WGS and amplicon sequencing

We used three approaches to sequence the segregant pools and quantify allele frequencies: (1) selective whole genome amplification combined with whole genome sequencing (sWGA-WGS), (2) direct whole genome sequencing (WGS) and (3) amplicon sequencing. The methods used were dependent on the level of host contamination and the total amount of DNA present in the samples (Table 1). We used multiple methods where possible to determine potential bias. We used the sWGA approach to enrich parasite DNA before WGS in samples with extensive host contamination, including the mosquito midgut and the FRG NOD human-chimeric mouse liver (Table 1, sample 2–5). With 0.2×106 copies of parasite genome as template, the sWGA-WGS approach yielded 0.6–1.4 μg of product after 3h of amplification, of which > 88% was from P. falciparum, for both mosquito and mouse samples. By sequencing pools to ~100× coverage, comparable results were obtained between samples prepared by the sWGA-WGS approach and the WGS approach (Fig 2A, Fig 3 and Fig 4). We used amplicon sequencing [16] to determine the frequencies of mtDNA from the two parents in those samples for which we used sWGA (Fig 2C and S2 Fig). This was necessary because our sWGA primers were specifically designed to minimize amplification of mtDNA, since we were concerned that sWGA with circular DNA would inundate autosomal sWGA products. For day 4 mosquito midgut samples, we only have amplicon sequencing data since there was insufficient parasite DNA for a successful sWGA.

Change in frequency of mitochondria and core genome at different infection stages.
Fig. 2. Change in frequency of mitochondria and core genome at different infection stages.
(A) Ridgeline plots showed genome-wide allele frequency distributions of NHP1337 through the Plasmodium life cycle. Each frequency distribution shows the frequency of genome-wide SNPs (9,462) found in progeny bulks at different time points of the parasite life cycle. * indicates Cohen’s d effect size > 0.5, and ** indicates effect size > 0.8. (B) We detected strong concordance between allele frequencies estimated from experimental replicates. (C) The allele frequency estimated from mitochondria and core genome showed the same pattern of skew with across the life cycle. (D) Natural log of the genotype ratio (NHP1337/MKK2835) plotted against asexual life cycles. The selection coefficient was estimated as the slope of the least-squares fit. Allele frequencies from day30 to day42 were used here. There was no significant difference between fitness costs estimated for the core genome and mitochondria (P = 0.363). Positive values of s indicate a selection disadvantage for NHP1337. MT, mitochondria; s, selection coefficients; R, correlation coefficient. X-axis in (A) and (C) indicated sample collecting days and corresponding parasite developmental stages.
Plot of allele frequencies across the genome through <i>Plasmodium falciparum</i> life cycle.
Fig. 3. Plot of allele frequencies across the genome through Plasmodium falciparum life cycle.
We divided the parasites into two replicates after two days’ in vitro culture (day 23). Orange and black indicated allele frequencies from these two parallel cultures. Red and blue lines showed tricube-smoothed allele frequencies. Black dashed lines indicated the average allele frequency across the genome. Sample collecting days were marked on the right. Day 10 shows allele frequencies of maturing oocysts, day 14 shows sporozoites, day 21.1 shows liver stage schizonts, day 21.2 shows transitioned blood stage parasites, and day 23–50 shows fifteen asexual cycles in blood stage culture.
Allele frequencies estimated before and after selective whole genome amplification (sWGA).
Fig. 4. Allele frequencies estimated before and after selective whole genome amplification (sWGA).
(A) Plot of allele frequencies across the genome. (B) Concordance between allele frequencies estimated before and after sWGA.

Evaluation of bias in allele frequency measurement

To evaluate the accuracy of allele frequencies estimated after sWGA, we sequenced blood samples using both the sWGA-WGS approach and the WGS approach. We plotted allele frequencies of the parent NHP1337 across the genome and tricube-smoothed the frequency with window size of 100kb to smooth out noise and estimate changes in adjacent regions. With 10 million 150 bp pair-end sequencing reads, there were fewer loci detected with coverage > 30× by the sWGA-WGS approach relative to WGS (5,024 loci by sWGA-WGS and 7,844 loci by direct WGS). The allele frequency trends, however, were highly consistent after smoothing (Fig 4A). The allele frequencies estimated before and after sWGA were strongly concordant (R2 = 0.985, Fig 4B, S6 Fig), which strongly supports the comparability of these two different methods.

Allele frequency changes in segregant pools

  • Mosquito stages: Plasmodium sexual blood stage infections differentiate into both male and female gametes and mate; consequently, selfed progeny, resulting from the fusion of gametes from the same parasite genotype, can occur (i.e., NHP1337 male gametes fertilizing NHP1337 female gametes and MKK2835 male gametes fertilizing MKK2835 female gametes). Selection towards selfed progeny is evident from skewing and shifting of whole genome allele frequencies. To investigate population composition at different infection stages, we plotted the allele frequency distribution of Plasmodium mitochondria and across the core genome (Fig 2). We observed a strong skew (>80%) towards alleles from the ART-R parent in the mosquito stages, which suggests that many selfed progeny from NHP1337 were present.

  • Liver stage: The allele frequency in the progeny parasite population shifted significantly towards the ART-R parent (NHP1337) at the liver stage. This is evident from comparisons of allele frequency distributions in the liver with those from sporozoites (Fig 2A, Cohen's d test, large effect size = 0.89). This skew observed in the liver stage is reduced in merozoites emerging from the liver (Cohen's d test, medium effect size = 0.61).

  • Blood stages: During in vitro culture, the allele frequency of NHP1337 (ART-R) dropped to 50%, between day 32 and day 40 (Cohen's d test, effect size = 1.65). We maintained replicate in vitro blood cultures from day 23 (corresponds to day 2 of in vitro blood stage culture). Highly repeatable skews were observed in allele frequencies across the genome in these two parallel cultures (Fig 2A and 2B, R2 = 0.985). Furthermore, we observed the same skews in both the mitochondria and across the core genome (Fig 2C), strongly suggesting that the selection was against NHP1337 selfed progeny. The NHP1337 selfed progeny were almost eliminated by day 42 and we thus estimated the selection coefficients against the NHP1337 selfed progeny. We observed strong selection against NHP1337 alleles, with s = 0.24±0.02 in the core autosomal genome and s = 0.22±0.01 in mitochondria (Fig 2D). There is no significant difference between these two estimates (p = 0.363, Least-Squares Means).

  • Gametocyte generation: The sexual commitment (gametocytogenesis) ratio of Plasmodium parasite is considered to be generally low (< 3%), but variable among different strains and under different conditions [3740]. Interestingly, we do not see specific allele frequency changes during the gametocyte enrichment experiment comparing to normal in vitro blood cultures (S3 Fig). These data suggest that progeny from this cross committed to gametocytes at similar rates.

Loci under selection

To pinpoint the loci that determine parasite fitness at each life cycle stage, we first plotted the whole genome allele frequencies throughout the life cycle (Fig 3). In addition to the whole genome skew described above, we also observed specific regions of the genome that showed distortion in allele frequency after day 32. The skews in allele frequencies were remarkably consistent between the two replicate blood stage cultures, suggesting pervasive selection at multiple loci across the genome. We calculated G’ values to measure the significance of allelic skews (Fig 5A, S4 Fig and S1 Table). Two strong quantitative trait loci (QTLs) were identified on chr 12 and 14, with genome-wide false discovery rate (FDR) < 0.01. We further used Δ (SNP-index) to determine the direction of the allele frequency changes (Fig 5B and S4 Fig). In both regions, alleles from NHP1337 (ART-R) were selected against. We then calculated selection coefficient (s) across the genome (Fig 5C). We observed particularly strong selection at these two QTL regions, with s = 0.12 on chr12 and s = 0.18 on chr14. In addition, there were a set of lower confidence QTLs with lower allele frequency changes and less impact on parasite fitness uncovered across the genome (Fig 5 and S1 Table).

Bulk segregant analysis.
Fig. 5. Bulk segregant analysis.
(A) QTLs were defined with G’ approach by comparing allele frequencies at each locus to the average allele frequency across the genome. Regions with FDR > 0.01 were taking as significant QTLs. (B) Δ(SNP-index) for day50 progeny pools. The Δ(SNP-index) was the difference between SNP-index of each locus and the genome-wide average SNP-index. A positive Δ(SNP-index) value indicates an increase in alleles from NHP1337. Red and blue lines showed the 95% and 99% confidential intervals that matched with the relevant window depth at each SNP. (C) Tricube-smoothed selection coefficients (s). Estimation of s was based on the changes of allele frequency from day25 to day50. The mean selection coefficient was adjusted to 0 to remove the influence of selfed progeny. Positive values of s indicate a disadvantage for alleles from NHP1337. Orange and black lines indicate experimental replicates.

Fine mapping of chr 12 and 14 QTLs

We calculated 95% confidence intervals to narrow down the genes driving selection within the two QTL regions. The QTL on chr 12 ranged from 1,102,148 to 1,327,968 (226 kb) and the QTL on chr 14 ranged from 2,378,002 to 2,541,869 (164 kb).

Chr 12: The QTL region contained 48 genes, with 27 genes bearing at least one non-synonymous mutation differentiating the two parents (Fig 6 and S2 Table). Among the candidate genes with functional annotation, the multidrug resistance-associated protein 2 gene (mrp2, PF3D7_1229100) was located at the peak of the chr 12 QTL (Fig 6A and S2 Table). The mrp2 allele from NHP1337 carries three indels (3–24 bp) within coding microsatellite sequences compared with that in MKK2835. These indels don’t interrupt the open reading frame.

Fig. 6.
Overview of the genes inside of QTL regions on chr 12 (A) and chr 14 (B). Black dashed vertical lines are boundaries of the 95% confidential intervals (CIs) of the QTL. The QTL on chr 12 spanned 226 kb and included 48 genes, and the QTL on ch14 spanned 164 kb and included 45 genes. 2D structure of MRP2 and ARPS10 are presented in boxes next to the G’ plot. The structure of MRP2 was adapted from Velga et al., 2014. There are 5 microindels in the coding region of Pfmrp2 gene (I-V, orange and green blocks). Four of the microindels (orange blocks, 1 SNP and 3 indels) are different between ART-S and ART-R parental strains. The changes in peptide length relative to P. falciparum 3D7 are indicated next to the microindels, as microindel: ART-S/ART-R. ART-S and ART-R parasites have the same amino acid insertion at microindel I, but the sequence includes a synonymous mutation. The structure of ARPS10 is predicted by I-TASSER. ART-R has two non-synonymous mutations in gene ARPS10, Val127Met and Asp128His (red star). TMD: transmembrane domain; NBD: nucleotide-binding domain; SP: signal peptides.

Chr 14: There are 45 genes located in this QTL and 13 contained non-synonymous mutations that distinguish the two parents (Fig 6 and S2 Table). The gene encoding apicoplast ribosomal protein S10 (arps10, PF3D7_1460900) was located at the peak of this QTL. There are two non-synonymous mutations (Val127Met and Asp128His) detected in arps10 from NHP1337 as compared to MKK2835. The Val127Met mutation was suggested to provide a permissive genetic background for artemisinin resistance-associated mutations in kelch13 in a genome-wide association analysis [31].

Discussion

Pervasive selection in a Plasmodium genetic cross

In this experiment, we observed both genome-wide selection against selfed progeny, and locus specific selection that resulted in skews in the frequency of particular parental alleles in progeny.

Genome-wide selection against selfed progeny

Initially, frequencies of alleles derived from the two parental parasites were strongly skewed (0.81 ± 0.08) towards the NHP1337 parent. This deviation from the expected 0.5 ratio for outcrossed progeny occurs because hermaphroditic malaria parasites produce both male and female gametocytes; fusion between male and female gametes of the same genotype (selfing) is possible. The simplest explanation for this observed skew is that an excess of selfed progeny were generated from the NHP1337 parent genotype compared to MKK2835 parent. We plotted the correlations between corrections between mitochondrial frequencies and frequencies of SNPs on different chromosomes across the experiment (S8 Fig). Strong correlations show that alleles are not segregated independently, supporting the selfing. We also cloned progeny collected on day 23 of the experiment, which confirmed our suspicion that selfing of NHP1337 leads to the skew in allele frequency (Button-Simmons et al. in preparation). Of 212 cloned genotyped progeny recovered, 144 (68%) are ART-R selfed, 5 (2%) are ART-S selfed, and 63 (30%) are recombinant progeny which represent 60 unique recombinants. Interestingly, our bulk sequencing data demonstrates 86% frequency of allele from ART-R parent. In the cloned progeny, we observed 68% + (30%/2) = 83% extremely close to our estimation from bulk data. The information from the cloned progeny strongly reinforces our conclusions from the bulk data. It is unclear whether the excess of selfed progeny from the NHP1337 parent relative to the MKK2835 parent results from an imbalance in gametocytes from these parental parasites when staging the cross, or from inherent differences in propensity to self in these two parasite clones.

Frequencies of NHP1337 remained high from day 10 (mature oocysts) until day 30 (after 10 days of in vitro blood culture). At this point, genome-wide frequencies of the NHP1337 parasite declined significantly from 0.85 to 0.54 on day 42. We observed a parallel decline of both mitochondrial and autosomal allele frequencies for the NHP1337 parasite. This is consistent with selection removing selfed NHP1337 genotypes from the progeny, otherwise we would expect selection on these two genomes to be decoupled. Selection was extremely strong (mitochondrial s = 0.22 ±0.01; autosomal s = 0.24 ±0.02) for both genomes. Furthermore, we observed the same patterns using whole genome sequencing and amplicon sequencing for measuring allele frequencies of mtDNA (S2 Fig), suggesting that our results are robust to any methodological biases. Analysis of further crosses will allow us to determine whether selection against selfed progeny is a general feature of crosses in malaria.

These data demonstrate systematic selection against genotypes generated by selfing of the ART-R parent. Reproduction by outcrossing is prevalent in nature, even in hermaphroditic species [41]. Inbreeding leads to reduced fitness of offspring (inbreeding depression), while outbreeding among genetically differentiated individuals improves the performance of the F1 generation (heterosis) [41, 42]. We observed strong selection against selfed NHP1337 genotypes which resulted in elimination of selfed progeny in six asexual cycles (day 30–42). Possible explanations for the lower fitness in selfed progeny are: (1) recombination allows removal of deleterious mutation in outcrossed progeny. Accumulation of deleterious mutations occurs during clonal expansion and in inbred parasite lineages. Both parental parasites used in this cross were isolated from Southeast Asia, an area of low parasite transmission intensity, where most infected patients harbor a single parasite genotype. As a consequence mosquito blood meals contain male and female sexual stages from the sample parasite clone, and therefore deleterious mutations can accumulate since inbreeding predominates [43, 44]. We speculate that recombinant genotypes generated by outcrossing between NHP1337 and MKK2835 parents have reduced numbers of deleterious alleles and therefore outcompete inbred parental genotypes. (2) In vitro culture, where the strongest selection was observed in this experiment, represents an ecological niche change for both parental genotypes. Recombinants generated by outbreeding may be more fit in these laboratory conditions. (3) The selfed NHP1337 parasites that predominate initially are ART-R which may carry a fitness cost relative to ART-S parasites [15, 16, 45]. We note that alleles at the two loci (chr 12 and 14) that are selected against (see discussion section “Locus Specific selection”) during days 30–50 are both derived from the ART-R parent.

There is an interesting shift in allele frequencies between sporozoites sampled from mosquito salivary glands and liver stage parasites recovered from infected mice on day 21 (Fig 2), with liver stage parasites carrying high frequencies of NHP1337 alleles (liver 0.89 vs sporozoites 0.79, with large Cohen’s d effect size [0.89]). The allele frequency of parasites from in vivo blood collected on the same day is 0.84, which is between those from sporozoites and liver stage parasites. During liver stage development, single sporozoites take up residency within hepatocytes and divide mitotically over the course of ~7 days (determined with laboratory strains of P. falciparum NF54 [46]) until liver schizonts burst releasing tens of thousands of merozoites into the blood. The simplest explanation of the observed allele frequency shift is a genotype-dependent variation in the duration of parasite liver stage development. We suggest that the selfed NHP1337 progeny remain in the liver longer and thus at the day 7 sampling, recombinant liver stage parasites have already transitioned to blood stage, generating the observed difference in allele frequencies. Further work is needed to directly determine the duration of liver stage development and if other liver stage parasite phenotypes (schizont size/merozoite numbers) differ among parasite genotypes.

Locus specific selection

We observe a progressive increase in the variance of allele frequencies of SNPs from day 30–50 (during blood stage culture) (Fig 2A). Several features of these data suggest that this is primarily driven by selection, rather than genetic drift. First, we note an extremely strong repeatability in allele frequency skews across the genome in the two replicate parasite cultures established from the humanized mouse infection. This is reflected in the high correlation between allele frequencies between these two replicates at the end of the experiment (Fig 3, day 50) when variance in allele frequency is at its maximum. The strong repeatability in patterns of skew observed suggest that there are multiple loci across the genome that influence parasite growth rate and competitive ability. Second, we see several regions of the genome that show extreme skew relative to the genome wide average. Two genome regions in particular (on chr 12 and on chr 14) show strong and significant skews that cannot be explained by drift. These allelic skews also increase progressively from 25–50 days, consistent with selection coefficients (s) of 0.18/48 hr asexual cycle for the chr 14 locus and (s) of 0.12/48 hr asexual generation for the chr 12 locus.

We observed strong selection against particular alleles segregating in this genetic cross (in the absence of drug pressure). How can such strongly disadvantageous alleles be maintained in natural parasite populations? We suggest three explanations. First, we think that the most likely explanation is that the fitness of these alleles may depend on genetic background [47] and reflect epistatic interactions. We note that of the two parental strains used in this study, MKK2835 (ART-S) was isolated in 2003, while NHP1337 (ART-R), was collected in 2013. In the 10 years between 2003–13, artemisinin-resistant parasites spread to high frequency on the Thailand-Myanmar border [25]. Intense drug selection in this 10-year interval has led to accumulation of additional genetic changes associated with ART-R, which may act epistatically together with other ART-R-associated genes [30]. It is certainly interesting that the chr 14 QTL contains arps10, which has been suggested to provide a permissive background for ART-R evolution [31]. Outcrossing between individuals with different adaptations can result in disruption of this selective advantage, resulting in a loss of fitness [48]. Further experimental work such as pairwise competition assays between recombinant progeny carrying different chr12 and chr14 haplotypes will be required to confirm the role of these loci, alone or in combination, in determining fitness. Second, there is a possibility that de novo deleterious mutations in these two QTL regions were fixed in the cloned NHP1337 parasites during the brief period of laboratory culture. We think this is unlikely because we also see pervasive selection at multiple genes outside these two major QTL regions, just with lower significance using G’ statistics. The overall recombination rate in this cross was 8.2 kb/cM (Button-Simons et al. in preparation). We further counted the recombination events between chromosome 12 and 14 loci (S5 Table). There were 35 recombination events observed in 60 unique recombinant progeny between the chromosome 12 and 14 segments. The recombination between these two loci was even, which indicated that the detecting of fitness traits of these two loci were independent. Third, we cannot discount the possibility that the strong selective disadvantages observed within these QTL regions reflects the artificial nature of this system with humanized mice and asexual culturing of parasites. During normal transmission in the field selection against these genes may not be present.

We note that similar bulk segregrant experiments examining fitness determinants in C. elegans [5] also detected QTLs with large effects sizes in several different genome regions. The QTLs regions identified corresponded to the location of known selfish elements, or co-localized with major eQTLs, consistent with the idea that epistatic interactions are important in fitness related traits. Similarly, analysis of C. elegans recombinant inbred lines generated in a 16-parent genetic cross revealed that ~40% of the variance of a key fitness trait (fertility) resulted from epistatic interactions between loci [49]. These C. elegans papers support the argument (above) that variation in fitness may be retained within natural populations due to epistasis among the genes involved.

We anticipate that intensity of competition among parasite clones within infected patients may closely parallel the patterns we observe within our genetic cross. The estimated occurrence rate of mixed infections ranges from 18% to 63% in African and Southeast Asia countries [43, 50]. Although there was likely more intense competition in this experimental cross, with millions of sporozoites infecting a single mouse, single cell sequencing has revealed seventeen unique clones in a single human infection [51], which suggests that similar competitive interactions will also occur in patients. We note that while the intensity of competition may be similar in humanized mice, in vitro parasite cultures or infected humans, the nature of selection may differ. In infected people, parasite genotypes that allow evasion of immunity or alter parasite cytoadhesion properties may be selected, while growth competition is likely to be the predominant selective force in immunosuppressed humanized mice or in vitro culture.

What drives QTL peaks on chr12 & chr14?

Inspection of the genes under the QTL peaks allows us to speculate about specific genes that may be driving the selection observed. Miotto et al. (2015) showed that four different non-kelch13 loci (ferredoxin, fd; apicoplast ribosomal protein S10, arps10; multidrug resistance protein 2, mdr2; chloroquine resistance transporter, crt) are associated with the resistance phenotype, but not directly responsible for resistance. They suggested that a suite of background mutations was a prerequisite for mutations in kelch13. In our experiment, arps10 falls near the peak of the strongly selected chr 14 loci (Fig 6), which could suggest a functional relationship. We examined the presence of the background mutations found in both parental strains. The ART-R parent, NHP1337, contains mutations in all of four of the genes described by Miotto et al. [31] (fd, mdr2, crt and arps10), while the ART-S parent, MKK2835, contains just three of these mutations (fd, mdr2 and crt), so only arps10 mutations are segregating in this cross. It will be interesting to test the role of the remaining three loci (fd, mdr2 and crt) by conducting additional experimental crosses.

The multidrug resistance-associated proteins (MRPs), belong to the C-family of ATP binding cassette (ABC) transport proteins that are well known for their role in multidrug resistance. Rodent malaria parasites encode one single MRP protein, whereas P. falciparum encodes two: MRP1 and MRP2 [52]. Several studies have shown that PfMRP1 is associated with P. falciparum’s response to multiple anti-malaria drugs and that disruption of PfMRP1 influences the fitness of parasites under normal culture conditions [5355]. The function of PfMRP2 is not as well understood. Transfection studies have shown that MRP2-deficient malaria parasites are not able to maintain a successful liver stage infection [52, 56]. In our study, mrp2 was found located at the peak QTL on chr12. We speculate that mrp2 may also play a role in parasite fitness during asexual parasite stages. However, we cannot exclude that other neighboring loci may drive the observed allele frequency changes. To confirm the roles of individual genes inside the QTL regions (on both chr12 and chr14), gene-editing studies will be required.

We further analyzed the ancestral/derived allelic state of genes inside QTLs regions of chr12 and chr14 (S1 File, S6 Table). For both pfmrp2 and pfarps10 genes, the selected alleles in the ART-S parent contained ancestral, rather than derived alleles. Deleterious derived alleles in the ART-R parasites may therefore explain the skews observed on chr12 and 14. We also analyzed the allele frequency distribution among Southeast Asia parasites. We used 678 single clone Southeast Asian parasite lines from Sanger pf3k project (ftp://ngs.sanger.ac.uk/production/pf3k/release5/). For pfmrp2, the selection was against minor alleles, while for pfarps10, selection was against major alleles.

No selection against the kelch13-C580Y allele conferring ART resistance

Interestingly, we do not see evidence for selection against the kelch13-C580Y allele (chr 13) that underlies resistance to ART treatment. We previously used CRISPR/Cas9 editing to add the C580Y substitution to a wild type parasite [16]. Head-to-head competition experiments revealed strong fitness costs (s = 0.15/asexual cycle) associated with this substitution. In agreement, Straimer et al [15] conducted similar experiments with Cambodian parasites: they showed that the addition of the C580Y resulted in strong fitness costs for some parasites, but had no fitness impact in recently isolated Cambodian parasites. These data also suggest that epistatic interactions with other loci may compensate and restore parasite fitness. We suspect that this may also be the case in our experiment. Technical considerations & caveats.

Maximizing statistical power

Our statistical power to detect QTLs is limited by the number of recombinants generated. In our experiment, the mouse was infected with 204 mosquitoes carrying on average three oocysts. Given that each oocyst is expected to contain sporozoites representing up to four different genotypes (i.e. a tetrad), the number of sporozoite genotypes is 204 × 3 × 4 = 2448 in this cross. We can increase the power of these experiments using mosquitoes with higher infection rates. We routinely obtain an average of 10 oocysts/mosquito, so can potentially increase numbers of recombinants by at least three-fold with the same number of mosquitoes. A second advantage of humanized mice over splenectomized chimpanzees as an infection model is that we can easily increase numbers of humanized mice used per cross. By using independent pools of mosquitoes to infect mice, we can multiply the numbers of recombinants generated, while also establishing true biological replicates of each experiment. A third advantage of the humanized mouse system is that we can stage independent crosses with different pairs of ART-R and ART-S parasites to determine if our conclusions are robust.

In this experiment, we found large numbers of inbred progeny generated by mating between male and female gametes of the same genotype. While we were able to use these to document selection against selfed progeny, this reduces the number of recombinant progeny and therefore limits statistical power for locating QTLs. For example, in our cross we estimated that 2448 sporozoite genotypes were initially used to infect the mouse. However, of these only 30% (estimated from dilution cloning of progeny) were recombinants, while the remaining 70% resulted from selfing (Button-Simmons et al. in preparation). Hence the number of independent recombinants used in this cross was 2448 × 30% = 734. A method that maximizes outcrossing would be particularly useful for future crosses. For example, aphidicolin treatment has been successfully used in rodent malaria systems to kill male gametes [57]. In yeast, BSA experiment use parental strains with different mating types to avoid inbreeding [4]. We are not yet able to do this with malaria P. falciparum crosses. The dynamics observed in our cross with selection against selfed progeny followed by allele specific selection reflects important differences between Plasmodium and yeast systems.

We expect that representation of individual parasite clones will be uneven within progeny pools. Elevated growth of particular “high fitness” clones can generate step-like changes in allele frequency at the recombination points. This has the potential to generate spurious peaks and to confuse the interpretation of BSA experiments [13, 58]. To identify such abrupt allele frequency jumps, we performed a jump-diffusion analysis as described by Abkallo et al [13]. This approach identified three allele frequency jump locations in the first experimental repeat of day 50 population, while no allele frequency jumps were found in the second experimental repeat (S3 Table and S5 Fig). One of the three jumps is located at the left end of chr 12 QTL, which indicates the possibility of chr 12 QTL generated by clonal growth. However, we detected no allele frequency jumps in the vicinity of the chr 14 QTL.

Combining BSA with cloning recombinant progeny to detect epistasis

BSA cannot be used to directly examine epistatic interactions, due to the lack of haplotype information. Fortunately, P. falciparum has a key advantage over rodent malaria systems because parasites can be grown in vitro and cloned by limiting dilution. Hence, BSA can be complemented by cloning progeny from the same genetic cross and directly examining haplotypes carrying different allele combinations. Furthermore, we can use BSA to directly test for interactions between genes. For example, we suspect that interactions between kelch13 mutations and arps10 may drive the skew observed at chr 14. This hypothesis can be directly tested by repeating the cross with parasites that have been edited to remove the kelch13 mutation or candidate arps10 mutations, to see if the skew on chr 14 disappears.

sWGA performance

The sWGA method efficiently enriches P. falciparum DNA from infected mosquito and mouse tissues, confirming the performance of this approach for enriching parasite DNA from dried blood spots [2023]. Our results further show that sWGA does not generate bias in allele frequency measurement (Fig 4). However, sWGA does have limitations with highly contaminated samples, such as early infected mosquitoes (four days post infection). DNA extracted from day 4 midguts typically contains > 99.99% mosquito DNA. Only 4.3% of sWGA products from these samples was Plasmodium DNA. In contrast, we were able to obtain > 88% of parasite DNA from sWGA, with starting material containing ≥ 1% P. falciparum DNA (Table 1).

Potential of BSA for examining selection in the mosquito stage

We did not observe allele frequency changes during mosquito infections in this experiment. We suggest two reasons for this. First, the Anopheles stephensi mosquito used is originally from urban India and widely spread across Southeast Asia, and therefore may show good compatibility with Southeast Asian parasites. Furthermore, this specific mosquito line has been long-term lab adapted, and is highly susceptible to infection with multiple parasite lines. Second, the infection period in mosquitoes in this experiment is relatively short, because we sacrificed all the mosquitos in two weeks. As a consequence, we can only detect very strong selection at this stage. However, hard selection resulting from incompatibility between parasites and mosquitoes should still be possible to detect and map in this system. We note that Molina-Cruz et al [59] were able to determine parasite QTLs for compatibility between P. falciparum and mosquitoes using parasite progeny derived from the original malaria crosses conducted in chimpanzees, providing proof-of-principal that this is possible.

Human malaria can now undergo liver stage development within humanized mice, while blood stage parasites can be grown in vitro in culture and cloned. The power of the BSA approach has been clearly demonstrated in rodent malaria, where it has been used to identify the genetic components controlling a broad range of selectable phenotypes, including virulence and immunity, growth rate and drug resistance [812]. However, human malaria parasites and rodent malaria parasites are genetically distant and human parasites show numerous unique biological features not found in rodent malaria parasites. Our approach can now be applied to directly study multiple selectable traits in the human parasite P. falciparum via genetic crosses. We anticipate that BSA will provide a powerful approach for the study of P. falciparum genetics.

Material and methods

Ethics approval and consent to participate

The study was performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health (NIH), USA. To this end, the Seattle Children’s Research Institute (SCRI) has an Assurance from the Public Health Service (PHS) through the Office of Laboratory Animal Welfare (OLAW) for work approved by its Institutional Animal Care and Use Committee (IACUC). All of the work carried out in this study was specifically reviewed and approved by the SCRI IACUC.

Preparation of genetic cross and sample collection

We generated the cross using FRG NOD huHep mice [60] with human chimeric livers and A. stephensi mosquitoes as described by Vaughan et al. [18] (see S1 File for details). We collected samples from infected mosquito midgut and salivary gland, mouse liver and in vivo blood, and in vitro blood cultures (Fig 1 and Table 1). We fed ~500 mosquitos with mixed gametocytes from each parent at equal ratio. This day was defined as day 0 for sample collecting. Forty-eight midguts were dissected at each oocyst collection time point (day 4 and day 10). The prevalence of infection was analyzed at day 10. Salivary gland were separated to collect sporozoites at day 14 after infection. Sporozoites from 204 mosquitoes were mixed together for infection into the mouse and for isolation of genomic DNA.

Six days after sporozoite injection (day 20), we injected mice intravenously with 400 μL of packed O+ huRBCs. The intravenous injection was repeated the next day (day 21). Four hours after the second huRBC injection, mice were sacrificed and blood was removed by cardiac puncture in order to recover P. falciparum–infected huRBCs. The mouse liver was dissected, immediately frozen in liquid nitrogen and then stored at −80°C.

The infected red blood cells were washed, mixed with equal volume of packed huRBCs, and resuspended in complete medium at 2% hematocrit. Two days after culture, the parasites were split equally into two wells (repeat A and repeat B) of a standard six-well plate. About 50 μL of freshly packed huRBCs were added every 2 days to each replicate. To maintain healthy cultures, serial dilutions of parasites were carried out once the parasitemia reached 4%. The cultures were maintained for 30 days in total (day21-day50), and 50ul packed red blood cells (RBCs) were collected and frozen down every 2–4 days. We also set up gametocyte enrichment cultures from day 32 progeny population with daily medium changes but no fresh huRBCs. Samples were collected 8 days (day 40) and 16 days (day 48) later.

Library preparation and sequencing

We extracted and purified genomic DNA using the Qiagen DNA mini kit, and quantified amounts using Qubit. We performed real-time quantitative PCR (qPCR) reactions to estimate the proportion of parasite genomes in each DNA sample (S1 File, S4 Table). We used selective whole genome amplification (sWGA) to enrich parasite DNA for samples obtained from infected mosquito and mouse tissues. We used selective whole genome amplification (sWGA, S1 File) to enrich parasite DNA for samples obtained from infected mosquito and mouse tissues. sWGA products were further quantified by qPCR (described above) to confirm that the majority of the products were from Plasmodium. We constructed next generation sequencing libraries using 50ng DNA or sWGA product following the KAPA HyperPlus Kit protocol with 3-cycle of PCR.

We used amplicon sequencing to trace the biases in mtDNA transmission, as sWGA with circular DNA may swamp out other sWGA products. We use at least 1000 copies of parasite genome as template for each reaction. Illumina adapters and index sequences were added to the PCR primers (S4 Table). Equal number of molecules were pooled from each reaction.

All libraries were sequenced to an average coverage of 100x using an Illumina NEXTseq 500 sequencer.

Genotype calling

We first genotyped the two parental strains. We mapped the whole-genome sequencing reads against the P. falciparum 3D7 reference genome (PlasmoDB, release32) using BWA mem (http://bio-bwa.sourceforge.net/) under the default parameters. To reduce false positives due to alignment errors, we excluded the high variable genome regions (subtelomeric repeats, hypervariable regions and centromeres) and only performed genotype calling in the 21 Mb core genome (defined in [34]). The resulting alignments were then converted to SAM format, sorted to BAM format, and deduplicated using picard tools v2.0.1 (http://broadinstitute.github.io/picard/). We used Genome Analysis Toolkit GATK v3.7 (https://software.broadinstitute.org/gatk/) to recalibrate the base quality score based on a set of verified known variants [34].

We called variants for each parent using HaplotypeCaller and then merged using GenotypeGVCFs with default parameters except for sample ploidy 1. We applied filters to the original GATK genotypes using standard filter methods described by McDew-White et al [61]. The recalibrated variant quality scores (VQSR) were calculated by comparing the raw variant distribution with the known and verified Plasmodium variant dataset. Loci with VQSR less than 1 were removed from further analysis.

We generated a “mock” genome using GATK FastaAlternateReferenceMaker from the genotype of parent NHP1337 (C580Y). The reads from bulk populations obtained at each stage of the lifecycle were mapped to this genome. Only loci with coverage > 30x were used for bulk segregant analysis. We counted reads with genotypes of each parent and calculated allele frequencies at each variable locus. Allele frequencies of NHP1337 were plotted across the genome, and outliers were removed following Hampel’s rule [62] with a window size of 100 loci (Fig 3).

Bulk segregant analysis

We performed the BSA analyses using the R package QTLseqr [63]. We first defined extreme-QTLs by looking for regions with false discovery rate (FDR) < 0.01 using the G’ approach [64]. We then calculated the Δ(SNP-index) to show the direction of the selection [65]. Once a QTL was detected, we calculated and approximate 95% confidence interval using Li’s method [66] to localize causative genes.

We also measured the fitness cost at each mutation by fitting a linear model between the natural log of the allele ratio (freq[allele1]/freq[allele2]) against time (measured in 48hr parasite asexual cycles). The slope provides a measure of the selection coefficient (s) driving each mutation [67]. The raw s values were tricube-smoothed with a window size of 100 kb to remove noise [68, 69]. A positive value of s indicates selection against alleles from the ART-R parent (NHP1337), while a negative value of s indicates selection for NHP1337 alleles.

Supporting information

S1 Fig [tif]
Number of SNPs between NHP1337 and MKK2835 in 100kb genome windows.

S2 Fig [red]

S3 Fig [tif]
Allele frequencies across the genome following gametocyte enrichment of progeny population.

S4 Fig [a]
Bulk segregant analysis by Δ (SNP-index) and G’ values.

S5 Fig [tif]
Sudden changes in allele frequency identified using a jump-diffusion model.

S6 Fig [tif]
Concordance between allele frequencies estimated before and after sWGA.

S7 Fig [tif]
PCA plot using 678 single clone infection samples from Southeast Asia.

S8 Fig [tif]
Plot of corrections between mitochondrial frequencies and frequencies of SNPs on different chromosomes across the experiment.

S1 Table [xlsx]
Summary of statistics from bulk segregant analyses.

S2 Table [xlsx]
Genes inside of QTL regions.

S3 Table [xlsx]
Sudden changes in allele frequency identified using a jump-diffusion model.

S4 Table [xlsx]
Primers used in this study.

S5 Table [xlsx]
Recombination events within the chromosome 12 and 14 QTL regions.

S6 Table [xlsx]
Ancestral/derived allelic state analysis under chr12 and chr14 QTL regions.

S1 File [docx]
Supplemental methods and materials.


Zdroje

1. Hopp CS, Chiou K, Ragheb DR, Salman AM, Khan SM, Liu AJ, et al. Longitudinal analysis of Plasmodium sporozoite motility in the dermis reveals component of blood vessel recognition. Elife. 2015;4:e07789.

2. Parts L, Cubillos F, Warringer J, Jain K, Salinas F, Bumpstead SJ, et al. Revealing the genetic structure of a trait by sequencing a population under selection. Genome research. 2011:gr. 116731.110.

3. Feng L, Jia H, Qin Y, Song Y, Liu Y, Tao S. Rapid identification of major QTLS associated with near-freezing temperature tolerance in Saccharomyces cerevisiae. Frontiers in microbiology. 2018;9:2110. doi: 10.3389/fmicb.2018.02110 30254614

4. Ehrenreich IM, Torabi N, Jia Y, Kent J, Martis S, Shapiro JA, et al. Dissection of genetically complex traits with extremely large pools of yeast segregants. Nature. 2010;464(7291):1039. doi: 10.1038/nature08923 20393561

5. Burga A, Ben-David E, Vergara TL, Boocock J, Kruglyak L. Fast genetic mapping of complex traits in C. elegans using millions of individuals in bulk. Nature communications. 2019;10(1):2680. doi: 10.1038/s41467-019-10636-9 31213597

6. Chevalier FD, Valentim CL, LoVerde PT, Anderson TJ. Efficient linkage mapping using exome capture and extreme QTL in schistosome parasites. BMC genomics. 2014;15(1):617.

7. Blake DP, Billington KJ, Copestake SL, Oakes RD, Quail MA, Wan K-L, et al. Genetic mapping identifies novel highly protective antigens for an apicomplexan parasite. PLoS pathogens. 2011;7(2):e1001279. doi: 10.1371/journal.ppat.1001279 21347348

8. Rosario V, Walliker D, Hall R, Beale G. Persistence of drug-resistant malaria parasites. The Lancet. 1978;311(8057):185–7.

9. Hunt P, Martinelli A, Modrzynska K, Borges S, Creasey A, Rodrigues L, et al. Experimental evolution, genetic analysis and genome re-sequencing reveal the mutation conferring artemisinin resistance in an isogenic lineage of malaria parasites. BMC genomics. 2010;11(1):499.

10. Martinelli A, Cheesman S, Hunt P, Culleton R, Raza A, Mackinnon M, et al. A genetic approach to the de novo identification of targets of strain-specific immunity in malaria parasites. Proceedings of the National Academy of Sciences. 2005;102(3):814–9.

11. Culleton R, Martinelli A, Hunt P, Carter R. Linkage group selection: rapid gene discovery in malaria parasites. Genome research. 2005;15(1):92–7. doi: 10.1101/gr.2866205 15632093

12. Pattaradilokrat S, Culleton RL, Cheesman SJ, Carter R. Gene encoding erythrocyte binding ligand linked to blood stage multiplication rate phenotype in Plasmodium yoelii yoelii. Proceedings of the National Academy of Sciences. 2009:pnas. 0811430106.

13. Abkallo HM, Martinelli A, Inoue M, Ramaprasad A, Xangsayarath P, Gitaka J, et al. Rapid identification of genes controlling virulence and immunity in malaria parasites. PLoS pathogens. 2017;13(7):e1006447. doi: 10.1371/journal.ppat.1006447 28704525

14. Petersen I, Gabryszewski SJ, Johnston GL, Dhingra SK, Ecker A, Lewis RE, et al. Balancing drug resistance and growth rates via compensatory mutations in the P lasmodium falciparum chloroquine resistance transporter. Molecular microbiology. 2015;97(2):381–95. doi: 10.1111/mmi.13035 25898991

15. Straimer J, Gnädig NF, Stokes BH, Ehrenberger M, Crane AA, Fidock DA. Plasmodium falciparum K13 mutations differentially impact ozonide susceptibility and parasite fitness in vitro. MBio. 2017;8(2):e00172–17. doi: 10.1128/mBio.00172-17 28400526

16. Nair S, Li X, Arya GA, McDew-White M, Ferrari M, Nosten F, et al. Fitness costs and the rapid spread of kelch13-C580Y substitutions conferring artemisinin resistance. Antimicrobial agents and chemotherapy. 2018;62(9):e00605–18. doi: 10.1128/AAC.00605-18 29914963

17. Walliker D, Hunt P, Babiker H. Fitness of drug-resistant malaria parasites. Acta tropica. 2005;94(3):251–9. doi: 10.1016/j.actatropica.2005.04.005 15845348

18. Vaughan AM, Pinapati RS, Cheeseman IH, Camargo N, Fishbaugher M, Checkley LA, et al. Plasmodium falciparum genetic crosses in a humanized mouse model. Nature methods. 2015;12(7):631. doi: 10.1038/nmeth.3432 26030447

19. Leichty AR, Brisson D. Selective whole genome amplification for re-sequencing target microbial species from complex natural samples. Genetics. 2014:genetics. 114.165498.

20. Guggisberg AM, Sundararaman SA, Lanaspa M, Moraleda C, González R, Mayor A, et al. Whole-genome sequencing to evaluate the resistance landscape following antimalarial treatment failure with fosmidomycin-clindamycin. The Journal of infectious diseases. 2016;214(7):1085–91. doi: 10.1093/infdis/jiw304 27443612

21. Sundararaman SA, Plenderleith LJ, Liu W, Loy DE, Learn GH, Li Y, et al. Genomes of cryptic chimpanzee Plasmodium species reveal key evolutionary events leading to human malaria. Nature communications. 2016;7:11078. doi: 10.1038/ncomms11078 27002652

22. Oyola SO, Ariani CV, Hamilton WL, Kekre M, Amenga-Etego LN, Ghansah A, et al. Whole genome sequencing of Plasmodium falciparum from dried blood spots using selective whole genome amplification. Malaria journal. 2016;15(1):597. doi: 10.1186/s12936-016-1641-7 27998271

23. Cowell AN, Loy DE, Sundararaman SA, Valdivia H, Fisch K, Lescano AG, et al. Selective whole-genome amplification is a robust method that enables scalable whole-genome sequencing of Plasmodium vivax from unprocessed clinical samples. MBio. 2017;8(1):e02257–16. doi: 10.1128/mBio.02257-16 28174312

24. Ariey F, Witkowski B, Amaratunga C, Beghain J, Langlois A-C, Khim N, et al. A molecular marker of artemisinin-resistant Plasmodium falciparum malaria. Nature. 2014;505(7481):50. doi: 10.1038/nature12876 24352242

25. Anderson TJ, Nair S, McDew-White M, Cheeseman IH, Nkhoma S, Bilgic F, et al. Population parameters underlying an ongoing soft sweep in southeast Asian malaria parasites. Molecular biology and evolution. 2016;34(1):131–44. doi: 10.1093/molbev/msw228 28025270

26. Fairhurst RM, Dondorp AM. Artemisinin-resistant Plasmodium falciparum malaria. Microbiology spectrum. 2016;4(3).

27. Project MPfC. Genomic epidemiology of artemisinin resistant malaria. elife. 2016;5:e08714. doi: 10.7554/eLife.08714 26943619

28. Takala-Harrison S, Jacob CG, Arze C, Cummings MP, Silva JC, Dondorp AM, et al. Independent emergence of artemisinin resistance mutations among Plasmodium falciparum in Southeast Asia. The Journal of infectious diseases. 2014;211(5):670–9. doi: 10.1093/infdis/jiu491 25180241

29. Imwong M, Suwannasin K, Kunasol C, Sutawong K, Mayxay M, Rekol H, et al. The spread of artemisinin-resistant Plasmodium falciparum in the Greater Mekong subregion: a molecular epidemiology observational study. The Lancet Infectious Diseases. 2017;17(5):491–7. doi: 10.1016/S1473-3099(17)30048-8 28161569

30. Cerqueira GC, Cheeseman IH, Schaffner SF, Nair S, McDew-White M, Phyo AP, et al. Longitudinal genomic surveillance of Plasmodium falciparum malaria parasites reveals complex genomic architecture of emerging artemisinin resistance. Genome biology. 2017;18(1):78. doi: 10.1186/s13059-017-1204-4 28454557

31. Miotto O, Amato R, Ashley EA, MacInnis B, Almagro-Garcia J, Amaratunga C, et al. Genetic architecture of artemisinin-resistant Plasmodium falciparum. Nature genetics. 2015;47(3):226. doi: 10.1038/ng.3189 25599401

32. Phyo AP, Nkhoma S, Stepniewska K, Ashley EA, Nair S, McGready R, et al. Emergence of artemisinin-resistant malaria on the western border of Thailand: a longitudinal study. The Lancet. 2012;379(9830):1960–6.

33. Ashley EA, Dhorda M, Fairhurst RM, Amaratunga C, Lim P, Suon S, et al. Spread of artemisinin resistance in Plasmodium falciparum malaria. New England Journal of Medicine. 2014;371(5):411–23. doi: 10.1056/NEJMoa1314981 25075834

34. Miles A, Iqbal Z, Vauterin P, Pearson R, Campino S, Theron M, et al. Indels, structural variation, and recombination drive genomic diversity in Plasmodium falciparum. Genome research. 2016;26(9):1288–99. doi: 10.1101/gr.203711.115 27531718

35. Miotto O, Almagro-Garcia J, Manske M, MacInnis B, Campino S, Rockett KA, et al. Multiple populations of artemisinin-resistant Plasmodium falciparum in Cambodia. Nature genetics. 2013;45(6):648. doi: 10.1038/ng.2624 23624527

36. Rosenberg R, Rungsiwongse J. The number of sporozoites produced by individual malaria oocysts. The American journal of tropical medicine and hygiene. 1991;45(5):574–7. doi: 10.4269/ajtmh.1991.45.574 1951866

37. Schneider P, Greischar MA, Birget PL, Repton C, Mideo N, Reece SE. Adaptive plasticity in the gametocyte conversion rate of malaria parasites. PLoS pathogens. 2018;14(11):e1007371. doi: 10.1371/journal.ppat.1007371 30427935

38. Ngotho P, Soares AB, Hentzschel F, Achcar F, Bertuccini L, Marti M. Revisiting gametocyte biology in malaria parasites. FEMS Microbiology Reviews. 2019.

39. Peatey CL, Dixon MW, Gardiner DL, Trenholme KR. Temporal evaluation of commitment to sexual development in Plasmodium falciparum. Malaria journal. 2013;12(1):134.

40. Talman AM, Domarle O, McKenzie FE, Ariey F, Robert V. Gametocytogenesis: the puberty of Plasmodium falciparum. Malaria journal. 2004;3(1):24.

41. Charlesworth D, Willis JH. The genetics of inbreeding depression. Nature reviews genetics. 2009;10(11):783. doi: 10.1038/nrg2664 19834483

42. Whitlock MC, Ingvarsson PK, Hatfield T. Local drift load and the heterosis of interconnected populations. Heredity. 2000;84(4):452.

43. Anderson TJ, Haubold B, Williams JT, Estrada-Franco § JG, Richardson L, Mollinedo R, et al. Microsatellite markers reveal a spectrum of population structures in the malaria parasite Plasmodium falciparum. Molecular biology and evolution. 2000;17(10):1467–82. doi: 10.1093/oxfordjournals.molbev.a026247 11018154

44. Nkhoma SC, Nair S, Al‐Saai S, Ashley E, McGready R, Phyo AP, et al. Population genetic correlates of declining transmission in a human pathogen. Molecular ecology. 2013;22(2):273–85. doi: 10.1111/mec.12099 23121253

45. Hott A, Casandra D, Sparks KN, Morton LC, Castanares G-G, Rutter A, et al. Artemisinin-resistant Plasmodium falciparum parasites exhibit altered patterns of development in infected erythrocytes. Antimicrobial agents and chemotherapy. 2015;59(6):3156–67. doi: 10.1128/AAC.00197-15 25779582

46. Vaughan AM, Mikolajczak SA, Camargo N, Lakshmanan V, Kennedy M, Lindner SE, et al. A transgenic Plasmodium falciparum NF54 strain that expresses GFP–luciferase throughout the parasite life cycle. Molecular and biochemical parasitology. 2012;186(2):143–7. doi: 10.1016/j.molbiopara.2012.10.004 23107927

47. Lynch M. The genetic interpretation of inbreeding depression and outbreeding depression. Evolution. 1991;45(3):622–9. doi: 10.1111/j.1558-5646.1991.tb04333.x 28568822

48. Coyne JA. I SPECIATJON. 2004.

49. Noble LM, Chelo I, Guzella T, Afonso B, Riccardi DD, Ammerman P, et al. Polygenicity and epistasis underlie fitness-proximal traits in the Caenorhabditis elegans multiparental experimental evolution (CeMEE) panel. 2017;207(4):1663–85. doi: 10.1534/genetics.117.300406 29066469

50. Zhu SJ, Hendry JA, Almagro-Garcia J, Pearson RD, Amato R, Miles A, et al. The origins and relatedness structure of mixed infections vary with local prevalence of P. falciparum malaria. bioRxiv. 2018:387266.

51. Nkhoma SC, Trevino SG, Gorena KM, Nair S, Khoswe S, Jett C, et al. Resolving within-host malaria parasite diversity using single-cell sequencing. bioRxiv. 2018:391268.

52. Rijpma SR, van der Velden M, González‐Pons M, Annoura T, van Schaijk BC, van Gemert GJ, et al. Multidrug ATP‐binding cassette transporters are essential for hepatic development of Plasmodium sporozoites. Cellular microbiology. 2016;18(3):369–83. doi: 10.1111/cmi.12517 26332724

53. Raj DK, Mu J, Jiang H, Kabat J, Singh S, Sullivan M, et al. Disruption of a Plasmodium falciparum multidrug resistance-associated protein (PfMRP) alters its fitness and transport of antimalarial drugs and glutathione. Journal of Biological Chemistry. 2009;284(12):7687–96. doi: 10.1074/jbc.M806944200 19117944

54. Mu J, Ferdig MT, Feng X, Joy DA, Duan J, Furuya T, et al. Multiple transporters associated with malaria parasite responses to chloroquine and quinine. Molecular microbiology. 2003;49(4):977–89. doi: 10.1046/j.1365-2958.2003.03627.x 12890022

55. Dahlström S, Veiga MI, Mårtensson A, Björkman A, Gil JP. Polymorphism in PfMRP1 (Plasmodium falciparum multidrug resistance protein 1) amino acid 1466 associated with resistance to sulfadoxine-pyrimethamine treatment. Antimicrobial agents and chemotherapy. 2009;53(6):2553–6. doi: 10.1128/AAC.00091-09 19364873

56. van der Velden M, Rijpma SR, Verweij V, van Gemert G-J, Chevalley-Maurel S, van de Vegte-Bolmer M, et al. Protective Efficacy Induced by Genetically Attenuated Mid-to-Late Liver-Stage Arresting Plasmodium berghei Δmrp2 Parasites. The American journal of tropical medicine and hygiene. 2016;95(2):378–82. doi: 10.4269/ajtmh.16-0226 27296385

57. Ramiro RS, Khan SM, Franke-Fayard B, Janse CJ, Obbard DJ, Reece SE. Hybridization and pre-zygotic reproductive barriers in Plasmodium. Proc R Soc B. 2015;282(1806):20143027. doi: 10.1098/rspb.2014.3027 25854886

58. Vázquez-García I, Salinas F, Li J, Fischer A, Barré B, Hallin J, et al. Clonal heterogeneity influences the fate of new adaptive mutations. Cell reports. 2017;21(3):732–44. doi: 10.1016/j.celrep.2017.09.046 29045840

59. Molina-Cruz A, Garver LS, Alabaster A, Bangiolo L, Haile A, Winikor J, et al. The human malaria parasite Pfs47 gene mediates evasion of the mosquito immune system. Science. 2013;340(6135):984–7. doi: 10.1126/science.1235264 23661646

60. Azuma H, Paulk N, Ranade A, Dorrell C, Al-Dhalimy M, Ellis E, et al. Robust expansion of human hepatocytes in Fah−/−/Rag2−/−/Il2rg−/− mice. Nature biotechnology. 2007;25(8):903. doi: 10.1038/nbt1326 17664939

61. McDew-White M, Li X, Nkhoma SC, Nair S, Cheeseman I, Anderson TJ. Mode and tempo of microsatellite length change in a malaria parasite mutation accumulation experiment. bioRxiv. 2019:560516.

62. Davies L, Gather U. The identification of multiple outliers. Journal of the American Statistical Association. 1993;88(423):782–92.

63. Mansfeld BN, Grumet R. QTLseqr: An R package for bulk segregant analysis with next-generation sequencing. The Plant Genome. 2018.

64. Magwene PM, Willis JH, Kelly JK. The statistics of bulk segregant analysis using next generation sequencing. PLoS computational biology. 2011;7(11):e1002255. doi: 10.1371/journal.pcbi.1002255 22072954

65. Takagi H, Abe A, Yoshida K, Kosugi S, Natsume S, Mitsuoka C, et al. QTL‐seq: rapid mapping of quantitative trait loci in rice by whole genome resequencing of DNA from two bulked populations. The Plant Journal. 2013;74(1):174–83. doi: 10.1111/tpj.12105 23289725

66. Li H. A quick method to calculate QTL confidence interval. Journal of genetics. 2011;90(2):355–60. 21869489

67. Dykhuizen D, Hartl DL. Selective neutrality of 6PGD allozymes in E. coli and the effects of genetic background. Genetics. 1980;96(4):801–17. 7021316

68. Nadaraya EA. On estimating regression. Theory of Probability & Its Applications. 1964;9(1):141–2.

69. Watson GS. Smooth regression analysis. Sankhyā: The Indian Journal of Statistics, Series A. 1964:359–72.

Štítky
Genetika Reprodukční medicína

Článek vyšel v časopise

PLOS Genetics


2019 Číslo 10
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy Podcasty Doporučená témata Časopisy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#