Plasmodium falciparum: genetic diversity and complexity of infections in an isolated village in western Thailand

Strong isolation by distance and evidence of population microstructure reflect ongoing Plasmodium falciparum transmission in Zanzibar

Background

The Zanzibar archipelago of Tanzania has become a low-transmission area for Plasmodium falciparum. Despite being considered an area of pre-elimination for years, achieving elimination has been difficult, likely due to a combination of imported infections from mainland Tanzania and continued local transmission.

Methods

To shed light on these sources of transmission, we applied highly multiplexed genotyping utilizing molecular inversion probes to characterize the genetic relatedness of 282 P. falciparum isolates collected across Zanzibar and in Bagamoyo district on the coastal mainland from 2016 to 2018.

Results

Overall, parasite populations on the coastal mainland and Zanzibar archipelago remain highly related. However, parasite isolates from Zanzibar exhibit population microstructure due to the rapid decay of parasite relatedness over very short distances. This, along with highly related pairs within shehias, suggests ongoing low-level local transmission. We also identified highly related parasites across shehias that reflect human mobility on the main island of Unguja and identified a cluster of highly related parasites, suggestive of an outbreak, in the Micheweni district on Pemba island. Parasites in asymptomatic infections demonstrated higher complexity of infection than those in symptomatic infections, but have similar core genomes.

Conclusions

Our data support importation as a main source of genetic diversity and contribution to the parasite population in Zanzibar, but they also show local outbreak clusters where targeted interventions are essential to block local transmission. These results highlight the need for preventive measures against imported malaria and enhanced control measures in areas that remain receptive to malaria reemergence due to susceptible hosts and competent vectors.

Funding

This research was funded by the National Institutes of Health, grants R01AI121558, R01AI137395, R01AI155730, F30AI143172, and K24AI134990. Funding was also contributed from the Swedish Research Council, Erling-Persson Family Foundation, and the Yang Fund. RV acknowledges funding from the MRC Centre for Global Infectious Disease Analysis (reference MR/R015600/1), jointly funded by the UK Medical Research Council (MRC) and the UK Foreign, Commonwealth & Development Office (FCDO), under the MRC/FCDO Concordat agreement and is also part of the EDCTP2 program supported by the European Union. RV also acknowledges funding by Community Jameel.

Plasmodium falciparum merozoite surface protein 1 as asexual blood stage malaria vaccine candidate

INTRODUCTION

Malaria represents a public health challenge in tropical and subtropical regions, and currently deployed control strategies are likely insufficient to drive elimination of malaria. Development and improvement of malaria vaccines might be key to reduce disease burden. Vaccines targeting asexual blood stages of the parasite have shown limited efficacy when studied in human trials conducted over the past decades.

AREAS COVERED

Vaccine candidates based on the merozoite surface protein 1 (MSP1) were initially envisioned as one of the most promising approaches to provide immune protection against asexual blood-stage malaria. Successful immunization studies in monkey involved the use of the full-length MSP1 (MSP1FL) as vaccine construct. Vaccines using MSP1FL for immunization have the potential benefit of including numerous conserved B-cell and T-cell epitopes. This could result in improved parasite strain-transcending, protective immunity in the field. We review outcomes of clinical trials that utilized a variety of MSP1 constructs and formulations, including MSP1FL, either alone or in combination with other antigens, in both animal models and humans.

EXPERT OPINION

Novel approaches to analyze breadth and magnitude of effector functions of MSP1-targeting antibodies in volunteers undergoing experimental vaccination and controlled human malaria infection will help to define correlates of protective immunity.

Structural organization and sequence diversity of the complete nucleotide sequence encoding the Plasmodium malariae merozoite surface protein-1

The merozoite surface protein-1 (MSP1) is a prime candidate for an asexual blood stage vaccine against malaria. However, polymorphism in this antigen could compromise the vaccine’s efficacy. Although the extent of sequence variation in MSP1 has been analyzed from various Plasmodium species, little is known about structural organization and diversity of this locus in Plasmodium malariae (PmMSP1). Herein, we have shown that PmMSP1 contained five conserved and four variable blocks based on analysis of the complete coding sequences. Variable blocks were characterized by short insertion and deletion variants (block II), polymorphic nonrepeat sequences (block IV), complex repeat structure with size variation (block VI) and degenerate octapeptide repeats (block VIII). Like other malarial MSP1s, evidences of intragenic recombination have been found in PmMSP1. The rate of nonsynonymous nucleotide substitutions significantly exceeded that of synonymous nucleotide substitutions in block IV, suggesting positive selection in this region. Codon-based analysis of deviation from neutrality has identified a codon under purifying selection located in close proximity to the homologous region of the 38 kDa/42 kDa cleavage site of P. falciparum MSP1. A number of predicted linear B-cell epitopes were identified across both conserved and variable blocks of the protein. However, polymorphism in repeat-containing blocks resulted in alteration of the predicted linear B-cell epitope scores across variants. Although a number of predicted HLA-class II-binding peptides were identified in PmMSP1, all variants of block IV seemed not to be recognized by common HLA-class II alleles among Thai population, suggesting that diversity in this positive selection region could probably affect host immune recognition. The data on structural diversity in PmMSP1 could be useful for further studies such as vaccine development and strain characterization of this neglected malaria parasite.

Host-Malaria Parasite Interactions and Impacts on Mutual Evolution

Malaria is the most deadly parasitic disease, affecting hundreds of millions of people worldwide. Malaria parasites have been associated with their hosts for millions of years. During the long history of host-parasite co-evolution, both parasites and hosts have applied pressure on each other through complex host-parasite molecular interactions. Whereas the hosts activate various immune mechanisms to remove parasites during an infection, the parasites attempt to evade host immunity by diversifying their genome and switching expression of targets of the host immune system. Human intervention to control the disease such as antimalarial drugs and vaccination can greatly alter parasite population dynamics and evolution, particularly the massive applications of antimalarial drugs in recent human history. Vaccination is likely the best method to prevent the disease; however, a partially protective vaccine may have unwanted consequences that require further investigation. Studies of host-parasite interactions and co-evolution will provide important information for designing safe and effective vaccines and for preventing drug resistance. In this essay, we will discuss some interesting molecules involved in host-parasite interactions, including important parasite antigens. We also discuss subjects relevant to drug and vaccine development and some approaches for studying host-parasite interactions.

Microarray analyses reveal strain-specific antibody responses to Plasmodium falciparum apical membrane antigen 1 variants following natural infection and vaccination

Vaccines based on Plasmodium falciparum apical membrane antigen 1 (AMA1) have failed due to extensive polymorphism in AMA1. To assess the strain-specificity of antibody responses to malaria infection and AMA1 vaccination, we designed protein and peptide microarrays representing hundreds of unique AMA1 variants. Following clinical malaria episodes, children had short-lived, sequence-independent increases in average whole-protein seroreactivity, as well as strain-specific responses to peptides representing diverse epitopes. Vaccination resulted in dramatically increased seroreactivity to all 263 AMA1 whole-protein variants. High-density peptide analysis revealed that vaccinated children had increases in seroreactivity to four distinct epitopes that exceeded responses to natural infection. A single amino acid change was critical to seroreactivity to peptides in a region of AMA1 associated with strain-specific vaccine efficacy. Antibody measurements using whole antigens may be biased towards conserved, immunodominant epitopes. Peptide microarrays may help to identify immunogenic epitopes, define correlates of vaccine protection, and measure strain-specific vaccine-induced antibodies.

An expanded global inventory of allelic variation in the most extremely polymorphic region of Plasmodium falciparum merozoite surface protein 1 provided by short read sequence data

Background

Within Plasmodium falciparum merozoite surface protein 1 (MSP1), the N-terminal block 2 region is a highly polymorphic target of naturally acquired antibody responses. The antigenic diversity is determined by complex repeat sequences as well as non-repeat sequences, grouping into three major allelic types that appear to be maintained within populations by natural selection. Within these major types, many distinct allelic sequences have been described in different studies, but the extent and significance of the diversity remains unresolved.

Methods

To survey the diversity more extensively, block 2 allelic sequences in the msp1 gene were characterized in 2400 P. falciparum infection isolates with whole genome short read sequence data available from the Pf3K project, and compared with the data from previous studies.

Results

Mapping the short read sequence data in the 2400 isolates to a reference library of msp1 block 2 allelic sequences yielded 3815 allele scores at the level of major allelic family types, with 46% of isolates containing two or more of these major types. Overall frequencies were similar to those previously reported in other samples with different methods, the K1-like allelic type being most common in Africa, MAD20-like most common in Southeast Asia, and RO33-like being the third most abundant type in each continent. The rare MR type, formed by recombination between MAD20-like and RO33-like alleles, was only seen in Africa and very rarely in the Indian subcontinent but not in Southeast Asia. A combination of mapped short read assembly approaches enabled 1522 complete msp1 block 2 sequences to be determined, among which there were 363 different allele sequences, of which 246 have not been described previously. In these data, the K1-like msp1 block 2 alleles are most diverse and encode 225 distinct amino acid sequences, compared with 123 different MAD20-like, 9 RO33-like and 6 MR type sequences. Within each of the major types, the different allelic sequences show highly skewed geographical distributions, with most of the more common sequences being detected in either Africa or Asia, but not in both.

Conclusions

Allelic sequences of this extremely polymorphic locus have been derived from whole genome short read sequence data by mapping to a reference library followed by assembly of mapped reads. The catalogue of sequence variation has been greatly expanded, so that there are now more than 500 different msp1 block 2 allelic sequences described. This provides an extensive reference for molecular epidemiological genotyping and sequencing studies, and potentially for design of a multi-allelic vaccine.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2475-2) contains supplementary material, which is available to authorized users.

Unraveling Haplotype Diversity of the Apical Membrane Antigen-1 Gene in Plasmodium falciparum Populations in Thailand

Development of an effective vaccine is critically needed for the prevention of malaria. One of the key antigens for malaria vaccines is the apical membrane antigen 1 (AMA-1) of the human malaria parasite Plasmodium falciparum, the surface protein for erythrocyte invasion of the parasite. The gene encoding AMA-1 has been sequenced from populations of P. falciparum worldwide, but the haplotype diversity of the gene in P. falciparum populations in the Greater Mekong Subregion (GMS), including Thailand, remains to be characterized. In the present study, the AMA-1 gene was PCR amplified and sequenced from the genomic DNA of 65 P. falciparum isolates from 5 endemic areas in Thailand. The nearly full-length 1,848 nucleotide sequence of AMA-1 was subjected to molecular analyses, including nucleotide sequence diversity, haplotype diversity and deduced amino acid sequence diversity and neutrality tests. Phylogenetic analysis and pairwise population differentiation (F_st indices) were performed to infer the population structure. The analyses identified 60 single nucleotide polymorphic loci, predominately located in domain I of AMA-1. A total of 31 unique AMA-1 haplotypes were identified, which included 11 novel ones. The phylogenetic tree of the AMA-1 haplotypes revealed multiple clades of AMA-1, each of which contained parasites of multiple geographical origins, consistent with the F_st indices indicating genetic homogeneity or gene flow among geographically distinct populations of P. falciparum in Thailand’s borders with Myanmar, Laos and Cambodia. In summary, the study revealed novel haplotypes and population structure needed for the further advancement of AMA-1-based malaria vaccines in the GMS.

Genetic Polymorphism of msp1 and msp2 in Plasmodium falciparum Isolates from Côte d'Ivoire versus Gabon

Introduction. The characterization of genetic profile of Plasmodium isolates from different areas could help in better strategies for malaria elimination. This study aimed to compare P. falciparum diversity in two African countries. Methods. Isolates collected from 100 and 73 falciparum malaria infections in sites of Côte d'Ivoire (West Africa) and Gabon (Central Africa), respectively, were analyzed by a nested PCR amplification of msp1 and msp2 genes. Results. The K1 allelic family was widespread in Côte d'Ivoire (64.6%) and in Gabon (56.6%). For msp2, the 3D7 alleles were more prevalent (>70% in both countries) compared to FC27 alleles. In Côte d'Ivoire, the frequencies of multiple infections with msp1 (45.1%) and msp2 (40.3%) were higher than those found for isolates from Gabon, that is, 30.2% with msp1 and 31.4% with msp2. The overall complexity of infection was 1.66 (SD = 0.79) in Côte d'Ivoire and 1.58 (SD = 0.83) in Gabon. It decreased with age in Côte d'Ivoire in contrast to Gabon. Conclusion. Differences observed in some allelic families and in complexity profile may suggest an impact of epidemiological facies as well as immunological response on genetic variability of P. falciparum.

Multiplicity of Infection and Disease Severity in Plasmodium vivax

Background

Multiplicity of infection (MOI) refers to the average number of distinct parasite genotypes concurrently infecting a patient. Although several studies have reported on MOI and the frequency of multiclonal infections in Plasmodium falciparum, there is limited data on Plasmodium vivax. Here, MOI and the frequency of multiclonal infections were studied in areas from South America where P. vivax and P. falciparum can be compared.

Methodology/Principal Findings

As part of a passive surveillance study, 1,328 positive malaria patients were recruited between 2011 and 2013 in low transmission areas from Colombia. Of those, there were only 38 P. vivax and 24 P. falciparum clinically complicated cases scattered throughout the time of the study. Samples from uncomplicated cases were matched in time and location with the complicated cases in order to compare the circulating genotypes for these two categories. A total of 92 P. vivax and 57 P. falciparum uncomplicated cases were randomly subsampled. All samples were genotyped by using neutral microsatellites. Plasmodium vivax showed more multiclonal infections (47.7%) than P. falciparum (14.8%). Population genetics and haplotype network analyses did not detect differences in the circulating genotypes between complicated and uncomplicated cases in each parasite. However, a Fisher exact test yielded a significant association between having multiclonal P. vivax infections and complicated malaria. No association was found for P. falciparum infections.

Conclusion

The association between multiclonal infections and disease severity in P. vivax is consistent with previous observations made in rodent malaria. The contrasting pattern between P. vivax and P. falciparum could be explained, at least in part, by the fact that P. vivax infections have lineages that were more distantly related among them than in the case of the P. falciparum multiclonal infections. Future research should address the possible role that acquired immunity and exposure may have on multiclonal infections and their association with disease severity.

Previous studies on rodent malarias and mathematical models have postulated a link between multiclonal infections and disease severity. This association has been tested in Plasmodium falciparum mostly in Africa with limited information on P. vivax. Furthermore, there is a paucity of information from areas with low transmission. Here, we used samples available from a passive surveillance carried out in Colombia, South America. We found an association between multiclonal infections and disease severity in P. vivax but not in P. falciparum. Although the number of complicated malaria cases is low, the contrasting pattern between these two species emphasizes their epidemiological differences. We discuss how this pattern could be the result of a higher divergence among the P. vivax lineages co-infecting a patient. We hypothesize that low levels of acquired immunity may play a role in the association between multiclonal infections and disease severity.

High-throughput quantitative proteomic analysis of dengue virus type 2 infected A549 cells

Disease caused by dengue virus is a global health concern with up to 390 million individuals infected annually worldwide. There are no vaccines or antiviral compounds available to either prevent or treat dengue disease which may be fatal. To increase our understanding of the interaction of dengue virus with the host cell, we analyzed changes in the proteome of human A549 cells in response to dengue virus type 2 infection using stable isotope labelling in cell culture (SILAC) in combination with high-throughput mass spectrometry (MS). Mock and infected A549 cells were fractionated into nuclear and cytoplasmic extracts before analysis to identify proteins that redistribute between cellular compartments during infection and reduce the complexity of the analysis. We identified and quantified 3098 and 2115 proteins in the cytoplasmic and nuclear fractions respectively. Proteins that showed a significant alteration in amount during infection were examined using gene enrichment, pathway and network analysis tools. The analyses revealed that dengue virus infection modulated the amounts of proteins involved in the interferon and unfolded protein responses, lipid metabolism and the cell cycle. The SILAC-MS results were validated for a select number of proteins over a time course of infection by Western blotting and immunofluorescence microscopy. Our study demonstrates for the first time the power of SILAC-MS for identifying and quantifying novel changes in cellular protein amounts in response to dengue virus infection.