Differential gene expression in the ookinete stage of the malaria parasite Plasmodium berghei.

Plasmodium, the malaria parasite, undergoes a complex developmental program in its mosquito vector. The ookinete is the parasite form which invades the mosquito midgut and is an important stage for genetic mixing. To identify genes expressed during ookinete development and mosquito midgut invasion, purified zygotes and ookinetes of the rodent parasite Plasmodium berghei were used to construct a suppression subtractive hybridization cDNA library, enriched in sequences expressed in the ookinete stage. In addition to four genes coding for previously described major ookinete-secreted proteins, we isolated ookinete-expressed sequences representing 18 predicted genes. Their gene products include proteins involved in signal transduction and regulatory processes. For six of these genes our analysis provides the first evidence for expression in the ookinete stage. A majority of the genes are not expressed in the zygote, the preceding developmental stage. Furthermore, four of the genes are also transcribed in sporozoites, and one of these in merozoites, suggesting that they code for proteins with a function common to Plasmodium invasive stages.

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 Drosophila melanogaster multidrug-resistance protein 1 (MRP1) homolog has a novel gene structure containing two variable internal exons.

Drosophila melanogaster has a gene very similar to human MRP1 that encodes a full ABC-transporter containing three membrane-spanning domains and two nucleotide-binding domains. This 19 exon insect gene, dMRP (FBgn0032456), spans slightly more than 22 kb. The cDNA SD07655 representing this gene was sequenced and found to contain sequences from 12 exons including single copies of two exons having multiple genomic copies. The gene contains two variant copies of exon 4 and seven of exon 8. While a cDNA contains only one version of each variable exon, all forms of these variable exons were detected in adult fly mRNA. These results predict that Drosophila could make 14 different MRP isoforms from a single gene by substituting different variable exons. This is the first report of any organism using differential splicing of alternative, internal exons, to produce such a large array of MRP isoforms having the same size, but with limited and defined internal variations. Defining the functional differences in the dMRP isoforms should provide clues to the structure/function relationships of the amino acids in these MRP domains, both for the insect enzyme and for those of other species.

Identification of the Anopheles gambiae ATP-binding cassette transporter superfamily genes.

The Anopheles gambiae genome sequence has been analyzed to find ATP-binding cassette protein genes based on deduced protein similarity to known family members. A nonredundant collection of 44 putative genes was identified including five genes not detected by the original Anopheles genome project machine annotation. These genes encode at least one member of all the human and Drosophila melanogaster ATP-binding protein subgroups. Like D. melanogaster, A. gambiae has subgroup ABCH genes encoding proteins different from the ABC proteins found in other complex organisms. The largest Anopheles subgroup is the ABCC genes which includes one member that can potentially encode ten different isoforms of the protein by differential splicing. As with Drosophila, the second largest Anopheles group is the ABCG subgroup with 12 genes compared to 15 genes in D. melanogaster, but only 5 genes in the human genome. In contrast, fewer ABCA and ABCB genes were identified in the mosquito genome than in the human or Drosophila genomes. Gene duplication is very evident in the Anopheles ABC genes with two groups of four genes, one group with three genes and three groups with two head to tail duplicated genes. These characteristics argue that the A. gambiae is actively using gene duplication as a mechanism to drive genetic variation in this important gene group.

SAGE analysis of mosquito salivary gland transcriptomes during Plasmodium invasion.

Invasion of the vector salivary glands by Plasmodium is a critical step for malaria transmission. To describe salivary gland cellular responses to sporozoite invasion, we have undertaken the analysis of Anopheles gambiae salivary gland transcriptome using Serial Analysis of Gene Expression (SAGE). Statistical analysis of the more than 160000 sequenced tags generated from four libraries, two from glands infected by Plasmodium berghei, two from glands of controls, revealed that at least 57 Anopheles genes are differentially expressed in infected salivary glands. Among the 37 immune-related genes identified by SAGE tags, four (Defensin1, GNBP, Serpin6 and Cecropin2) were found to be upregulated during salivary gland invasion, while five genes encoding small secreted proteins display induction patterns strongly reminiscent of that of Cecropin2. Invasion by Plasmodium has also an impact on the expression of genes involved in transport, lipid and energy metabolism, suggesting that the sporozoite may exploit the metabolism of its host. In contrast, protein composition of saliva is predicted to be only slightly modified after infection. This study, which is the first transcriptome analysis of the salivary gland response to Plasmodium infection, provides a basis for a better understanding of Plasmodium/Anopheles salivary gland interactions.

Pilot Anopheles gambiae full-length cDNA study: sequencing and initial characterization of 35,575 clones.

We describe the preliminary analysis of over 35,000 clones from a full-length enriched cDNA library from the malaria mosquito vector Anopheles gambiae. The clones define nearly 3,700 genes, of which around 2,600 significantly improve current gene definitions. An additional 17% of the genes were not previously annotated, suggesting that an equal percentage may be missing from the current Anopheles genome annotation.

Assessing the Drosophila melanogaster and Anopheles gambiae genome annotations using genome-wide sequence comparisons.

We performed genome-wide sequence comparisons at the protein coding level between the genome sequences of Drosophila melanogaster and Anopheles gambiae. Such comparisons detect evolutionarily conserved regions (ecores) that can be used for a qualitative and quantitative evaluation of the available annotations of both genomes. They also provide novel candidate features for annotation. The percentage of ecores mapping outside annotations in the A. gambiae genome is about fourfold higher than in D. melanogaster. The A. gambiae genome assembly also contains a high proportion of duplicated ecores, possibly resulting from artefactual sequence duplications in the genome assembly. The occurrence of 4063 ecores in the D. melanogaster genome outside annotations suggests that some genes are not yet or only partially annotated. The present work illustrates the power of comparative genomics approaches towards an exhaustive and accurate establishment of gene models and gene catalogues in insect genomes.