https://orcid.org/0000-0002-7667-1025
Figures
Abstract
Individuals acquire immunity to clinical malaria after repeated Plasmodium falciparum infections. Immunity to disease is thought to reflect the acquisition of a repertoire of responses to multiple alleles in diverse parasite antigens. In previous studies, we identified polymorphic sites within individual antigens that are associated with parasite immune evasion by examining antigen allele dynamics in individuals followed longitudinally. Here we expand this approach by analyzing genome-wide polymorphisms using whole genome sequence data from 140 parasite isolates representing malaria cases from a longitudinal study in Malawi and identify 25 genes that encode possible targets of naturally acquired immunity that should be validated immunologically and further characterized for their potential as vaccine candidates.
Author summary
Each year, there are hundreds of millions of cases of malaria resulting in several hundred thousand deaths. In malaria-endemic areas with high transmission, individuals experience malaria illness multiple times during their lifetimes, and after many infections, develop immunity that prevents symptoms. The proteins targeted by acquired immunity are not fully known. Understanding which proteins are targets of protective immune responses may aid in development of a malaria vaccine. Here, we used whole-genome sequence data from malaria parasites collected from symptomatic individuals in Malawi to identify targets of malaria immunity based on the frequency of different sequence variants observed in people with different levels of immunity and in individuals over time. Using the combined results of these approaches, we identified genes encoding 25 proteins that may be targets of clinical immunity to malaria that will be further investigated in future studies for their potential as vaccine candidate antigens.
Citation: Shah Z, Naung MT, Moser KA, Adams M, Buchwald AG, Dwivedi A, et al. (2021) Whole-genome analysis of Malawian Plasmodium falciparum isolates identifies possible targets of allele-specific immunity to clinical malaria. PLoS Genet 17(5): e1009576. https://doi.org/10.1371/journal.pgen.1009576
Editor: Taane G. Clark, London School of Hygiene and Tropical Medicine, UNITED KINGDOM
Received: November 6, 2020; Accepted: May 4, 2021; Published: May 25, 2021
Copyright: © 2021 Shah et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: Sequencing data are available at NCBI under the project numbers and accession numbers listed in S6 Table.
Funding: This work was supported by funding from the following awards granted by the National Institutes of Health: R01AI101713 (ST-H), R01AI125579 (ST-H), R01AI141900 (JCS), U19AI110820 (PI: Fraser, Project Lead: JCS), K24AI114996 (MKL), and the Malawi International Center of Excellence for Malaria Research U19AI089683 (DPM), as well as an award from the National Health and Medical Research Council: NHMRC APP1161066 (AEB, ST-H). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Despite recent progress in reducing the burden of malaria, this disease remains a leading cause of mortality worldwide, resulting in an estimated 409,000 deaths in 2019 [1]. In areas with high transmission of Plasmodium falciparum, individuals develop immunity to malaria [2]. This immunity does not provide sterile protection against all infections, but decreases the risk of clinical disease, and increases with age as individuals are repeatedly exposed to the parasite [3,4]. This age-related pattern of immunity to disease is thought to reflect the need for a repertoire of immune responses to multiple alleles in diverse parasite antigens [2,3,5,6]. The extensive genetic diversity in P. falciparum surface antigens is thought to have evolved over millennia as a means of parasite immune evasion [7]. Allele-specific immune responses have been demonstrated for several parasite antigens [8–15]. In previous work, we examined parasite alleles in repeated infections occurring in individuals followed longitudinally and identified specific polymorphic sites within parasite surface antigens (i.e. AMA1 and MSP1) where amino acid changes were associated with immune escape and increased risk of disease, consistent with allele-specific acquisition of immunity to these antigens [16,17]. Furthermore, malaria subunit vaccines based on a single antigen allele have displayed greater efficacy against parasites with alleles matching the vaccine strain compared to the diverse alleles observed in natural parasite populations [17–20]. Such allele-specific vaccine efficacy could lead to poor overall vaccine efficacy when the vaccine target allele is at low frequency in the parasite population and could result in selection of non-vaccine alleles capable of vaccine escape [18]. Overcoming this scenario may require the design of a multivalent malaria vaccine [18]. However, the design of such a vaccine is hampered by an incomplete knowledge of which parasite proteins are targets of acquired natural immunity.
The advent of technologies allowing whole genome sequencing at epidemiological scales has allowed investigators to transition from investigation of single antigens to performing genome-wide screens to identify loci likely to be involved in the acquisition of protective immunity to malaria, including uncharacterized genes encoding products of unknown function [21,22]. Although there have been previous genome-wide studies in P. falciparum to identify immune targets based on genomic signatures of balancing or diversifying selection at a population level [23–27], these studies do not relate identified signatures with clinical outcomes in individuals in endemic areas, making it difficult to directly link such signatures to immune selection. In addition, identification of correlates of protection based strictly on immunological approaches is challenging, as it can be difficult to differentiate responses reflective of exposure versus protection [28]. Expanding on our previous approach where we examined the dynamics of vaccine antigen alleles in individuals’ repeated infections over time in relation to the development of symptoms[16,17], we compared whole genome sequence data generated from P. falciparum infections collected from participants in a longitudinal cohort study conducted in Malawi to identify targets of allele-specific immunity to malaria. Specifically, we compared the frequency of parasite alleles in symptomatic infections occurring in individuals with different levels of malaria immunity to identify significantly differentiated sites, and also compared alleles in repeated infections within an individual versus between individuals to identify polymorphic sites that vary most within individuals. Genes identified using both approaches were considered possible immune targets and were further examined for their potential as vaccine candidates. As a proof of concept of the utility of our approach in identifying targets of allele-specific immune responses, we compared the frequency of alleles in one of the identified antigens (previously considered as a potential vaccine candidate) in individuals with different levels of malaria immunity to test the hypothesis that individuals who are more immune become ill when infected with a parasite having rarer alleles to which they have not yet developed immunity.
Results
Participant/Infection characteristics and definition of immune status
To identify targets of allele-specific immunity to malaria, we generated whole-genome sequence data from 140 parasite isolates collected from symptomatic infections occurring in participants in a longitudinal cohort study in Malawi [29].
Although age is often used as a proxy for immune status in high transmission areas [14,30–32], this metric does not account for heterogeneous exposure to infectious mosquito bites, which has been observed in endemic areas [33–38]. To better account for heterogeneous exposure at an individual level, we used the proportion of total infections that were symptomatic over the two-year study period to categorize immune status, using the median as a cutoff to define high and low immunity groups. Although these groups consisted of individuals with a range of ages, the median age (13.2 years) of individuals in the group with high immunity was significantly greater than the median age (7.3 years) of individuals in the low immunity group (p-value = 0.0002, Wilcoxon-Mann-Whitney test) (Fig 1 and Table 1), as would be expected in a high malaria transmission setting such as Malawi.
Scatterplot, including linear regression line (blue), shows the relationship between the proportion of symptomatic infections per individual over the course of the study and age of the individual at enrollment. The dashed line shows median proportion of symptomatic infections, which was used to define the high and low immunity groups.
Only one infection from each individual was included in comparisons between the high and low immunity groups, with samples selected in a manner to reduce temporal variability between infections (see Methods). DEploid-IBD [39] was used to estimate the proportion of each parasite clone within an infection. Infections without a predominant clone (i.e., where the majority clone had a frequency <60% within the infection) were defined as complex infections. Although the median frequency of the majority clone was not significantly different between infections in the two immunity groups (S1 Fig, p-value = 0.34, Wilcoxon-Mann-Whitney test), the high immunity group had a greater number of complex infections (n = 7) than the low immunity group (n = 2). These nine complex infections were excluded from further analysis to avoid confounding by infection complexity and misclassification of alleles likely contributing to clinical illness.
The median parasite density of infections in the high immunity group was significantly lower than in the low immunity group (Table 1, p-value = 7.74 x 10−06, Wilcoxon-Mann-Whitney test). However, there was no significant difference in the percentage of the parasite genome with at least 20-fold coverage (Table 1, p-value = 0.27, Wilcoxon-Mann-Whitney test), or in the median average depth of coverage (Table 1, p-value = 0.25, Wilcoxon-Mann-Whitney test) between whole genome sequence data generated from infections in the two groups.
Differentiated loci between groups with different levels of immunity to clinical malaria
We hypothesized that, because of allele-specific immune responses, individuals with greater protective immunity to malaria would experience disease when infected with parasite antigen alleles that are rarer in the parasite population, having already developed immunity to more common alleles circulating in the population. Thus, we would expect significant genetic differentiation between antigen alleles in individuals with high versus low immunity at loci that are targets of allele-specific immunity.
To test this hypothesis, Wright’s fixation index (FST), a measure of genetic differentiation between two populations [40], was estimated per non-synonymous single nucleotide polymorphism (SNP) to identify genetically differentiated sites between parasites from the high and low immunity groups with the significance threshold determined through resampling (10,000 permutations). We identified 160 sites (in 145 genes) in the parasite genome that were significantly differentiated between the two immunity groups (p-value ≤ 0.009, Fig 2A and S1 Table). Fifty-five of the genes containing significantly differentiated sites (38%) encode proteins of unknown function that are not associated with any computed or curated molecular function or biological process based on Gene Ontology (S1 Table). These 145 gene products included some proteins previously identified as potential vaccine candidates, including AMA1 (apical membrane antigen 1) [17,20], ASP (apical sushi protein, PF3D7_0405900) [41], CLAG8 (cytoadherence linked asexual protein 8, PF3D7_0831600) [42,43], SLARP (sporozoite and liver asparagine-rich protein, PF3D7_1147000) [44] and a conserved protein of unknown function (PF3D7_1359000) [45].
a) Genome-wide genetic differentiation (FST) between parasites from individuals with higher immunity vs. lower immunity. Each point represents a variable, non-synonymous site. Results are plotted as–log10 p-values on the y-axis. The color of each point represents the FST value, with darker points indicating higher FST values. The dashed line denotes statistical significance (p-value = 0.0095), with p-value determined by permutation. b) Nucleotide diversity for significantly differentiated SNPs in parasites from individuals with higher immunity and lower immunity. c) Box-plot of mean allele frequency per individual based on SNPs in the whole genome sequences, and which are significantly differentiated SNPs from (a). Red indicates the high immunity group and blue color indicates the low immunity group.
To further test the hypothesis that individuals who are more immune become symptomatic when infected with parasites having antigen alleles that are rarer in the parasite population, we estimated nucleotide diversity at significantly differentiated non-synonymous sites in parasites from the two immunity groups and observed a significantly greater median nucleotide diversity in parasites from the high immunity group compared to parasites from the low immunity group (Fig 2B, p-value = 4.74 x 10−05, Wilcoxon-Mann-Whitney test). In addition, we estimated the average frequency of alleles in each infection at both significantly differentiated sites and genome-wide variable sites. The median frequency of alleles at genome-wide variable sites was not significantly different between immunity groups (Fig 2C, p-value = 0.88, Wilcoxon-Mann-Whitney test); in contrast, at the differentiated sites, the median frequency of alleles was significantly lower in the high immunity group compared to the low immunity group (Fig 2C, p-value = 0.007, Wilcoxon-Mann-Whitney test). These results are consistent with the scenario that individuals in the high immunity group are infected with parasites having different lower-frequency alleles compared to individuals in the low immunity group who are infected with parasites sharing more common alleles.
Within 23 polyclonal infections, we also compared the proportion of mismatched alleles between the predominant and minor clones at both significantly differentiated sites and genome-wide variable sites in order to assess whether these clones differ at sites thought to be relevant for immunity. We observed a significantly greater median proportion of mismatches between the major and minor clones within an infection at the differentiated sites compared to genome-wide variable sites (S2 Fig, p-value = 8 x 10−05, Wilcoxon-Mann-Whitney test), consistent with the hypothesis that the predominant clone represents a breakthrough infection that has escaped allele-specific immune responses that maintain minor clones at a subclinical level.
Loci that differ more within individuals than between individuals
We expected that allele-specific immune responses would result in a greater proportion of genetic differences in parasites causing multiple symptomatic infections within an individual compared to parasites causing infection in different individuals, at antigenic loci that are targets of immunity. To identify regions of the parasite genome that are most different in infections occurring within an individual versus between individuals, we compared the allelic states at each variable non-synonymous site across the genome between pairs of isolates sampled within an individual and between individuals and estimated the proportion of mismatches per site in each group (Fig 3A). The distribution of the difference in the proportion of mismatches across all non-synonymous variable sites between the two groups (within minus between) is shown in Fig 3B (see also S3A Fig). The difference in the proportion of mismatches at each site had a median and mode equal to zero, indicating that the proportion of mismatches was not different within and between individuals at most sites. There was no significant correlation between the number of days between infections in a pair and the proportion of mismatches (S3B Fig, Pearson’s correlation r = 0.065, p-value = 0.3), suggesting that time was not a significant confounding factor in the analysis. We further examined the top 1% of sites that differed most within individuals compared to between individuals, which included 223 SNPs, located in 173 genes (S2 Table). Sixty-eight (39%) of these 173 genes encode proteins of unknown function that are not associated with any computed or curated molecular function or biological process based on Gene Ontology, and at least 15 (8.7%) encode for proteins that have been previously identified as potential vaccine candidates.
a) Illustration of analysis to identify regions of the genome that vary more in parasites causing illness within the same individual over time (within individuals) compared to random pairs of parasites in the population (between individuals). b) Distribution of differences in the proportion of mismatched alleles in the within group and the between group. The difference was calculated as the proportion of mismatches at each non-synonymous SNP in the within group minus the proportion of mismatches at each non-synonymous SNP in the between group. The dashed black line indicates the threshold for the top 1% most different SNPs in the within group compared to the between group.
Loci identified as possible targets of immunity by both analytical approaches
Twenty-five genes were identified by both analytical approaches, of which 11 (44%) encode proteins of unknown function (Table 2). Based on publicly available data, 20 of the 25 genes have a moderate to high level of expression in the erythrocytic stage of the parasite life cycle [46] (S3 Table). Eight genes have a mutagenesis index score near zero, suggesting that they are likely essential [47], and at least 12 genes have either a known or a predicted transmembrane domain or a signal peptide [48] (S4 Table). Among these 25 gene products are proteins whose putative function make them biologically plausible vaccine candidates, including SURFIN4.2 (thought to be involved in formation of the moving junction during erythrocyte invasion) [49], as well as members of the clag and phist multigene families (both thought to have a role in parasite cytoadherence) [50,51].
As a proof of concept of the utility of our approach for identifying targets of allele-specific immunity, we used publicly available sequence data from Pf3K [52] combined with 156 sequences from Papua New Guinea (PNG) [53,54] from the MalariaGEN P. falciparum Community Project [55] to examine global diversity in clag8, one of the 25 identified genes that has been considered previously as a potential vaccine candidate, and found high nucleotide diversity, Tajima’s D, and number of segregating sites in the C-terminal region of clag8. This pattern was consistent when comparing data from parasites collected in Malawi to data from parasites from other regions of the world (Fig 4A). A haplotype network generated using clag8 sequences from multiple geographic areas showed no evidence of regional adaptation (Fig 4B), including low FST values between clag8 sequences from Africa, Asia and PNG (S5 Table). Further analysis of CLAG8 protein sequences show regions with high protein disorder and B-cell epitope sites, especially in the C-terminal region (S4 Fig). To test our initial hypothesis that immune individuals become ill when infected with parasite antigen alleles that are rarer in the parasite population, we compared the frequency of clag8 haplotypes based on the C-terminal region of the gene in individuals with different levels of immunity and observed that parasites causing illness in individuals in the high immunity group tended to have clag8 haplotypes that were lower in frequency compared to individuals in the low immunity group, although significance was borderline (S5 Fig, p-value = 0.09, Wilcoxon-Mann-Whitney test).
a) Tajima’s D along the clag8 gene in samples from Thailand, PNG, Malawi, and Ghana. The dotted red lines represent a positive Tajima’s D value (≥ 2), suggestive of balancing selection. b) Haplotype network of clag8 using P. falciparum sequence data.
Discussion
Previous studies have shown evidence of allele-specific acquisition of immunity to P. falciparum in single genes or proteins [8–10,12,16,17] that have been identified as potential vaccine antigens based on traditional vaccinology approaches that empirically identify immunogenic proteins. Other studies have performed genome-wide screens to identify antigens based on genomic signatures of balancing selection but have not related these signatures to clinical outcomes that might link specific signatures with protective immune responses [23–27]. In this study, we conducted a genome-wide, individual-level analysis to identify targets of allele-specific immunity to clinical malaria using P. falciparum whole genome sequence data by identifying parasite genes that are genetically differentiated between individuals with different levels of immunity to malaria and genomic regions that are most different in parasites causing illness within the same individual versus between individuals. Twenty-five genes were identified using both analytical approaches and encode possible targets of allele-specific acquired immunity to clinical malaria, including genes thought to be involved in erythrocyte invasion and cytoadherence, among other functions. Examination of global diversity in clag8, a gene encoding a protein that has previously been considered as a vaccine candidate, provided evidence of immune selection in the C-terminal region of the gene and did not indicate geographical differences that might be indicative of local adaptation or genetic drift. Consistent with the hypothesis of allele-specific acquisition of immunity, clag8 haplotypes of parasites causing illness in individuals with greater immunity were lower in frequency compared to those causing illness in less immune individuals. These findings support the utility of our approach for identification of targets of allele-specific immunity and further investigation of these 25 loci as potential vaccine candidates.
In this study, we estimated the genetic complexity of infections and observed that individuals in the group with higher immunity generally had more complex infections than individuals with lower immunity. Although not statistically significant, this pattern is broadly consistent with results from other studies that have reported associations between infection complexity and either age and/or risk of disease [16,56]. The greater complexity of infections in individuals with higher immunity is consistent with the idea that immune individuals are capable of maintaining some parasite clones at a subclinical level and could suggest a possible role of polyclonal infections in maintaining immunity through continuous exposure to different strains [57,58]. Alternatively, the difference in complexity of infections between the two groups could also be reflective of transmission levels within the community, rather than just reflective of immunity. In polyclonal infections from our data set, we observed that the predominant clone was more likely to have different alleles than the minor clone at sites identified as potential targets of immunity in our analysis. This finding is consistent with the hypothesis that minority clones are maintained at a subclinical level by acquired immunity, while the predominant clone is able to escape the immune response (because it has unrecognized alleles), resulting in a symptomatic infection. Prior to downstream analysis, infections lacking a predominant clone were excluded to avoid confounding by infection complexity and misclassification of clones likely responsible for disease symptoms. The exclusion of infections lacking a predominant clone led to the removal of a greater proportion of infections from the high immunity group than the low immunity group, which may have resulted in an underestimation of parasite diversity in the high immunity group. However, we do not believe this underestimation impacted our conclusions, as infections within the high immunity group were significantly more diverse at differentiated loci even after exclusion of these highly complex infections.
We hypothesized that individuals who have a higher degree of immunity to malaria would have a greater risk of malaria symptoms when infected with parasites having protein variants that are of lower frequency in the local parasite population, owing to their having already acquired immunity to variants encoded by more common alleles. When we estimated nucleotide diversity at differentiated loci in both immunity groups, we found significantly greater diversity in the group with higher immunity compared to the group with lower immunity. This difference in diversity may reflect the fact that individuals with higher immunity have symptomatic infections with parasite alleles that are rarer in the population and therefore more likely to be different from one another at differentiated loci. This scenario is also supported by our finding that the median frequency of the infecting allele was significantly lower in more immune individuals compared to less immune individuals when examining differentiated loci. These results also demonstrate the importance of accounting for rare alleles in vaccine design, as they may also lead to escape from vaccine-induced immunity.
Genome-wide screens for signatures of balancing selection have previously identified some of the 25 genes that we identified in this study, including PF3D7_0710200, clag8 and other genes from the surfin and phist multigene families [23–25,27]. PF3D7_0710200, a gene encoding a conserved protein of unknown function, was identified by four separate studies as a potential immune target [23–25,27]; however, little is known about the function of this protein. Indeed, 11 of the 25 genes (44%) identified in this study encode for proteins of unknown function, highlighting the importance of further genetic screens to determine the role of such genes, which make up ~35% of the parasite genome [59]. A recent study has suggested that most genes involved in host-parasite interactions, including many known antigens, are non-essential [47]; however, eight of the 25 genes identified in this study would be considered essential genes, based on their low mutagenesis index score [47]. Such essential genes may be attractive vaccine targets as there may be less redundancy in function that would allow the parasite to adapt and escape vaccine-induced inhibition.
As a proof of concept of the utility of our approach for identifying targets of allele-specific immunity, we further examined the global diversity of one of the 25 identified genes, clag8, which has previously been considered as a malaria vaccine candidate. clag8 belongs to the clag multigene family and is one of the least studied genes in the clag family. clag8 is highly expressed during the first few hours following erythrocyte invasion (early ring stage), displays decreased expression from 5 hours to 35 hours post invasion, and then increased expression during the schizont stage [46]. It is thought to be part of the RhopH complex, formed by the members of rhoph1/clag gene families. The RhopH complex is an erythrocyte-binding protein complex inside the rhoptry and has been suggested to play an important role in establishment of the parasitophorous vacuole [43,60,61]. Previous studies have reported evidence of positive diversifying selection in this gene, with high nucleotide diversity and a high proportion of non-synonymous substitutions per site (dN) [43]. Our global analysis of clag8 diversity displayed high nucleotide diversity and Tajima’s D values in the C-terminal region of the gene, as well as evidence of high protein disorder and predicted B-cell epitopes in the C-terminal region of the protein, suggesting that it is likely to be immunogenic [62]. These results were consistent in all the countries included in the dataset, which represented parasite isolates from three major malaria endemic regions [52]. Additionally, clag8 haplotypes displayed no evidence of geographical adaptation, with major haplotypes being observed in all geographic areas, and seemed to cluster into three main groups. Further studies to identify functional epitopes within the protein and potential cross-reactivity are necessary to determine whether related haplotypes can be grouped into serotypes for the purpose of designing a broadly protective vaccine. At an individual level, and in support of our overarching hypothesis, individuals with higher immunity to malaria were infected with clag8 haplotypes that were less frequent compared to the clag8 haplotypes infecting individuals with lower immunity.
It is noteworthy that our combined list of genes based on both analytical approaches did not include leading blood stage vaccine candidates such as AMA1 and the MSPs, although some of these genes were identified in a single approach. Using a protein microarray, Crompton et al. found that antibody responses to leading vaccine candidates, such as AMA1, MSP1, and MSP2, did not distinguish individuals who were protected from clinical infection versus those who were not in a cohort of individuals from Mali [15], suggesting the possibility that responses to these proteins are not the primary drivers of clinical immunity. Other studies [17,63,64] have supported the hypothesis that responses to antigens such as AMA1 contribute to allele-specific clinical immunity but may be short-lived. In our analyses, infection pairs within an individual were not necessarily consecutive infections, with the time between infections ranging from 27 to 696 days. Although we did not see a significant correlation between the proportion of allele mismatches and the number of days between infections in a pair, it is possible that we could have failed to identify antigens involved in allele-specific acquisition of immunity that elicit short-lived immune responses. Studies with larger sample size would likely be required to distinguish antigens that have different antibody kinetics.
In addition, analyses in this study included only the core genome, owing to our use of selective whole genome amplification to enrich for parasite DNA prior to sequencing, which has been shown to result in poor sequencing coverage in the telomeric and centromeric regions [65,66]. Limiting analysis to the core genome could have prevented us from identifying members of multigene families that may be important for development of immunity to clinical malaria, as many of these genes are located in these low-coverage regions of the genome. However, because multigene families have been implicated in severe malaria, it is possible that immunity to these diverse antigens may be more relevant to preventing severe rather than uncomplicated malaria.
Other limitations of this study include the inability to examine asymptomatic infections, owing to the difficulty in obtaining whole genome sequence data from these low parasitemia infections, as well as the focus on SNP data versus other types of variants, such as indels or structural variants, that could also contribute to allele-specific immunity. Specifically, our approach would not have detected allelic dimorphism observed in repetitive regions of some known malaria vaccine antigens, such as MSP1 and MSP2 [67]. Such low-complexity regions are difficult to resolve using short-read sequencing data as reads spanning these regions do not map uniquely to the reference genome. Targeted sequencing of specific genes with long-read sequencing platforms may allow sequencing of asymptomatic infections and examination of the contribution of additional types of variants to allele-specific immunity.
Here, we describe a promising genome-wide approach to identify potential targets of allele-specific immunity to clinical malaria. This approach identifies individual-level associations between antigen allele dynamics and patient clinical outcomes. Inferences are based on changes in parasite allele frequencies in individuals over time, potentially allowing identification of subdominant, but protective, epitopes that might otherwise be difficult to detect in immunological studies, because they are masked by responses to immunodominant, but not protective, loci. Using this approach, we identified 25 genes, many of unknown function, that encode proteins that can be further characterized for their potential as candidates for a multivalent subunit malaria vaccine through immunological epitope mapping and functional studies of antibodies elicited to these proteins.
Methods
Ethics statement
Clinical samples were collected under protocols approved by the ethics committees at the College of Medicine in Blantyre, Malawi, and the University of Maryland, Baltimore. Written informed consent was provided by study participants or their guardians.
Study design and samples
Parasite isolates were collected from participants in a longitudinal cohort study conducted in the Chikhwawa district of southern Malawi. Details about the participants and study procedures have been described previously by Buchwald et al [29]. Briefly, 120 children and adults reporting to the Mfera Health Centre with uncomplicated malaria between June 2014 and March 2015 were followed monthly over two years. Blood samples were collected at each monthly visit and all unscheduled visits where individuals reported to the Health Centre with symptoms of malaria. For each visit, parasitemia was diagnosed by both microscopy and PCR. The data analyzed in this study were generated from red blood cell pellets collected from symptomatic, uncomplicated malaria infections identified during passive follow up. The median parasitemia of the sampled infections as determined by microscopy was 21,960 parasites/μL and ranged from 0 parasites/μL (but positive by a rapid diagnostic test) to 241,260 parasites/μL. All samples were confirmed to be positive for P. falciparum by PCR. To ensure only independent infections were included in the analysis, infections within an individual separated by <14 days were excluded. DNA from red blood cell pellets was extracted using the method of Zainabadi et al [68]. Extracted DNA was enriched for parasite DNA using an optimized selective whole genome amplification approach described by Shah et al [65].
Whole genome sequencing
Genomic DNA libraries were constructed for sequencing using the KAPA Library Preparation Kit (Kapa Biosystems, Woburn, MA). DNA (≥ 200 ng) was fragmented with the Covaris E210 to ~200 bp. Libraries were prepared using a modified version of the manufacturer’s protocol. The DNA was purified between enzymatic reactions and library size selection was performed with AMPure XT beads. Libraries were assessed for concentration and fragment size using the DNA High Sensitivity Assay on the LabChip GX (Perkin Elmer, Waltham, MA). Library concentrations were also assessed by qPCR using the KAPA Library Quantification Kit. Libraries were pooled and subsequently sequenced on an Illumina HiSeq 4000 (Illumina, San Diego, CA) to generate 150 bp paired-end reads.
Read mapping and SNP calling
Sequencing data were analyzed by mapping raw fastq files to the 3D7 reference genome using Bowtie2 [69]. Binary Alignment Map (BAM) files were processed following the GATK Best Practices workflow to obtain analysis-ready reads [70,71]. Bedtools [72] was used to generate coverage and depth estimates from the processed reads, and the GATK Best Practices workflow was followed for variant calling [70,71]. Haplotype Caller was used to create genomic variant call format (GVCF) files for each sample and joint SNP Calling was performed (GATK v3.7). Variants were removed if they met the following filtering criteria: variant confidence/quality by depth (QD) < 2.0, strand bias (FS) > 60.0, root mean square of the mapping quality (MQ) < 40.0, mapping quality rank sum (MQRankSum) < -12.5, read position rank sum (ReadPosRankSum) < -8.0, quality (QUAL) < 50. Variant sites with >20% missing genotypes and samples with >30% missing data were additionally removed using vcftools. Variants were also removed if the minor allele was not present in at least two samples. Only the core genome was used for further analysis, which has been previously defined by exclusion of the highly variable telomeric and centromeric regions of the genome [73]. The median percentage of the genome covered by ≥ 20 reads was 88.9% [65]. After applying quality-control filters, 55,970 SNPs were called in the core genome, including 22,177 non-synonymous SNPs, with an average 11.6 variants called per gene.
Definition of immune status
The degree of immunity to clinical malaria was defined based on the proportion of symptomatic infections out of all P. falciparum infections experienced by each study participant over the course of the two-year study. To account for exposure, individuals with less than five total infections, including symptomatic and asymptomatic infections, were excluded from the analysis. The median proportion of symptomatic infections was used as the cutoff to categorize individuals into higher and lower immunity groups. The limited sample size of our study did not allow us to categorize immune status as an ordinal variable.
Complexity of infection and genetic differentiation
Only one infection from each individual was included in comparisons between high a low immunity groups. Infections were selected based on proximity to the median of the distribution of sampling dates to reduce temporal variability. DEploid-IBD [39] was used to estimate the proportion of each clone within an infection. Infections without a predominant clone (i.e., where the majority clone had a frequency <60% within the infection) were defined as complex infections and were excluded from downstream analysis. For the remaining samples, the major allele was called at heterozygous positions if the allele was supported by ≥70% of reads; otherwise, the genotype was coded as missing. A Wilcoxon-Mann-Whitney test was used to assess differences in the frequency of the majority clone in infections from the two immunity groups.
Vcftools [74] was used to estimate Weir and Cockerham FST in variable non-synonymous, bi-allelic sites. Significance was determined using 10,000 permutations, where the observed population was resampled without replacement. To determine the impact of performing the analysis based on the predominant allele at biallelic sites, we also performed the analysis using multiallelic sites and all alleles within an infection. Although FST values were generally higher in the analysis with multiple alleles compared to the analysis with a single major allele, sites that were significantly differentiated in the analysis based on the major allele were also significantly differentiated in the analysis where minor alleles were also included. Nucleotide diversity at significantly differentiated sites was estimated using vcftools [74]. PlasmoDB (v44) [22] was used to identify genes containing differentiated SNPs.
In all polyclonal infections, the major and minor clones (defined by clone frequencies obtained from DEploid-IBD [39]) were compared, provided clone frequency was less than 80% and greater than 10% (n = 23). At each non-synonymous site, the proportion of samples with mismatched alleles from major and minor clones was estimated. The proportion of mismatches was then compared between significantly differentiated sites and all remaining variable sites from the genome. The p-value was estimated by conducting a Wilcoxon-Mann-Whitney test to determine if there is a significant difference in mismatches between clones at different sites versus remaining genome-wide variable sites.
Paired infection analysis
Individuals with parasite whole genome sequence data from at least two symptomatic infections occurring at least 14 days apart were included in the comparison of infections occurring within the same host to infections occurring in different hosts. Multi-allelic sites were included in the analysis of paired infections, in contrast to analyses of genetic differentiation. The ‘within’ group included all pairs of parasites collected at different time points from the same individual. The ‘between group’, included all pairs of parasites from different individuals. A total of 116 samples were included in this analysis. The within group contained 124 pairs of samples and the between group contained 6546 pairs of samples. For all pairs, the allelic state was compared at each site and the proportion of pairs with non-matching allelic states was estimated by site (illustrated in Fig 3). The difference between the within group and the between group was calculated by subtracting the proportion of pairs with non-matching allelic states for each site. The p-value was estimated by conducting a one-sided z-test using the difference in proportion of mismatched alleles between the two groups. PlasmoDB [22] was used to identify genes containing the SNPs of interest.
Global diversity in clag8
The MalariaGEN Pf3K project release 5.1 data [52] was used to estimate global diversity in these genes identified in this study. The Pf3K dataset includes whole genome sequencing data from 2,512 samples collected in multiple locations in Asia and Africa. Data [53,54] from 156 additional isolates from Papua New Guinea were also included in the analysis. VaxPack (https://github.com/BarryLab01/vaxpack) was used for global population genetic analysis. GATKv4.0 was used for variant calling. Samples containing ambiguous bases were removed. Singleton SNPs were converted back to reference to prevent false positive variants. Nucleotide diversity and Tajima’s D were calculated for all polymorphic sites separately for every country that had a sample size greater than 50. Templeton, Crandall, and Sing (TCS) [75] method on PopArt [76] was used to construct the haplotype network using non-synonymous SNPs. Protein disorder region and B-cell epitope regions were predicted using PlasmoSIP [62]. The haplotype frequencies of the C-terminal region in Malawian isolates from different immunity groups were estimated for non-synonymous sites using DnaSP v6 [77].
Supporting information
S1 Fig. Proportion of the infection comprised by the majority clone in infections from the high (n = 35) and low (n = 35) immunity groups.
Samples without a predominant clone (≥ 0.6), indicated by the dashed line, were defined as complex infections and were removed from downstream analyses.
https://doi.org/10.1371/journal.pgen.1009576.s001
(TIF)
S2 Fig. Proportion of mismatches between the major and minor clones within 23 polyclonal infections.
The median proportion of mismatches was significantly greater at differentiated sites thought to be targets of allele-specific immunity compared to genome-wide sites (p-value = 8x10-05, Wilcoxon-Mann-Whitney test), consistent with the hypothesis that the predominant clone represents a breakthrough infection that has escaped allele-specific immune responses that maintain minor clones at a subclinical level.
https://doi.org/10.1371/journal.pgen.1009576.s002
(TIF)
S3 Fig.
(a) Proportion of mismatched alleles within individuals vs. between individuals. Each point is the proportion of mismatches at a non-synonymous SNP. Black points represent the top 1% most mismatched alleles within individuals. (b) Correlation between proportion of mismatched SNPs per pair within individuals (y-axis) and time between infections (x-axis). The blue line represents the linear regression line with 95% confidence region shown by the shaded region.
https://doi.org/10.1371/journal.pgen.1009576.s003
(TIF)
S4 Fig. Predicted protein disorder and B-cell epitopes in clag8.
The orange line and blocks show the linear B-cell epitope mapping score and predicted B-cell epitope sites, respectively. The blue line and blocks show the protein disorder score and highly disordered region, respectively. The asterisks along the bottom represent known SNPs.
https://doi.org/10.1371/journal.pgen.1009576.s004
(TIF)
S5 Fig. Haplotype frequency of clag8 c-terminal region (amino acid position > 1000), from Malawian whole-genome sequences used in this study, in individuals from both the immunity groups.
https://doi.org/10.1371/journal.pgen.1009576.s005
(TIF)
S1 Table. Significantly differentiated SNPs in parasites infecting individuals with different levels of immunity to clinical malaria.
https://doi.org/10.1371/journal.pgen.1009576.s006
(DOCX)
S2 Table. Top 1% most different SNPs within individuals compared to between individuals.
https://doi.org/10.1371/journal.pgen.1009576.s007
(DOCX)
S3 Table. Heatmap of Fragments Per Kilobase of transcript per Million mapped reads (FPKM) values to show expression of 25 identified genes during the intraerythrocytic stages of the parasite life-cycle.
https://doi.org/10.1371/journal.pgen.1009576.s008
(DOCX)
S4 Table. Protein/Gene features for genes identified with both analytical approaches.
https://doi.org/10.1371/journal.pgen.1009576.s009
(DOCX)
S5 Table. Pairwise genetic differentiation (FST) between clag8 sequences from Africa, Asia and Papua New Guinea (PNG).
https://doi.org/10.1371/journal.pgen.1009576.s010
(DOCX)
Acknowledgments
We thank the participants in the Mfera Cohort Study. We would like to thank Biraj Shrestha and Gillian Mbambo for their assistance with sample processing. We would also like to acknowledge Timothy D. O’Connor, Michael P. Cummings and Alexis Boleda for their valuable input and suggestions, and Terrie Taylor for her role in leading the Malawi International Center of Excellence for Malaria Research.
References
- 1.
WHO | World malaria report 2018 [Internet]. WHO. [cited 2019 May 30]. Available from: http://www.who.int/malaria/publications/world-malaria-report-2018/report/en/
- 2. Marsh K, Kinyanjui S. Immune effector mechanisms in malaria. Parasite Immunol. 2006;28(1–2):51–60. pmid:16438676
- 3. Cowman AF, Healer J, Marapana D, Marsh K. Malaria: Biology and Disease. Cell. 2016 Oct 20;167(3):610–24. pmid:27768886
- 4. Ryg-Cornejo V, Ly A, Hansen DS. Immunological processes underlying the slow acquisition of humoral immunity to malaria. Parasitology. 2016 Feb;143(2):199–207. pmid:26743747
- 5. Portugal S, Pierce SK, Crompton PD. Young Lives Lost as B Cells Falter: What We Are Learning About Antibody Responses in Malaria. J Immunol. 2013 Apr 1;190(7):3039–46. pmid:23526829
- 6. Buchwald AG, Sorkin JD, Sixpence A, Chimenya M, Damson M, Wilson ML, et al. Association Between Age and Plasmodium falciparum Infection Dynamics. Am J Epidemiol. 2019 Jan 1;188(1):169–76. pmid:30252032
- 7. Weedall GD, Conway DJ. Detecting signatures of balancing selection to identify targets of anti-parasite immunity. Trends Parasitol. 2010 Jul 1;26(7):363–9. pmid:20466591
- 8. Polley SD, Tetteh KKA, Lloyd JM, Akpogheneta OJ, Greenwood BM, Bojang KA, et al. Plasmodium falciparum Merozoite Surface Protein 3 Is a Target of Allele-Specific Immunity and Alleles Are Maintained by Natural Selection. J Infect Dis. 2007 Jan 15;195(2):279–87. pmid:17191173
- 9. Osier FHA, Fegan G, Polley SD, Murungi L, Verra F, Tetteh KKA, et al. Breadth and magnitude of antibody responses to multiple Plasmodium falciparum merozoite antigens are associated with protection from clinical malaria. Infect Immun. 2008 May;76(5):2240–8. pmid:18316390
- 10. Cortés A, Mellombo M, Masciantonio R, Murphy VJ, Reeder JC, Anders RF. Allele specificity of naturally acquired antibody responses against Plasmodium falciparum apical membrane antigen 1. Infect Immun. 2005 Jan;73(1):422–30. pmid:15618180
- 11. Barry AE, Trieu A, Fowkes FJI, Pablo J, Kalantari-Dehaghi M, Jasinskas A, et al. The Stability and Complexity of Antibody Responses to the Major Surface Antigen of Plasmodium falciparum Are Associated with Age in a Malaria Endemic Area. Mol Cell Proteomics [Internet]. 2011 Nov 1 [cited 2020 Apr 14];10(11). Available from: https://www.mcponline.org/content/10/11/M111.008326
- 12. Early AM, Lievens M, MacInnis BL, Ockenhouse CF, Volkman SK, Adjei S, et al. Host-mediated selection impacts the diversity of Plasmodium falciparum antigens within infections. Nat Commun [Internet]. 2018 Apr 11 [cited 2018 Jul 24];9. Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5895824/ pmid:29643376
- 13. Dutta S, Lee SY, Batchelor AH, Lanar DE. Structural basis of antigenic escape of a malaria vaccine candidate. Proc Natl Acad Sci. 2007 Jul 24;104(30):12488–93. pmid:17636123
- 14. Tran TM, Ongoiba A, Coursen J, Crosnier C, Diouf A, Huang C-Y, et al. Naturally Acquired Antibodies Specific for Plasmodium falciparum Reticulocyte-Binding Protein Homologue 5 Inhibit Parasite Growth and Predict Protection From Malaria. J Infect Dis. 2014 Mar 1;209(5):789–98. pmid:24133188
- 15. Crompton PD, Kayala MA, Traore B, Kayentao K, Ongoiba A, Weiss GE, et al. A prospective analysis of the Ab response to Plasmodium falciparum before and after a malaria season by protein microarray. Proc Natl Acad Sci. 2010 Apr 13;107(15):6958–63. pmid:20351286
- 16. Takala SL, Coulibaly D, Thera MA, Dicko A, Smith DL, Guindo AB, et al. Dynamics of polymorphism in a malaria vaccine antigen at a vaccine-testing site in Mali. PLoS Med. 2007 Mar;4(3):e93. pmid:17355170
- 17. Takala SL, Coulibaly D, Thera MA, Batchelor AH, Cummings MP, Escalante AA, et al. Extreme Polymorphism in a Vaccine Antigen and Risk of Clinical Malaria: Implications for Vaccine Development. Sci Transl Med. 2009 Oct 14;1(2):2ra5–2ra5. pmid:20165550
- 18. Ouattara A, Barry AE, Dutta S, Remarque EJ, Beeson JG, Plowe CV. Designing malaria vaccines to circumvent antigen variability. Vaccine. 2015 Dec 22;33(52):7506–12. pmid:26475447
- 19. Ouattara A, Takala-Harrison S, Thera MA, Coulibaly D, Niangaly A, Saye R, et al. Molecular basis of allele-specific efficacy of a blood-stage malaria vaccine: vaccine development implications. J Infect Dis. 2013 Feb 1;207(3):511–9. pmid:23204168
- 20. Thera MA, Doumbo OK, Coulibaly D, Laurens MB, Ouattara A, Kone AK, et al. A field trial to assess a blood-stage malaria vaccine. N Engl J Med. 2011 Sep 15;365(11):1004–13. pmid:21916638
- 21. Gardner MJ, Hall N, Fung E, White O, Berriman M, Hyman RW, et al. Genome sequence of the human malaria parasite Plasmodium falciparum. Nature. 2002 Oct;419(6906):498–511. pmid:12368864
- 22.
PlasmoDB: The Plasmodium Genomics Resource [Internet]. [cited 2018 May 1]. Available from: http://plasmodb.org/plasmo/
- 23. Amambua-Ngwa A, Tetteh KKA, Manske M, Gomez-Escobar N, Stewart LB, Deerhake ME, et al. Population Genomic Scan for Candidate Signatures of Balancing Selection to Guide Antigen Characterization in Malaria Parasites. PLoS Genet [Internet]. 2012 Nov 1 [cited 2019 Aug 15];8(11). Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3486833/ pmid:23133397
- 24. Mobegi VA, Duffy CW, Amambua-Ngwa A, Loua KM, Laman E, Nwakanma DC, et al. Genome-wide analysis of selection on the malaria parasite Plasmodium falciparum in West African populations of differing infection endemicity. Mol Biol Evol. 2014 Jun;31(6):1490–9. pmid:24644299
- 25. Mu J, Awadalla P, Duan J, McGee KM, Keebler J, Seydel K, et al. Genome-wide variation and identification of vaccine targets in the Plasmodium falciparum genome. Nat Genet. 2007 Jan;39(1):126–30. pmid:17159981
- 26. Ravenhall M, Benavente ED, Mipando M, Jensen ATR, Sutherland CJ, Roper C, et al. Characterizing the impact of sustained sulfadoxine/pyrimethamine use upon the Plasmodium falciparum population in Malawi. Malar J. 2016 Nov 29;15(1):575. pmid:27899115
- 27. Ocholla H, Preston MD, Mipando M, Jensen ATR, Campino S, MacInnis B, et al. Whole-Genome Scans Provide Evidence of Adaptive Evolution in Malawian Plasmodium falciparum Isolates. J Infect Dis. 2014 Dec 15;210(12):1991–2000. pmid:24948693
- 28. Hoffman SL, Oster CN, Plowe CV, Woollett GR, Beier JC, Chulay JD, et al. Naturally acquired antibodies to sporozoites do not prevent malaria: vaccine development implications. Science. 1987 Aug 7;237(4815):639–42. pmid:3299709
- 29. Buchwald AG, Sixpence A, Chimenya M, Damson M, Sorkin JD, Wilson ML, et al. Clinical Implications of Asymptomatic Plasmodium falciparum Infections in Malawi. Clin Infect Dis Off Publ Infect Dis Soc Am. 2019 01;68(1):106–12. pmid:29788054
- 30. Metzger WG, Okenu DMN, Cavanagh DR, Robinson JV, Bojang KA, Weiss HA, et al. Serum IgG3 to the Plasmodium falciparum merozoite surface protein 2 is strongly associated with a reduced prospective risk of malaria. Parasite Immunol. 2003 Jun;25(6):307–12. pmid:14507328
- 31. Egan AF, Morris J, Barnish G, Allen S, Greenwood BM, Kaslow DC, et al. Clinical immunity to Plasmodium falciparum malaria is associated with serum antibodies to the 19-kDa C-terminal fragment of the merozoite surface antigen, PfMSP-1. J Infect Dis. 1996 Mar;173(3):765–9. pmid:8627050
- 32. Perraut R, Varela M-L, Joos C, Diouf B, Sokhna C, Mbengue B, et al. Association of antibodies to Plasmodium falciparum merozoite surface protein-4 with protection against clinical malaria. Vaccine. 2017 04;35(48 Pt B):6720–6. pmid:29042203
- 33. Kang SY, Battle KE, Gibson HS, Cooper LV, Maxwell K, Kamya M, et al. Heterogeneous exposure and hotspots for malaria vectors at three study sites in Uganda. Gates Open Res [Internet]. 2018 Nov 13 [cited 2020 Mar 4];2. Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6350504/ pmid:30706054
- 34. Bejon P, Williams TN, Liljander A, Noor AM, Wambua J, Ogada E, et al. Stable and unstable malaria hotspots in longitudinal cohort studies in Kenya. PLoS Med. 2010 Jul 6;7(7):e1000304. pmid:20625549
- 35. Bousema T, Drakeley C, Gesase S, Hashim R, Magesa S, Mosha F, et al. Identification of Hot Spots of Malaria Transmission for Targeted Malaria Control. J Infect Dis. 2010 Jun 1;201(11):1764–74. pmid:20415536
- 36. Clark TD, Greenhouse B, Njama-Meya D, Nzarubara B, Maiteki-Sebuguzi C, Staedke SG, et al. Factors determining the heterogeneity of malaria incidence in children in Kampala, Uganda. J Infect Dis. 2008 Aug 1;198(3):393–400. pmid:18522503
- 37. Elissa N, Migot-Nabias F, Luty A, Renaut A, Touré F, Vaillant M, et al. Relationship between entomological inoculation rate, Plasmodium falciparum prevalence rate, and incidence of malaria attack in rural Gabon. Acta Trop. 2003 Mar;85(3):355–61. pmid:12659973
- 38. Gaudart J, Poudiougou B, Dicko A, Ranque S, Toure O, Sagara I, et al. Space-time clustering of childhood malaria at the household level: a dynamic cohort in a Mali village. BMC Public Health. 2006 Nov 21;6:286. pmid:17118176
- 39. Zhu SJ, Almagro-Garcia J, McVean G. Deconvolution of multiple infections in Plasmodium falciparum from high throughput sequencing data. Bioinforma Oxf Engl. 2017 Aug 22;
- 40. Bhatia G, Patterson N, Sankararaman S, Price AL. Estimating and interpreting FST: The impact of rare variants. Genome Res. 2013 Sep 1;23(9):1514–21. pmid:23861382
- 41. Vanegas M, Bermúdez A, Guerrero YA, Cortes-Vecino JA, Curtidor H, Patarroyo ME, et al. Protecting capacity against malaria of chemically defined tetramer forms based on the Plasmodium falciparum apical sushi protein as potential vaccine components. Biochem Biophys Res Commun. 2014 Aug 15;451(1):15–23. pmid:25063026
- 42. Gupta A, Thiruvengadam G, Desai SA. The conserved clag multigene family of malaria parasites: essential roles in host-pathogen interaction. Drug Resist Updat Rev Comment Antimicrob Anticancer Chemother. 2015 Jan;18:47–54. pmid:25467627
- 43. Iriko H, Kaneko O, Otsuki H, Tsuboi T, Su X, Tanabe K, et al. Diversity and evolution of the rhoph1/clag multigene family of Plasmodium falciparum. Mol Biochem Parasitol. 2008 Mar 1;158(1):11–21. pmid:18155305
- 44. van Schaijk BCL, Ploemen IHJ, Annoura T, Vos MW, Foquet L, van Gemert G-J, et al. A genetically attenuated malaria vaccine candidate based on P. falciparum b9/slarp gene-deficient sporozoites. eLife. 2014 Nov 19;3.
- 45. Krzyczmonik K, Świtnicki M, Kaczanowski S. Analysis of immunogenicity of different protein groups from malaria parasite Plasmodium falciparum. Infect Genet Evol. 2012 Dec 1;12(8):1911–6. pmid:22986003
- 46. Toenhake CG, Fraschka SA-K, Vijayabaskar MS, Westhead DR, van Heeringen SJ, Bártfai R. Chromatin Accessibility-Based Characterization of the Gene Regulatory Network Underlying Plasmodium falciparum Blood-Stage Development. Cell Host Microbe. 2018 Apr 11;23(4):557–569.e9. pmid:29649445
- 47. Zhang M, Wang C, Otto TD, Oberstaller J, Liao X, Adapa SR, et al. Uncovering the essential genes of the human malaria parasite Plasmodium falciparum by saturation mutagenesis. Science [Internet]. 2018 May 4 [cited 2020 Feb 18];360(6388). Available from: https://science.sciencemag.org/content/360/6388/eaap7847 pmid:29724925
- 48.
PlasmoDB: a functional genomic database for malaria parasites.—PubMed—NCBI [Internet]. [cited 2020 Feb 18]. Available from: https://www.ncbi.nlm.nih.gov/pubmed?cmd=search&term=18957442
- 49. Quintana MDP, Ch’ng J-H, Zandian A, Imam M, Hultenby K, Theisen M, et al. SURGE complex of Plasmodium falciparum in the rhoptry-neck (SURFIN4.2-RON4-GLURP) contributes to merozoite invasion. PloS One. 2018;13(8):e0201669. pmid:30092030
- 50. Proellocks NI, Herrmann S, Buckingham DW, Hanssen E, Hodges EK, Elsworth B, et al. A lysine-rich membrane-associated PHISTb protein involved in alteration of the cytoadhesive properties of Plasmodium falciparum-infected red blood cells. FASEB J Off Publ Fed Am Soc Exp Biol. 2014 Jul;28(7):3103–13. pmid:24706359
- 51. Holt DC, Gardiner DL, Thomas EA, Mayo M, Bourke PF, Sutherland CJ, et al. The cytoadherence linked asexual gene family of Plasmodium falciparum: are there roles other than cytoadherence? Int J Parasitol. 1999 Jun;29(6):939–44. pmid:10480731
- 52.
Pf3k pilot data release 5 | MalariaGEN [Internet]. [cited 2019 Dec 23]. Available from: https://www.malariagen.net/data/pf3K-5
- 53. Tessema SK, Nakajima R, Jasinskas A, Monk SL, Lekieffre L, Lin E, et al. Protective Immunity against Severe Malaria in Children Is Associated with a Limited Repertoire of Antibodies to Conserved PfEMP1 Variants. Cell Host Microbe. 2019 Nov 13;26(5):579–590.e5. pmid:31726028
- 54. Valencia-Hernandez AM, Ng WY, Ghazanfari N, Ghilas S, Menezes MN de, Holz LE, et al. A Natural Peptide Antigen within the Plasmodium Ribosomal Protein RPL6 Confers Liver TRM Cell-Mediated Immunity against Malaria in Mice. Cell Host Microbe. 2020 Jun 10;27(6):950–962.e7. pmid:32396839
- 55. MalariaGEN Plasmodium falciparum Community Project. Genomic epidemiology of artemisinin resistant malaria. Neher RA, editor. eLife. 2016 Mar 4;5:e08714. pmid:26943619
- 56. Färnert A, Rooth I, Svensson Å, Snounou G, Björkman A. Complexity of Plasmodium falciparum Infections Is Consistent over Time and Protects against Clinical Disease in Tanzanian Children. J Infect Dis. 1999 Apr 1;179(4):989–95. pmid:10068596
- 57. Bereczky S, Liljander A, Rooth I, Faraja L, Granath F, Montgomery SM, et al. Multiclonal asymptomatic Plasmodium falciparum infections predict a reduced risk of malaria disease in a Tanzanian population. Microbes Infect. 2007 Jan 1;9(1):103–10. pmid:17194613
- 58. Sondén K, Doumbo S, Hammar U, Vafa Homann M, Ongoiba A, Traoré B, et al. Asymptomatic Multiclonal Plasmodium falciparum Infections Carried Through the Dry Season Predict Protection Against Subsequent Clinical Malaria. J Infect Dis. 2015 Aug 15;212(4):608–16. pmid:25712968
- 59. Sexton AE, Doerig C, Creek DJ, Carvalho TG. Post-Genomic Approaches to Understanding Malaria Parasite Biology: Linking Genes to Biological Functions. ACS Infect Dis. 2019 09;5(8):1269–78. pmid:31243988
- 60. Sam-Yellowe TY, Perkins ME. Interaction of the 140/130/110 kDa rhoptry protein complex of Plasmodium falciparum with the erythrocyte membrane and liposomes. Exp Parasitol. 1991 Aug 1;73(2):161–71. pmid:1889471
- 61. Kaneko O, Lim BYSY, Iriko H, Ling IT, Otsuki H, Grainger M, et al. Apical expression of three RhopH1/Clag proteins as components of the Plasmodium falciparum RhopH complex. Mol Biochem Parasitol. 2005 Sep 1;143(1):20–8. pmid:15953647
- 62. Guy AJ, Irani V, MacRaild CA, Anders RF, Norton RS, Beeson JG, et al. Insights into the Immunological Properties of Intrinsically Disordered Malaria Proteins Using Proteome Scale Predictions. PLOS ONE. 2015 Oct 29;10(10):e0141729. pmid:26513658
- 63. Kinyanjui SM, Conway DJ, Lanar DE, Marsh K. IgG antibody responses to Plasmodium falciparum merozoite antigens in Kenyan children have a short half-life. Malar J. 2007 Jun 28;6:82. pmid:17598897
- 64. Akpogheneta OJ, Duah NO, Tetteh KKA, Dunyo S, Lanar DE, Pinder M, et al. Duration of Naturally Acquired Antibody Responses to Blood-Stage Plasmodium falciparum Is Age Dependent and Antigen Specific. Infect Immun. 2008 Apr 1;76(4):1748–55. pmid:18212081
- 65. Shah Z, Adams M, Moser KA, Shrestha B, Stucke EM, Laufer MK, et al. Optimization of parasite DNA enrichment approaches to generate whole genome sequencing data for Plasmodium falciparum from low parasitaemia samples. Malar J. 2020 Mar 30;19(1):135. pmid:32228559
- 66. 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. Malar J [Internet]. 2016 Dec 20 [cited 2018 Mar 28];15. Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5175302/ pmid:27998271
- 67. Roy SW, Ferreira MU, Hartl DL. Evolution of allelic dimorphism in malarial surface antigens. Heredity. 2008 Feb;100(2):103–10. pmid:17021615
- 68. Zainabadi K, Adams M, Han ZY, Lwin HW, Han KT, Ouattara A, et al. A novel method for extracting nucleic acids from dried blood spots for ultrasensitive detection of low-density Plasmodium falciparum and Plasmodium vivax infections. Malar J. 2017 Sep 18;16(1):377. pmid:28923054
- 69. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012 Apr;9(4):357–9.
- 70. DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011 May;43(5):491–8. pmid:21478889
- 71. Van der Auwera GA, Carneiro MO, Hartl C, Poplin R, Del Angel G, Levy-Moonshine A, et al. From FastQ data to high confidence variant calls: the Genome Analysis Toolkit best practices pipeline. Curr Protoc Bioinforma. 2013;43:11.10.1–33. pmid:25431634
- 72. Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010 Mar 15;26(6):841–2. pmid:20110278
- 73. Manske M, Miotto O, Campino S, Auburn S, Almagro-Garcia J, Maslen G, et al. Analysis of Plasmodium falciparum diversity in natural infections by deep sequencing. Nature. 2012 Jul;487(7407):375–9. pmid:22722859
- 74. Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, et al. The variant call format and VCFtools. Bioinforma Oxf Engl. 2011 Aug 1;27(15):2156–8. pmid:21653522
- 75.
Clement M, Snell Q, Walke P, Posada D, Crandall K. TCS: estimating gene genealogies. In: Proceedings 16th International Parallel and Distributed Processing Symposium [Internet]. Ft. Lauderdale, FL: IEEE; 2002 [cited 2020 Feb 4]. p. 7 pp. Available from: http://ieeexplore.ieee.org/document/1016585/
- 76.
popart: full-feature software for haplotype network construction—Leigh—2015—Methods in Ecology and Evolution—Wiley Online Library [Internet]. [cited 2020 Feb 4]. Available from: https://besjournals.onlinelibrary.wiley.com/doi/full/10.1111/2041-210X.12410
- 77. Rozas J, Ferrer-Mata A, Sánchez-DelBarrio JC, Guirao-Rico S, Librado P, Ramos-Onsins SE, et al. DnaSP 6: DNA Sequence Polymorphism Analysis of Large Data Sets. Mol Biol Evol. 2017 Dec 1;34(12):3299–302. pmid:29029172