Profiling the malaria genome: a gene survey of three species of malaria parasite with comparison to other apicomplexan species

Molecular surveillance of chloroquine resistance in Plasmodium vivax isolates from malaria cases in Yunnan Province of China using pvcrt-o gene polymorphisms

BACKGROUND

The efficacy of chloroquine treatment for vivax malaria has been rarely evaluated due to a lack of an appropriate testing method. The objective of this study was to conduct molecular monitoring of chloroquine resistance in Plasmodium vivax strains from vivax malaria patients in Yunnan Province, focusing on the analysis of polymorphism in the P. vivax chloroquine resistance transporter protein orthologous gene (pvcrt-o).

METHODS

In accordance with the principles of a cohort study, blood samples were collected from malaria cases diagnosed with a P. vivax mono-infection in Yunnan Province from 2020 to 2022. Segmental PCR was used to amplify the whole pvcrt-o gene in the blood samples and their products were subsequently sequenced. The sequencing data were arranged to obtain the full coding DNA sequence (CDS) as well as the gene's promoter region sequences. The CDSs were aligned with the reference sequence (XM_001613407.1) of the P. vivax SalI isolate to identify the mutant loci.

RESULTS

From a total of 375 blood samples taken from vivax malaria cases, 272 both whole gene CDSs (1272-1275 bp) and promoter DNA sequences (707 bp) of pvcrt-o gene were obtained. Among the whole CDSs, there were 7 single nucleotide polymorphic sites in which c.7 A>G was the minor allele frequency (MAF) site with 4.4% (12/272) detection rate. The mutation detection rate showed a significant decrease from 9.8% (10/102) in 2020 to 1.1% (1/92) in 2021 and 1.3% (1/78) in 2022, indicating statistical significance (χ2 = 11.256, P < 0.05). Among the identified 12 haplotypes, the majority of which were wild type (75.7%; 206/272). These four mutant haplotypes (Hap_3, Hap_5, Hap_9, and Hap_10) were classified as "K10 insertion type" and accounted for 12.1% (33/272). The detection rate of Hap_3 increased from 1.0% (1/102) in 2020 to 13.0% (12/92) in 2021 and 14.1% (11/78) in 2022, indicating statistical significance. A total of 23.8% (65/272) of the samples exhibited 14 bp (bp) deletions in the promoter region, occurring most frequently in the wild type haplotype (Hap_1) samples at a rate of 28.6% (59/206).

CONCLUSIONS

In recent years in Yunnan Province, a notable proportion of vivax malaria patients are infected by P. vivax strains with a "K10 insertion" and partial sequence deletions in the promoter region of the pvcrt-o gene, necessitating vigilance.

Directional gene expression and antisense transcripts in sexual and asexual stages of Plasmodium falciparum

BACKGROUND

It has been shown that nearly a quarter of the initial predicted gene models in the Plasmodium falciparum genome contain errors. Although there have been efforts to obtain complete cDNA sequences to correct the errors, the coverage of cDNA sequences on the predicted genes is still incomplete, and many gene models for those expressed in sexual or mosquito stages have not been validated. Antisense transcripts have widely been reported in P. falciparum; however, the extent and pattern of antisense transcripts in different developmental stages remain largely unknown.

RESULTS

We have sequenced seven bidirectional libraries from ring, early and late trophozoite, schizont, gametocyte II, gametocyte V, and ookinete, and four strand-specific libraries from late trophozoite, schizont, gametocyte II, and gametocyte V of the 3D7 parasites. Alignment of the cDNA sequences to the 3D7 reference genome revealed stage-specific antisense transcripts and novel intron-exon splicing junctions. Sequencing of strand-specific cDNA libraries suggested that more genes are expressed in one direction in gametocyte than in schizont. Alternatively spliced genes, antisense transcripts, and stage-specific expressed genes were also characterized.

CONCLUSIONS

It is necessary to continue to sequence cDNA from different developmental stages, particularly those of non-erythrocytic stages. The presence of antisense transcripts in some gametocyte and ookinete genes suggests that these antisense RNA may play an important role in gene expression regulation and parasite development. Future gene expression studies should make use of directional cDNA libraries. Antisense transcripts may partly explain the observed discrepancy between levels of mRNA and protein expression.

Expressed sequence tag analysis of the erythrocytic stage of Plasmodium berghei

Rodent malaria parasites, such as Plasmodium berghei, are practical and useful model organisms for human malaria research because of their analogies to the human malaria in terms of structure, physiology, and life cycle. Exploiting the available genetic sequence information, we constructed a cDNA library from the erythrocytic stages of P. berghei and analyzed the expressed sequence tag (EST). A total of 10,040 ESTs were generated and assembled into 2,462 clusters. These EST clusters were compared against public protein databases and 48 putative new transcripts, most of which were hypothetical proteins with unknown function, were identified. Genes encoding ribosomal or membrane proteins and purine nucleotide phosphorylases were highly abundant clusters in P. berghei. Protein domain analyses and the Gene Ontology functional categorization revealed translation/protein folding, metabolism, protein degradation, and multiple family of variant antigens to be mainly prevalent. The presently-collected ESTs and its bioinformatic analysis will be useful resources to identify for drug target and vaccine candidates and validate gene predictions of P. berghei.

Oil palm (Elaeis guineensis Jacq.) tissue culture ESTs: identifying genes associated with callogenesis and embryogenesis

BACKGROUND

Oil palm (Elaeis guineensis Jacq.) is one of the most important oil bearing crops in the world. However, genetic improvement of oil palm through conventional breeding is extremely slow and costly, as the breeding cycle can take up to 10 years. This has brought about interest in vegetative propagation of oil palm. Since the introduction of oil palm tissue culture in the 1970s, clonal propagation has proven to be useful, not only in producing uniform planting materials, but also in the development of the genetic engineering programme. Despite considerable progress in improving the tissue culture techniques, the callusing and embryogenesis rates from proliferating callus cultures remain very low. Thus, understanding the gene diversity and expression profiles in oil palm tissue culture is critical in increasing the efficiency of these processes.

RESULTS

A total of 12 standard cDNA libraries, representing three main developmental stages in oil palm tissue culture, were generated in this study. Random sequencing of clones from these cDNA libraries generated 17,599 expressed sequence tags (ESTs). The ESTs were analysed, annotated and assembled to generate 9,584 putative unigenes distributed in 3,268 consensi and 6,316 singletons. These unigenes were assigned putative functions based on similarity and gene ontology annotations. Cluster analysis, which surveyed the relatedness of each library based on the abundance of ESTs in each consensus, revealed that lipid transfer proteins were highly expressed in embryogenic tissues. A glutathione S-transferase was found to be highly expressed in non-embryogenic callus. Further analysis of the unigenes identified 648 non-redundant simple sequence repeats and 211 putative full-length open reading frames.

CONCLUSION

This study has provided an overview of genes expressed during oil palm tissue culture. Candidate genes with expression that are modulated during tissue culture were identified. However, in order to confirm whether these genes are suitable as early markers for embryogenesis, the genes need to be tested on earlier stages of tissue culture and a wider range of genotypes. This collection of ESTs is an important resource for genetic and genome analyses of the oil palm, particularly during tissue culture development.

cDNA sequences reveal considerable gene prediction inaccuracy in the Plasmodium falciparum genome

Background

The completion of the Plasmodium falciparum genome represents a milestone in malaria research. The genome sequence allows for the development of genome-wide approaches such as microarray and proteomics that will greatly facilitate our understanding of the parasite biology and accelerate new drug and vaccine development. Designing and application of these genome-wide assays, however, requires accurate information on gene prediction and genome annotation. Unfortunately, the genes in the parasite genome databases were mostly identified using computer software that could make some erroneous predictions.

Results

We aimed to obtain cDNA sequences to examine the accuracy of gene prediction in silico. We constructed cDNA libraries from mixed blood stages of P. falciparum parasite using the SMART cDNA library construction technique and generated 17332 high-quality expressed sequence tags (EST), including 2198 from primer-walking experiments. Assembly of our sequence tags produced 2548 contigs and 2671 singletons versus 5220 contigs and 5910 singletons when our EST were assembled with EST in public databases. Comparison of all the assembled EST/contigs with predicted CDS and genomic sequences in the PlasmoDB database identified 356 genes with predicted coding sequences fully covered by EST, including 85 genes (23.6%) with introns incorrectly predicted. Careful automatic software and manual alignments found an additional 308 genes that have introns different from those predicted, with 152 new introns discovered and 182 introns with sizes or locations different from those predicted. Alternative spliced and antisense transcripts were also detected. Matching cDNA to predicted genes also revealed silent chromosomal regions, mostly at subtelomere regions.

Conclusion

Our data indicated that approximately 24% of the genes in the current databases were predicted incorrectly, although some of these inaccuracies could represent alternatively spliced transcripts, and that more genes than currently predicted have one or more additional introns. It is therefore necessary to annotate the parasite genome with experimental data, although obtaining complete cDNA sequences from this parasite will be a formidable task due to the high AT nature of the genome. This study provides valuable information for genome annotation that will be critical for functional analyses.

Plasmodium cysteine repeat modular proteins 1-4: complex proteins with roles throughout the malaria parasite life cycle

The Cysteine Repeat Modular Proteins (PCRMP1-4) of Plasmodium, are encoded by a small gene family that is conserved in malaria and other Apicomplexan parasites. They are very large, predicted surface proteins with multipass transmembrane domains containing motifs that are conserved within families of cysteine-rich, predicted surface proteins in a range of unicellular eukaryotes, and a unique combination of protein-binding motifs, including a >100 kDa cysteine-rich modular region, an epidermal growth factor-like domain and a Kringle domain. PCRMP1 and 2 are expressed in life cycle stages in both the mosquito and vertebrate. They colocalize with PfEMP1 (P. falciparum Erythrocyte Membrane Antigen-1) during its export from P. falciparum blood-stage parasites and are exposed on the surface of haemolymph- and salivary gland-sporozoites in the mosquito, consistent with a role in host tissue targeting and invasion. Gene disruption of pcrmp1 and 2 in the rodent malaria model, P. berghei, demonstrated that both are essential for transmission of the parasite from the mosquito to the mouse and has established their discrete and important roles in sporozoite targeting to the mosquito salivary gland. The unprecedented expression pattern and structural features of the PCRMPs thus suggest a variety of roles mediating host-parasite interactions throughout the parasite life cycle.

Plasmodium vivax genetic diversity: microsatellite length matters

The Plasmodium vivax genome is very diverse but has a relatively low abundance of microsatellites. Leclerc et al. had shown that these di-nucleotide repeats have a low level of polymorphism, suggesting a recent bottleneck event in the evolutionary history of P. vivax. By contrast, in a recent paper, Imwong et al. show that there is a very high level of microsatellite diversity. The difference in these results is probably due to the set array lengths chosen by each group. Longer arrays are more diverse than are shorter ones because slippage mutations become exponentially more common with an increase in array length. These studies highlight the need to consider carefully the application and design of studies involving microsatellites.

Immunomics: discovering new targets for vaccines and therapeutics

T-cell-epitope mapping has emerged as one of the most powerful new drug discovery tools for a range of biomedical applications. Initially, T-cell-epitope discovery was applied to the development of vaccines for infectious diseases and cancer. T-cell-epitope-mapping applications have now expanded to include reengineering of protein therapeutics (a process now called deimmunization), as well as the fields of autoimmunity, endocrinology, allergy, transplantation and diagnostics. Research employing T-cell-epitope mapping falls within the realm of immunomics, a new field that addresses the interface between host and (pathogen) proteome, bridging informatics, genomics, proteomics, immunology and clinical medicine. This review highlights aspects of recent immunomics research that are related to the discovery of the T-cell immunome.

Gene discovery in Plasmodium vivax through sequencing of ESTs from mixed blood stages

Despite the significance of Plasmodium vivax as the most widespread human malaria parasite and a major public health problem, gene expression in this parasite is poorly understood. To accelerate gene discovery and facilitate the annotation phase of the P. vivax genome project, we have undertaken a transcriptome approach to study gene expression in the mixed blood stages of a P. vivax field isolate. Using a cDNA library constructed from purified blood stages, we have obtained single-pass sequences for approximately 21,500 expressed sequence tags (ESTs), the largest number of transcript tags obtained so far for this species. Cluster analysis revealed that the library is highly redundant, resulting in 5407 clusters. Clustered ESTs were searched against public protein databases for functional annotation, and more than one-third showed a significant match, the majority of these to Plasmodium falciparum proteins. The most abundant clusters were to genes encoding ribosomal proteins and proteins involved in metabolism, consistent with the predominance of trophozoites in the field isolate sample. In spite of the scarcity of other parasite stages in the field isolate, we could identify genes that are expressed in rings, schizonts and gametocytes. This study should facilitate our understanding of the gene expression in P. vivax asexual stages and provide valuable data for gene prediction and annotation of the P. vivax genome sequence.

Genome Sequencing and Comparative Genomics of Tropical Disease Pathogens

The sequencing of eukaryotic genomes has lagged behind sequencing of organisms in the other domains of life, archae and bacteria, primarily due to their greater size and complexity. With recent advances in high-throughput technologies such as robotics and improved computational resources, the number of eukaryotic genome sequencing projects has increased significantly. Among these are a number of sequencing projects of tropical pathogens of medical and veterinary importance, many of which are responsible for causing widespread morbidity and mortality in peoples of developing countries. Uncovering the complete gene complement of these organisms is proving to be of immense value in the development of novel methods of parasite control, such as antiparasitic drugs and vaccines, as well as the development of new diagnostic tools. Combining pathogen genome sequences with the host and vector genome sequences is promising to be a robust method for the identification of host-pathogen interactions. Finally, comparative sequencing of related species, especially of organisms used as model systems in the study of the disease, is beginning to realize its potential in the identification of genes, and the evolutionary forces that shape the genes, that are involved in evasion of the host immune response.

Integrative analysis of intraerythrocytic differentially expressed transcripts yields novel insights into the biology of Plasmodium falciparum

Background

The intraerythrocytic development of Plasmodium falciparum, the most virulent human malaria parasite involves asexual and gametocyte stages. There has been a significant increase in disparate datasets derived from genomic and post-genomic analysis of the parasite that necessitates delivery of integrated analysis from which biological processes important to the survival of the parasite can be determined.

Methods

In order to resolve genes associated with stage differentially expressed transcripts, we have developed and implemented an integrative approach that combines evidence from P. falciparum expressed sequence tags (ESTs), genomic, microarray, proteomic and gene ontology data.

Results

A total of 143 gametocyte-overexpressed and 51 asexual-overexpressed transcripts were identified. A subset of 74 genes associated with these transcripts showed evidence of stage-correlated protein expression, of which 53 have not been experimentally characterised. Our study has revealed (1) possible regulatory mechanisms in malaria parasites' gametocyte maturation, (2) correlation between EST and microarray data for a P. falciparum gene family to present unique EST-derived information, (3) candidate drug and antigenic targets on which computational and experimental studies can be performed, and (4) the need for more empirical studies on gene and protein expression in malaria parasites.

Conclusion

Applying different domains of data to the same underlying gene set has yielded novel insights into the biology of the parasite and presents an approach to appraise critically the data quality of post-genomic datasets from malaria parasites.

Malaria parasites lacking eef1a have a normal S/M phase yet grow more slowly due to a longer G1 phase

Eukaryotic elongation factor 1A (eEF1A) plays a central role in protein synthesis, cell growth and morphology. Malaria parasites possess two identical genes encoding eEF1A (eef1aa and eef1ab). Using pbeef1a-Plasmodium berghei mutants that lack an eEF1a gene, we demonstrate that the level of eEF1A production affects the proliferation of blood stages and parasite fitness. Pbeef1a- parasites can complete the vertebrate and mosquito phases of the life cycle, but the growth phase of the asexual blood stages is extended by up to 20%. Analysis of the cell cycle by flow cytometry as well as transcriptional analyses revealed that the duration of the S and M phases and the number of daughter cells produced were not detectably affected, but that the G1 phase is elongated. Thus, as in budding yeast, a growth threshold must be achieved by blood-stage Plasmodium parasites to permit transition from G1 into S/M phase. Initial analyses indicate that transcriptional events associated with gametocyte development were not remarkably retarded. Insight into protein synthesis and its influence on cell proliferation might be used to generate slow-growing (attenuated) parasites.

Plasmepsin 4, the food vacuole aspartic proteinase found in all Plasmodium spp. infecting man

Plasmepsins are aspartic proteinases of the malaria parasite, and seven groups of plasmepsins have been identified by comparing genomic sequence data available for the genes encoding these enzymes from Plasmodium falciparum, Plasmodium vivax, Plasmodium knowlesi, Plasmodium berghei, and Plasmodium yoelii. The food vacuole plasmepsins typified by plasmepsin 4 from P. falciparum (PfPM4) constitute one of these groups. Genes encoding the ortholog of PfPM4 have been cloned from Plasmodium ovale, Plasmodium malariae, and P. vivax. In addition, P. falciparum contains three paralagous food vacuole plasmepsins or plasmepsin-like enzymes that appear to have arisen by gene duplication, plasmepsins 1 (PfPM1), 2 (PfPM2) and HAP, and all four were localized to purified food vacuole preparations by two-dimensional gel electrophoresis and mass spectroscopic analysis. The three paralogs of PfPM4 do not have counterparts in the six other Plasmodium spp. examined by genomic DNA blot analysis and by review of available genomic sequence data. The presence of these paralogs among the food vacuole plasmepsins in P. falciparum as compared with the other three species causing malaria in man will impact efforts to rationally design antimalarials targeting the food vacuole plasmepsins.

Pilot survey of expressed sequence tags (ESTs) from the asexual blood stages of Plasmodium vivax in human patients

BACKGROUND

Plasmodium vivax is the most widely distributed human malaria, responsible for 70-80 million clinical cases each year and large socio-economical burdens for countries such as Brazil where it is the most prevalent species. Unfortunately, due to the impossibility of growing this parasite in continuous in vitro culture, research on P. vivax remains largely neglected.

METHODS

A pilot survey of expressed sequence tags (ESTs) from the asexual blood stages of P. vivax was performed. To do so, 1,184 clones from a cDNA library constructed with parasites obtained from 10 different human patients in the Brazilian Amazon were sequenced. Sequences were automatedly processed to remove contaminants and low quality reads. A total of 806 sequences with an average length of 586 bp met such criteria and their clustering revealed 666 distinct events. The consensus sequence of each cluster and the unique sequences of the singlets were used in similarity searches against different databases that included P. vivax, Plasmodium falciparum, Plasmodium yoelii, Plasmodium knowlesi, Apicomplexa and the GenBank non-redundant database. An E-value of <10(-30) was used to define a significant database match. ESTs were manually assigned a gene ontology (GO) terminology

RESULTS

A total of 769 ESTs could be assigned a putative identity based upon sequence similarity to known proteins in GenBank. Moreover, 292 ESTs were annotated and a GO terminology was assigned to 164 of them.

CONCLUSION

These are the first ESTs reported for P. vivax and, as such, they represent a valuable resource to assist in the annotation of the P. vivax genome currently being sequenced. Moreover, since the GC-content of the P. vivax genome is strikingly different from that of P. falciparum, these ESTs will help in the validation of gene predictions for P. vivax and to create a gene index of this malaria parasite.

Single-nucleotide polymorphisms and genome diversity in Plasmodium vivax

The study of genetic variation in malaria parasites has practical significance for developing strategies to control the disease. Vaccines based on highly polymorphic antigens may be confounded by allelic restriction of the host immune response. In response to drug pressure, a highly plastic genome may generate resistant mutants more easily than a monomorphic one. Additionally, the study of the distribution of genomic polymorphisms may provide information leading to the identification of genes associated with traits such as parasite development and drug resistance. Indeed, the age and diversity of the human malaria parasite Plasmodium falciparum has been the subject of recent debate, because an ancient parasite with a complex genome is expected to present greater challenges for drug and vaccine development. The genome diversity of the important human pathogen Plasmodium vivax, however, remains essentially unknown. Here we analyze an approximately 100-kb contiguous chromosome segment from five isolates, revealing 191 single-nucleotide polymorphisms (SNPs) and 44 size polymorphisms. The SNPs are not evenly distributed across the segment with blocks of high and low diversity. Whereas the majority (approximately 63%) of the SNPs are in intergenic regions, introns contain significantly less SNPs than intergenic sequences. Polymorphic tandem repeats are abundant and are more uniformly distributed at a frequency of about one polymorphic tandem repeat per 3 kb. These data show that P. vivax has a highly diverse genome, and provide useful information for further understanding the genome diversity of the parasite.

The Plasmodium vivax genome sequencing project

With the successful completion of the project to sequence the Plasmodium falciparum genome, researchers are now turning their attention to other malaria parasite species. Here, an update on the Plasmodium vivax genome sequencing project is presented, as part of the Trends in Parasitology series of reviews expanding on various aspects of P. vivax research.

Toward a catalog for the transcripts and proteins (sialome) from the salivary gland of the malaria vectorAnopheles gambiae

Hundreds of Anopheles gambiae salivary gland cDNA library clones have been sequenced. A cluster analysis based on sequence similarity at e-60 grouped the 691 sequences into 251 different clusters that code for proteins with putative secretory, housekeeping, or unknown functions. Among the housekeeping cDNAs, we found sequences predicted to code for novel thioredoxin, tetraspanin, hemopexin, heat shock protein, and TRIO and MBF proteins. Among secreted cDNAs, we found 21 novel A. gambiaesalivary sequences including those predicted to encode amylase, calreticulin,selenoprotein, mucin-like protein and 30-kDa allergen, in addition to antigen 5- and D7-related proteins, three novel salivary gland (SG)-like proteins and eight unique putative secreted proteins (Hypothetical Proteins, HP). The electronic version of this paper contains hyperlinks to FASTA-formatted files for each cluster with the best match to the nonredundant (NR) and conserved domain databases (CDD) in addition to CLUSTAL alignments of each cluster. The N terminus of 12 proteins (SG-1, SG-1-like 2, SG-6, HP 8, HP 9-like, 5′nucleotidase, 30-kDa protein, antigen 5- and four D7-related proteins) has been identified by Edman degradation of PVDF-transferred, SDS/PAGE-separated salivary gland proteins. Therefore, we contribute to the generation of a catalog of A. gambiae salivary transcripts and proteins. These data are freely available and will eventually become an invaluable tool to study the role of salivary molecules in parasite-host/vector interactions.

Plasmodium, human and Anopheles genomics and malaria

The Plasmodium spp. parasites that cause malaria are transmitted to humans by Anopheles spp. mosquitoes. Scientists have now amassed a great body of knowledge about the parasite, its mosquito vector and human host. Yet this year there will be 300-500 million new malaria infections and 1-3 million deaths caused by the disease. We believe that integrated analyses of genome sequence, DNA polymorphisms, and messenger RNA and protein expression profiles will lead to greater understanding of the molecular basis of vector-human and host-parasite interactions and provide strategies to build upon these insights to develop interventions to mitigate human morbidity and mortality from malaria.

A status report on the sequencing and annotation of the P. falciparum genome

Almost 5 years ago, an international consortium of sequencing centers and funding agencies was formed to sequence the genome of the human malaria parasite Plasmodium falciparum. A novel chromosome by chromosome shotgun strategy was devised to sequence this very AT-rich genome. Two of the 14 chromosomes have been completed and the remaining chromosomes are in the final stages of gap closure. The consortium recently developed plans for the annotation and analysis of the complete genome sequence and its publication in 2002.