Skip to main content

Genetic diversity and natural selection on the thrombospondin-related adhesive protein (TRAP) gene of Plasmodium falciparum on Bioko Island, Equatorial Guinea and global comparative analysis

Abstract

Background

Thrombospondin-related adhesive protein (TRAP) is a transmembrane protein that plays a crucial role during the invasion of Plasmodium falciparum into liver cells. As a potential malaria vaccine candidate, the genetic diversity and natural selection of PfTRAP was assessed and the global PfTRAP polymorphism pattern was described.

Methods

153 blood spot samples from Bioko malaria patients were collected during 2016–2018 and the target TRAP gene was amplified. Together with the sequences from database, nucleotide diversity and natural selection analysis, and the structural prediction were preformed using bioinformatical tools.

Results

A total of 119 Bioko PfTRAP sequences were amplified successfully. On Bioko Island, PfTRAP shows its high degree of genetic diversity and heterogeneity, with π value for 0.01046 and Hd for 0.99. The value of dN–dS (6.2231, p < 0.05) hinted at natural selection of PfTRAP on Bioko Island. Globally, the African PfTRAPs showed more diverse than the Asian ones, and significant genetic differentiation was discovered by the fixation index between African and Asian countries (Fst > 0.15, p < 0.05). 667 Asian isolates clustered in 136 haplotypes and 739 African isolates clustered in 528 haplotypes by network analysis. The mutations I116T, L221I, Y128F, G228V and P299S were predicted as probably damaging by PolyPhen online service, while mutations L49V, R285G, R285S, P299S and K421N would lead to a significant increase of free energy difference (ΔΔG > 1) indicated a destabilization of protein structure.

Conclusions

Evidences in the present investigation supported that PfTRAP gene from Bioko Island and other malaria endemic countries is highly polymorphic (especially at T cell epitopes), which provided the genetic information background for developing an PfTRAP-based universal effective vaccine. Moreover, some mutations have been shown to be detrimental to the protein structure or function and deserve further study and continuous monitoring.

Background

Malaria is a major public health threat in many parts of the globe and is responsible for half a million deaths annually. According to the World Malaria Report 2019, an estimated 228 million persons suffered from malaria worldwide, with 435,000 malaria deaths in 2018 [1]. Plasmodium falciparum is the parasite that cause 99.7% malaria cases in African regions [1].

The sub-Saharan African country, Equatorial Guinea has a total population of 1.31 million (2018). Bioko, an island of Equatorial Guinea off the coast of Cameroon, with historically high malaria transmission. The Bioko Island Malaria Elimination Project (BIMEP) is a fusion of two long-standing anti-malaria programmes (Bioko Island Malaria Control Project, BIMCP and Equatorial Guinea Malaria Vaccine Initiative, EGMVI) in Equatorial Guinea, with Medical Care Development Institution (MCDI) as the lead implementing partner [2,3,4]. Recently evaluations indicated that malaria prevalence had dropped considerably from 43.3 to 10.5% between 2004 and 2016, resulting in a 13.3% reduction of moderate to severe anaemia in children aged 1–5 years. Despite considerable success in reducing the burden on the island, malaria is still a major public health concern in Bioko Island.

Malaria infection in humans starts when infected female Anopheles mosquitoes take their blood meal and inject sporozoites into the host (human) skin. Sporozoites quickly pass into the liver, where they infect hepatocytes. Thrombospondin-related adhesive protein (TRAP) is the most extensively studied Plasmodium transmembrane protein and it would be released on the surface when the sporozoite into contact with the host cell [5]. Plasmodium TRAP also assists the sporozoite in several pivotal functions, such as sporozoite gliding motility, hepatocyte invasion and establishment of infection in the vertebrate host [6,7,8]. The P. falciparum TRAP (PfTRAP) extracellular domain (ECD) consists of three domains/motifs that include the A-domain (similar to the A-or I-domain which is found in integrins), the TSP (thrombospondin repeat motif, a heparin-binding module, also called the RII region) and a proline-rich segment at the C-terminus [9]. The previous studies [10] revealed a higher frequency of nonsynonymous to synonymous single nucleotide polymorphisms (SNPs) of TRAP within the P. falciparum population of Gambian and Thailand. In their studies, McDonald–Kreitman test showed that the ratio of the number of nonsynonymous to synonymous SNPs within P. falciparum was significantly higher than that of the number of nonsynonymous to synonymous fixed sites between P. falciparum and Plasmodium reichenowi. Furthermore, the value of Tajima’D test and Fu and Li’s test also suggested that the TRAP gene is under diversifying selection in the P. falciparum population in Gambian and Thailand [10]. However, no evidence for balancing selection was reported for parasites obtained from sub-Saharan African countries, except The Gambia. Thus, it remains unclear whether positive selection acts on the TRAP gene of P. falciparum in other geographic regions of sub-Saharan Africa.

The aims of the present study were to investigate whether the TRAP gene is under diversifying selection in Bioko P. falciparum and to elucidate how TRAP gene is differentiated among P. falciparum populations. This research would be helpful not only for understanding the molecular evolution of the TRAP gene in P. falciparum but also for the improvement of peptide vaccines based on the TRAP antigen.

Methods

Study area

This study was carried out in Malabo Regional Hospital, in private diagnostic clinics (PDCs) in different regions of Bioko Island, and in the clinic of the Chinese medical aid team to the Republic of Equatorial Guinea. Bioko is an island 32 km off the west coast of Africa, and the northernmost part of Equatorial Guinea. The island has a population of 334,463 inhabitants (2015 census, of which approximately 90% live in Malabo, the capital city of Equatorial Guinea) and a humid tropical environment. The launch of the BIMCP has had a marked impact on malaria transmission, with the decrease of parasite prevalence from over 45% in 2004 to 8.5% in 2016 and the reduction of entomological inoculation rate (EIR) from more than 1000 before 2004 to 14 in 2015 (www.mcdinternational.org) in Bioko Island.

Samples collection

A total of 153 blood spot samples were collected from the patients with uncomplicated malaria during January 2016–December 2018. Consents were obtained from all participating subjects or their parents and the ethical approval was obtained from the Ethics Committee of Malabo Regional Hospital. Included patients were aged between 4 months and 80 years, were residents on Bioko Island. Malaria patients were classified into uncomplicated malaria states according to World Health Organization (WHO) criteria, which were defined as positive smear for P. falciparum and presence of fever (≥ 37.5 °C). Dried blood spots were collected on day zero of enrollment through finger prick bleeding spotted onto Whatman 903® filter paper (GE Healthcare, Pittsburgh, USA) for future use. Laboratory screening for malaria was done using rapid diagnostic test (RDT) and confirmed using microscopic examination of blood smears. For quality control, archived malaria-positive microslides were re-examined, and parasite density was recorded; The Plasmodium species were identified by a real-time PCR followed by high-resolution melting (HRM) [11]. The pGEM-T standard plasmids of the four human Plasmodium species, P. falciparum, Plasmodium ovale, Plasmodium malariae and Plasmodium vivax were used as control. In this study, although only clear non-superimposed signals of the sequences were included from electropherograms for further analysis, these sequences could derive from either monoclonal infections or the most predominant clones in these isolates.

Genomic DNA extraction

Parasite genomic DNA was extracted from dried filter blood spots by Chelex-100 extraction method described in a previous article [12]. The DNA products were collected in sterile tubes and stored at − 80 °C in reserve.

Sequencing analysis of the entire TRAP gene

The full-length TRAP gene (NCBI Gene ID: 814170) was divided into two segments and amplified by nested PCR. The primers designed for nested PCR are presented in Table 1. For the first round PCR, l μl of genomic DNA was amplified with 12.5 μl 2 × Master Mix (DNA Polymerase, dNTP Mixture, PCR buffer), 1 μl 10 nM outer forward primer, 1 μl 10 nM outer reverse primer, and sterile ultra-pure water to a final volume of 25 μl. Thermal cycling parameters for PCR were as follows: initial denaturation at 94 °C for 3 min; 30 cycles of 94 °C for 30 s, annealing at 55 °C for 30 s and extension at 72 °C for 2 min; followed by a final extension step at 72 °C for 5 min. For the second round PCR, 2 μl of the primary PCR product was amplified in a 50 μl reaction volume comprised of 25 μl 2 × Master Mix (DNA Polymerase, dNTP Mixture, PCR buffer), 2 μl 10 nM inner forward primer, 2 μl 10 nM inner reverse primer, and sterile ultrapure water to a final volume of 50 μl. All PCR products were analysed using 1.2% agarose gel electrophoresis, and then, they were purified and sequenced by using an ABI 3730 × L automated sequencer (Shanghai Yingjun Biotechnology Co., LTD, Guangzhou branch). To ensure the accuracy of the sequencing, each isolate was sequenced at least twice. Sequencing primers were the reverse primers of the second round PCR. All the sequences were analysed and integrated by MEGA6 software [13]. These nucleotide sequences have been deposited at NCBI under Accession Numbers (MK981410–MK981530).

Table 1 Primer design information of PfTRAP gene

Sequence polymorphism and natural selection analysis

The nucleotide and deduced amino acid sequences of PfTRAP were analysed using EditSeq and SeqMan in the DNASTAR package (DNASTAR, Madison, WI, USA). The PfTRAP sequence of the laboratory-adapted P. falciparum strain 3D7 (NC_004331.3) was included in the alignment for comparison as a reference sequence. The values of segregating sites (S), the number of haplotypes (H), haplotype diversity (Hd), and nucleotide diversity (π) were calculated using DnaSP version 5.10.00 [14]. In order to test the null hypothesis of neutrality of PfTRAP, the rates of synonymous (dS) and nonsynonymous (dN) substitutions were estimated and were compared using the Z-test (P < 0.05) in MEGA6 program [13] using Nei and Gojobori’s method [15] with the Jukes and Cantor (JC) correction of 1000 bootstrap replications. Tajima’s D test [16] were performed using DnaSP ver. 5.10.00 in order to evaluate the neutral theory of natural selection. The probability of recombination between adjacent sites per generation (Rb), and the minimum number of recombination events (Rm) were analysed using DnaSP ver. 5.10.00.

Global PfTRAP sequences acquisition and comparative analysis

In order to analyse the genetic diversities and natural selection of PfTRAP among global P. falciparum isolates, a total of 1287 sequences from 13 malaria endemic regions (Bangladesh, Cambodia, Congo, Gambia, Ghana, Guinea, Laos, Malawi, Mali, Myanmar, Senegal, Thailand and Vietnam) were acquired by mining the MalariaGEN Pf3k Project database (release 5) [16, 17] using samtools [18] and vcftools [19]. Genetic polymorphism and tests of neutrality were calculated for each population using DnaSP ver. 5.10.00 and MEGA6 as described above. In order to investigate the genetic relationships among global PfTRAP haplotypes, the haplotype network for 1406 full-length sequences of PfTRAP from Bioko and other countries listed above was constructed by Network 10.1.0.0 program using Median-Joining method [20].

Prediction of impact of amino acid change upon protein structure

The 3-dimensional structure of PfTRAP (3D7 isolate) was predicted and modeled by I-TASSER online server (https://zhanglab.ccmb.med.umich.edu/I-TASSER/) [21] and presented using YASARA software [22]. PolyPhen-2 [23] online serve was used to predict potential impact of amino acid substitutions on the structure or function. A prediction of probably damaging indicates that the query substitution is predicted to be damaging with high confidence, while a prediction of benign indicates that the query substitution is predicted to be benign with high confidence. A prediction of possibly damaging indicates that the query substitution is predicted to be damaging, but with low confidence. Using FOLDX plugin [24] in YASARA [22] to predict the changes in free energy before and after the mutations: ΔΔG(change) = ΔG(mutation)  − ΔG(wild-type). As a rule of thumb for use: ΔΔG (change) > 0: the mutation is destabilizing; ΔΔG (change) < 0: the mutation is stabilizing.

Results

Genetic polymorphism and natural selection of PfTRAP on Bioko Island

Of the 153 blood samples extracted from Bioko Island, 121 yielded suitable PfTRAP amplicons for sequencing and deposited to NCBI Genbank (MK981410–MK981530). Finally, 119 monoclonal sequences were applied in the further analysis while 2 polyclonal sequences (MK981439, MK981530) were excluded.

As the result of genetic polymorphism and natural selection analysis shown in Table 2, on Bioko Island, a total of 87 haplotypes were found among 119 samples, with the haplotype diversity (Hd) for 0.99. Furthermore, nucleotide diversity (π) of Bioko Island was detected as 0.01046, which was a relatively high value among the 14 countries in this analysis. As for the parameters related to natural selection, the value of dN–dS was 6.2231 (p < 0.05), which hinted at natural selection on Bioko Island. The value of Tajima’s D of Bioko Island was detected as − 0.41438, which was not significant based on Tajima’s D significance levels [16]. To analyse the recombination degree of PfTRAP on Bioko Island, the minimum number of recombination event (Rm) and the probability of recombination between adjacent sites (Rb) was detected as 22 and 119, respectively, which shows relatively high level among the 14 countries or areas included in this analysis.

Table 2 Result of genetic diversity, natural selection and recombination among global PfTRAP sequences

Global PfTRAP comparative analysis

Not only Bioko Island, but also 7 other African countries and 6 Asian countries were included in the genetic diversity and natural selection analysis of PfTRAP. 1287 global sequences mined from Pf3K database (23 for Bangladesh, 393 for Cambodia, 54 for Congo, 34 for Gambia, 228 for Ghana, 49 for Guinea, 42 for Laos, 106 for Malawi, 39 for Mali, 45 for Myanmar, 110 for Senegal, 104 for Thailand, 60 for Vietnam, shown in Additional file 4) and 119 Bioko sequences were applied in the global comparative analysis. According to the results, a total of 161 SNPs were detected from the 1406 PfTRAP sequences (7 SNPs for DomainI, 41 for DomainII, 23 for DomainIII, 83 for DomainIV, 3 for DomainV and 4 for DomainVI, respectively), among them there were 42 SNPs appeared only once. Due to the uncertainty, these single SNPs were not included in the following analysis (original SNP information was shown in Additional file 1). There were 7 SNPs (S91N, K130R, S179N, L287P, N317K, S419P and A489G) popular globally, with percentage over 80%, and notably, there are 60 SNPs unique to Africa, while there were 12 SNPs unique to Asia. Not surprisingly, no mutations were found in some important motifs including WSPCSVTCG in A-domain, metal ion-dependent adhesion site (MIDAS), IQQ motif, A glycoprotein with a cellular recognition function (RGD) and two important points (aa131 and aa162).

As shown in Table 2, both Hd and π show that the polymorphism in Africa is greater than in Asia (except Bangladesh). All the 14 countries and areas in this analysis shown evidences of natural selection with the positive value of dN–dS (p < 0.05). Except for Bioko, Congo, Guinea, Myanmar, Bangladesh, the Tajima’s D of other 9 countries shown positive value, which signified low levels of both low and high frequency polymorphisms, indicating a decrease in population size and/or balancing selection (Table 2). In the sight of the Tajima’s D of the full-length PfTRAP, on Bioko, the value of DomainII is positive and obviously deviated from 0, while most of the rest region got the negative value (Fig. 1). The overall trend is similar in Bioko and Africa, while in Asia, higher positive Tajima’s D was found in N-terminal of DomainII and C-terminal of DomainIV (Fig. 1). When turns to the recombination effect, the parameters (Rm and Rb) revealed that more frequent recombination events were taking place in Africa rather than Asia (Table 2).

Fig. 1
figure 1

Trend chart of Tajima’s D value of PfTRAP from Africa, Asia and Bioko Island. Red dot indicated p value < 0.05; Yellow dot indicated p value < 0.1; Blue dot indicated p value > 0.1. DomainI to Domain VI were presented in different color

For further exploration of the relationship of haplotypes, a haplotype network was constructed based on global PfTRAP sequences (excluded repeat region) using NETWORK software. A total of 661 haplotypes were detected among 1406 isolates (Fig. 2), among them, 667 Asian isolates clustered in 136 haplotypes and 739 African isolates clustered in 528 haplotypes, which indicated the higher heterogeneity found among African P. falciparum population. Generally, the vast majority of haplotypes were limited to one continent (Africa or Asia) and only several haplotypes (Hap_8, Hap_18, Hap_69, Hap_104 and Hap_370) were shared between two continents isolates (Detailed information about haplotypes was shown in Additional files 2, 3 and 5). In the meantime, it was found that the haplotypes of Bioko Island isolates were scattered. Obviously, haplotypes from Asian isolates were distributed relatively concentrated while the haplotypes from African isolates were in greater dispersion with long branch, which indicated that the recent mutations have been more active in Africa than Asia. In order to analyse the genetic differentiation among PfTRAPs from 14 malaria endemic regions, Fixation Index (FST) was calculated and shown in Table 3. According to Sewall Wright rules, FST range from 0 to 0.05 reflects little genetic differentiation; FST range from 0.05 to 0.15 means moderate genetic differentiation; FST range from 0.15 to 0.25 for great genetic differentiation, while FST over 0.25 means extremely high genetic differentiation. According to Table 3, it seems like the TRAP gene from Bioko Island have not so much genetic differentiation with other 7 African countries in this analysis according to the FST between them shows little or moderate level of genetic differentiation (p < 0.05). The FST between Bioko and Asian countries were over 0.15 (p < 0.05), which indicated a greater genetic differentiation. Generally, the differentiation of PfTRAP among African countries is less pronounced, and the same phenomenon was found among Asian countries, but the high degree of genetic differentiation between African and Asian countries cannot be ignored (p < 0.05).

Fig. 2
figure 2

Haplotype network of global PfTRAP sequences. Each circle stands for one haplotype. The size of circle indicated the sample size of the haplotype. The length of the line connecting two haplotypes indicated the genetic distance. Blue circles represent sequences from Asia; Orange circles represent sequences from Africa; Green circles represent sequences from Bioko Island

Table 3 Result of Fixation Index (FST) of PfTRAP between populations of 14 countries and 3D7 isolate

Effect prediction of mutations located at immune epitopes

In this study, all the proven T cell epitopes, B cell epitopes, as well as nonsynonymous substitutions were marked and presented in Fig. 3. As it shown, 98 amino acid substitutions were located at the B cell epitopes while 25 substitutions were at T cell epitopes. 24 substitutions located at the overlap region of T cell epitope and B cell epitope. Since previous reports shown that cellular immunity is more dominant than humoral immunity in the PfTRAP-based vaccine clinical trials, the following mutation effect prediction analysis was focus on the variations at T cell epitopes. In Table 4, among the 30 mutations in T cell epitopes, 13 of them were predicted as benign and 12 for possibly damaging. Not surprisingly, the high frequency (> 90%) mutations were all predicted as benign. It is worth noting that there were 5 mutations (I116T, L221I, Y128F, G228V and P299S) predicted as probably damaging and interestingly, all these 5 mutations were mostly or totally distributed in Africa region, with the occurrence frequency ranging from 0 to 28% worldwide.

Fig. 3
figure 3

CD8+ T cell and linear B cell Epitopes and amino acids substitution distribution of PfTRAP gene. The full-length TRAP sequences (1–559aa) was referred to the 3D7 isolates (NCBI Gene ID: 814170). Sequences with blue block above is linear B cell epitope; Sequences with orange block above is CD8+ T cell epitope (Detail information of CD8+ T cell epitopes was shown in Additional File 6); Numbers in the color blocks are the corresponding epitope ID of IEDB database. Sequences in gray block is the repeat region. Amino acid in red is for mutant

Table 4 Information of mutations at T cell epitopes and mutation effect prediction

Furthermore, the three dimensions (3D) structure of full-length 3D7 PfTRAP was predicted and modeled by using I-TASSER server and displayed in YASARA application (Fig. 4). Based on the predicted structure of TRAP, the changes in free energy difference before and after the mutations (ΔΔG) were calculated, and the result shows that these mutations L49V, R285G, R285S, P299S and K421N would lead to the destabilization because of the obviously increased free energy difference (ΔΔG > 1). Furthermore, the spatial conformational changes of the twelve relatively common (> 10% globally) mutations at T cell epitopes were predicted and presented in Fig. 5. According to it, 7 amino acid substitutions (E46Q, I116S, K119R, S123N, Y128F, L287P and R290W) were identified as would lead to the breakage of hydrogen bonding, which would further impair the stability of the protein structure after these mutations (Fig. 5).

Fig. 4
figure 4

Predicted 3-dimensional structure of TRAP protein. a Molecular surface of predicted TRAP structure. Yellow for DomainI (Signal sequences); Cyan for DomainII (A domain); Green for Domain III (TSR region); Blue for Domain IV (Proline-rich region, including repeat region); Red for Domain V (transmembrane domain); Orange for Domain VI (cytoplasmic tail domain). b Predicted secondary structure of TRAP protein. Helix, Strand and Coil were shown in different color and form

Fig. 5
figure 5

Structure changes before and after relatively common (> 10% globally) mutations which located in the T cell epitopes. The primary structures of the target amino acids were marked as yellow. The green dotted line indicated hydrogen bonding

Discussion

As PfTRAP was the potential candidate for anti-malarial vaccine and there are several PfTRAP-based vaccines undergoing the clinical trials [25,26,27], the worldwide information of its polymorphism was necessary and important for design and improvement of an effective vaccine. In this study, PfTRAP gene data of Bioko Island, Equatorial Guinea was presented and submitted to public database, which had improved the global malarial database. Overall, with the polymorphism analysis, Bioko PfTRAP exhibited the high polymorphism, which is consistent with PfTRAP from other African countries but significantly distinct with the relatively low polymorphic Asian PfTRAP. These phenomena consistence with the previous report about P. falciparum circumsporozoite protein (PfCSP) gene [28], which also indicated that the similar polymorphic pattern between Bioko and Africa mainland countries. Moreover, the exploration and analysis of other vaccine candidate genes were also analysed previously, including pre-erythrocytic stage CSP gene, erythrocytic phase merozoite surface protein (MSP-1/2) and asexual blood stage apical membrane antigen-1 (AMA-1) gene [28,29,30]. These results shown that the polymorphism of candidate genes associated with malaria vaccines in Africa has been at a high level for a long time.

As for the natural selection analysis of Bioko Island PfTRAP, the results obtained by the two different analytical methods are not consistent. The result of dN–dS is 6.2231 (p < 0.05), which is significantly hinted at a natural selection, while the value of Tajima’s D is − 0.41438 (p > 0.05), which is not significant based on Tajima’s D significance levels. For these two results are different, it cannot be concluded with certainty that the PfTRAP on Bioko island has been affected by natural selection, but continued monitoring is still recommended.

In the sight of the haplotype distribution and relation network, haplotypes of Bioko and other African countries distributed scattered while Asian ones are tended to cluster. This phenomenon is in line with the polymorphism result, which shows more mutations appear in African area and thereby lead to the abundant haplotypes. Moreover, the vast majority of these haplotypes are presented as singletons and the high prevalence of singleton is probably associated with the intensity of transmission or rapid expansion of population. Even so, the genetic differentiation between African countries shows limited and so do Asian countries, but an obvious genetic differentiation is found between countries which come from different continents. It might be explained with the population segregation of Anopheles, the host of P. falciparum, caused by the geographical separation. But interestingly, though Bioko Island is an island separated from African mainland by the Atlantic Ocean, its PfTRAP gene shows low genetic differentiation with other African mainland countries. This result might be explained by the work of Guerra et al. [31, 32], which reported that the strong connection of human movement between Bioko and the mainland Equatorial Guinea (EG), determine a high vulnerability of Bioko to malaria importation; these studies reported that the odds of malaria infection in travellers who had been to mainland EG were more than three times the rest of the population, which confirmed that the majority malaria cases are actively imported by off-island travelers to mainland EG [31, 32]. In general, the non-negligible geographical characteristic might provide a new insight for the development of universal PfTRAP-based vaccine.

Antigen polymorphism has been a major obstacle in the way of developing effective vaccines. Recent studies have highlighted the importance of protective roles of CD8+ T-cell and memory T-cell responses to PfTRAP from clinical malaria cases [33, 34]. Mutations on the surface of the antigen make it more difficult for the host immune system to recognize the antigen, and even lead to immune avoidance and reduce the immune effect. The present study had found a large number of substitutions located at the antigen epitopes, and the destructive prediction was presented. The destructive prediction result in this study reflect the potential destructive effect on the protein structure and further impaction of TRAP function. As it is known that TRAP is a protein that plays an important role in the sporozoite gliding motility, hepatocyte invasion and establishment of infection in the vertebrate host [6,7,8], thus ‘damaging mutations’ might affect these functions but it would not be a threat to survival, which could explain that why the parasites with 'damaging mutations' still alive and be detected. Furthermore, TRAP as a vaccine candidate gene of great potential, the ‘damaging mutations’ probably affect the effectiveness of the vaccine, for the protein surface structure were changed. Several mutations were predicted as damaging (I116T, L122I, Y128F, G228V and P299S), and Y128F is noteworthy with global frequency of 21%. Some ‘possibly damaging’ mutations such as E46Q, I116S and T134S were also found in high frequency globally (24%, 13% and 37%, respectively). As these mutations were in such high frequency and were predicted as damaging to TRAP function or structure, the follow-up continuous monitoring project is worth to pay special attention to.

Nowadays, several malaria vaccines encoding the pre-erythrocytic antigen ME-TRAP, which often coupled with different adjuvants such as Chimpanzee adenovirus (ChAd63) or modified vaccinia Ankara (MVA), had been developed and issued by researchers have entered the clinical trial stage [26, 27]. Still, these effects are not ideal for the goal of developing a globally effective vaccine. Some researchers had put forward a new insight that CSP-TRAP fusion antigens (TRAP N-terminal domains fused to circumsporozoite protein C-terminus with or without repeat region) could produce an effective host immune and it has been validated in vitro cell experiments [35]. In this global analysis, the TRAP N-terminal domains (including Signal Sequences, A-domain and TSR) showed active genetic mutation, especially in A-domain, and there were 8 mutations predicted to do harm to protein structure or function. It presents the suggestion that the application of TRAP N-terminal domains in the vaccine components might deserve more observation and in-depth assessment. Moreover, as the previous research shown, there are numerous epitopes were predicted in the conserved regions that signifies the possibility of the development of a universal TRAP-based malaria vaccine design [36]. As TRAP showed its great potential in inducing immunity, the systematic statistics and analysis of the polymorphism based on the genetic data from global malaria endemic regions shows a more important role in the development of TRAP-related vaccines.

Conclusion

The overall trend of Bioko PfTRAP and African PfTRAP shows no significant difference. The PfTRAPs from the same continent shows more homogeneity while the ones from different continents shows more heterogeneity. Taken together, results of the present investigation showed that global PfTRAP gene is highly polymorphic (especially at T cell epitopes) and some mutations shows destructive to the protein structure or function, which may lead to a change in the affinity of these epitopes to immune molecules. Moreover, evidences in this analysis supported that the global PfTRAP gene is under natural selection. Although TRAP related vaccine clinical trials have not been deployed on Bioko Island, in this paper, polymorphisms and natural selection of PfTRAP are analysed to provide relevant genetic background information for future deployment of clinical trials.

Availability of data and materials

All data generated or analysed during this study are included in this published article [and its additional information files].

Abbreviations

PfTRAP:

Plasmodium falciparum Thrombospondin-related adhesive protein

BIMEP:

Bioko Island Malaria Elimination Project

BIMCP:

Bioko Island Malaria Control Project

EGMVI:

Equatorial Guinea Malaria Vaccine Initiative

SNP:

Single nucleotide polymorphism

PDCs:

Private diagnostic clinics

EIR:

Entomological inoculation rate

RDT:

Rapid diagnostic test

HRM:

High-resolution melting

PCR:

Polymerase chain reaction

Rb:

Recombination between adjacent sites

Rm:

The minimum number of recombination events

MIDAS:

Metal ion-dependent adhesion site

FST:

Fixation Index

AMA-1:

Apical membrane antigen-1

CSP:

Circumsporozoite protein

MSP-1/2:

Merozoite surface proteins

EG:

Equatorial Guinea

References

  1. WHO. World Malaria Report 2019. Geneva: World Health Organization; 2019. https://www.who.int/malaria/publications/world_malaria_report/en/.

  2. Bioko Island Malaria Elimination Project. https://www.mcdinternational.org/bimep.

  3. Bioko Island Malaria Control Project. https://www.mcdinternational.org/bimcp.

  4. Equatorial Guinea Malaria Vaccine Initiative. https://www.mcdinternational.org/egmvi/en/

  5. Muller HM, Scarselli E, Crisanti A. Thrombospondin related anonymous protein (TRAP) of Plasmodium falciparum in parasite-host cell interactions. Parassitologia. 1993;35(Suppl):69–72.

    PubMed  Google Scholar 

  6. Muller HM, Reckmann I, Hollingdale MR, Bujard H, Robson KJ, Crisanti A. Thrombospondin related anonymous protein (TRAP) of Plasmodium falciparum binds specifically to sulfated glycoconjugates and to HepG2 hepatoma cells suggesting a role for this molecule in sporozoite invasion of hepatocytes. EMBO J. 1993;12:2881–9.

    Article  CAS  Google Scholar 

  7. Robson KJ, Frevert U, Reckmann I, Cowan G, Beier J, Scragg IG, et al. Thrombospondin-related adhesive protein (TRAP) of Plasmodium falciparum: expression during sporozoite ontogeny and binding to human hepatocytes. EMBO J. 1995;14:388–94.

    Article  Google Scholar 

  8. Sultan AA, Thathy V, Frevert U, Robson KJ, Crisanti A, Nussenzweig V, et al. TRAP is necessary for gliding motility and infectivity of Plasmodium sporozoites. Cell. 1997;90:511–22.

    Article  CAS  Google Scholar 

  9. Akhouri RR, Bhattacharyya A, Pattnaik P, Malhotra P, Sharma A. Structural and functional dissection of the adhesive domains of Plasmodium falciparum thrombospondin-related anonymous protein (TRAP). Biochem J. 2004;379:815–22.

    Article  CAS  Google Scholar 

  10. Ohashi J, Suzuki Y, Naka I, Hananantachai H, Patarapotikul J. Diversifying selection on the thrombospondin-related adhesive protein (TRAP) gene of Plasmodium falciparum in Thailand. PLoS ONE. 2014;9:e90522.

    Article  Google Scholar 

  11. Wang SQ, Zhou HY, Li Z, Liu YB, Fu XF, Zhu JJ, et al. Quantitative detection and species identificaton of human Plasmodium spp. by using SYBR Green I based real-time PCR. Zhongguo Xue Xi Chong Bing Fang Zhi Za Zhi. 2011;23:677–81 (in Chinese).

    CAS  PubMed  Google Scholar 

  12. Li J, Chen J, Xie D, Eyi UM, Matesa RA, Obono MMO, et al. Molecular mutation profile of Pfcrt and Pfmdr1 in Plasmodium falciparum isolates from Bioko Island. Equatorial Guinea Infect Genet Evol. 2015;36:552–6.

    Article  CAS  Google Scholar 

  13. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.

    Article  CAS  Google Scholar 

  14. Rozas J, Ferrer-Mata A, Sanchez-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;34:3299–302.

    Article  CAS  Google Scholar 

  15. Ina Y. New methods for estimating the numbers of synonymous and nonsynonymous substitutions. J Mol Evol. 1995;40:190–226.

    Article  CAS  Google Scholar 

  16. Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123:585–95.

    Article  CAS  Google Scholar 

  17. Pringle JC, Carpi G, Almagro-Garcia J, Zhu SJ, Kobayashi T, Mulenga M, et al. RTS, S/AS01 malaria vaccine mismatch observed among Plasmodium falciparum isolates from southern and central Africa and globally. Sci Rep. 2018;8:6622.

    Article  Google Scholar 

  18. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25:2078–9.

    Article  Google Scholar 

  19. Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, et al. The variant call format and VCF tools. Bioinformatics. 2011;27:2156–8.

    Article  CAS  Google Scholar 

  20. Bandelt HJ, Forster P, Rohl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999;16:37–48.

    Article  CAS  Google Scholar 

  21. Roy A, Kucukural A, Zhang Y. I-TASSER: a unified platform for automated protein structure and function prediction. Nat Protoc. 2010;5:725–38.

    Article  CAS  Google Scholar 

  22. Krieger E, Vriend G. YASARA View - molecular graphics for all devices—from smartphones to workstations. Bioinformatics. 2014;30:2981–2.

    Article  CAS  Google Scholar 

  23. Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, et al. A method and server for predicting damaging missense mutations. Nat Methods. 2010;7:248–9.

    Article  CAS  Google Scholar 

  24. Schymkowitz J, Borg J, Stricher F, Nys R, Rousseau F, Serrano L. The FoldX web server: an online force field. Nucleic Acids Res. 2005;33:W382-388.

    Article  CAS  Google Scholar 

  25. Moorthy VS, Imoukhuede EB, Milligan P, Bojang K, Keating S, Kaye P, et al. A randomised, double-blind, controlled vaccine efficacy trial of DNA/MVA ME-TRAP against malaria infection in Gambian adults. PLoS Med. 2004;1:e33.

    Article  Google Scholar 

  26. Bliss CM, Drammeh A, Bowyer G, Sanou GS, Jagne YJ, Ouedraogo O, et al. Viral vector malaria vaccines induce high-level T cell and antibody responses in West African children and infants. Mol Ther. 2017;25:547–59.

    Article  CAS  Google Scholar 

  27. Mensah VA, Gueye A, Ndiaye M, Edwards NJ, Wright D, Anagnostou NA, et al. Safety, immunogenicity and efficacy of prime-boost vaccination with ChAd63 and MVA encoding ME-TRAP against Plasmodium falciparum infection in adults in Senegal. PLoS ONE. 2016;11:e0167951.

    Article  Google Scholar 

  28. Huang HY, Liang XY, Lin LY, Chen JT, Ehapo CS, Eyi UM, et al. Genetic polymorphism of Plasmodium falciparum circumsporozoite protein on Bioko Island, Equatorial Guinea and global comparative analysis. Malar J. 2020;19:245.

    Article  CAS  Google Scholar 

  29. Chen JT, Li J, Zha GC, Huang G, Huang ZX, Xie DD, et al. Genetic diversity and allele frequencies of Plasmodium falciparum msp1 and msp2 in parasite isolates from Bioko Island Equatorial Guinea. Malar J. 2018;17:458.

    Article  CAS  Google Scholar 

  30. Wang YN, Lin M, Liang XY, Chen JT, Xie DD, Wang YL, et al. Natural selection and genetic diversity of domain I of Plasmodium falciparum apical membrane antigen-1 on Bioko Island. Malar J. 2019;18:317.

    Article  Google Scholar 

  31. Guerra CA, Citron DT, Garcia GA, Smith DL. Characterising malaria connectivity using malaria indicator survey data. Malar J. 2019;18:440.

    Article  Google Scholar 

  32. Guerra CA, Kang SY, Citron DT, Hergott DEB, Perry M, Smith J, et al. Human mobility patterns and malaria importation on Bioko Island. Nat Commun. 2019;10:2332.

    Article  Google Scholar 

  33. Ewer KJ, O’Hara GA, Duncan CJ, Collins KA, Sheehy SH, Reyes-Sandoval A, et al. Protective CD8+ T-cell immunity to human malaria induced by chimpanzee adenovirus-MVA immunisation. Nat Commun. 2013;4:2836.

    Article  Google Scholar 

  34. Todryk SM, Bejon P, Mwangi T, Plebanski M, Urban B, Marsh K, et al. Correlation of memory T cell responses against TRAP with protection from clinical malaria, and CD4 CD25 high T cells with susceptibility in Kenyans. PLoS ONE. 2008;3:e2027.

    Article  Google Scholar 

  35. Lu C, Song G, Beale K, Yan J, Garst E, Feng J, et al. Design and assessment of TRAP-CSP fusion antigens as effective malaria vaccines. PLoS ONE. 2020;15:e0216260.

    Article  CAS  Google Scholar 

  36. Mehrizi AA, Jafari Zadeh A, Zakeri S, Djadid ND. Population genetic structure analysis of thrombospondin-related adhesive protein (TRAP) as a vaccine candidate antigen in worldwide Plasmodium falciparum isolates. Infect Genet Evol. 2020;80:104197.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

The authors thank the Department of Health of Guangdong Province and the Department of Aid to Foreign Countries of the Ministry of Commerce of the People’s Republic of China for their help. The authors also thank Santiago-m Monte-Nguba for his technical help during the samples collection and diagnosis.

Funding

This study was supported by the Key scientific research projects of Guangdong Provincial Department of Education (Grant No. 2019-GDXK-0031); Guangxi Key Research and Development Foundation (Grant No. 2019JJD140052 and AB18221029); Special technology program of Chaozhou for novel coronavirus infection control (Grant No. CZK20200602 and CZK20200603); Special research project of Guangdong for Novel coronavirus pneumonia epidemic prevention and control (Grant No. 2020KZDZX1146).

Author information

Authors and Affiliations

Authors

Contributions

Field work was performed on Bioko Island, EG. ML and LYL conceived and designed the experiments. JTC, DDX, YLW, CSE and UME contributed the blood sample collection and diagnosis. Laboratory work was conducted at Hanshan Normal University and Chaozhou People’s Hospital Affiliated to Shantou University Medical College, LYL, HYH, XYL, WZC, XZL, YZZ, GCZ, HTM and XYC carried out molecular studies and performed statistical analysis; LYL, HYH, XYL, HGW, WYH and JLW collated data results and making tables and charts. ML, LYL and HHY wrote the draft of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Min Lin.

Ethics declarations

Ethics approval and consent to participate

Participants in the clinical study provided written informed consent before their enrolment, and the study was approved by the institutional ethics committee of Malabo Regional Hospital, Bioko, Equatorial Guinea. All participants received adequate anti-malarial treatment.

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.

Global PfTRAP SNP detail information.

Additional file 2.

Geographical distribution of haplotypes.

Additional file 3.

Haplotype information.

Additional file 4.

1287 global PfTRAP sequences from Pf3K Database.

Additional file 5.

Original file for haplotype network.

Additional file 6.

Detail information of CD8+ T cell epitopes of PfTRAP.

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

Lin, LY., Huang, HY., Liang, XY. et al. Genetic diversity and natural selection on the thrombospondin-related adhesive protein (TRAP) gene of Plasmodium falciparum on Bioko Island, Equatorial Guinea and global comparative analysis. Malar J 20, 124 (2021). https://0-doi-org.brum.beds.ac.uk/10.1186/s12936-021-03664-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12936-021-03664-8

Keywords