Multiple genetic loci define Ca++ utilization by bloodstream malaria parasites
BMC Genomics volume 20, Article number: 47 (2019)
Bloodstream malaria parasites require Ca++ for their development, but the sites and mechanisms of Ca++ utilization are not well understood. We hypothesized that there may be differences in Ca++ uptake or utilization by genetically distinct lines of P. falciparum. These differences, if identified, may provide insights into molecular mechanisms.
Dose response studies with the Ca++ chelator EGTA (ethylene glycol-bis(β-aminoethyl ether)-N,N,N′,N′-tetraacetic acid) revealed stable differences in Ca++ requirement for six geographically divergent parasite lines used in previous genetic crosses, with the largest difference seen between the parents of the HB3 x Dd2 cross. Genetic mapping of Ca++ requirement yielded complex inheritance in 34 progeny clones with a single significant locus on chromosome 7 and possible contributions from other loci. Although encoded by a gene in the significant locus and a proposed Ca++ target, PfCRT (P. falciparum chloroquine resistance transporter), the primary determinant of clinical resistance to the antimalarial drug chloroquine, does not appear to contribute to this quantitative trait. Stage-specific application of extracellular EGTA also excluded determinants associated with merozoite egress and erythrocyte reinvasion.
We have identified differences in Ca++ utilization amongst P. falciparum lines. These differences are under genetic regulation, segregating as a complex trait in genetic cross progeny. Ca++ uptake and utilization throughout the bloodstream asexual cycle of malaria parasites represents an unexplored target for therapeutic intervention.
Development of Plasmodium parasites within erythrocytes accounts for most of the clinical sequelae of malaria, which remains a leading infectious cause of morbidity and mortality worldwide. In the virulent human pathogen, P. falciparum, intraerythrocytic growth protects the parasite from host immune responses, but also limits the pathogen from access to nutrients and critical ions in host plasma.
Ca++ is one such critical ion because it is present at relatively high concentrations in plasma (~ 1.1 mM ionized and 1.4 mM bound to proteins ), but available at much lower levels within the erythrocyte (< 100 nM, [2, 3]). Multiple studies have determined that the intracellular parasite requires Ca++ [4,5,6] and that this pathogen may fulfil this need by increasing Ca++ permeability at the host membrane [7,8,9,10,11]. Interestingly, the plasmodial surface anion channel (PSAC), a broad selectivity nutrient and ion channel induced on the host erythrocyte membrane [12, 13], is not responsible for the increased Ca++ permeability . In addition to roles in parasite development, Ca++ is also needed for egress from the host cell as determined through selective loading of membrane-permeable Ca++ chelators . A separate requirement for the subsequent invasion of erythrocytes has also been identified [4, 15, 16].
Despite evidence for multiple distinct roles of Ca++ in the bloodstream malaria parasites, the molecular targets and mechanisms of calcium acquisition remain poorly characterized. Computational analyses of the parasite genome database have identified multiple genes encoding Ca++-binding proteins via either the EF-hand or C2 domains . A single calmodulin gene, two calcineurin subunits, several centrins, and Ca++ dependent protein kinases have been identified as EF-hand domain containing proteins in P. falciparum; several C2 domain containing proteins have also been identified . While these genes have been characterized to varying extents, few have been studies through biochemical studies to confirm regulated response to Ca++ binding. As these domains were identified in higher organisms [19, 20], it is possible that Plasmodium spp. use additional, uncharacterized protein motifs and Ca++ regulatory mechanisms. Thus, the molecular basis of Ca++ requirement and the roles of Ca++ signaling in malaria parasite biology remain unclear.
To explore the global effects of Ca++ and overcome the limitations of computational approaches, we reasoned that quantitative trait locus (QTL) mapping in parasite genetic crosses may provide an unbiased method of identifying the major Ca++ targets. This method has successfully identified parasite molecules involved in drug resistance, invasion, and nutrient uptake [21,22,23,24]. Here, we identified differences in Ca++ requirement for laboratory lines of P. falciparum. We then tracked inheritance of a significant difference in the Dd2 x HB3 genetic cross and used QTL mapping to identify possible Ca++ targets. Our studies implicate complex roles for Ca++ throughout the parasite cycle and suggest a genomic locus with a significant effect.
P. falciparum lines exhibit differing susceptibilities to the Ca++ chelator EGTA
Previous studies used addition of EGTA to standard culture medium to determine that P. falciparum lines require extracellular Ca++ for propagation [4,5,6]; EGTA toxicity from mechanisms other than chelation of divalent cations was excluded because equimolar addition of CaCl2 fully restores parasite growth. Here, we hypothesized that there may be reproducible differences in Ca++ requirement for parasite clones; such differences could be used to explore possible molecular targets and the precise roles served by Ca++. We therefore surveyed several laboratory clones and measured parasite growth inhibition by a range of EGTA concentrations (Fig. 1a). Under our experimental conditions, 0.45 mM EGTA effectively abolished expansion of cultures for each of the examined parasite lines. This concentration is consistent with stochiometric chelation of Ca++, which is present at a nominal 0.42 mM concentration in standard RPMI 1640 medium; addition of lipid-rich bovine albumin preparations, as required for parasite cultivation, may affect the free Ca++ available for EGTA chelation. Lower EGTA concentrations yielded survival and growth that differed significantly depending on parasite genotype (P < 0.01 for one-way ANOVA comparisons at 0.33, 0.36, and 0.39 mM EGTA, n = 3 to 7 trials each). At concentrations that produced incomplete growth inhibition, we observed that Pf803 and Dd2 were consistently the least and most resistant to EGTA addition, respectively.
While these growth inhibition experiments utilized relatively small changes in EGTA concentration, the free Ca++ concentration in the culture medium is expected to change dramatically and nonlinearly with incremental addition of chelator . We therefore used a Ca++-sensitive electrode to estimate free Ca++ in culture media supplemented with EGTA. At concentrations up to 0.4 mM, EGTA addition reduced free Ca++ stoichiometrically; higher concentrations produced logarithmic reductions in free Ca++, but precise estimation was limited by electrode sensitivity (Fig. 1b and c). Because the Ca++ affinity of EGTA depends on temperature, pH and other factors, we performed these electrode measurements at 37 °C and attempted to simulate the conditions encountered in parasite culture. Our estimates may nevertheless be adversely affected by changes in pH under culture conditions due to metabolic acids produced by parasites and the presence of Ca++-chelating agents such as intracellular proteins and incompletely characterized metabolic byproducts. With these caveats in mind, the IC50 values for the parasite lines we examined, between 0.34 and 0.40 mM EGTA, correspond to half-maximal requirement (EC50) between 11 and 40 μM free extracellular Ca++ for the parasite lines we examined.
Each of the parasite clones in Fig. 1a has been used to generate one or more P. falciparum genetic crosses . Thus, the differences in EGTA susceptibility in Fig. 1 could be examined in cross progeny to identify genetic loci involved in Ca++ requirement. Amongst the four available genetic crosses, the Dd2 x HB3 pair of parental lines exhibited the greatest difference in EGTA susceptibility and was, therefore, selected for further study.
Inheritance in the Dd2 x HB3 cross
To search for genetic loci that may determine in vitro susceptibility, we performed EGTA growth inhibition studies with each of the 34 available independent progeny clones from the Dd2 x HB3 genetic cross. Each clone was examined in dose response experiments using 0.36, 0.39, and 0.40 mM EGTA and evaluated in 7–17 independent trials. For each trial, growth was normalized to 100% for matched controls grown without EGTA to correct for known differences in propagation rates for these lines . Throughout these studies, we continued to evaluate EGTA sensitivities of Dd2 and HB3 in parallel with the progeny to evaluate stability of the parental phenotypes on continued cultivation and as controls for measurements using progeny clones. These studies confirmed stable differences between the parental lines (P < 10− 3 at 0.36, 0.39, and 0.4 mM EGTA, n = 49–53 independent trials) and revealed a range of EGTA sensitivities for the progeny clones (Fig. 2a), consistent with one or more genetic loci that determine in vitro Ca++ requirement. Many progeny clones were more sensitive to EGTA than either parent; some exhibited intermediate values. Because the progeny were ordered in Fig. 2a according to increasing tolerance to 0.36 mM EGTA (black bars), this plot also reveals that the relative sensitivities to 0.39 and 0.40 mM EGTA could not be confidently predicted by measurements at a single, lower concentration. For example, the QC101 and 3BA6 daughters expanded equally well with 0.36 mM EGTA, but 3BA6 tolerated the two higher EGTA concentrations better than QC101.
These observed differences in dose response steepness have several explanations . Here, a conservative hypothesis is that Ca++ serves multiple roles over the course of the parasite cycle such as facilitating intracellular development, parasite egress and re-invasion. Each of these target sites may have a distinct dissociation constant, Kd, for Ca++ binding. The overall shape of the EGTA dose response is then determined by summing the growth inhibitory effects of Ca++ removal at each target site based on the corresponding Kd values. Distinct parasite genomic loci presumably define these target sites and their Kd values. Thus, the observed range of EGTA susceptibilities and dose response profiles in the HB3 x Dd2 cross progeny are consistent with non-Mendelian or complex inheritance.
We then used QTL analysis to search for inherited genetic loci associated with growth at each EGTA concentration. We sought correlations between inheritance of mean growth at each EGTA concentration and each of 448 microsatellite markers from the 14 parasite chromosomes. This statistical approach identified a single locus that reached weak statistical significance at 0.4 mM EGTA (Fig. 2b, P = 0.05 for the 3E7 microsatellite marker on chromosome 7 at 0.4 mM EGTA). It also found a few minor peaks that did not reach statistical significance, including the 3E7 marker when growth was evaluated at 0.39 mM EGTA.
No obvious candidate genes or association with haplotypes of PfCRT, the chloroquine resistance transporter
The chromosome 7 locus associated with the positive 3E7 microsatellite marker is flanked by two negative markers, HSP86 and 7A11. Because the polymorphism(s) that mediate differences in EGTA sensitivity could be situated anywhere between these negative markers, we examined this entire 45 kB region and found it contains 12 protein-coding genes and two putative small nucleolar RNAs. None of these candidates has conventional Ca++-binding domains or annotated roles that could readily account for changes in EGTA sensitivity.
Interestingly, this chromosome locus has previously been implicated in inheritance of resistance to chloroquine, once the mainstay of antimalarial chemotherapy. Molecular studies have implicated mutations in PfCRT, a polytopic membrane protein at the parasite’s digestive vacuole . While it is clear that these mutations confer chloroquine resistance, whether this protein is a transporter and its physiological role remain debated; in addition to direct efflux of chloroquine via PfCRT, hypotheses include roles in transport of H+ , amino acids , iron , glutathione and other organic solutes [33, 34], Cl− , and Ca++ . Our observation that tolerance to EGTA may also be partially determined by this locus suggests a role in Ca++ transport for PfCRT or other gene products in the mapped locus.
To test this possibility, we performed EGTA growth inhibition dose response experiments with the C1GC03, C2GC03, C4Dd2, and C67G8 transfectant lines . These lines were generated in the GC03 chloroquine sensitive progeny clone from the Dd2 x HB3 cross and designed to confer resistant PfCRT haplotypes (C4Dd2 and C67G8) or serve as transfection controls (C1GC03 and C2GC03) for possible effects of transfection or the associated hDHFR and BSD selectable markers (human dihydrofolate reductase and blasticidin S deaminase, respectively). EGTA dose response studies with these lines revealed no significant differences between these four transfectants (Fig. 3); there were also insignificant effects of transfection or selectable markers (P > 0.1 at each EGTA concentration for one-way ANOVA comparisons amongst the four transfectant lines and in comparisons of the two transfection controls with wild-type GC03, n = 10 or 11 trials for each transfectant line). These experiments suggest that PfCRT haplotype does not influence Ca++ utilization in the intracellular parasite. As other genomic loci may interact with this locus and these transfections may not yeild the maximally effective haplotypes, additional studies including transfections that do not use multiple selectable markers may be required to determine whether PfCRT or other genes in this locus contribute to parasite Ca++ requirement.
How might this locus determine variations in Ca++ requirement in the HB3 x Dd2 genetic cross? Although the absence of Ca++-binding domains in any of the 12 protein-coding genes in this locus may be discouraging, we recognize that Plasmodium spp. may use proteins without clear orthologs in other systems to serve critical roles. A good example is the clag gene family , which is conserved in all examined Plasmodium spp. but absent from other genera. Although the encoded CLAG proteins do not resemble ion channel proteins, genetic mapping and DNA transfection implicated a primary role for this gene family in the plasmodial surface anion channel (PSAC), a nutrient uptake channel at the host erythrocyte membrane . As used in that study, an agnostic approach to identifying the gene(s) that govern Ca++ requirement would involve allelic exchange DNA transfection for each of the 12 candidates in the 3E7 marker-associated locus; EGTA dose response studies as we have performed here could then uncover the responsible gene(s). Ongoing improvements in parasite DNA transfection technology may make this a less effort- and time-consuming approach.
Ca++ requirements during merozoite egress and invasion
Two of the most critical events that require Ca++ are the egress of daughter merozoites at the end of the intracellular cycle and the subsequent reinvasion of new erythrocytes to form ring-stage parasites. Egress is dependent on intracellular Ca++ and does not appear to require extracellular Ca++ [14, 38]; in contrast, reinvasion is thought to require primarily extracellular Ca++ . We therefore examined whether Ca++ requirements during egress and reinvasion differ between Dd2 and HB3 parasites. Experiments using a short ~ 5 h application of EGTA to tightly synchronized late-stage schizonts followed by microscopic examination and quantification of ring-stage parasitemias with flow cytometry revealed dose-dependent reduction in successful progression to rings for both Dd2 and HB3 (Fig. 4), with inhibition plateaus at 65–70% of those seen with the positive control cytochalasin D. The EGTA dose responses for Dd2 and HB3 were indistinguishable, with IC50 values of 0.57 mM ± 0.04 for both lines (P 0.31, n = 3 independent trials); these values were higher that observed in our standard growth inhibition experiments, as might have been predicted given the shorter duration of EGTA application. Thus, differing Ca++ requirements for egress or reinvasion do not account for the difference in EGTA sensitivities of these parasite lines.
Calcium serves many roles in cells, including control of cell cycle regulation, DNA replication, modulation of enzymes, and signaling . Over the course of their asexual bloodstream cycle, malaria parasites undergo each of these activities. Nevertheless, the roles served by Ca++ in these and other cellular processes unique to these pathogens such as host cell egress and invasion or ion homeostasis within specialized organelles have been difficult to study for several reasons. The multiple compartments within infected cells make precise quantification of Ca++ at target sites difficult. Manipulation of free Ca++ concentrations at individual target sites is also currently not feasible. The low Ca++ concentration within the host cell, maintained by an Ca++ ATPase extrusion pump at the erythrocyte surface, also complicates basic studies aimed at mechanistic insights. Finally, transient changes in Ca++ concentrations, as observed and tracked in other cells [40, 41], are more difficult to detect within infected erythrocytes because of their small size and quenching by host cell hemoglobin. While Ca++ transients and spatial gradients have been well-studied in larger muscle cells and neurons [42, 43], the small compartments within infected cells require measurements near or beyond the limit of spatial resolution in conventional confocal microscopy.
To overcome these limitations, our study used the Ca++ chelator EGTA to determine that P. falciparum lines exhibit stable differences in extracellular Ca++ requirement for in vitro cultivation. We then tracked this difference in the progeny clones of the Dd2 x HB3 genetic cross; this approach revealed stable, heritable differences in parasite Ca++ requirement that may be defined by multiple genetic loci. A key limitation we encountered was the relatively small range of EGTA concentrations over which differences between Dd2 and HB3 are observed. This narrow range is determined primarily by exponential decreases in free Ca++ concentrations resulting from further addition of EGTA once the total Ca++ and EGTA concentrations are approximately equal . The steep reduction in parasite survival with small increases in EGTA may also reflect synergistic killing due to action on multiple parasite Ca++-dependent targets. While these differences in EGTA tolerance were modest, their statistical significance provides a framework that may eventually lead to identification of key molecular determinants. For example, a modest difference between parental lines for PSAC inhibition by furosemide in the 3D7xHB3 cross failed to detect significant QTL . Although unsuccessful in finding the genes, that study stimulated high-throughput screens to find a Dd2-specific inhibitor and enabled successful genetic mapping to implicate clag3 genes in PSAC formation [12, 23].
Complex inheritance is consistent with multiple targets for Ca++, which may act at various developmental points in the 72 h growth inhibition assays we used. The observed clone-dependent differences in EGTA susceptibility presumably reflect polymorphisms in key target proteins. Polymorphisms in transporters or Ca++-dependent enzymes may directly lead to altered Ca++ affinities or may indirectly affect growth if the polymorphisms affect enzyme kinetics. Heritable differences could also reflect polymorphisms in proteins that do not bind Ca++; one such possibility is a Ca++-independent transcription factor activated by Ca++-dependent phosphatases to then control cell cycle progression . The stable differences in progeny clone phenotypes shown in Fig. 2 could also reflect stable epigenetic marks under in vitro culture conditions rather than heritable genome-level differences, as exemplified by epigenetic regulation of var and clag3 genes in malaria [46, 47].
The potential locus we identified on chromosome 7 is particularly interesting as it encodes PfCRT, a putative transporter implicated in clinical resistance to chloroquine and other antimalarial drugs . Notably, several studies have proposed Ca++ affects chloroquine action and digestive vacuolar physiology, based primarily on the effects of Ca++ channel inhibitors and calmodulin antagonists [48, 49]. Our growth inhibition studies with transfectant lines carrying modified pfcrt alleles did not alter EGTA sensitivity. Although this result does not itself support a link between Ca++ availability and PfCRT activity, it is possible that additional transfections using different genetic backgrounds and other polymorphisms may confirm the proposed roles of Ca++ in chloroquine susceptibility. Alternatively, one or more of the other genes within the mapped locus may encode Ca++-dependent activities. As linkage analysis produced only a marginally significant QTL, it will be important to have greater confidence in this locus, possibly though identification and characterization of a larger number of progeny clones. Studies that examine how Ca++ concentration affects digestive vacuolar physiology may also be complicated by Ca++ indicator dye import via PfMDR1 (P. falciparum multidrug resistance gene product), another transporter on this vacuole .
Our studies have identified differences in Ca++ utilization by P. falciparum parasite lines. These differences may reflect variable requirement amongst parasite clones in human infections to confer survival advantage under stress conditions, such as hypocalcemia associated with malnutrition . Alternatively, they may have arisen through selective pressures during adaptation to in vitro culture. In either case, the observed differences and unbiased studies such as these linkage analyses have the potential to uncover Ca++-dependent processes as targets for therapy development.
Plasmodium falciparum laboratory lines were cultivated separately in RPMI 1640 with L-glutamine (Thermo Fisher Scientific, Waltham, MA) supplemented with 50 mg/L hypoxanthine, 31 mM NaHCO3, 25 mM HEPES, 10 μg/mL gentamicin, and 0.5% NZ microbiological BSA (MP Biomedicals, Santa Ana, CA) using O+ human erythrocytes obtained from anonymous donors (VA Blood Services, Richmond, VA). Cultures were maintained at 5% hematocrit under 5% O2, 5% CO2, 90% N2 at 37 °C. Parasite genotypes were confirmed with molecular studies.
Growth inhibition assays
Parasite Ca++ requirement for in vitro propagation was evaluated by addition of EGTA to culture media; experiments were identically performed for all laboratory parasite lines and progeny clones. A standard 72 h growth assay and the SYBR Green I dye were used to measure parasite nucleic acid production, as described previously . Parasite cultures were synchronized with 5% sorbitol prior to each trial. These ring-stage cultures were seeded in 96-well microplates at 2% hematocrit and 0.5–1.0% parasitemia in the above culture medium with indicated EGTA concentrations. pH measurements were used to ensure that EGTA addition did not change the pH of the medium. Microplate cultures were then incubated at 37 °C under 5% O2, 5% CO2, 90% N2 without medium changes. Expansion of cultures during this incubation was evaluated by addition of lysis buffer (20 mM Tris, 10 mM EDTA, 0.016% saponin, 1.6% triton X100, pH 7.5) with SYBR Green I nucleic acid gel stain at a 5000x dilution (Invitrogen, Carlsbad CA). After a 45 min incubation without ambient light, fluorescence measurements (excitation/emission wavelengths of 485 nm and 528 nm) were used to quantify parasite DNA. Growth was calculated as the mean of triplicate wells after subtracting background fluorescence from matched cultures seeded with 20 μM chloroquine to abolish parasite expansion. These values were then normalized to 100% for matched cultures grown without chelator or inhibitor. Control experiments with stoichiometric addition of CaCl2 to EGTA confirmed that growth is fully restored, excluding nonspecific EGTA toxicity.
Free [ca++] measurements
Free Ca++ concentrations in culture media with and without EGTA addition were measured with an ion sensitive electrode (Cole-Parmer, Vernon Hills, IL). Media were gassed and measured at 37 °C to match pH and temperature encountered by parasites under in vitro culture. Freshly prepared Ca++ standards were also measured at 37 °C and used to estimate an electrode slope of 26.7 mV/decade, indicating electrode specificity for Ca++. Media and standards were measured after addition of 4 M KCl as an ionic strength adjuster, according to the electrode manufacturer.
QTL analysis to search for parasite genomic loci involved in Ca++ utilization was performed as described . We used 448 previously selected microsatellite markers that distinguish the Dd2 and HB3 parasite lines and define the available progeny clones from this genetic cross ; the full list of markers and their sequences are available at https://www.ncbi.nlm.nih.gov/probe/?term=10558988[pmid]. QTL analysis was performed by the multiple imputation method  as implemented in the R/qtl software (freely available at http://www.rqtl.org/) . A P = 0.05 significance threshold, shown as a dashed line in each QTL genome scan, was determined through analysis of 1000 permutations.
Merozoite egress and reinvasion assays
EGTA dose response studies were used to examine the Ca++ requirement during egress and reinvasion at the end of the asexual parasite lifecycle using modifications to a previously described method . Dd2 and HB3 parasite lines were tightly synchronized with Percoll/sorbitol enrichment of schizont-stage cultures, cultivation with fresh erythrocytes for 6 h, and a 5% sorbitol treatment. Subsequently, serial Giemsa-stained smears were performed until a majority of mature schizonts were observed, typically after 38–40 h. Synchronized schizonts were again enriched with Percoll/sorbitol, added to fresh erythrocytes, and resuspended in 96-well microplates at a 1% hematocrit and 2% parasitemia with complete medium and indicated EGTA concentrations. These parasites were cultivated at 37 °C with 5% CO2 and regular monitoring until the ring-stage parasitemia reached 2%, approximately 5 h. 20 μL of each suspension culture was then transferred to a fresh microplate with 80 μL of 1X SYBR Green I in 1x PBS, incubated for 30 min at room temperature, and washed twice in 1x PBS. The stained cells were resuspended in 100 μL of 1x PBS and the number of ring-infected erythrocytes counted using an Accuri C6 flow cytometer (BD Biosciences, Franklin Lakes, NJ). The data were normalized to percent inhibition of invasion using two matched controls (no inhibitor and 5 μM cytochalasin D) present in each microplate. IC50 values were calculated in SigmaPlot 10 (Systat Software, San Jose, CA).
Statistical significance was evaluated using Student’s t tests or one-way analysis of variance (ANOVA) as indicated, using SigmaPlot 10 or Prism 7 software.
- BSD :
Blasticidin S deaminase gene
- clag and CLAG:
Cytoadherence linked antigen gene and protein, respectively
Ethylene glycol-bis(β-aminoethyl ether)-N,N,N′,N′-tetraacetic acid)
- hDHFR :
Human dihydrofolate reductase gene
- IC 50 :
Half-maximal inhibitory concentration
- K d :
P. falciparum chloroquine resistance transporter
P. falciparum multidrug resistance gene product
Plasmodial surface anion channel
Quantitative trait locus
Goldstein DA, Calcium S. In: Walker HK, Hall WD, Hurst JW, editors. Clinical methods: the history, physical, and laboratory examinations. Boston: Butterworths; 1990. p. 677–9.
Murphy E, Levy L, Berkowitz LR, Orringer EP, Gabel SA, London RE. Nuclear magnetic resonance measurement of cytosolic free calcium levels in human red blood cells. Am J Phys. 1986;251(4 Pt 1):C496–504.
Simonsen LO, Gomme J, Lew VL. Uniform ionophore A23187 distribution and cytoplasmic calcium buffering in intact human red cells. Biochim Biophys Acta. 1982;692(3):431–40.
Wasserman M, Alarcon C, Mendoza PM. Effects of ca++ depletion on the asexual cell cycle of Plasmodium falciparum. Am J Trop Med Hyg. 1982;31(4):711–7.
McCallum-Deighton N, Holder AA. The role of calcium in the invasion of human erythrocytes by Plasmodium falciparum. Mol Biochem Parasitol. 1992;50(2):317–23.
Zipprer EM, Neggers M, Kushwaha A, Rayavara K, Desai SA. A kinetic fluorescence assay reveals unusual features of ca++ uptake in Plasmodium falciparum-infected erythrocytes. Malar J. 2014;13:184.
Tanabe K, Mikkelsen RB, Wallach DF. Calcium transport of Plasmodium chabaudi-infected erythrocytes. J Cell Biol. 1982;93(3):680–4.
Krungkrai J, Yuthavong Y. Enhanced Ca2+ uptake by mouse erythrocytes in malarial (Plasmodium berghei) infection. Mol Biochem Parasitol. 1983;7(3):227–35.
Desai SA, McCleskey EW, Schlesinger PH, Krogstad DJ. A novel pathway for ca++ entry into Plasmodium falciparum-infected blood cells. Am J Trop Med Hyg. 1996;54(5):464–70.
Tiffert T, Staines HM, Ellory JC, Lew VL. Functional state of the plasma membrane Ca2+ pump in Plasmodium falciparum-infected human red blood cells. J Physiol. 2000;525(Pt 1):125–34.
Kushwaha AK, Apolis L, Ito D, Desai SA. Increased ca++ uptake by erythrocytes infected with malaria parasites: evidence for exported proteins and novel inhibitors. Cell Microbiol. 2018;e12853.
Nguitragool W, Bokhari AA, Pillai AD, Rayavara K, Sharma P, Turpin B, Aravind L, Desai SA. Malaria parasite clag3 genes determine channel-mediated nutrient uptake by infected red blood cells. Cell. 2011;145(5):665–77.
Ito D, Schureck MA, Desai SA. An essential dual-function complex mediates erythrocyte invasion and channel-mediated nutrient uptake in malaria parasites. Elife. 2017;6:e23485.
Glushakova S, Lizunov V, Blank PS, Melikov K, Humphrey G, Zimmerberg J. Cytoplasmic free Ca2+ is essential for multiple steps in malaria parasite egress from infected erythrocytes. Malar J. 2013;12:41.
Johnson JG, Epstein N, Shiroishi T, Miller LH. Factors affecting the ability of isolated Plasmodium knowlesi merozoites to attach to and invade erythrocytes. Parasitology. 1980;80(3):539–50.
Gao X, Gunalan K, Yap SS, Preiser PR. Triggers of key calcium signals during erythrocyte invasion by Plasmodium falciparum. Nat Commun. 2013;4:2862.
Lourido S, Moreno SN. The calcium signaling toolkit of the apicomplexan parasites Toxoplasma gondii and Plasmodium spp. Cell Calcium. 2015;57(3):186–93.
Brochet M, Billker O. Calcium signalling in malaria parasites. Mol Microbiol. 2016;100(3):397–408.
Ikura M. Calcium binding and conformational response in EF-hand proteins. Trends Biochem Sci. 1996;21(1):14–7.
Rizo J, Sudhof TC. C2-domains, structure and function of a universal Ca2+-binding domain. J Biol Chem. 1998;273(26):15879–82.
Wellems TE, Walker-Jonah A, Panton LJ. Genetic mapping of the chloroquine-resistance locus on Plasmodium falciparum chromosome 7. Proc Natl Acad Sci U S A. 1991;88(8):3382–6.
Hayton K, Gaur D, Liu A, Takahashi J, Henschen B, Singh S, Lambert L, Furuya T, Bouttenot R, Doll M, Nawaz F, Mu J, Jiang L, Miller LH, Wellems TE. Erythrocyte binding protein PfRH5 polymorphisms determine species-specific pathways of Plasmodium falciparum invasion. Cell Host Microbe. 2008;4(1):40–51.
Pillai AD, Nguitragool W, Lyko B, Dolinta K, Butler MM, Nguyen ST, Peet NP, Bowlin TL, Desai SA. Solute restriction reveals an essential role for clag3-associated channels in malaria parasite nutrient acquisition. Mol Pharmacol. 2012;82(6):1104–14.
Nguitragool W, Rayavara K, Desai SA. Proteolysis at a specific extracellular residue implicates integral membrane CLAG3 in malaria parasite nutrient channels. PLoS One. 2014;9(4):e93759.
Patton C, Thompson S, Epel D. Some precautions in using chelators to buffer metals in biological solutions. Cell Calcium. 2004;35(5):427–31.
Figan CE, Sa JM, Mu J, Melendez-Muniz VA, Liu CH, Wellems TE. A set of microsatellite markers to differentiate Plasmodium falciparum progeny of four genetic crosses. Malar J. 2018;17(1):60.
Reilly Ayala HB, Wacker MA, Siwo G, Ferdig MT. Quantitative trait loci mapping reveals candidate pathways regulating cell cycle duration in Plasmodium falciparum. BMC Genomics. 2010;11:577.
Prinz H. Hill coefficients, dose-response curves and allosteric mechanisms. J Chem Biol. 2010;3:37–44.
Sidhu AB, Verdier-Pinard D, Fidock DA. Chloroquine resistance in Plasmodium falciparum malaria parasites conferred by pfcrt mutations. Science. 2002;298(5591):210–3.
Bennett TN, Kosar AD, Ursos LM, Dzekunov S, Singh Sidhu AB, Fidock DA, Roepe PD. Drug resistance-associated pfCRT mutations confer decreased Plasmodium falciparum digestive vacuolar pH. Mol Biochem Parasitol. 2004;133(1):99–114.
Juge N, Moriyama S, Miyaji T, Kawakami M, Iwai H, Fukui T, Nelson N, Omote H, Moriyama Y. Plasmodium falciparum chloroquine resistance transporter is a H+-coupled polyspecific nutrient and drug exporter. Proc Natl Acad Sci U S A. 2015;112(11):3356–61.
Bakouh N, Bellanca S, Nyboer B, Moliner CS, Karim Z, Sanchez CP, et al. Iron is a substrate of the Plasmodium falciparum chloroquine resistance transporter PfCRT in Xenopus oocytes. J Biol Chem. 2017;292:16109–21.
Patzewitz EM, Salcedo-Sora JE, Wong EH, Sethia S, Stocks PA, Maughan SC, Murray JA, Krishna S, Bray PG, Ward SA, Muller S. Glutathione transport: a new role for PfCRT in chloroquine resistance. Antioxid Redox Signal. 2013;19(7):683–95.
van Schalkwyk DA, Nash MN, Shafik SH, Summers RL, Lehane AM, Smith PJ, Martin RE. Verapamil-sensitive transport of quinacrine and methylene blue via the Plasmodium falciparum chloroquine resistance transporter reduces the parasite's susceptibility to these tricyclic drugs. J Infect Dis. 2016;213(5):800–10.
Tan W, Gou DM, Tai E, Zhao YZ, Chow LM. Functional reconstitution of purified chloroquine resistance membrane transporter expressed in yeast. Arch Biochem Biophys. 2006;452(2):119–28.
Scheibel LW, Colombani PM, Hess AD, Aikawa M, Atkinson CT, Igarashi I, Matsumoto Y, Milhous WK. Calcium/calmodulin functions in P. falciparum in vitro--implications for antiprotozoal drug design. Prog Clin Parasitol. 1989;1:21–56.
Gupta A, Thiruvengadam G, Desai SA. The conserved clag multigene family of malaria parasites: essential roles in host-pathogen interaction. Drug Resist Updates. 2015;18:47–54.
Garg S, Agarwal S, Kumar S, Yazdani SS, Chitnis CE, Singh S. Calcium-dependent permeabilization of erythrocytes by a perforin-like protein during egress of malaria parasites. Nat Commun. 2013;4:1736.
Weiss GE, Gilson PR, Taechalertpaisarn T, Tham WH, de Jong NW, Harvey KL, Fowkes FJ, Barlow PN, Rayner JC, Wright GJ, Cowman AF, Crabb BS. Revealing the sequence and resulting cellular morphology of receptor-ligand interactions during Plasmodium falciparum invasion of erythrocytes. PLoS Pathog. 2015;11(2):e1004670.
Berridge MJ, Lipp P, Bootman MD. The versatility and universality of calcium signalling. Nat Rev Mol Cell Biol. 2000;1(1):11–21.
Kao, J.P., Alderton, J.M., Tsien, R.Y., & Steinhardt, R.A. (1990). Active involvement of Ca2+ in mitotic progression of Swiss 3T3 fibroblasts. J Cell Biol , 111(1), 183–196.
Hernandez-Cruz A, Sala F, Adams PR. Subcellular calcium transients visualized by confocal microscopy in a voltage-clamped vertebrate neuron. Science. 1990;247:858–62.
Baylor SM, Hollingworth S. Simulation of Ca2+ movements within the sarcomere of fast-twitch mouse fibers stimulated by action potentials. J Gen Physiol. 2007;130:283–302.
Alkhalil A, Pillai AD, Bokhari AA, Vaidya AB, Desai SA. Complex inheritance of the plasmodial surface anion channel in a Plasmodium falciparum genetic cross. Mol Microbiol. 2009;72(2):459–69.
Hogan PG, Chen L, Nardone J, Rao A. Transcriptional regulation by calcium, calcineurin. and NFAT Genes Dev. 2003;17(18):2205–32.
Voss TS, Bozdech Z, Bartfai R. Epigenetic memory takes center stage in the survival strategy of malaria parasites. Curr Opin Microbiol. 2014;20:88–95.
Rovira-Graells N, Crowley VM, Bancells C, Mira-Martinez S, Ribas de PL, Cortes A. Deciphering the principles that govern mutually exclusive expression of Plasmodium falciparum clag3 genes. Nucleic Acids Res. 2015;43(17):8243–57.
Krogstad DJ, Gluzman IY, Kyle DE, Oduola AM, Martin SK, Milhous WK, Schlesinger PH. Efflux of chloroquine from Plasmodium falciparum: mechanism of chloroquine resistance. Science. 1987;238(4831):1283–5.
Ginsburg H. Effect of calcium antagonists on malaria susceptibility to chloroquine. Parasitol Today. 1988;4(8):209–11.
Rohrbach P, Sanchez CP, Hayton K, Friedrich O, Patel J, Sidhu AB, Ferdig MT, Fidock DA, Lanzer M. Genetic linkage of pfmdr1 with food vacuolar solute import in Plasmodium falciparum. EMBO J. 2006;25(13):3000–11.
Teaema FH, Al AK. Nineteen cases of symptomatic neonatal hypocalcemia secondary to vitamin D deficiency: a 2-year study. J Trop Pediatr. 2010;56(2):108–10.
Su X, Ferdig MT, Huang Y, Huynh CQ, Liu A, You J, Wootton JC, Wellems TE. A genetic map and recombination parameters of the human malaria parasite Plasmodium falciparum. Science. 1999;286:1351–3.
Sen S, Churchill GA. A statistical framework for quantitative trait mapping. Genetics. 2001;159(1):371–87.
Broman KW, Wu H, Sen S, Churchill GA. R/QTL: QTL mapping in experimental crosses. Bioinformatics. 2003;19(7):889–90.
Srinivasan P, Yasgar A, Luci DK, Beatty WL, Hu X, Andersen J, et al. Disrupting malaria parasite AMA1-RON2 interaction with a small molecule prevents erythrocyte invasion. Nat Commun. 2013;4:2261.
LA and JO were supported by the Intramural NIAID Research Opportunties (INRO) program. This study was supported by the Intramural Research Program of National Institutes of Health, National Institute of Allergy and Infectious Diseases. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
All data generated or analyzed during this study are included in this published article. Raw experimental data including those from parasite growth studies and genetic analysis are available from the corresponding author upon request.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Apolis, L., Olivas, J., Srinivasan, P. et al. Multiple genetic loci define Ca++ utilization by bloodstream malaria parasites. BMC Genomics 20, 47 (2019). https://doi.org/10.1186/s12864-018-5418-y