Genetic diversity and natural selection of three blood-stage 6-Cys proteins in Plasmodium vivax populations from the China-Myanmar endemic border

Comparison of high throughput RNA sequences between Babesia bigemina and Babesia bovis revealed consistent differential gene expression that is required for the Babesia life cycle in the vertebrate and invertebrate hosts

Bovine babesiosis caused by Babesia bigemina and Babesia bovis is an economically important disease that affects cattle worldwide. Both B. bigemina and B. bovis are transovarially transmitted by Rhipicephalus ticks. However, little is known regarding parasite gene expression during infection of the tick vector or mammalian host, which has limited the development of effective control strategies to alleviate the losses to the cattle industry. To understand Babesia gene regulation during tick and mammalian host infection, we performed high throughput RNA-sequencing using samples collected from calves and Rhipicephalus microplus ticks infected with B. bigemina. We evaluated gene expression between B. bigemina blood-stages and kinetes and compared them with previous B. bovis RNA-seq data. The results revealed similar patterns of gene regulation between these two tick-borne transovarially transmitted Babesia parasites. Like B. bovis, the transcription of several B. bigemina genes in kinetes exceeded a 1,000-fold change while a few of these genes had a >20,000-fold increase. To identify genes that may have important roles in B. bigemina and B. bovis transovarial transmission, we searched for genes upregulated in B. bigemina kinetes in the genomic datasets of B. bovis and non-transovarially transmitted parasites, Theileria spp. and Babesia microti. Using this approach, we identify genes that may be potential markers for transovarial transmission by B. bigemina and B. bovis. The findings presented herein demonstrate common Babesia genes linked to infection of the vector or mammalian host and may contribute to elucidating strategies used by the parasite to complete their life cycle.

The contributions and achievements on malaria control and forthcoming elimination in China over the past 70 years by NIPD-CTDR

Although the past decades have seen a remarkable decrease in malaria-caused mortality and morbidity, the infection remains a significant challenge to global health. In the battle against malaria, China has gained notable feat and achievement since the 1940s through the efforts of several generations. Notably, China has not recorded a single indigenous malaria case since August 2016. The National Institute of Parasitic Diseases of the Chinese Center for Disease Control and Prevention (NIPD), as the only specialized institution for parasitic disease at the national level, has played a significant role in the malaria control, prevention, and elimination in China in the different historical periods. In order to transfer Chinese experiences on malaria control and elimination to other Low and Middle Income Countries (LMICs) and to improve global health collaboration, we have summarized and reviewed the contributions and achievements by the NIPD over the past 70 years, covering the epidemic situation; field investigation and laboratory experimental research on both parasite and vector; research and development on diagnostics, drugs, and insecticides; surveillance and response; technical and international. Support and cooperation. In addition, we also focus in particular on malaria retransmission risk, strategies on management of imported malaria cases and mobile populations, surveillance and response capacity to be maintained in post-elimination stage, challenges on diagnosis, drug resistance, and insecticide resistance as future concerns.

The Plasmodium knowlesi Pk41 surface protein diversity, natural selection, sub population and geographical clustering: a 6-cysteine protein family member

Introduction

The zoonotic malaria parasite Plasmodium knowlesi has currently become the most dominant form of infection in humans in Malaysia and is an emerging infectious disease in most Southeast Asian countries. The P41 is a merozoite surface protein belonging to the 6-cysteine family and is a well-characterized vaccine candidate in P. vivax and P. falciparum; however, no study has been done in the orthologous gene of P. knowlesi. This study investigates the level of polymorphism, haplotypes and natural selection of pk41 genes in clinical isolates from Malaysia.

Method

Thirty-five full-length pk41 sequences from clinical isolates of Malaysia along with four laboratory lines (along with H-strain) were downloaded from public databases. For comparative analysis between species, orthologous P41 genes from P. falciparum, P. vivax, P. coatneyi and P. cynomolgi were also downloaded. Genetic diversity, polymorphism, haplotype and natural selection were determined using DnaSP 5.10 software. Phylogenetic relationships between Pk41 genes were determined using MEGA 5.0 software.

Results

Analysis of 39 full-length pk41 sequences along with the H-strain identified 36 SNPs (20 non-synonymous and 16 synonymous substitutions) resulting in 31 haplotypes. Nucleotide diversity across the full-length gene was low and was similar to its ortholog in P. vivax; pv41. Domain-wise amino acid analysis of the two s48/45 domains indicated low level of polymorphisms for both the domains, and the glutamic acid rich region had extensive size variations. In the central domain, upstream to the glutamate rich region, a unique two to six (K-E)_n repeat region was identified within the clinical isolates. Overall, the pk41 genes were indicative of negative/purifying selection due to functional constraints. Domain-wise analysis of the s48/45 domains also indicated purifying selection. However, analysis of Tajima’s D across the genes identified non-synonymous SNPs in the s48/45 domain II with high positive values indicating possible epitope binding regions. All the 6-cysteine residues within the s48/45 domains were conserved within the clinical isolates indicating functional conservation of these regions. Phylogenetic analysis of full-length pk41 genes indicated geographical clustering and identified three subpopulations of P. knowlesi; one originating in the laboratory lines and two originating from Sarawak, Malaysian Borneo.

Conclusion

This is the first study to report on the polymorphism and natural selection of pk41 genes from clinical isolates of Malaysia. The results reveal that there is low level of polymorphism in both s48/45 domains, indicating that this antigen could be a potential vaccine target. However, genetic and molecular immunology studies involving higher number of samples from various parts of Malaysia would be necessary to validate this antigen’s candidacy as a vaccine target for P. knowlesi.

Whole-genome sequencing and analysis of Plasmodium falciparum isolates from China-Myanmar border area

BACKGROUND

China has made progress in malaria control and aims to eliminate malaria nationwide, but implementing effective interventions along the border regions remain a huge task. The Plasmodium falciparum cases imported from Southeast Asia has frequently reported especially in the China-Myanmar border (CMB) area. Though, information is scant on P. falciparum genetic variability in this area.

METHODS

This study reported P. falciparum isolates genome sequence of six clinical isolates in the CMB area. Furthermore, we estimated the nucleotide diversity, Watterson's estimator and Tajima's D value for the whole genome mutation rate in slide window.

RESULTS

Our data were aligned onto 96.05-98.61% of the reference 3D7 genome in high fold coverages. Principal component analysis result showed that P. falciparum clustered generally according to their geographic origin. A total of 91 genes were identified as positive selection with Ka/Ks ratio significantly higher than 1, and most of them were multigene families encoding variant surface antigens (VSAs) such as var, rif and stevor. The enrichment of the positive selection on VSA genes implied that the environment complexity subjected CMB's P. falciparum to more pressure for survival.

CONCLUSIONS

Our research suggests that greater genetic diversity in CMB area and the positive selection signals in VSA genes, which allow P. falciparum to fit the host immune system well and aggravate the difficulty of treatment. Meanwhile, results obtained from this study will provide the fundamental basis for P. falciparum population genomic research in CMB area.

On the Evolution and Function of Plasmodium vivax Reticulocyte Binding Surface Antigen (pvrbsa)

The RBSA protein is encoded by a gene described in Plasmodium species having tropism for reticulocytes. Since this protein is antigenic in natural infections and can bind to target cells, it has been proposed as a potential candidate for an anti-Plasmodium vivax vaccine. However, genetic diversity (a challenge which must be overcome for ensuring fully effective vaccine design) has not been described at this locus. Likewise, the minimum regions mediating specific parasite-host interaction have not been determined. This is why the rbsa gene’s evolutionary history is being here described, as well as the P. vivax rbsa (pvrbsa) genetic diversity and the specific regions mediating parasite adhesion to reticulocytes. Unlike what has previously been reported, rbsa was also present in several parasite species belonging to the monkey-malaria clade; paralogs were also found in Plasmodium parasites invading reticulocytes. The pvrbsa locus had less diversity than other merozoite surface proteins where natural selection and recombination were the main evolutionary forces involved in causing the observed polymorphism. The N-terminal end (PvRBSA-A) was conserved and under functional constraint; consequently, it was expressed as recombinant protein for binding assays. This protein fragment bound to reticulocytes whilst the C-terminus, included in recombinant PvRBSA-B (which was not under functional constraint), did not. Interestingly, two PvRBSA-A-derived peptides were able to inhibit protein binding to reticulocytes. Specific conserved and functionally important peptides within PvRBSA-A could thus be considered when designing a fully-effective vaccine against P. vivax.

Study roadmap for high-throughput development of easy to use and affordable biomarkers as diagnostics for tropical diseases: a focus on malaria and schistosomiasis

BACKGROUND

Interventions are currently being used against 'infectious diseases of poverty', which remain highly debilitating and deadly in most endemic countries, especially malaria, schistosomiasis, echinococcosis and African sleeping sickness. However, major limitations of current 'traditional' methods for diagnosis are neither simple nor convenient for population surveillance, and showed low sensitivity and specificity. Access to novel technologies for the development of adequate and reliable tools are expressly needed. A collaborative project between African Network for Drugs and Diagnostics Innovation and partner institutions in Africa and China aims to screen suitable serological biomarkers for diagnostic pipelines against these 'diseases of the poor'.

METHODS

Parasite-specific exposed versus unexposed individuals were screened and sera or urine/stools were collected through case-control studies in China and African countries. Target genes/open reading frames were selected, then will be cloned and cell-free expressed, quantified and immuno-detected. Target antigens/epitopes will be probed and screened with sera from exposed or unexposed individuals using a high-throughput antigen screening platform as the study progresses. The specificity and sensitivity of highly immunoreactive biomarkers will be evaluated as well, using enzyme-linked immunosorbent assays or dipsticks.

DISCUSSION

This roadmap explicitly unfolds the integrated operating procedures with focus on malaria and schistosomiasis, for the identification of suitable biomarkers that will aid the prioritization of diagnostics for population use. However, there is need to further validate any new diagnostic through comparison with standard methods in field deployable tests for each region. Our expectations for the future are to seek regulatory approval and promote the use of diagnostics in endemic areas.

Genetic Diversity and Natural Selection in 42 kDa Region of Plasmodium vivax Merozoite Surface Protein-1 from China-Myanmar Endemic Border

Plasmodium vivax merozoite surface protein-1 (PvMSP1) gene codes for a major malaria vaccine candidate antigen. However, its polymorphic nature represents an obstacle to the design of a protective vaccine. In this study, we analyzed the genetic polymorphism and natural selection of the C-terminal 42 kDa fragment within PvMSP1 gene (Pv MSP142) from 77 P. vivax isolates, collected from imported cases of China-Myanmar border (CMB) areas in Yunnan province and the inland cases from Anhui, Yunnan, and Zhejiang province in China during 2009-2012. Totally, 41 haplotypes were identified and 30 of them were new haplotypes. The differences between the rates of non-synonymous and synonymous mutations suggest that PvMSP142 has evolved under natural selection, and a high selective pressure preferentially acted on regions identified of PvMSP133. Our results also demonstrated that PvMSP142 of P. vivax isolates collected on China-Myanmar border areas display higher genetic polymorphisms than those collected from inland of China. Such results have significant implications for understanding the dynamic of the P. vivax population and may be useful information towards China malaria elimination campaign strategies.

PvGAMA reticulocyte binding activity: predicting conserved functional regions by natural selection analysis

Background

Adhesin proteins are used by Plasmodium parasites to bind and invade target cells. Hence, characterising molecules that participate in reticulocyte interaction is key to understanding the molecular basis of Plasmodium vivax invasion. This study focused on predicting functionally restricted regions of the P. vivax GPI-anchored micronemal antigen (PvGAMA) and characterising their reticulocyte binding activity.

Results

The pvgama gene was initially found in P. vivax VCG-I strain schizonts. According to the genetic diversity analysis, PvGAMA displayed a size polymorphism very common for antigenic P. vivax proteins. Two regions along the antigen sequence were highly conserved among species, having a negative natural selection signal. Interestingly, these regions revealed a functional role regarding preferential target cell adhesion.

Conclusions

To our knowledge, this study describes PvGAMA reticulocyte binding properties for the first time. Conserved functional regions were predicted according to natural selection analysis and their binding ability was confirmed. These findings support the notion that PvGAMA may have an important role in P. vivax merozoite adhesion to its target cells.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-017-2183-8) contains supplementary material, which is available to authorized users.

Whole-genome sequencing of a Plasmodium vivax clinical isolate exhibits geographical characteristics and high genetic variation in China-Myanmar border area

Background

Currently in China, the trend of Plasmodium vivax cases imported from Southeast Asia was increased especially in the China-Myanmar border area. Driven by the increase in P. vivax cases and stronger need for vaccine and drug development, several P. vivax isolates genome sequencing projects are underway. However, little is known about the genetic variability in this area until now.

Results

The sequencing of the first P. vivax isolate from China-Myanmar border area (CMB-1) generated 120 million paired-end reads. A percentage of 10.6 of the quality-evaluated reads were aligned onto 99.9% of the reference strain Sal I genome in 62-fold coverage with an average of 4.8 SNPs per kb. We present a 539-SNP marker data set for P. vivax that can identify different parasites from different geographic origins with high sensitivity. We also identified exceptionally high levels of genetic variability in members of multigene families such as RBP, SERA, vir, MSP3 and AP2. The de-novo assembly yielded a database composed of 8,409 contigs with N50 lengths of 6.6 kb and revealed 661 novel predicted genes including 78 vir genes, suggesting a greater functional variation in P. vivax from this area.

Conclusion

Our result contributes to a better understanding of P. vivax genetic variation, and provides a fundamental basis for the geographic differentiation of vivax malaria from China-Myanmar border area using a direct sequencing approach without leukocyte depletion. This novel sequencing method can be used as an essential tool for the genomic research of P. vivax in the near future.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3523-y) contains supplementary material, which is available to authorized users.

The s48/45 six-cysteine proteins: mediators of interaction throughout the Plasmodium life cycle

During their life cycle Plasmodium parasites rely upon an arsenal of proteins that establish key interactions with the host and vector, and between the parasite sexual stages, with the purpose of ensuring infection, reproduction and proliferation. Among these is a group of secreted or membrane-anchored proteins known as the six-cysteine (6-cys) family. This is a small but important family with only 14 members thus far identified, each stage-specifically expressed during the parasite life cycle. 6-cys proteins often localise at the parasite surface or interface with the host and vector, and are conserved in different Plasmodium species. The unifying feature of the family is the s48/45 domain, presumably involved in adhesion and structurally related to Ephrins, the ligands of Eph receptors. The most prominent s48/45 members are currently under functional investigation and are being pursued as vaccine candidates. In this review, we examine what is known about the 6-cys family, their structure and function, and discuss future research directions.

The Structure of Plasmodium falciparum Blood-Stage 6-Cys Protein Pf41 Reveals an Unexpected Intra-Domain Insertion Required for Pf12 Coordination

Plasmodium falciparum is an apicomplexan parasite and the etiological agent of severe human malaria. The complex P. falciparum life cycle is supported by a diverse repertoire of surface proteins including the family of 6-Cys s48/45 antigens. Of these, Pf41 is localized to the surface of the blood-stage merozoite through its interaction with the glycophosphatidylinositol-anchored Pf12. Our recent structural characterization of Pf12 revealed two juxtaposed 6-Cys domains (D1 and D2). Pf41, however, contains an additional segment of 120 residues predicted to form a large spacer separating its two 6-Cys domains. To gain insight into the assembly mechanism and overall architecture of the Pf12-Pf41 complex, we first determined the 2.45 Å resolution crystal structure of Pf41 using zinc single-wavelength anomalous dispersion. Structural analysis revealed an unexpected domain organization where the Pf41 6-Cys domains are, in fact, intimately associated and the additional residues instead map predominately to an inserted domain-like region (ID) located between two β-strands in D1. Notably, the ID is largely proteolyzed in the final structure suggesting inherent flexibility. To assess the contribution of the ID to complex formation, we engineered a form of Pf41 where the ID was replaced by a short glycine-serine linker and showed by isothermal titration calorimetry that binding to Pf12 was abrogated. Finally, protease protection assays showed that the proteolytic susceptibility of the ID was significantly reduced in the complex, consistent with the Pf41 ID directly engaging Pf12. Collectively, these data establish the architectural organization of Pf41 and define an essential role for the Pf41 ID in promoting assembly of the Pf12-Pf41 heterodimeric complex.