The conserved genome organisation of non-falciparum malaria species: the need to know more

The current knowledge on genomes of non-falciparum malaria species and the potential of model malaria parasites for functional analyses are reviewed and compared with those of the most pathogenic human parasite, Plasmodium falciparum. There are remarkable similarities in overall genome composition among the different species at the level of chromosome organisation and chromosome number, conserved order of individual genes, and even conserved functions of specific gene domains and regulatory control elements. With the initiative taken to sequence the genome of P. falciparum, a wealth of information is already becoming available to the scientific community. In order to exploit the biological information content of a complete genome sequence, simple storage of the bulk of sequence data will be inadequate. The requirement for functional analyses to determine the biological role of the open reading frames is commonly accepted and knowledge of the genomes of the animal model malaria species will facilitate these analyses. Detailed comparative genome information and sequencing of additional Plasmodium genomes will provide a deeper insight into the evolutionary history of the species, the biology of the parasite, and its interactions with the mammalian host and mosquito vector. Therefore, an extended and integrated approach will enhance our knowledge of malaria and will ultimately lead to a more rational approach that identifies and evaluates new targets for anti-malarial drug and vaccine development.     

Comparative analysis of BAC and whole genome shotgun sequences from an Anopheles gambiae region related to Plasmodium encapsulation

The only natural mechanism of malaria transmission in sub-Saharan Africa is the mosquito, generally Anopheles gambiae. Blocking malaria parasite transmission by stopping the development of Plasmodium in the insect vector would provide a useful alternative to the current methods of malaria control. Toward this end, it is important to understand the molecular basis of the malaria parasite refractory phenotype in An. gambiae mosquito strains. We have selected and sequenced six bacterial artificial chromosome (BAC) clones from the Pen-1 region that is the major quantitative trait locus involved in Plasmodium encapsulation. The sequence and the annotation of five overlapping BAC clones plus one adjacent, but not contiguous clone, totaling 585kb of genomic sequence from the centromeric end of the Pen-1 region of the PEST strain were compared to that of the genome sequence of the same strain produced by the whole genome shotgun technique. This project identified 23 putative mosquito genes plus putative copies of the retrotransposable elements BEL12 and TRANSIBN1_AG in the six BAC clones. Nineteen of the predicted genes are most similar to their Drosophila melanogaster homologs while one is more closely related to vertebrate genes. Comparison of these new BAC sequences plus previously published BAC sequences to the cognate region of the assembled genome sequence identified three retrotransposons present in one sequence version but not the other. One of these elements, Indy, has not been previously described. These observations provide evidence for the recent active transposition of these elements and demonstrate the plasticity of the Anopheles genome. The BAC sequences strongly support the public whole genome shotgun assembly and automatic annotation while also demonstrating the benefit of complementary genome sequences and of human curation. Importantly, the data demonstrate the differences in the genome sequence of an individual mosquito compared to that of a hypothetical, average genome sequence generated by whole genome shotgun assembly.     

The malaria genome sequencing project: complete sequence of Plasmodium falciparum chromosome 2

An international consortium has been formed to sequence the entire genome of the human malaria parasite Plasmodium falciparum. We sequenced chromosome 2 of clone 3D7 using a shotgun sequencing strategy. Chromosome 2 is 947 kb in length, has a base composition of 80.2% A + T, and contains 210 predicted genes. In comparison to the Saccharomyces cerevisiae genome, chromosome 2 has a lower gene density, a greater proportion of genes containing introns, and nearly twice as many proteins containing predicted non-globular domains. A group of putative surface proteins was identified, rifins, which are encoded by a gene family comprising up to 7% of the protein-encoding gene in the genome. The rifins exhibit considerable sequence diversity and may play an important role in antigenic variation. Sixteen genes encoded on chromosome 2 showed signs of a plastid or mitochondrial origin, including several genes involved in fatty acid biosynthesis. Completion of the chromosome 2 sequence demonstrated that the A + T-rich genome of P. falciparum can be sequenced by the shotgun approach. Within 2-3 years, the sequence of almost all P. falciparum genes will have been determined, paving the way for genetic, biochemical, and immunological research aimed at developing new drugs and vaccines against malaria.     

Fine-scale genetic characterization of Plasmodium falciparum chromosome 7 encompassing the antigenic var and the drug-resistant pfcrt genes

The fact that malaria is still an uncontrolled disease is reflected by the genetic organization of the parasite genome. Efforts to curb malaria should begin with proper understanding of the mechanism by which the parasites evade human immune system and evolve resistance to different antimalarial drugs. We have initiated such a study and presented herewith the results from the in silico understanding of a seventh chromosomal region of the malarial parasite Plasmodium falciparum encompassing the antigenic var genes (coding pfemp1) and the drug-resistant gene pfcrt located at a specified region of the chromosome 7. We found 60 genes of various functions and lengths, majority (61.67%) of them were performing known functions. Almost all the genes have orthologs in other four species of Plasmodium, of which P. chabaudi seems to be the closest to P. falciparum. However, only two genes were found to be paralogous. Interestingly, the drug-resistant gene, pfcrt was found to be surrounded by seven genes coding for several CG proteins out of which six were reported to be responsible for providing drug resistance to P. vivax. The intergenic regions, in this specified region were generally large in size, majority (73%) of them were of more than 500 nucleotide bp length. We also designed primers for amplification of 21 noncoding DNA fragments in the whole region for estimating genetic diversity and inferring the evolutionary history of this region of P. falciparum genome.     

Genome data: what do we learn?

Genome sequence information has continued to accumulate at a spectacular pace during the past year. Details of the sequence and gene content of human chromosome 22 were published. The sequencing and annotation of the first two Arabidopsis thaliana chromosomes was completed. The sequence of chromosome 3 from Plasmodium falciparum, the second sequenced malaria chromosome, was reported, as was that of chromosome 1 from Leishmania major. The complete genomic sequences of five microbes were reported. Approaches to using data from completely sequenced microbial genomes in phylogenetic studies are being explored, as is the application of microarrays to whole genome expression analysis.     

The physics of DNA and the annotation of the Plasmodium falciparum genome

A gene identification procedure is formulated, based on large-scale structural analyses of genomic sequences. The structural property is the physical - thermal - stability of the DNA double-helix, as described by the classical helix-coil model. The analyses are detailed for the Plasmodium falciparum genome, which represents one of the most difficult cases for the gene identification problem (notably because of the extreme AT-richness of the genome). In this genome, the coding domains (either uninterrupted genes or exons in split genes) are accurately identified as regions of high thermal stability. The conclusion is based on the study of the available cloned genes, of which 17 examples are described in detail. These examples demonstrate that the physical criterion is valid for the detection of coding regions whose lengths extend from a few base pairs up to several thousand base pairs. Accordingly, the structural analyses can provide a powerful and convenient tool for the identification of complex genes in the P. falciparum genome. The limits of such a scheme are discussed. The gene identification procedure is applied to the completely sequenced chromosomes (2 and 3), and the results are compared with the database annotations. The structural analyses suggest more or less extensive revision to the annotations, and also allow new putative genes to be identified in the chromosome sequences. Several examples of such new genes are described in detail.     

Interpolated Markov models for eukaryotic gene finding.

Computational gene finding research has emphasized the development of gene finders for bacterial and human DNA. This has left genome projects for some small eukaryotes without a system that addresses their needs. This paper reports on a new system, GlimmerM, that was developed to find genes in the malaria parasite Plasmodium falciparum. Because the gene density in P. falciparum is relatively high, the system design was based on a successful bacterial gene finder, Glimmer. The system was augmented with specially trained modules to find splice sites and was trained on all available data from the P. falciparum genome. Although a precise evaluation of its accuracy is impossible at this time, laboratory tests (using RT-PCR) on a small selection of predicted genes confirmed all of those predictions. With the rapid progress in sequencing the genome of P. falciparum, the availability of this new gene finder will greatly facilitate the annotation process.     

Plasmodium falciparum: an intervening sequence in the GBP 130/96 tR gene

The 96 tR antigen is a heat stable protein produced during the late stages of the intraerythrocytic development of the malaria parasite Plasmodium falciparum and is released into the culture supernatant or the sera of infected patients at the time of schizont rupture. This antigen, identified as a putative protective antigen, was shown to be identical to the glycophorin-binding protein GBP 130 (Perkins 1988, Bonnefoy et al. 1988). We report here that the gene contains a small undescribed intervening sequence located immediately after the sequence coding for the signal sequence. This shows that in P. falciparum, all the genes described so far coding for proteins exported outside the parasitophorous vacuole share a common organization.     

The gene encoding DNA polymerase alpha from Plasmodium falciparum.

The gene encoding DNA polymerase alpha from the human malaria parasite Plasmodium falciparum has been sequenced and characterised. The deduced amino acid sequence possesses the seven sequence motifs which characterise eukaryotic replicative DNA polymerases (I-VII) and four of five motifs (A-E) identified in alpha DNA polymerases. The predicted protein also contains sequences which are reminiscent of Plasmodium proteins but absent from other DNA polymerases. These include four blocks of additional amino acids interspersed with the conserved motifs of the DNA polymerases, four asparagine rich sequences and a novel carboxy-terminal extension. Repetitive sequences similar to those found in other malarial proteins are also present. cDNA-directed PCR was used to establish the presence of these features in the approximately 7kb mRNA. The coding sequence contains a single intron. The gene for DNAPol alpha is located on chromosome 4 and is transcribed in both asexual and sexual erythrocytic stages of the parasite.     

The serine repeat antigen (SERA) gene family phylogeny in Plasmodium: the impact of GC content and reconciliation of gene and species trees

Plasmodium falciparum is the parasite responsible for the most acute form of malaria in humans. Recently, the serine repeat antigen (SERA) in P. falciparum has attracted attention as a potential vaccine and drug target, and it has been shown to be a member of a large gene family. To clarify the relationships among the numerous P. falciparum SERAs and to identify orthologs to SERA5 and SERA6 in Plasmodium species affecting rodents, gene trees were inferred from nucleotide and amino acid sequence data for 33 putative SERA homologs in seven different species. (A distance method for nucleotide sequences that is specifically designed to accommodate differing GC content yielded results that were largely compatible with the amino acid tree. Standard-distance and maximum-likelihood methods for nucleotide sequences, on the other hand, yielded gene trees that differed in important respects.) To infer the pattern of duplication, speciation, and gene loss events in the SERA gene family history, the resulting gene trees were then "reconciled" with two competing Plasmodium species tree topologies that have been identified by previous phylogenetic studies. Parsimony of reconciliation was used as a criterion for selecting a gene tree/species tree pair and provided (1) support for one of the two species trees and for the core topology of the amino acid-derived gene tree, (2) a basis for critiquing fine detail in a poorly resolved region of the gene tree, (3) a set of predicted "missing genes" in some species, (4) clarification of the relationship among the P. falciparum SERA, and (5) some information about SERA5 and SERA6 orthologs in the rodent malaria parasites. Parsimony of reconciliation and a second criterion--implied mutational pattern at two key active sites in the SERA proteins-were also seen to be useful supplements to standard "bootstrap" analysis for inferred topologies.     

Circumsporozoite protein gene from Plasmodium brasilianum. Animal reservoirs for human malaria parasites?

We describe here the sequence of the circumsporozoite protein gene of the monkey malaria parasite Plasmodium brasilianum and show that the immunodominant repeat domain is the same as that of the human malaria parasite, Plasmodium malariae. The immunodominant epitope on the surface of sporozoites of a third species of human malaria parasite has, therefore, been identified. This genetic based data and the biological similarities between P. brasilianum and P. malariae support their putative zoonotic/anthroponotic relationship. We also show that an ape malaria parasite, Plasmodium reichenowi, and the human malaria parasite, Plasmodium falciparum, have a similar relationship. The implications of these observations are discussed with respect to vaccine development.     

A Plasmodium whole-genome synteny map: indels and synteny breakpoints as foci for species-specific genes

Whole-genome comparisons are highly informative regarding genome evolution and can reveal the conservation of genome organization and gene content, gene regulatory elements, and presence of species-specific genes. Initial comparative genome analyses of the human malaria parasite Plasmodium falciparum and rodent malaria parasites (RMPs) revealed a core set of 4,500 Plasmodium orthologs located in the highly syntenic central regions of the chromosomes that sharply defined the boundaries of the variable subtelomeric regions. We used composite RMP contigs, based on partial DNA sequences of three RMPs, to generate a whole-genome synteny map of P. falciparum and the RMPs. The core regions of the 14 chromosomes of P. falciparum and the RMPs are organized in 36 synteny blocks, representing groups of genes that have been stably inherited since these malaria species diverged, but whose relative organization has altered as a result of a predicted minimum of 15 recombination events. P. falciparum-specific genes and gene families are found in the variable subtelomeric regions (575 genes), at synteny breakpoints (42 genes), and as intrasyntenic indels (126 genes). Of the 168 non-subtelomeric P. falciparum genes, including two newly discovered gene families, 68% are predicted to be exported to the surface of the blood stage parasite or infected erythrocyte. Chromosomal rearrangements are implicated in the generation and dispersal of P. falciparum-specific gene families, including one encoding receptor-associated protein kinases. The data show that both synteny breakpoints and intrasyntenic indels can be foci for species-specific genes with a predicted role in host-parasite interactions and suggest that, besides rearrangements in the subtelomeric regions, chromosomal rearrangements may also be involved in the generation of species-specific gene families. A majority of these genes are expressed in blood stages, suggesting that the vertebrate host exerts a greater selective pressure than the mosquito vector, resulting in the acquisition of diversity.     

Curation of the Plasmodium falciparum genome

The malaria genome has proved invaluable to researchers worldwide in the continuing fight against malaria by stimulating and underpinning molecular approaches in gene expression studies, vaccine and drug discovery research, and by providing data to facilitate hypothesis-driven research. The combination of in silico and experimental investigations has already yielded dividends by strengthening our understanding of the many facets of the malaria parasite Plasmodium falciparum. The recently initiated curation of the genome resource is a vital investment for maintaining and enhancing the use of this genomic information in the post-genomic era.     

Mammalian keratin gene families: organisation of genes coding for the B2 high-sulphur proteins of sheep wool.

We have isolated two genomic clones containing three B2 high-sulphur keratin genes from a sheep genomic library constructed in Charon 4A. These genes do not contain intervening sequences. Two genes, encoding the B2A and B2D proteins are closely linked in the genome, being separated by 1.9 kb, and are transcribed in the same direction. Although there is extensive sequence conservation in the 5' non-coding and coding regions, the 3' non-coding regions diverge both in length and sequence. Within the 5' non-coding region adjacent to the initiating AUG there is a highly conserved 18 bp sequence which is also present in another gene coding for a member of a different, unrelated high-sulphur keratin family. In the B2A-B2D intergene region, tightly linked to the B2D gene, there is a putative, divergently transcribed gene.