Skip to main content

Genetic polymorphism of Plasmodium falciparum circumsporozoite protein on Bioko Island, Equatorial Guinea and global comparative analysis

Abstract

Background

Plasmodium falciparum circumsporozoite protein (PfCSP) is a potential malaria vaccine candidate, but various polymorphisms of the pfcsp gene among global P. falciparum population become the major barrier to the effectiveness of vaccines. This study aimed to investigate the genetic polymorphisms and natural selection of pfcsp in Bioko and the comparison among global P. falciparum population.

Methods

From January 2011 to December 2018, 148 blood samples were collected from P. falciparum infected Bioko patients and 96 monoclonal sequences of them were successfully acquired and analysed with 2200 global pfcsp sequences mined from MalariaGEN Pf3k Database and NCBI.

Results

In Bioko, the N-terminus of pfcsp showed limited genetic variations and the numbers of repetitive sequences (NANP/NVDP) were mainly found as 40 (35%) and 41 (34%) in central region. Most polymorphic characters were found in Th2R/Th3R region, where natural selection (p > 0.05) and recombination occurred. The overall pattern of Bioko pfcsp gene had no obvious deviation from African mainland pfcsp (Fst = 0.00878, p < 0.05). The comparative analysis of Bioko and global pfcsp displayed the various mutation patterns and obvious geographic differentiation among populations from four continents (p < 0.05). The global pfcsp C-terminal sequences were clustered into 138 different haplotypes (H_1 to H_138). Only 3.35% of sequences matched 3D7 strain haplotype (H_1).

Conclusions

The genetic polymorphism phenomena of pfcsp were found universal in Bioko and global isolates and the majority mutations located at T cell epitopes. Global genetic polymorphism and geographical characteristics were recommended to be considered for future improvement of malaria vaccine design.

Background

Malaria, caused by Plasmodium spp. infections, is one of the most significant life-threatening infectious diseases to humans worldwide. According to the World Malaria Report 2019 [1], an estimated 228 million (95% confidence interval [CI] 206–258 million) persons suffered from malaria infections worldwide, with 405,000 malaria deaths in 2018. Twenty countries accounted for 85% of global malaria cases in 2018; all these countries are in sub-Saharan Africa, except for India. Resistance to anti-malarial drugs and insecticides, coupled with the lack of availability of an effective vaccine, is the leading factors behind the parasite’s continuing burden. Apart from its complex life cycle, which alternates between the human and the mosquito host, the malaria parasite also exhibits stages characterized by extensive genetic and antigenic diversity which may present adverse obstacles to anti-malarial control measures.

Currently, there are many efforts and studies have been performed in order to develop effective vaccines, several potential vaccine candidates targeted against pre-erythrocytic, erythrocytic and sexual stages of Plasmodium falciparum are under various stages of clinical development [2, 3]. RTS, S/AS01 vaccine is a pre-erythrocytic stage vaccine based on the P. falciparum circumsporozoite protein (PfCSP) [4, 5]. In 2015, the European Medicines Agency for the immunization of children against malaria approved the RTS, S/AS01 vaccine [6] and the phase 3 clinical trials conducted in various sites in Africa showed that the RTS, S/AS01 vaccine has a protective efficacy of 45% in children in the first twenty months after vaccination [7, 8]. In 2018, the World Health Organization through a large-scale pilot malaria vaccine implementation program (MVIP) aimed to introduce this vaccine in three sub-Saharan countries (Ghana, Kenya, Malawi) [6]. Besides of RTS, S/AS01, a live attenuated Plasmodium falciparum whole sporozoite (SPZ) vaccine is also regarded as a great potential malarial vaccine. Sanaria® PfSPZ Vaccine had conducted a clinical trial on Bioko Island where 70% vaccinees developed antibodies to P. falciparum circumsporozoite protein, which was the first clinical trial conducted in Equatorial Guinea [9]. It is not hard to see that pfcsp is a very important gene for the host immune response to the P. falciparum invasion.

PfCSP is predominantly distributed on the surface of the sporozoites with a molecular mass of about 58 kDa. PfCSP is GPI-anchored on the sporozoite surface and plays a critical role in sporozoite development, motility and hepatocyte invasion [10, 11]. The structure of PfCSP can be divided into three distinct regions: a highly variable central repeat region flanked by a conserved N-terminal region and a C-terminal non-repeat region [12]. The central repeat region, which has been recognized as a major target for antibody-mediated neutralization, is rich in Asn-Ala-Asn-Pro (NANP) tandem repeats and contains a small number of Asn-Val-Asp-Pro (NVDP) motifs [12], constitutes immunodominant B cell epitopes. The C-terminal non-repeat region includes two polymorphic sub-regions, Th2R and Th3R, where T cell epitopes were identified.

The previous studies revealed higher single nucleotide polymorphisms (SNPs) of pfcsp within the P. falciparum population from different geographic regions [13]. Indeed, most P. falciparum vaccine candidate gene including pfcsp have been found to show various genetic and antigenic polymorphisms in global parasites, which might obstruct or reduce the efficacy of vaccines [14, 15].

Understanding the genetic nature of vaccine candidate antigens is critical for designing an effective vaccine. The aims of the present study are to investigate the polymorphism pattern of pfcsp gene and its diversifying selection of P. falciparum on Bioko Island, and to elucidate how pfcsp gene is differentiated among global P. falciparum populations. This study will fill in the blank of Bioko Island pfcsp data, as well as be helpful not only for understanding the molecular evolution of the pfcsp gene in P. falciparum, but also for designing peptide-based vaccines for the PfCSP antigen.

Methods

Study area

The study was carried out in Malabo Regional Hospital and 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 located in the northernmost part of Equatorial Guinea. The island has a population of 334,463 (2015 census), of which approximately 90% live in Malabo (the capital city of Equatorial Guinea) in a humid tropical environment. Malaria due to P. falciparum is a major public health problem on the island [16]. Since the Bioko Island Malaria Control Project (BIMCP) has launched at 2004, the parasite prevalence on Bioko decreased from over 45% prevalence in 2004 to 8.5% in 2016, and the reduction of entomological inoculation rate from more than 1000 before 2004 to 14 in 2015 (www.mcdinternational.org).

Ethical approval

Verbal informed consent was obtained from all participating subjects or their parents, and this study, as well as the consent process, was approved by the Ethics Committee of Malabo Regional Hospital. The Ethical approval letter had been shown as Additional files 1 and 2.

Samples collection

A total of 148 blood spot samples were collected from the patients with uncomplicated malaria during January 2011–December 2018 in Bioko Island. Included patients were residents on Bioko Island aged between 4 months and 80 years. Malaria patients were classified into uncomplicated malaria states according to the 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 tests (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 was identified by a real-time PCR followed by high-resolution melting (HRM) [17]. The pGEM-T standard plasmids of four human Plasmodium species including P. falciparum, Plasmodium ovale, Plasmodium malariae and Plasmodium vivax, which were kindly provided by Dr. Cao (Jiangsu Institute of Parasitic Diseases, Wuxi, Jiangsu Province, China), were used as control.

Genomic DNA extraction

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

Amplification of the entire pfcsp gene

The entire pfcsp gene (NCBI Gene ID: 814364) was amplified by nested PCR. For the first round PCR, 2μl of genomic DNA was amplified with 0.25μl 2× HotStart DNA Polymerase, 2μl dNTP Mixture, 5μl 5× PCR buffer, 1μ1 10 mol/L forward primer (5′-CCGGTCATAAATTCTGAATTATCAA-3′), 1μl 10 mol/L reverse primer (5′-CTACAATTAATCGCAAACGTA-3′), and sterile ultra-pure water to a final volume of 25μl. Thermal cycling parameters for PCR were as follows: initial denaturation at 95 °C for 3 min; 30 cycles of 98 °C for 10 s and 68 °C for 90 s. For the second round PCR, 3μl of the primary PCR product was amplified in a 50μl reaction volume comprised of 0.4μl HotStart DNA Polymerase, 3.2μl dNTP Mixture, 8μl 5 × PCR buffer), 1.6 μl 10 mol/L forward primer (5′-CGTGTAAAAATAAGTAGAAA CCACG-3′), 1.6 μl 10 mol/L reverse primer (5′-GTACAACTCAAACTAAG ATGTGTTC-3′), and sterile ultra-pure water to a final volume of 50μl. PCR procedure was as follows: initial denaturation at 95 °C for 3 min; 30 cycles of 98 °C for 10 s and 68 °C for 90 s. All PCR products were analysed using 1.2% agarose gel electrophoresis, and then, they were purified and sequenced by using ABI 3730×L automated sequencer (Shanghai Yingjun Biotechnology Co., LTD, Guangzhou branch). To ensure the accuracy of the sequencing, at least two clones for each isolate were sequenced. Sequencing primers were the reverse primers of the second round PCR; all the sequences were analysed and integrated by MEGA 6.0 software [19].

Sequences analysis

The pfcsp sequence of the laboratory-adapted P. falciparum strain 3D7 (NCBI Gene ID: 814364) was included in the alignment for comparison as a reference sequence. The values of segregating sites (S), number of Haplotypes (H), haplotype diversity (Hd), and observed average pairwise nucleotide diversity (π) were calculated using DnaSP version 6.12.01 [20]. The π was also calculated on a sliding window plot of 10 bases with a step size of 5 bp in order to estimate the stepwise diversity across the sequences. In order to test the null hypothesis of neutrality of pfcsp, the rates of synonymous (dS) and nonsynonymous (dN) substitutions were estimated and were compared by MEGA 6.0 program using Nei and Gojobori’s method [21] with the Jukes and Cantor (JC) correction of 1000 bootstrap replications. Tajima’s D test [22], Fu and Li’s D and F statistics analysis [23] were performed using DnaSP 6.12.01 in order to evaluate the neutral theory of natural selection (Table 1). The recombination parameter (R), which included the effective population size and probability of recombination between adjacent nucleotides per generation, and the minimum number of recombination events (Rm) were analysed using DnaSP 6.12.01 (Table 1).

Table 1 Genetic diversity and natural selection test and recombination analysis of global PfCSP C-terminus

Sequence acquisition and global analysis

The genetic diversities of pfcsp among global P. falciparum isolates were analysed. A total of the 2200 pfcsp sequences from 24 countries or areas were acquired as follows: (i) 1747 monoclonal sequences of Bangladesh, Cambodia, Congo, Gambia, Ghana, Guinea, Laos, Malawi, Mali, Myanmar, Nigeria, Senegal, Thailand and Vietnam were extracted successfully by mining the MalariaGEN Pf3k Project (release 5) [13] using samtools [24] and vcftools [25]; (ii) 453 sequences of Philippines, Iran, India, Papua New Guinea (PNG), Vanuatu, Solomon Islands, Cameroon, Tanzania, Venezuela and Brazil were obtained from NCBI database (Additional file 3). Genetic polymorphism and tests of neutrality were calculated for each population using DnaSP 6.12.01 and MEGA 6.0 as described above. A logo plot was constructed for each pfcsp population using the WebLogo program (https://weblogo.berkeley.edu/logo.cgi). In order to investigate the genetic relationships among global pfcsp haplotypes, the haplotype network for C-terminal of pfcsp from Bioko and other 24 countries and areas listed above was constructed by Popart program (http://popart.otago.ac.nz) using Median-Joining method [26].

Prediction of impact of amino acid change upon protein structure

The crystallized structure of PfCSP C-terminus, PDBID 3VDK [27] was applied in analysis. PolyPhen-2 [28] and SIFT [29] online serve was used to predict potential impact of amino acid substitutions on the structure or function. Using FOLDX plugin [30] in YASARA [31] to predict the changes in free energy before and after the mutations: ΔΔG(change) = ΔG(mutation) − ΔG(wild-type). As a ‘rule of thumb’: ΔΔG (change) > 0: the mutation is destabilizing; ΔΔG (change) < 0: the mutation is stabilizing.

Results

Amplification of Bioko pfcsp

Of the 148 blood samples extracted from the collections in Bioko Island, 118 yielded suitable pfcsp amplicons for sequencing. Finally, 96 full-length monoclonal pfcsp were analysed in this study and 22 polyclonal pfcsp were excluded. As expected, size variations were observed in the amplified pfcsp sequences. The approximate sizes of amplified products varied from 1.1 to 1.2 kb, which was mainly caused by differences in the number of tandem repeats in the central repeat region. These nucleotide sequences have been deposited at GenBank under Accession Numbers (MN623126–MN623221).

Genetic polymorphisms of N-terminal region of Bioko and global pfcsp

The N-terminal non-repeat region was relatively conserved in Bioko pfcsp. Compared with the 3D7 reference sequence (XM_001351086), five variations were found in pfcsp N-terminal region of Bioko parasites including L5F (2.08%, 2/96), R70K (1.04%, 1/96), D82N (1.04%, 1/96), A98G (24%, 23/96) and a 57 bp (encoding 19 amino acids of 80NNGDNGREGKDEDKRDGNN81) insertion (50%, 48/96). A comparative analysis of the N-terminal non-repeat region in global pfcsp also showed that the region is relatively well-conserved in global parasites. As shown in Fig. 1a, the 19 amino acids length insertion and A98G were two major variations observed in global pfcsp. Almost all Asian and Oceanian countries showed a high frequency of insertion and A98G (ranging from 80 to 100%), but lower in African and American isolates (ranging from 15 to 79%). Meanwhile, some variations showed uneven geographic distributions and in relatively low frequencies. As shown in Fig. 1a, D99G and G100D were only detected from about 50% of Indian and Iranian parasites.

Fig. 1
figure1

Global sequence polymorphism of pfcsp N-terminus and central repeat region. a Mutation type and frequency of N-terminal region of global pfcsp. Mutations with frequency < 10% were not marked. b Frequencies composition of NANP/NVDP repeat numbers in central repeat region among global pfcsp

Genetic polymorphisms of central repeat region of Bioko and global pfcsp

A total of 7 haplotypes of Bioko pfcsp central region was found at amino acid levels (Fig. 1b). The number of NANP/NVDP repeats were analysed and compared among Bioko and global isolates. In Bioko pfcsp, the number of repetitive sequences (NANP/NVDP) were mainly found as 40 (35%, 34/96) and 41 (34%, 33/96). Globally, the number of NANP/NVDP repeat were differed by geographic location. As shown in Fig. 1b, repeat number of majority global isolates in this study were ranging from 40 to 43, while the patterns of Philippines, India and Iran were more polymorphic than others.

Genetic polymorphisms and natural selection of the C-terminal non-repeat region in Bioko and global pfcsp

Nucleotide diversity (π) of the C-terminal non-repeat region was analysed in Bioko and global pfcsp (Fig. 2). Both Th2R (314KHIKEYLNKIQNSL327) and Th3R (352NKPKDELDYAND363) region, the proven T cell epitopes, are in high nucleotide diversity, while the connecting region between Th2R and Th3R was conserved. The pattern of nucleotide diversity in Bioko pfcsp was perfectly matched with other African countries ones. Compared to patterns of Asia, Africa and America, the one of Oceania was in relatively low diversity, especially in Th2R region, which nearly shows no nucleotide diversity (Fig. 2).

Fig. 2
figure2

Global nucleotide diversity (π) of C-terminal (311–363) region. The values of nucleotide diversity (π) were calculated using DNASP version 6.12.01 with the sliding window length of 10 bp and step size of 5 bp

The parameters associated with nucleotide diversity and natural selection were also evaluated on C-terminus non-repeat region (311–363) of Bioko and global pfcsp (Table 1). The average number of nucleotide diversity (K) of Bioko pfcsp was 5.775 and the overall haplotype diversity (Hd) was 0.962 ± 0.008. The estimated value of dN-dS in Bioko pfcsp was found to be 0.0166 (Table 1). For further analysis of natural selection in the C-terminus of Bioko pfcsp, Tajima’s test and Fu and Li’s test were performed and the result was shown in Table 1. Both Tajima’s D (− 0.68556, p > 0.1) and Fu and Li’s F and D (− 1.23926, p > 0.1 and − 1.22255, p > 0.1, respectively) values were found to be negative.

As for globally situation, Hd of African countries were generally higher than others (Hd > 0.9), which verified the higher level of genetic diversity on African pfcsp. The global dN-dS were shown as positive except Nigeria, and global Tajima’s D values were deviation from 0 in different extents. Recombination events were also evaluated among both Bioko and global pfcsp. As shown in Table 1, relative high recombination parameters were shown in all African countries and Philippines, Bangladesh and Venezuela, while lower recombination parameters in other countries.

In terms of amino acid, the mutation types and its frequencies in C-terminus (311–363) were briefly presented in Fig. 3. There were totally 26 logos generated, one for 3D7 reference isolate and 25 for isolates from different countries and areas. As for Bioko pfcsp, mutations were detected at twelve positions (314, 317, 318, 321, 322, 324, 327, 352, 356, 357, 359, 361). All these positions were situated at two T-cell epitopes (Th2R and Th3R). The overall pattern of Bioko is similar to those of African countries. Relatively, more kinds of mutations existed in African isolates, as well as in Philippine and Venezuelan isolates. In contrast, the Oceanian mutation patterns were tended to more uncomplicated. Rare mutation L320I was only found in Philippines while S326A was only found in Venezuela. The high frequency mutation, A361E, existed in all 25 countries, while its wild type (A361) was mainly found in Africa. Notably, the wild type residues of 317, 318, and 321 positions were rarely seen in global isolates, instead, K317E, E318K, E318Q, N321K were mainly found in these positions (Fig. 3).

Fig. 3
figure3

Non-synonymous mutation in C-terminus (311–363) of global pfcsp. Each logo consists of stacks of symbols, one stack for each position in the sequence. The height of the amino acid abbreviation indicates its relative frequency at that position. As for the pattern of 3D7 isolate, all positions were marked in coloured; as for global pfcsp patterns, mutation sites are marked in colour, while the conserved sites are in gray

Mutation distribution and C-terminus point mutation effect prediction

By analysing with global data, a total of 66 amino acid substitutions were found in the full-length pfcsp sequences. In order to know about the distribution of T cell epitopes of pfcsp, the proven epitopes (CD8+ and CD4+) were searched from IEDB database [32,33,34,35,36,37,38,39].

As shown in Fig. 4, 54 mutations were distributed in T cell epitopes. Majority mutations (74%) were located at the C-terminus of pfcsp, as well as the CD8+ T cell epitopes. Notably, there were 28 variances found in the TSR region (including Th2R and Th3R), which also is the overlap of CD4+ and CD8+ T cell epitopes. Furthermore, mutation effect prediction was conducted among these 28 variances. As shown in Table 2, the mutations K322I, N325Y and S326A were predicted to be deleterious using SIFT program (SIFT < 0.05). According to Humdiv score predicted by PolyPhen 2.0 program, 13 mutants were predicted as benign, 4 mutants were possibly damaging and 11 for probably damaging. Among these probably damaging mutants, the protein structures of K317T, K317A, L327I, N352G, P354S and A361I were tending to destabilize (ΔΔG > 0). Some high frequency mutations such as K317E (84.32%), N321K (84.76%) and A361E (72.43%), were predicted as benign. Some extremely low frequency but predicted damaging mutations like K317A (0.17%), S326A (0.09%), G349D (0.13%) and D356G (0.09%), were lack of persuasion (Table 2).

Fig. 4
figure4

Mutations distribution and T cell epitopes map of pfcsp (3D7 isolate). Capital letters in black are amino acid sequences of 3D7 isolate; The red capital letters under the black ones are for mutants. Sequences with black solid line below indicated CD8+ T cell epitopes, sequences with blue dotted line above indicated CD4+ T cell epitopes. Repeat region is in gray shadow; Th2R region is in orange shadow; Th3R region is in green shadow

Table 2 Global PfCSP C-terminal Mutant types and effect prediction

Population differentiation analysis of pfcsp C-terminus among global P. falciparum isolates

A haplotype network was constructed using 96 samples from Bioko in addition to 2200 global pfcsp C-terminal monoclonal sequences mining from the Pf3k database and NCBI (Fig. 5). The 2296 pfcsp C-terminal sequences were clustered into 138 unique haplotypes (H_1 to H_138). Detailed information of haplotypes was presented in Additional file 4. Fifty-eight haplotypes were shared by pfcsp sequences from at least two different countries; 70 haplotypes were limited to singleton (only composed by 1 sequence). And as for the H_1, which belongs to the 3D7 standard isolate, as well as the component of RTS,S malaria vaccine, only hold 2.08% (2/96) in Bioko isolates and 3.35% (77/2296) in the worldwide isolates, among which 74 isolates were found in Africa. Only H_62 was composed of samples from four continents (Africa, Asia, America and Oceania) but in a low prevalence (24/2296). Interestingly, the isolates from Africa and America shared the same haplotypes or the related ones (H_54, H_131), while the haplotypes of Oceanian isolates (H_35, H_134) have closer relationship with Asian’s. These phenomena correspond to the Fst index results shown in Table 3. As the Table 3 shown, Fst between Bioko Island and African mainland showed no significant population differentiation (Fst = 0.00878, p < 0.05). Meanwhile, clear population differentiation was identified between American, Asian, Oceanian and African parasite population (p < 0.05). Relatively closer genetic relationships were found in African & American parasite population and Asian & Oceanian parasite population (Fst = 0.19194, p < 0.05 and Fst = 0.06564, p < 0.05, respectively).

Fig. 5
figure5

Haplotype network of C-terminal region among global pfcsp. Isolates from four continents and Bioko Island were marked in five different series colours, blue series for Africa, red series for Asia, khaki series for Oceania, green series for America, and yellow for Bioko Island

Table 3 Population pairwise fixation index (Fst) result

Discussion

Bioko Island, Equatorial Guinea, is a historically high malaria transmission region [16, 40]. Though BIMCP had launched in Bioko Island since 2004 and achieved a remarkable result, malaria is still a major health problem in this region. The genetic diversity and natural selection were analysed in Bioko pfcsp and global pfcsp. In general, the polymorphism patterns between Bioko pfcsp and African mainland pfcsp have no obvious differentiation, although the geographic location of Bioko Island was relatively isolated. This result might be explained by the work of Guerra et al., 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 travellers to mainland EG [41, 42]. Furthermore, it is worth mention that the PfSPZ vaccine had been tested in Malabo and a series of clinical trials are undergoing, which might likely to affect the genetic background of the malaria parasites in this region [9]. According to the report [9], PfSPZ vaccine could induced the immune response to PfCSP, which might influence the genetic diversity and natural selection of pfcsp in Malabo. The natural selection analysis revealed that Bioko pfcsp might under a selection effect although there is no statistical significance (p > 0.1). These findings were in line with the prior studies about P. falciparum merozoite surface protein-1/2 (PfMSP-1/2) and P. falciparum apical membrane antigen-1 (PfAMA-1) genes in Bioko Island [43, 44].

N-terminal region of PfCSP plays an important role in the procedure of sporozoite invades to the hepatocytes [45]. In Bioko and global pfcsp, the genetic polymorphism of N-terminus was in a relatively low level. 19 amino acids length insertion and A98G were universally popular while several novel mutations were found with low frequency. Some scientists verified previously that the antibodies against to N-terminal region could be produced by host immune system and could evoke a partial inhibition of sporozoite invasion of hepatocytes in vitro [46]. Now the evidences of relatively conservative N-terminus might raise the possibility that whether the N-terminus has the potential to be a component of anti-malarial vaccine.

Central repeat region is an immunodominant epitope of PfCSP, and it had been applied to the component of RTS,S malaria vaccine [47]. Different numbers of tetrapeptide repeat was an important cause of pfcsp polymorphism. As expected, this study revealed the diversity of the number of tetrapeptide repeat (NANP/NVNP). Through the analysis among global different geographic regions, it was found that majority of samples possessed the tetrapeptide repeat ranging from 39 to 44 times. Though some scientists hold the view that the various number of tetrapeptide repeat make no significant impacts on RTS,S vaccine efficacy [14], it was known to correlated with the stability of CS protein structure [48]. However, the mechanism and effect of this variation is still unclear. For the universality of this variation, deeper research towards to this region is still necessary.

In the analysis of C-terminus of pfcsp, there were abundant polymorphisms found, especially in the TSR region (including Th2R and Th3R), the proven T cell immunogenic epitopes. The C-terminus of African, Asian, American and Oceanian samples presented their own distinctive diversity patterns. Not surprisingly, more polymorphisms were performed in the two larger-size parasite population (African and Asian) compared to those of America and Oceania. Because of the geographical isolation effect, some mutations showed the regional difference, for example the mutant at 325 position (N325Y) was only occurred in Asian countries; S326A was only found in Venezuela; wild type A361 was mainly observed in Africa, and so on. These phenomena indicated us that continuous monitor to these regional characteristic mutations, and exploration on their association with regional malaria epidemic situation are necessary.

In terms of C-terminal haplotypes analysis, 29 of 34 Bioko pfcsp haplotypes were shared with African continent samples while only 5 were limited to singleton, which implied that Bioko pfcsp was not completely independent of African continent. An obvious phenomenon was found that haplotypes from Oceanian pfcsp have closer genetic relationship with Asian haplotypes. Additionally, the same phenomenon happened among the parasites from America and Africa. It reflects that worldwide genotype of pfcsp C-terminus might divide into two major groups (Africa & America and Asia & Oceania), which probably caused by the frequent communication due to geographical advantages. It provides an insight of the vaccine design based on PfCSP that the regional differentiation might be took into consideration.

The absence of 3D7-matched pfcsp was not the uncommon finding anymore [13, 49]. Unsurprisingly, in Bioko Island, only 2% 3D7-matched pfcsp were found. A study about genetic diversity and protective efficacy of the RTS,S/AS01 malaria vaccine stated that the 3D7-mismatched malaria might probably weaken the efficacy of vaccine, especially the mutations at 299, 301, 317, 354, 356, 359 and 361 amino acid position [14]. In this research, the polymorphism situation of these loci showed different degrees. It is worth mentioning that mutation rate of position 317 reached 91% and mutation rate of position 361 reached 73%. As these mutations are so common and probably affect the vaccine effect, a question raised that whether these high-frequency alleles instead of the wild-type ones could be applied in the vaccine component.

In terms of the distribution of mutations, all the 66 mutations found from global sequences were located at CD8+ T cell epitopes, while 28 of them were located at the overlap of CD8 + and CD4+ T cell epitopes. It is well known that CD8+ and CD4+ T cell are thought to play a role in natural and sporozoite vaccine induced immunity in P. falciparum malaria [50, 51]. This raises the question of whether these mutations affect host immunity. Mutation-effect prediction of these 28 mutations showed that more than half of them were predicted as damaging (15 of 28). Notably, when mutations located at some specific positions (including the probably harmful position 317 and 354) [14], great changes have taken place on the free energy difference, which would result in destabilization on CS protein structure in difference extent. However, the specific mechanism of whether and how these mutations do harm to vaccine efficacy are still not clear. Therefore, continuous monitoring on these mutations and deeper exploration on the mechanism is still necessary.

According to this study, there are several new insights might be considered in the design and improvement of PfCSP-based vaccine: (1) The globally high frequency alleles instead of the wild-type ones of C-terminus might be used for composing vaccine. (2) The immunogenic and conservative N-terminus might be applied in the composition of vaccine. (3) The regional differences should be considered in the improvement of universal malaria vaccine, mainly divided as Asia-Oceania region and Africa-America region.

Conclusions

In this study, the genetic diversity of Bioko and global pfcsp was analysed. The genetic polymorphism of pfcsp was found to be universal. Besides this, significant geographical differentiation of pfcsp were found around the world, which could mainly be divided into Asia & Oceania group and Africa & America group. Meanwhile, the 3D7 isolate was rare to found worldwide. Some mutations which are located at T-cell epitopes might impair the PfCSP-based vaccine efficacy by using prediction tools. Findings in this study filled in missing data of Bioko pfcsp. A holistic view of global pfcsp polymorphism was presented in the article and provides more insight for the improvement of malaria vaccine design.

Availability of supporting data

The datasets supporting the conclusions of this article are included with in the article.

Abbreviations

BIMCP:

Bioko Island Malaria Control Project

dS:

Synonymous

dN:

Nonsynonymous

EG:

Equatorial Guinea

HRM:

High-resolution melting

PfCSP:

Plasmodium falciparum circumsporozoite protein

PfMSP-1/2:

Plasmodium falciparum merozoite surface protein-1/2

PfAMA-1:

Plasmodium falciparum apical membrane antigen-1

PNG:

Papua New Guinea

RDT:

Rapid diagnostic tests

R:

Recombination parameter

Rm:

Minimum number of recombination events

SNPs:

Single nucleotide polymorphisms

References

  1. 1.

    WHO. World malaria report. Geneva, World Health Orgnization, 2019. [https://www.who.int/news-room/feature-stories/detail/world-malaria-report-2019].

  2. 2.

    Mac-Daniel L, Menard R. Live vaccines against Plasmodium preerythrocytic stages. Methods Mol Biol. 2019;2013:189–98.

    CAS  PubMed  Google Scholar 

  3. 3.

    Draper SJ, Sack BK, King CR, Nielsen CM, Rayner JC, Higgins MK, et al. Malaria vaccines: recent advances and new horizons. Cell Host Microbe. 2018;24:43–56.

    CAS  PubMed  PubMed Central  Google Scholar 

  4. 4.

    Casares S, Brumeanu TD, Richie TL. The RTS, S malaria vaccine. Vaccine. 2010;28:4880–94.

    CAS  PubMed  Google Scholar 

  5. 5.

    Zeeshan M, Alam MT, Vinayak S, Bora H, Tyagi RK, Alam MS, et al. Genetic variation in the Plasmodium falciparum circumsporozoite protein in India and its relevance to RTS, S malaria vaccine. PLoS ONE. 2012;7:e43430.

    CAS  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Dimala CA, Kika BT, Kadia BM, Blencowe H. Current challenges and proposed solutions to the effective implementation of the RTS, S/AS01 malaria vaccine program in sub-Saharan Africa: a systematic review. PLoS ONE. 2018;13:e0209744.

    CAS  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Rts SCTP. Efficacy and safety of the RTS, S/AS01 malaria vaccine during 18 months after vaccination: a phase 3 randomized, controlled trial in children and young infants at 11 African sites. PLoS Med. 2014;11:e1001685.

    Google Scholar 

  8. 8.

    Olotu A, Fegan G, Wambua J, Nyangweso G, Leach A, Lievens M, et al. Seven-year efficacy of RTS, S/AS01 malaria vaccine among young African children. N Engl J Med. 2016;374:2519–29.

    CAS  PubMed  PubMed Central  Google Scholar 

  9. 9.

    Olotu A, Urbano V, Hamad A, Eka M, Chemba M, Nyakarungu E, et al. Advancing global health through development and clinical trials partnerships: a randomized, placebo-controlled, double-blind assessment of safety, tolerability, and immunogenicity of PfSPZ vaccine for malaria in healthy Equatoguinean men. Am J Trop Med Hyg. 2018;98:308–18.

    PubMed  Google Scholar 

  10. 10.

    Nussenzweig V, Nussenzweig RS. Circumsporozoite proteins of malaria parasites. Cell. 1985;42:401–3.

    CAS  PubMed  Google Scholar 

  11. 11.

    Plassmeyer ML, Reiter K, Shimp RL Jr, Kotova S, Smith PD, et al. Structure of the Plasmodium falciparum circumsporozoite protein, a leading malaria vaccine candidate. J Biol Chem. 2009;284:26951–63.

    CAS  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Enea V, Ellis J, Zavala F, Arnot DE, Asavanich A, Masuda A, et al. DNA cloning of Plasmodium falciparum circumsporozoite gene: amino acid sequence of repetitive epitope. Science. 1984;225:628–30.

    CAS  PubMed  Google Scholar 

  13. 13.

    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.

    PubMed  PubMed Central  Google Scholar 

  14. 14.

    Neafsey DE, Juraska M, Bedford T, Benkeser D, Valim C, Griggs A, et al. Genetic diversity and protective efficacy of the RTS, S/AS01 malaria vaccine. N Engl J Med. 2015;373:2025–37.

    CAS  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Dobano C, Ubillos I, Jairoce C, Gyan B, Vidal M, Jimenez A, et al. RTS, S/AS01E immunization increases antibody responses to vaccine-unrelated Plasmodium falciparum antigens associated with protection against clinical malaria in African children: a case–control study. BMC Med. 2019;17:157.

    PubMed  PubMed Central  Google Scholar 

  16. 16.

    Cook J, Hergott D, Phiri W, Rivas MR, Bradley J, Segura L, et al. Trends in parasite prevalence following 13 years of malaria interventions on Bioko island, Equatorial Guinea: 2004–2016. Malar J. 2018;17:62.

    PubMed  PubMed Central  Google Scholar 

  17. 17.

    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(in Chinese). Zhongguo Xue Xi Chong Bing Fang Zhi Za Zhi. 2011;23:677–81.

    CAS  PubMed  Google Scholar 

  18. 18.

    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.

    CAS  PubMed  Google Scholar 

  19. 19.

    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.

    CAS  PubMed  PubMed Central  Google Scholar 

  20. 20.

    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.

    CAS  PubMed  Google Scholar 

  21. 21.

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

    CAS  PubMed  Google Scholar 

  22. 22.

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

    CAS  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Fu YX, Li WH. Statistical tests of neutrality of mutations. Genetics. 1993;133:693–709.

    CAS  PubMed  PubMed Central  Google Scholar 

  24. 24.

    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.

    PubMed  PubMed Central  Google Scholar 

  25. 25.

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

    CAS  PubMed  PubMed Central  Google Scholar 

  26. 26.

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

    CAS  PubMed  Google Scholar 

  27. 27.

    Doud MB, Koksal AC, Mi LZ, Song G, Lu C, Springer TA. Unexpected fold in the circumsporozoite protein target of malaria vaccines. Proc Natl Acad Sci USA. 2012;109:7817–22.

    CAS  PubMed  Google Scholar 

  28. 28.

    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.

    CAS  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Sim NL, Kumar P, Hu J, Henikoff S, Schneider G, Ng PC. SIFT web server: predicting effects of amino acid substitutions on proteins. Nucleic Acids Res. 2012;40:W452–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  30. 30.

    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–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  31. 31.

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

    CAS  PubMed  PubMed Central  Google Scholar 

  32. 32.

    Kaba SA, McCoy ME, Doll TA, Brando C, Guo Q, Dasgupta D, et al. Protective antibody and CD8 + T-cell responses to the Plasmodium falciparum circumsporozoite protein induced by a nanoparticle vaccine. PLoS ONE. 2012;7:e48304.

    CAS  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Kumar A, Kumar S, Le TP, Southwood S, Sidney J, Cohen J, et al. HLA-A*01-restricted cytotoxic T-lymphocyte epitope from the Plasmodium falciparum circumsporozoite protein. Infect Immun. 2001;69:2766–71.

    CAS  PubMed  PubMed Central  Google Scholar 

  34. 34.

    Malik A, Egan JE, Houghten RA, Sadoff JC, Hoffman SL. Human cytotoxic T lymphocytes against the Plasmodium falciparum circumsporozoite protein. Proc Natl Acad Sci USA. 1991;88:3300–4.

    CAS  PubMed  Google Scholar 

  35. 35.

    Pinder M, Reece WH, Plebanski M, Akinwunmi P, Flanagan KL, Lee EA, et al. Cellular immunity induced by the recombinant Plasmodium falciparum malaria vaccine, RTS, S/AS02, in semi-immune adults in The Gambia. Clin Exp Immunol. 2004;135:286–93.

    CAS  PubMed  PubMed Central  Google Scholar 

  36. 36.

    Sedegah M, Kim Y, Ganeshan H, Huang J, Belmonte M, Abot E, et al. Identification of minimal human MHC-restricted CD8+ T-cell epitopes within the Plasmodium falciparum circumsporozoite protein (CSP). Malar J. 2013;12:185.

    CAS  PubMed  PubMed Central  Google Scholar 

  37. 37.

    Vita R, Mahajan S, Overton JA, Dhanda SK, Martini S, Cantrell JR, et al. The immune epitope database (IEDB): 2018 update. Nucleic Acids Res. 2019;47:D339–43.

    CAS  PubMed  Google Scholar 

  38. 38.

    Wang R, Doolan DL, Le TP, Hedstrom RC, Coonan KM, Charoenvit Y, et al. Induction of antigen-specific cytotoxic T lymphocytes in humans by a malaria DNA vaccine. Science. 1998;282:476–80.

    CAS  PubMed  Google Scholar 

  39. 39.

    Wang R, Epstein J, Baraceros FM, Gorak EJ, Charoenvit Y, Carucci DJ, et al. Induction of CD4(+) T cell-dependent CD8(+) type 1 responses in humans by a malaria DNA vaccine. Proc Natl Acad Sci USA. 2001;98:10817–22.

    CAS  PubMed  Google Scholar 

  40. 40.

    Cano J, Berzosa PJ, Roche J, Rubio JM, Moyano E, Guerra-Neira A, et al. Malaria vectors in the Bioko Island (Equatorial Guinea): estimation of vector dynamics and transmission intensities. J Med Entomol. 2004;41:158–61.

    CAS  PubMed  Google Scholar 

  41. 41.

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

    PubMed  PubMed Central  Google Scholar 

  42. 42.

    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.

    PubMed  PubMed Central  Google Scholar 

  43. 43.

    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.

    PubMed  PubMed Central  Google Scholar 

  44. 44.

    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.

    CAS  PubMed  PubMed Central  Google Scholar 

  45. 45.

    Rathore D, Sacci JB, deVega P, McCutchan TF. Binding and invasion of liver cells by Plasmodium falciparum sporozoites. Essential involvement of the amino terminus of circumsporozoite protein. J Biol Chem. 2002;277:7092–8.

    CAS  PubMed  Google Scholar 

  46. 46.

    Bongfen SE, Ntsama PM, Offner S, Smith T, Felger I, Tanner M, et al. The N-terminal domain of Plasmodium falciparum circumsporozoite protein represents a target of protective immunity. Vaccine. 2009;27:328–35.

    CAS  PubMed  Google Scholar 

  47. 47.

    Gordon DM, McGovern TW, Krzych U, Cohen JC, Schneider I, LaChance R, et al. Safety, immunogenicity, and efficacy of a recombinantly produced Plasmodium falciparum circumsporozoite protein-hepatitis B surface antigen subunit vaccine. J Infect Dis. 1995;171:1576–85.

    CAS  PubMed  Google Scholar 

  48. 48.

    Escalante AA, Grebert HM, Isea R, Goldman IF, Basco L, Magris M, et al. A study of genetic diversity in the gene encoding the circumsporozoite protein (CSP) of Plasmodium falciparum from different transmission areas–XVI. Asembo Bay cohort project. Mol Biochem Parasitol. 2002;125:83–90.

    CAS  PubMed  Google Scholar 

  49. 49.

    Le HG, Kang JM, Moe M, Jun H, Thai TL, Lee J, et al. Genetic polymorphism and natural selection of circumsporozoite surface protein in Plasmodium falciparum field isolates from Myanmar. Malar J. 2018;17:361.

    CAS  PubMed  PubMed Central  Google Scholar 

  50. 50.

    Rathore D, McCutchan TF. The cytotoxic T-lymphocyte epitope of the Plasmodium falciparum circumsporozoite protein also modulates the efficiency of receptor-ligand interaction with hepatocytes. Infect Immun. 2000;68:740–3.

    CAS  PubMed  PubMed Central  Google Scholar 

  51. 51.

    Kurup SP, Butler NS, Harty JT. T cell-mediated immunity to malaria. Nat Rev Immunol. 2019;19:457–71.

    CAS  PubMed  PubMed Central  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 Natural Science Foundation of Guangdong Province to Jiang-Tao Chen (Grant No. 2016A03031311).

Author information

Affiliations

Authors

Contributions

Field work was performed on Bioko Island, EG. ML and JTC 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, HYH, XYL, LYL, WZC, XZL, YZZ, GCZ HTM and XYC carried out molecular studies and performed statistical analyses; HHY and ML wrote the draft of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Min Lin.

Ethics declarations

Ethical 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.

Ethical approval letter (Spanish version)

Additional file 2.

Ethical approval letter (Chinese version).

Additional file 3.

Global pfcsp sequences acquired from NCBI.

Additional file 4.

Detail information of haplotypes.

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

Verify currency and authenticity via CrossMark

Cite this article

Huang, H., Liang, X., Lin, L. et al. Genetic polymorphism of Plasmodium falciparum circumsporozoite protein on Bioko Island, Equatorial Guinea and global comparative analysis. Malar J 19, 245 (2020). https://doi.org/10.1186/s12936-020-03315-4

Download citation

Keywords

  • Malaria
  • Plasmodium falciparum
  • Circumsporozoite protein
  • Genetic polymorphism
  • Bioko Island