Skip to main content

Recent increase in low complexity polygenomic infections and sialic acid-independent invasion pathways in Plasmodium falciparum from Western Gambia

Abstract

Background

The malaria parasite Plasmodium falciparum utilizes multiple alternative receptor-ligand interactions for the invasion of human erythrocytes. While some P. falciparum clones make use of sialic acid (SA) residues on the surface of the human glycophorin receptors to invade the erythrocyte, others use alternative receptors independent of sialic acid residues. We hypothesized that over the years, intensified malaria control interventions and declining prevalence in The Gambia have resulted in a selection of parasites with a dominant invasion pathways and ligand expression profiles.

Methods

Blood samples were collected from 65 malaria-infected participants with uncomplicated malaria across 3 years (2015, 2016, and 2021). Genetic diversity was determined by genotyping the merozoite surface protein 2 (msp2) polymorphic gene of P. falciparum. Erythrocyte invasion phenotypes were determined using neuraminidase, trypsin, and chymotrypsin enzymes, known to cleave different receptors from the surface of the erythrocyte. Schizont-stage transcript levels were obtained for a panel of 6 P. falciparum invasion ligand genes (eba175, eba181, Rh2b, Rh4, Rh5, and clag2) using 48 successfully cultured isolates.

Results

Though the allelic heterozygosity of msp2 repeat region decreased as expected with reduced transmission, there was an increase in infections with more than a single msp2 allelotype from 2015 to 2021. The invasion phenotypes of these isolates were mostly SA independent with a continuous increase from 2015 to 2021. Isolates from 2021 were highly inhibited by chymotrypsin treatment compared to isolates from 2015 and 2016. Higher invasion inhibition for 2021 isolates was further obtained following erythrocyte treatment with a combination of chymotrypsin and trypsin. The transcript levels of invasion ligand genes varied across years. However, levels of clag2, a rhoptry-associated protein, were higher in 2015 and 2016 isolates than in 2021 isolates, while Rh5 levels were higher in 2021 compared to other years.

Conclusions

Overall, these findings suggest increasing mixed infections with an increase in the use of sialic-acid independent invasion pathways by P. falciparum clinical isolates in the Western part of Gambia.

Graphical Abstract

Background

Plasmodium falciparum malaria parasite invade and multiply in human red blood cells (RBCs) during the entire erythrocytic phase of its life cycle. Erythrocyte invasion by P. falciparum merozoite stages is a crucial, complex, and multistep process requiring multiple alternative receptor-ligand interactions [1]. Plasmodium falciparum ligands are made up of two protein families: the erythrocyte binding antigen (EBA-175, EBA-140, EBA-165, EBL1, and EBA-181) [2,3,4,5] and the reticulocyte binding protein-like homologs (Rh1, RH2A, Rh2B, Rh4, and Rh5) [6,7,8,9]. These ligands determine P. falciparum invasion phenotypes, which have mostly been characterized by the treatment of erythrocytes with enzymes known to cleave part of the receptor repertoire on which they bind. Commonly used enzymes include neuraminidase, which cleaves sialic acids from the glycophorin (Gly) receptors (GlyA, B, and C), trypsin cleaves glycophorins A, C, and complement receptor 1 (CR1), while chymotrypsin cleaves GlyB, band 3, and other invasion-related proteins [10, 11]. Thus, based on the sialic acid (SA) residues of the glycophorin receptors, P. falciparum invasion pathways are described as either SA dependent or SA independent. These ligand-receptor interactions are variably used by P. falciparum isolates in different malaria populations, and some parasite lines activate alternative pathways or switch invasion phenotypes during the in vitro life cycle [12,13,14,15].

Previous studies in The Gambia have reported the predominant use of the SA-dependent pathways during a period of relatively higher malaria transmission [16, 17]. Over the last 20 years, malaria transmission has been drastically reduced because of intensified and sustained malaria control interventions. Currently, The Gambia is moving towards malaria elimination (pre-elimination phase) with an overall prevalence of 0.2% in the country and incidence of 25–100 cases per 1000 population in the Western region [18]. Hence, how this transition from high to low malaria transmission has affected erythrocyte invasion mechanisms, which involves ligands targeted for vaccine development, is not known.

In this study, we investigated the invasion phenotypes and transcript levels of six invasion ligand genes (eba175, eba181, Rh2b, Rh4, Rh5, clag2) of P. falciparum field isolates from three different years in The Gambia (2015, 2016, and 2021). We also used the highly polymorphic P. falciparum surface protein, merozoite surface protein 2 (msp2), to assess the complexity of infections from these different years. We hypothesize that intensified malaria control interventions and declining prevalence in The Gambia have resulted in a selection of parasites with the most dominant invasion pathways driven by expression of specific ligands.

Methods

Study site and malaria sampling

This study was conducted in the Western Region of The Gambia. Uncomplicated malaria patients attending two health facilities (Brikama and Fajikunda health centres) who tested positive for P. falciparum by immunochromatic rapid diagnostic testing (RDT) and had reported not taking antimalarial drugs during the preceding 3 days were invited to participate in the study.

Patients were recruited in 2015 (n = 15), 2016 (n = 26), and 2021 (n = 24) during the annual peak of malaria transmission (September to December). Isolates from 2015 and 2016 were collected from children ≤ 14 years old while those from 2021 were collected from all age groups. The inclusion criteria were patients with axillary temperatures of > 37.5 °C, or a history of fever in the previous 48 h, and a positive RDT test for malaria. For each consenting qualified case, approximately 2 ml venous blood was collected into EDTA anticoagulant tubes and further confirmed for malaria by light microscopy. Plasma was removed and infected erythrocytes separated from lymphocytes and buffy coat by NycoPrep gradient centrifugation. The infected erythrocytes of samples collected in 2015 were cryopreserved in glycerolyte and stored in liquid nitrogen until use while those collected in 2016 and 2021 were directly resuspended at 2% hematocrit in Roswell Park Memorial Institute 1640 complete medium (RPMI 1640) (incorporating 25 mmol/l HEPES, 2 mmol/l, l glutamine, 25 mmol/l glucose, 25 mg of Gentamicin per litre, 10 mg of hypoxanthine per litre, and 10% human AB serum or 5000 mg/l of Albumax) and used for invasion assays.

DNA extraction and msp2 genotyping

For all patients included in the study, blood was spotted on filter papers (Whatman 3 mm, GE Healthcare, USA) prior to culture (day zero samples), labelled, air-dried, and stored in sealed plastic bags at 4 °C. Plasmodium falciparum DNA was subsequently extracted using the QIAamp DNA Mini Kit (Qiagen, Germany) according to the Qiagen-DNA purification from dried blood spot protocol. Merozoite surface protein 2 (msp2) repeat genotyping was done using specific primers as previously described [19]. Briefly, the primary PCR targeted the polymorphic block 3 region of msp2 while the nested PCR primer sets were specific to msp2 allelic families (3D7 and FC27). The primary PCR was carried out in a final volume of 15 μl containing 2 μl gDNA, 0.076 U/μl Taq polymerase (New England Biolabs), 250 nM forward and reverse primers, 125 μM dNTPs, 1× (1.5 μl) Thermopol buffer (New England Biolabs), and 10.35 μl nuclease-free water. Nested PCR was done at a final volume of 20 μl containing 1 μl of the primary PCR product, 0.05 U/μl Taq polymerase (New England Biolabs), 125 nM forward and reverse primers, 100 μM dNTPs, 1× (2 μl) Thermopol buffer (New England Biolabs) and 16.1 μl nuclease-free water. Cycling conditions used were as previously described [20] with positive (3D7 and Dd2) and negative (nuclease-free water) controls incorporated in each PCR run. Amplificons from the allelic families were analysed using QIAxcel (Qiagen) and fragment sizes were determined using QIAxcel ScreenGel Software.

Plasmodium falciparum ligand gene expression assays

Isolates from samples with adequate volumes and at least 0.5% parasitaemia were cultured to the schizont stage and frozen at − 80 °C in three parts of TRIzol reagent (Ambion/Life Technologies, USA) as previously described [17]. RNA was extracted with the phenol-chloroform RNA extraction protocol according to the instructions of the Tri Reagent® manufacturer and concentrations were estimated using the Qubit fluorometer (Life Technologies, UK). RNA from each isolate was reverse transcribed using the First Strand cDNA Synthesis Kit (New England Biolabs) and the transcript levels for eba175, eba181, clag2, Rh2b, Rh4, and Rh5 genes were determined using gene-specific primer/probe sets [17, 21, 22]. Reverse transcriptase quantitative PCRs were performed using the TaqMan universal PCR Master Mix (ThermoFisher) in 15 μl volumes with 330 nM concentrations of each primer and 160 nM concentrations of each probe at 95 °C for 10 min followed by 40 cycles of 95 °C for 30 s and 50 °C for 1:30 s. Each run included controls and 3D7 genomic DNA standards, with standard curves generated for each run. The threshold fluorescence value for each run was determined automatically by the BioRad CFX96 software.

Erythrocyte invasion assays

Erythrocytes used for invasion assays were from blood group O+ uninfected human volunteers who had not taken antimalarial in the previous 1 month and were not on any medication. Fresh target erythrocytes were treated weekly with neuraminidase (Sigma-Aldrich, UK), trypsin (Sigma-Aldrich, UK), and chymotrypsin (Sigma-Aldrich, UK) in different combinations and concentrations as previously described [23] with some modifications. Briefly, the target erythrocytes were treated with 66.7 mU/ml of neuraminidase (NM), 66.7 µg/ml for low trypsin (LT), 1.0 mg/ml for high trypsin (HT); chymotrypsin and low trypsin (CHY_LT) with 1.0 mg/ml and 66.7 µg/ml, respectively, and 1.0 mg/ml chymotrypsin (CHY). Complete RPMI media were used as the negative control (no invasion inhibition) and a combination of all enzymes was used as the positive control (100% inhibition). All enzyme treatments of erythrocytes were performed at 37 °C for 1 h, followed by two washes with 1× PBS and once with RPMI 1640, and suspension at 2% hematocrit. The treated erythrocytes at 2% hematocrit were fluorescently labelled with Cell Trace(R) Far Red (CTFR) Proliferation Kit (Invitrogen, C34564) as previously described [17, 23].

The invasion phenotypes of P. falciparum clinical isolates across 3 years, 2015 (n = 15), 2016 (n = 26), and 2021 (n = 24) were determined through their ability to invade receptor-depleted CTFR-labelled target erythrocytes. Assays were set up in triplicate in flat-bottom 96-well plates. For each assay, ring-stage parasite isolates from infected patients were added to CTFR-stained erythrocytes in a 1:1 ratio at 2% hematocrit in flat-bottomed 96-well culture plates for a total culture volume of 100 µl per well. Each isolate was tested in triplicate per enzyme-treated cell and plates were incubated at 37 °C in a gas atmosphere of 1% O2, 3% CO2, and 96% N2. After the first complete invasion cycle (48 h), cells were labelled with SYBR Green I nucleic acid gel stain (Invitrogen) to differentiate infected from uninfected erythrocytes. Invasion levels were determined using BD Accuri™ C6 Plus flow cytometer (BD Biosciences), and the percentage of cells stained with both CTFR and SYBR Green I was recorded as newly re-invaded.

Statistical analyses

For msp2 genotyping, the multiplicity of infection (MOI) was calculated by dividing the total number of alleles detected by the total number of samples [24] while the expected heterozygosity (a measure of genetic diversity) was used to assess population structure of parasites. Heterozygosity (He) was calculated using the formula

$${\text{He}}\, = \,\left[ {n/\left( {n - {1}} \right)} \right]\left[ {\left( {{1} - \Sigma {\text{Pi2}}} \right)} \right],$$

where n is the sample size and Pi is the allele frequency as previously described [25].

Gene transcript expression and invasion data analyses were performed with R version 4.2.0 software. The Gaussian distribution of data was assessed by the Shapiro-Wilk normality test, and variables that passed the test for normality were analysed using parametric methods; otherwise, non-parametric procedures were used. Enzyme inhibition rates and relative gene expressions were compared across sampling years using the Kruskal-Wallis test, and post hoc pairwise comparisons were performed using the Mann-Whitney U test with Bonferroni correction of multiple pairwise comparisons. Correlations between the distributions of the different enzyme inhibition phenotypes were tested by Spearman’s rho. P values < 0.05 were considered statistically significant for all analyses.

Results

Participant’s characteristics

A total of 65 malaria-infected people were retained with the lowest proportion from 2015 (23.1%) and similar proportions from 2016 and 2021 (40% and 36.9%, respectively). While the mean age was similar between the 2015 and 2016 participants (9 and 8.21, respectively), it was higher for 2021 participants (20.48). The proportion of male-to-female participants in the study was not different (31 and 30, respectively), and four participants did not indicate their gender. Samples collected in 2015 had higher parasitaemia ranging from 1 to 7% with a median of 1.85% compared to 2016 and 2021 with parasitaemia ranging from 0.21 to 5.2% (median 1.24%) and 0.4 to 3.6% (median 1%), respectively (Table 1).

Table 1 Characteristics of study participants from Western Gambia across three different years

Genetic diversity of P. falciparum infections

There was a drastic reduction of malaria prevalence in The Gambia from 4% in 2010 to 0.2% in 2017 [26]. However, since 2017 till present the overall countrywide prevalence still stands at 0.2% despite continuous control efforts. This is probably due to a stall or a slight rebound of infections in some parts of the country as indicated by the world malaria reports in recent years [18, 27, 28]. This could affect the genetic diversity of infections and patterns of invasion. Therefore, we genotyped the P. falciparum isolates for msp2 gene repeat polymorphisms to investigate whether parasites collected in the 3 years were genetically different. Alleles were classified according to the size of amplified fragments. Thirty-two individual alleles were detected with 22 alleles for 3D7 (fragment range 230–700 bp) (Fig. 1A) and 10 alleles for FC27 (fragment range 250–450 bp) (Fig. 1B). The frequencies of all alleles detected in the 3D7 family was < 20% and the 350 bp allele was the most prevalent in both 3D7 and FC27 families (18% and 39%, respectively).

Fig. 1
figure 1

Prevalence and allelic frequencies of msp2. A Frequencies of 3D7 allelic family, B frequencies of FC27 allelic family, C prevalence of msp2 repeat polymorphism allelic families of Plasmodium falciparum in clinical isolates collected in 2015, 2016, and 2021 from Western Gambia. bp base pairs

The frequency of likely monogenomic infections with single 3D7 allele was higher in 2016 (61.5%) compared to 2015 (40%) and 2021 (40.9%). However, these differences were not statistically significant (P = 0.37). Furthermore, there was an increase in polygenomic infections of the same 3D7 allelotype in 2021 (13.6%) compared to 2015 (6.6%) and 2016 (7.7%), but these differences were also not statistically significant (P = 0.4). Somewhat similarly, isolates with the FC27 allelotype had a higher frequency (33.3%) of likely monogenomic infections in 2015, while polygenomic FC27 types were detected in 2016 and 2021. Complex (polygenomic) infections with both the 3D7/FC27 allelic types were more prevalent in the 2021 isolates with a frequency of 36.6% compared to 2015 (20%) and 2016 (15.4%) (Fig. 1C).

Overall, the presence of multiple clones in a single infection as defined by the multiplicity of infection (MOI) index was higher in 2021 isolates (1.55) than in 2015 (1.33) and 2016 (1.27). This was contrary to an overall decrease in the heterozygosity (He) of the msp2 gene from 2015 (0.29) through 2016 (0.2) to 2021 (0.06) (Table 2).

Table 2 MOI and heterozygosity of msp2 gene of Plasmodium falciparum across 3 years

Transcript analysis of invasion ligands

Specific merozoite ligand genes are predominately expressed by P. falciparum isolates during the invasion process and could determine the pathways used. To determine variation in transcript levels of six invasion ligands (CLAG-2, EBA-175, EBA-181, RH-2B, RH4, and RH5), each isolate was cultured ex vivo to predominantly schizont stages and bulk RNA extracted. A total of 48 of the 65 isolates genotyped were successfully cultured for transcript analysis [2015 (n = 12), 2016 (n = 15), and 2021 (n = 21)]. For the other 17 isolates, data were not available due to either poor growth in culture or too low RNA yield after extraction to allow for reliable transcript quantification. The combined transcript profiles across all ligands grouped the isolates into four main clusters (A, B, C, and D) (Fig. 2A). Transcript level profiles for isolates clustered RH4 and RH-2B as well as RH5 and RH-2B together, which are associated with sialic acid-independent invasion pathway, while EBA-181 and EBA-175 genes associated with sialic acid-dependent invasion pathway also clustered together (Fig. 2A). Thus, significant positive correlations were observed between EBA-175 and EBA-181 (Spearman’s r = 0.5, P = 2e−04), RH-2B and RH4 (r = 0.5, P = 0.001), and RH5 and RH4 (Spearman’s r = 0.7, P = 1e−06) (Fig. 2B).

Fig. 2
figure 2

Relative gene expression and correlation between gene expression levels. A Hierarchical clustering heatmap of relative expression levels of six Plasmodium falciparum ligand genes (rows) between clinical isolates (columns). White to light red colour signifies lower expression while brown colour indicates higher expression. B Correlation between gene expression levels; numbers out of the brackets in square boxes represent the correlation coefficients while those in brackets represent the P values. Significant P values are indicated in red

Overall, CLAG2 (mean 2.61) and RH5 (mean 2.93) showed the relative highest copy numbers of transcripts across the 3 years with no significant temporal differences between their expression (P = 0.41), followed by EBA-181 (mean 1.05) and EBA-175 (mean 1.03) (Additional file 1: Fig. S1). In addition, analysis of the invasion efficiency of isolates according to the ligand genes expression clusters in Fig. 2A showed no statistically significant differences in invasion efficiency of all enzyme treatments across all clusters (Additional file 1: Fig. S2).

Fig. 3
figure 3

Invasion phenotypes of P. falciparum isolates across 3 years in The Gambia. NM neuraminidase, LT low trypsin, HT high trypsin, CHY_LT chymotrypsin and low trypsin, and CHY chymotrypsin. Each point represents an isolate while the vertical line of each box represents the median. P < 0.05 was considered statistically significant

Invasion inhibition of P. falciparum clinical isolates across three years

Following depletion of receptors using various enzymes, we observed very low levels (mean 27.3%) of invasion inhibition into neuraminidase (NM)-treated erythrocytes by isolates from all 3 years. Higher levels of inhibition were observed for low trypsin (LT; 45.6%, P = 6.1 e−6), high trypsin (HT; 86.6%, P = 2 e−16), and chymotrypsin (CHY; 84.9%, P = 2 e−16). Combining chymotrypsin and trypsin (CHY_LT) resulted in higher invasion inhibition (mean 88.9%) than when the enzymes were used separately but these differences were not statistically significant (CHY_LT and CHY P = 0.058, CHY_LT and HT P = 0.065 (Fig. 3).

Comparison of invasion phenotypes across the 3 years showed a trend towards increased neuraminidase inhibition between 2015 to 2021 (mean: 2015 = 22.21%, 2016 = 26.69%, and 2021 = 31.6%) but the differences were not statistically significant (Kruskal-Wallis test, P = 0.3) (Fig. 4A). Low trypsin-treated erythrocytes resulted in a significant increase (P = 0.004 between 2015 and 2016 and 0.0045 between 2016 and 2021) in invasion inhibition between years with a mean of 37% in 2015 and 41.1% in 2021, and isolates collected in 2016 showed a higher (mean 53.65%) inhibition rate compared to the preceding year 2015 isolates (uncorrected P = 0.004, P corrected for multiple pairwise comparisons = 0.012) (Fig. 4B). Furthermore, inhibition by chymotrypsin was significantly higher in 2016 isolates (mean 87.38%) compared to 2015 (Mann-Whitney uncorrected P = 0.0023, corrected P = 0.007 (Fig. 4C) and even more in 2021 (mean 90.6%) compared to 2015 isolates (Mann-Whitney uncorrected P = 0.00037, corrected P = 0.001). However, inhibition by this enzyme combination (chymotrypsin and low trypsin) was more effective across all years (mean: 2015 = 84.57%, 2016 = 91.65%, and 2021 = 88.35) (Fig. 4D).

Fig. 4
figure 4

Comparison of erythrocyte invasion inhibition phenotypes of Plasmodium falciparum isolates across 3 years in The Gambia. A Neuraminidase, B low trypsin, C chymotrypsin, and D chymotrypsin and low trypsin. Each box shows the invasion efficiencies in enzyme-treated erythrocytes relative to untreated erythrocytes, and the horizontal black lines in each box represent the median. P values < 0.05 were considered statistically significant

Following the observed invasion inhibition phenotypes and variations across the years, we tested whether the phenotypes also varied by msp2 allelic families. For the most dominant allelic type (350 bp) for both families (3D7 and FC27), invasion inhibition was significantly higher in 3D7 allelic family only in the case of erythrocyte treatment with high trypsin (P = 0.027) (Additional file 1: Fig. S3). For neuraminidase-treated cells the 3D7 type isolates were also less inhibited, even though the difference with FC27 type isolates was not significant.

Correlation between invasion inhibition by enzyme treatments and age

The magnitude and pattern of inhibition of invasion by CHY, CHY_LT, and NM were positively correlated with age but negatively correlated with parasitaemia while HT was negatively correlated with age but positively correlated with parasitaemia. These correlations among age, parasitaemia, and enzyme treatments were not statistically significant, although marginally significant for HT inhibition and age (0.05) (Fig. 5). Further grouping of the study participants into three age groups as young children (0–5 years), older children (6–17 years), and adults (18 years and above) resulted in significant differences among age groups and enzyme treatments (Additional file 1: Fig. S4). Inhibition by low trypsin and chymotrypsin significantly correlated negatively with parasitaemia in younger children (Spearman’s r = − 0.76 and r = − 0.82, respectively) (Additional file 1: Fig. S4).

Fig. 5
figure 5

Correlation among enzyme inhibition, age, and parasitaemia. Numbers out of the brackets represent the correlation coefficients while those in brackets indicate the P-values with significant P-values highlighted in red. PCT parasitaemia

Discussion

This study characterized and compared the genetic diversity, ligand genes expression, and invasion phenotypes of P. falciparum clinical isolates across 3 years in Western Gambia, where enhanced control interventions for over a decade have led to a significant decline in malaria transmission and clinical incidence. This overall decline in malaria transmission was reflected by lower heterozygosity in msp2 polymorphic repeat alleles in this study compared to previous reports. However, world malaria reports in the recent years covered by the study indicate a slight rebound of infections in The Gambia [18, 27, 28]. Thus, we observed an increase in polygenomic infections in the infections sampled from 2021 compared to 2015 and 2016, despite the lower overall heterozygosity. There were more infections in 2021 with more than a single clone despite the fewer overall number of clones in circulation. This could be attributed to the slight increase in transmission. Such a phenomenon had been described between populations in Cameroon, where high multiplicity of infections was observed despite the low number of circulating genotypes [29]. This observation is unexpected since The Gambia is in the pre-elimination phase of malaria and MOI is expected to be decreasing over the years. Moreover, it is generally expected that an increase in MOI should lead to increasing heterozygosity but that was not the case in this study [30,31,32]. It is probable that despite the low transmission in the Western region of The Gambia, multiple clones with few polymorphic loci have been maintained or are circulating in the population or imported in this most urban region of the country. These results should be interpreted with caution because of the small sample size, and sampling only from Western Gambia does not reflect the entire country. Moreover, only the msp2 gene was genotyped, and the pattern could differ if diversity at other polymorphic loci such as microsatellites and genome-wide single-nucleotide polymorphisms (SNPs) was evaluated.

Parasitaemia is established in malaria following several merozoite ligand-erythrocyte receptor interactions, merozoite invasion and replication in erythrocytes [10]. Erythrocyte invasion can occur through several alternative pathways, which could be determined by the expression levels and type or use of different ligand combinations [33,34,35]. Analyses of transcript levels relative to apical membrane antigen 1 (AMA1) gene for selected invasion ligands in our clinical isolates showed significant differences in the rhoptry protein, RH5 [36, 37], relative to EBA-181 and EBA-175. RH5 is a host tropism determinant and an unavoidable ligand in the invasion process, binding to basigin (CD147) and P113 [38, 39]. It is one of the leading blood-stage P. falciparum malaria vaccine antigens, and its higher expression relative to EBA-175 is in contrast to previous studies in Ghana and The Gambia, where dominance of EBA-175 transcripts and SA-dependent pathways was reported [17, 40]. Unlike isolates from these previous years, most of those we characterised here were SA independent. RH5 expression correlated strongly with that of RH4 and RH2b, two ligands also important for SA-independent invasion, supporting the switch from SA-dependent invasion in The Gambia. Parasites may vary invasion pathways to non-EBA-dependent and SA-independent pathways, given immune pressure against the EBA family of proteins as demonstrated from disruption of EBA-175 gene [41]. Malaria has been in a decline in The Gambia, with expected reduction in population level immunity [42, 43]. While this should ideally reduce the pressure on EBA-175, and allow for sustained SA-dependent invasion, the reverse seen here points to immune pressure not being the only factor in pathway switch. Indeed, pathway switch to more efficient mechanisms may be conditioned by the environment as seen during in vitro culture adaptation [44]. High expression was also detected for CLAG-2, a gene that is associated with the rhoptry bulb, and a member of a multigene family of rhoptry proteins, implicated in cytoadherence, infected cell permeability, and invasion [45, 46]. These results together further indicate that the level of antigen expression alone cannot explain the variation in patterns of invasion pathways, as these correlate with other infection matrices. For example, there was a positive correlation between RH5 expression and parasitaemia as previously reported from The Gambia [17], Ghana [40], Mali [47], and Papua New Guinea [48]. RH5 is an essential invasion antigen, a potent vaccine candidate in advanced development. Its protective role against infections and disease needs to be further evaluated as malaria parasite populations evolve against intense interventions with changing epidemiology.

The increased use of SA-independent receptors and neuraminidase resistance from 2015 to 2021 in the Western region of The Gambia may be because of adaptation to reduced transmission pressure. Clinical isolates of P. falciparum from other regions in Africa have also shown similar levels of resistance to neuraminidase treatment [23, 40, 49]. Most of these were also isolates from recent populations, in contrast to previous findings of field isolates collected over a decade ago in The Gambia, where the parasites’ ability to invade neuraminidase-treated erythrocytes was much higher [16, 17]. Malaria transmission during the earlier studies was more intense in Western Gambia, and the variation observed here could also be due to some of the precious specimens tested from severe malaria cases compared to the mostly low parasitaemia uncomplicated cases in this study [17]. There is also the possibility of variations in the activities of commercial neuraminidase enzyme batches used. Overall, the use of different invasion pathways remains heterogeneous, and the molecular patterns determining this remain complicated, needing further evaluation as parasite ligands are targeted for vaccine development. This study is limited by the small number of specimens from each year, a factor in other invasion studies most likely due to the significant capacity and resources needed for parasite culturing and cellular assays.

Conclusion

In conclusion, this study showed a small recent increase in low complexity polygenomic infections in the Western region of The Gambia and suggests that recent clinical isolates of P. falciparum from uncomplicated malaria are increasingly using sialic acid-independent invasion pathways against the sialic acid-dependent alternative. With The Gambia accelerating interventions towards pre-elimination, and possibility of new elimination tools (e.g., vaccines), continuous molecular and phenotypic surveillance of the parasite population is necessary.

Availability of data and materials

Data supporting the conclusion of this article are included within the article.

Abbreviations

Rh5:

Reticulocyte binding protein homolog 5

Rh4:

Reticulocyte binding protein homolog 4

Rh2b:

Reticulocyte binding protein homolog 2b

CR1:

Complement receptor 1

CLAG-2:

Cytoadherence linked asexual protein 2

EBA-181:

Erythrocyte binding antigen 181

EBA-175:

Erythrocyte binding antigen 175

Gly:

Glycophorin

RDT:

Rapid diagnostic test

EDTA:

Ethylenediaminetetraacetic acid

SNPs:

Single nucleotide polymorphisms

AMA1:

Apical membrane antigen 1

CTFR:

Cell trace far red

SA:

Sialic acid

References

  1. Gaur D, Mayer DCG, Miller LH. Parasite ligand-host receptor interactions during invasion of erythrocytes by Plasmodium merozoites. Int J Parasitol. 2004;34:1413–29.

    Article  CAS  PubMed  Google Scholar 

  2. Gilberger T-W, Thompson JK, Triglia T, Good RT, Duraisingh MT, Cowman AF. A novel erythrocyte binding antigen-175 paralogue from Plasmodium falciparum defines a new trypsin-resistant receptor on human erythrocytes. J Biol Chem. 2003;278:14480–6.

    Article  CAS  PubMed  Google Scholar 

  3. Mayer DCG, Cofie J, Jiang L, Hartl DL, Tracy E, Kabat J, et al. Glycophorin B is the erythrocyte receptor of Plasmodium falciparum erythrocyte-binding ligand, EBL-1. Proc Natl Acad Sci. 2009;106:5348–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Narum DL, Fuhrmann SR, Luu T, Sim BKL. A novel Plasmodium falciparum erythrocyte binding protein-2 (EBP2/BAEBL) involved in erythrocyte receptor binding. Mol Biochem Parasitol. 2002;119:159–68.

    Article  CAS  PubMed  Google Scholar 

  5. Sim BK, Orlandi PA, Haynes JD, Klotz FW, Carter JM, Camus D, et al. Primary structure of the 175K Plasmodium falciparum erythrocyte binding antigen and identification of a peptide which elicits antibodies that inhibit malaria merozoite invasion. J Cell Biol. 1990;111:1877–84.

    Article  CAS  PubMed  Google Scholar 

  6. DeSimone TM, Bei AK, Jennings CV, Duraisingh MT. Genetic analysis of the cytoplasmic domain of the PfRh2b merozoite invasion protein of Plasmodium falciparum. Int J Parasitol. 2009;39:399–405.

    Article  CAS  PubMed  Google Scholar 

  7. Rayner JC, Vargas-Serrato E, Huber CS, Galinski MR, Barnwell JW. A Plasmodium falciparum homologue of Plasmodium vivax reticulocyte binding protein (PvRBP1) defines a trypsin-resistant erythrocyte invasion pathway. J Exp Med. 2001;194:1571–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Tham W-H, Wilson DW, Reiling L, Chen L, Beeson JG, Cowman AF. Antibodies to reticulocyte binding protein-like homologue 4 inhibit invasion of Plasmodium falciparum into human erythrocytes. Infect Immun. 2009;77:2427–35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Triglia T, Thompson J, Caruana SR, Delorenzi M, Speed T, Cowman AF. Identification of proteins from Plasmodium falciparum that are homologous to reticulocyte binding proteins in Plasmodium vivax. Infect Immun. 2001;69:1084–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Cowman AF, Crabb BS. Invasion of red blood cells by malaria parasites. Cell. 2006;124:755–66.

    Article  CAS  PubMed  Google Scholar 

  11. Triglia T, Chen L, Lopaticki S, Dekiwadia C, Riglar DT, Hodder AN, et al. Plasmodium falciparum merozoite invasion is inhibited by antibodies that target the PfRh2a and b binding domains. PLoS Pathog. 2011;7:e1002075.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Hadley TJ, Klotz FW, Pasvol G, Haynes JD, McGinniss MH, Okubo Y, et al. Falciparum malaria parasites invade erythrocytes that lack glycophorin A and B (MkMk). Strain differences indicate receptor heterogeneity and two pathways for invasion. J Clin Invest. 1987;80:1190–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Mitchell G, Hadley T, McGinniss M, Klotz F, Miller L. Invasion of erythrocytes by Plasmodium falciparum malaria parasites: evidence for receptor heterogeneity and two receptors. Blood. 1986;67:1519–21.

    Article  CAS  PubMed  Google Scholar 

  14. Okoyeh JN, Pillai CR, Chitnis CE. Plasmodium falciparum field isolates commonly use erythrocyte invasion pathways that are independent of sialic acid residues of glycophorin A. Infect Immun. 1999;67:5784–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Perkins ME, Holt EH. Erythrocyte receptor recognition varies in Plasmodium falciparum isolates. Mol Biochem Parasitol. 1988;27:23–34.

    Article  CAS  PubMed  Google Scholar 

  16. Baum J, Pinder M, Conway DJ. Erythrocyte invasion phenotypes of Plasmodium falciparum in The Gambia. Infect Immun. 2003;71:1856–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Gomez-Escobar N, Amambua-Ngwa A, Walther M, Okebe J, Ebonyi A, Conway DJ. Erythrocyte invasion and merozoite ligand gene expression in severe and mild Plasmodium falciparum malaria. J Infect Dis. 2010;201:444–52.

    Article  CAS  PubMed  Google Scholar 

  18. WHO. World malaria report 2022. Geneva: World Health Organization; 2022. https://www.who.int/publications/i/item/9789240064898. Accessed 21 Feb 2023.

  19. Snounou G, Zhu X, Siripoon N, Jarra W, Thaithong S, Brown KN, et al. Biased distribution of msp1 and msp2 allelic variants in Plasmodium falciparum populations in Thailand. Trans R Soc Trop Med Hyg. 1999;93:369–74.

    Article  CAS  PubMed  Google Scholar 

  20. Oboh MA, Ndiaye T, Diongue K, Ndiaye YD, Sy M, Deme AB, et al. Allelic diversity of MSP1 and MSP2 repeat loci correlate with levels of malaria endemicity in Senegal and Nigerian populations. Malar J. 2021;20:38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Blair PL. Transcripts of developmentally regulated Plasmodium falciparum genes quantified by real-time RT-PCR. Nucleic Acids Res. 2002;30:2224–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Nery S, Deans A-M, Mosobo M, Marsh K, Rowe JA, Conway DJ. Expression of Plasmodium falciparum genes involved in erythrocyte invasion varies among isolates cultured directly from patients. Mol Biochem Parasitol. 2006;149:208–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Bowyer PW, Stewart LB, Aspeling-Jones H, Mensah-Brown HE, Ahouidi AD, Amambua-Ngwa A, et al. Variation in Plasmodium falciparum erythrocyte invasion phenotypes and merozoite ligand gene expression across different populations in areas of malaria endemicity. Infect Immun. 2015;83:2575–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Mayengue PI, Ndounga M, Davy MM, Tandou N, Ntoumi F. In vivo chloroquine resistance and prevalence of the pfcrt codon 76 mutation in Plasmodium falciparum isolates from the Republic of Congo. Acta Trop. 2005;95:219–25.

    Article  CAS  PubMed  Google Scholar 

  25. Nei M. Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics. 1978;89:583–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Mwesigwa J, Achan J, Di Tanna GL, Affara M, Jawara M, Worwui A, et al. Residual malaria transmission dynamics varies across The Gambia despite high coverage of control interventions. PLoS ONE. 2017;12:e0187059.

    Article  PubMed  PubMed Central  Google Scholar 

  27. WHO. World malaria report 2020: 20 years of global progress and challenges. Geneva: World Health Organization; 2020. https://apps.who.int/iris/handle/10665/337660. Accessed 17 May 2022.

  28. WHO. World malaria report 2021. Geneva: World Health Organization; 2021. https://apps.who.int/iris/handle/10665/350147. Accessed 17 May 2022.

  29. Atuh NI, Anong DN, Jerome F-C, Oriero E, Mohammed NI, D’Alessandro U, et al. High genetic complexity but low relatedness in Plasmodium falciparum infections from Western Savannah Highlands and coastal equatorial Lowlands of Cameroon. Pathog Glob Health. 2022;116:428–37.

    Article  CAS  PubMed  Google Scholar 

  30. Fola AA, Abby Harrison GL, Hazairin MH, Barnadas C, Hetzel MW, Iga J, et al. Higher complexity of infection and genetic diversity of Plasmodium vivax than Plasmodium falciparum across all malaria transmission zones of Papua New Guinea. Am J Trop Med Hyg. 2017;96:630.

    PubMed  PubMed Central  Google Scholar 

  31. Roh ME, Tessema SK, Murphy M, Nhlabathi N, Mkhonta N, Vilakati S, et al. High genetic diversity of Plasmodium falciparum in the low-transmission setting of the Kingdom of Eswatini. J Infect Dis. 2019;220:1346–54.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Touray AO, Mobegi VA, Wamunyokoli F, Herren JK. Diversity and Multiplicity of P. falciparum infections among asymptomatic school children in Mbita, Western Kenya. Sci Rep. 2020;10:5924.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Dolan SA, Miller LH, Wellems TE. Evidence for a switching mechanism in the invasion of erythrocytes by Plasmodium falciparum. J Clin Invest. 1990;86:618–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Duraisingh MT, Maier AG, Triglia T, Cowman AF. Erythrocyte-binding antigen 175 mediates invasion in Plasmodium falciparum utilizing sialic acid-dependent and -independent pathways. Proc Natl Acad Sci. 2003;100:4796–801.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Taylor HM, Grainger M, Holder AA. Variation in the expression of a Plasmodium falciparum protein family implicated in erythrocyte invasion. Infect Immun. 2002;70:5779–89.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Douglas AD, Williams AR, Illingworth JJ, Kamuyu G, Biswas S, Goodman AL, et al. The blood-stage malaria antigen PfRH5 is susceptible to vaccine-inducible cross-strain neutralizing antibody. Nat Commun. 2011;2:601.

    Article  PubMed  Google Scholar 

  37. Illingworth JJ, Alanine DG, Brown R, Marshall JM, Bartlett HE, Silk SE, et al. Functional comparison of blood-stage Plasmodium falciparum malaria vaccine candidate antigens. Front Immunol. 2019;10:1254.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Wanaguru M, Liu W, Hahn BH, Rayner JC, Wright GJ. RH5–Basigin interaction plays a major role in the host tropism of Plasmodium falciparum. Proc Natl Acad Sci. 2013;110:20735–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Campeotto I, Galaway F, Mehmood S, Barfod LK, Quinkert D, Kotraiah V, et al. The structure of the cysteine-rich domain of Plasmodium falciparum P113 identifies the location of the RH5 binding site. MBio. 2020;11:e01566-20.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Mensah-Brown HE, Amoako N, Abugri J, Stewart LB, Agongo G, Dickson EK, et al. Analysis of erythrocyte invasion mechanisms of Plasmodium falciparum clinical isolates across 3 malaria-endemic areas in Ghana. J Infect Dis. 2015;212:1288–97.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Stubbs J, Simpson KM, Triglia T, Plouffe D, Tonkin CJ, Duraisingh MT, et al. Molecular mechanism for switching of P. falciparum invasion pathways into human erythrocytes. Science. 2005;309:1384–7.

    Article  CAS  PubMed  Google Scholar 

  42. Ceesay SJ, Casals-Pascual C, Erskine J, Anya SE, Duah NO, Fulford AJ, et al. Changes in malaria indices between 1999 and 2007 in The Gambia: a retrospective analysis. The Lancet. 2008;372:1545–54.

    Article  Google Scholar 

  43. Ceesay SJ, Casals-Pascual C, Nwakanma DC, Walther M, Gomez-Escobar N, Fulford AJC, et al. Continued decline of malaria in The Gambia with implications for elimination. PLoS ONE. 2010;5:e12242.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Awandare GA, Nyarko PB, Aniweh Y, Ayivor-Djanie R, Stoute JA. Plasmodium falciparum strains spontaneously switch invasion phenotype in suspension culture. Sci Rep. 2018;8:5782.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Nguitragool W, Bokhari AAB, Pillai AD, Rayavara K, Sharma P, Turpin B, et al. Malaria parasite clag3 genes determine channel-mediated nutrient uptake by infected red blood cells. Cell. 2011;145:665–77.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Trenholme KR, Gardiner DL, Holt DC, Thomas EA, Cowman AF, Kemp DJ. clag9: A cytoadherence gene in Plasmodium falciparum essential for binding of parasitized erythrocytes to CD36. Proc Natl Acad Sci. 2000;97:4029–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. 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;209:789–98.

    Article  CAS  PubMed  Google Scholar 

  48. Chiu CYH, Healer J, Thompson JK, Chen L, Kaul A, Savergave L, et al. Association of antibodies to Plasmodium falciparum reticulocyte binding protein homolog 5 with protection from clinical malaria. Front Microbiol. 2014. https://doi.org/10.3389/fmicb.2014.00314/abstract.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Deans A-M, Nery S, Conway DJ, Kai O, Marsh K, Rowe JA. Invasion pathways and malaria severity in Kenyan Plasmodium falciparum clinical isolates. Infect Immun. 2007;75:3014–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We are grateful to all the malaria patients for their participation in the study and the field workers who facilitated the sample collection at Brikama and Fajikunda Health Centres. We appreciate the support of many colleagues at the Medical Research Council Unit The Gambia at London School of Hygiene and Tropical Medicine (MRCG@LSHTM), especially blood group O volunteer donors.

Funding

This work was supported by the Pan-African Malaria Genetic Epidemiology Network (PAMGEN) PhD fellowship. PAMGEN is a team of African scientists studying the genetic interactions between human populations and malaria parasites in different environmental settings across Africa with funding from the Wellcome Trust through the African Academy of Science H3Africa Program (H3A/18/002).

Author information

Authors and Affiliations

Authors

Contributions

Funding acquisition: AAN, LNA, UD, and GAA. Experimental design and conception: AAN and LNA. Sample collection: ECO, AAN, FB, and NNN, Experiments: NNN, FB, AO, ASJ, NFD, HM, and SC. Data analysis: NNN and AAN, Manuscript drafting: NNN, Critical review, and editing: AAN, LNA, ECO, AO, and HM. Read and approved the final version: all authors.

Corresponding author

Correspondence to Alfred Amambua-Ngwa.

Ethics declarations

Ethics approval and consent to participate

The study protocol was approved by The Gambia's government and MRC Joint Ethics Committee. Written informed consent was obtained from each participant or parent/guardian of minors before enrolment.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

: Figure S1. Fold gene expression of genes encoding six invasion ligands in ex vivo schizont-stage cultures of Plasmodium falciparum clinical isolates across 3 years in The Gambia. (A) combined expression of isolates across all years (n = 48), (B) 2015 isolates (n = 12), (C) 2016 isolates (n = 15), and (D) 2021 isolates (n = 11). Horizontal lines indicate the mean for each transcript across all isolates sampled in each year, and each dot denotes the transcript for each gene in a single clinical isolate. Figure S2. Invasion inhibition phenotypes of Plasmodium falciparum clinical isolates according to gene expression clusters. Figure S3. Invasion inhibition phenotypes of Plasmodium falciparum clinical isolates by dominant allele per msp2 gene allelic family. Neuraminidase (NM), low trypsin (LT), high trypsin (HT), chymotrypsin and low trypsin (CHY_LT), and chymotrypsin (CHY). Figure S4. Scatter matrix plot of enzyme treatment and parasitaemia (x-axes and y-axes) between different age groups. Within each panel, each bullet point represents the % invasion inhibition of each enzyme. The asterisk (*) represents groups with significant P-values. NM = neuraminidase, LT = low trypsin, HT = high trypsin, CHY_LT = chymotrypsin/low trypsin, PCT = parasitaemia, and CHY = chymotrypsin.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Nganyewo, N.N., Bojang, F., Oriero, E.C. et al. Recent increase in low complexity polygenomic infections and sialic acid-independent invasion pathways in Plasmodium falciparum from Western Gambia. Parasites Vectors 16, 309 (2023). https://doi.org/10.1186/s13071-023-05929-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-023-05929-4

Keywords

  • Erythrocyte
  • Invasion
  • P. falciparum
  • Ligand genes
  • Expression
  • Malaria
  • Inhibition