# Serial Analysis of Gene Expression in Plasmodium berghei salivary gland sporozoites.

Isabelle Rosinski-Chupin, Thomas Chertemps, Bertrand Boisson, Sylvie Perrot, Emmanuel Bischoff, Jérôme Briolay, Pierre Couble, Robert Ménard, Paul Brey, Patricia Baldacci

► To cite this version:

Isabelle Rosinski-Chupin, Thomas Chertemps, Bertrand Boisson, Sylvie Perrot, Emmanuel Bischoff, et al.. Serial Analysis of Gene Expression in Plasmodium berghei salivary gland sporozoites.. BMC Genomics, 2007, 8, pp.466. 10.1186/1471-2164-8-466. hal-00294035

HAL Id: hal-00294035

https://hal.science/hal-00294035v1

Submitted on 5 Jun 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Research article

Open Access

# Serial Analysis of Gene Expression in *Plasmodium berghei* salivary gland sporozoites

Isabelle Rosinski-Chupin*¹, Thomas Chertemps¹, Bertrand Boisson², Sylvie Perrot¹, Emmanuel Bischoff³, Jérôme Briolay⁴, Pierre Double⁵, Robert Ménard², Paul Brey¹ and Patricia Baldacci*²

Address: ¹Biochimie et Biologie Moléculaire des Insectes, Institut Pasteur, 28 rue du Dr Roux, 75724, Paris cedex 15, France, ²Biologie et Génétique du Paludisme, Institut Pasteur, 28 rue du Dr Roux, 75724, Paris cedex 15, France, ³Plateforme 2: Puces à ADN, Institut Pasteur, 28 rue du Dr Roux, 75724, Paris cedex 15, France, ⁴DTAMB, Université de Lyon; IFR 41; 43 boulevard du 11 Novembre 1918, 69622 Villeurbanne cedex, France and ⁵Centre de Génétique Moléculaire et Cellulaire, UMR 5534, C.N.R.S.-Université de Lyon, 69622 Villeurbanne cedex, France

Email: Isabelle Rosinski-Chupin* - ichupin@pasteur.fr; Thomas Chertemps - tcher@pasteur.fr; Bertrand Boisson - boisson@pasteur.fr; Sylvie Perrot - sperrot@pasteur.fr; Emmanuel Bischoff - bischoff@pasteur.fr; Jérôme Briolay - briolay@biomserv.univ-lyon1.fr; Pierre Double - Pierre.Couple@univ-lyon1.fr; Robert Ménard - rmenard@pasteur.fr; Paul Brey - pbrey@pasteur.fr; Patricia Baldacci* - baldacci@pasteur.fr

Published: 19 December 2007

BMC Genomics 2007, 8:466 doi:10.1186/1471-2164-8-466

This article is available from: http://www.biomedcentral.com/1471-2164/8/466

© 2007 Rosinski-Chupin et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: 18 April 2007

Accepted: 19 December 2007

## Abstract

**Background:** The invasion of *Anopheles* salivary glands by *Plasmodium* sporozoites is an essential step for transmission of the parasite to the vertebrate host. Salivary gland sporozoites undergo a developmental programme to express genes required for their journey from the site of the mosquito bite to the liver and subsequent invasion of, and development within, hepatocytes. A Serial Analysis of Gene Expression was performed on *Anopheles gambiae* salivary glands infected or not with *Plasmodium berghei* and we report here the analysis of the *Plasmodium* sporozoite transcriptome.

**Results:** Annotation of 530 tag sequences homologous to *Plasmodium berghei* genomic sequences identified 123 genes expressed in salivary gland sporozoites and these genes were classified according to their transcript abundance. A subset of these genes was further studied by quantitative PCR to determine their expression profiles. This revealed that sporozoites modulate their RNA amounts not only between the midgut and salivary glands, but also during their storage within the latter. Among the 123 genes, the expression of 66 is described for the first time in sporozoites of rodent *Plasmodium* species.

**Conclusion:** These novel sporozoite expressed genes, especially those expressed at high levels in salivary gland sporozoites, are likely to play a role in *Plasmodium* infectivity in the mammalian host.

## Background

The sporozoite is the stage of the malaria parasite transmitted from the mosquito vector to the mammalian host. The success of infection depends on the sporozoites leav-

ing their site of inoculation in the dermis, rapidly attaining the liver, invading hepatocytes, and developing therein [1-4]. This results in the release of thousands of merozoites from the infected hepatocytes that subse-

<footnote>* Corresponding authors</footnote>
quently invade red blood cells, causing the malaria disease [5]. Sporozoites are formed within oocysts of the mosquito midgut and are initially poorly infectious. They migrate to mosquito salivary glands (SG) and must undergo a developmental programme, with concomitant changes in gene expression, in order to become highly infectious to the mammalian host. These SG sporozoites exhibit a circular gliding movement and in certain conditions can elicit a strong protective immune response in the mammalian host [6-8].

Recent technological advances have improved our understanding, at the molecular level, of the *Plasmodium* parasite including the less well known sporozoite stage. The completion of the *Plasmodium falciparum* (*P. falciparum*) genome sequence allowed the analysis of gene expression at different stages of the parasite life cycle with microarrays [9]. The data showed that approximately 2000 genes are expressed in SG sporozoites, 500 of which are expressed at high levels and over a 100 of these are not expressed at significant levels in blood stages. In addition, proteomic analyses detected 274 proteins in *P. falciparum* sporozoites [10,11].

*Plasmodium* species infecting rodents are powerful laboratory models as they are more easily amenable to genetic studies and their genomes have also been sequenced [11,12]. However, the transcriptome of sporozoites from these species is less well known and has been obtained mainly from *Plasmodium yoelii* (*P. yoelii*) cDNA libraries [13] and Suppressive Subtractive Hybridisation studies [14,15]. In addition, 119 proteins have been identified in *Plasmodium berghei* (*P. berghei*) sporozoites, of which 34 are sporozoite specific (PlasmoDB release 5.1).

We recently reported a SAGE aimed at characterising *Anopheles gambiae* (*A. gambiae*) SG genes that are differentially expressed in response to infection with *P. berghei* sporozoites [16]. In that study 530 tag sequences, found exclusively in libraries from infected mosquitoes, were identical to *P. berghei* genomic sequences, and 41 of these were readily annotated. Considering the data presently available for the transcriptome of *P. berghei* sporozoites, we decided that it would be worthwhile to annotate the remaining tags. SAGE allows gene-expression-profiling based on the quantification of short 14–15 nucleotide (nt) sequence tags, each sequence being, in theory, associated with the transcript of a single gene. It provides an overall estimation of the abundance of transcripts, while requiring no *a priori* information about the sequence of the transcripts to be studied [17,18]. However, only short cDNA sequences, usually located in the 3'UTR, are obtained and consequently the attribution of a tag sequence to a gene is highly dependent on the quality of annotation of the genome of the organism studied and

available cDNA or EST studies [19,20]. On the other hand, positioning tags on the genomic sequence can provide information on both the orientation of a transcript and the length of the 3' UTR.

The tags from our SAGE were annotated by combining information on their position in the *P. berghei* genomic sequences, predicted gene models and ESTs from *P. berghei* and *P. yoelii*. We unambiguously identified 123 genes expressed in *P. berghei* sporozoites, of which 66 are detected for the first time in rodent *Plasmodium* species. The hierarchical classification of these transcripts according to the abundance of the tags was confirmed by qPCR and the characterisation of the gene structure and/or gene expression was undertaken for some. Finally, our results provide evidence that mRNA levels may vary not only between midgut and SG sporozoites but also during storage of the sporozoites in the SG.

## Results

### Annotation of Plasmodium berghei tags: rationale

SAGE libraries were constructed from four different SG RNA preparations: SG isolated from *A. gambiae* mosquitoes 14, or 18, days after a blood meal on *P. berghei* infected mice and SG isolated 14, or 15–18, days after feeding on uninfected mice [16]. Sequence analysis showed that among the tags identified in the infected SG libraries, 17181 were absent from uninfected SG libraries and were thus considered as potential *P. berghei* sequences (see Figure 1). These sequences were then screened using the following criteria: 1) to be found more than once in the cumulative libraries, 2) to give a single hit in the *P. berghei* genome, 3) to be derived from the most 3' *Nla*III site predicted in the annotated gene sequence in the sense orientation ("primary" tags in Figure 1) or to be found within 500 nt of the stop codon of an annotated gene in the sense orientation.

We found 2154 tag sequences at least twice in the cumulative libraries and 530 of these gave one or more perfect matches by BLASTN comparison on *P. berghei* genomic sequences and annotated genes. These tag sequences and their annotations are provided in Additional file 1. We discarded from further analysis 162 tags that gave multiple hits in the genome due to the impossibility of assigning them to a unique gene. Of the 368 tags that gave a single hit in the genome, 166 were found within annotated genes and 62 of these were derived from the most 3' *Nla*III predicted site.

Sampling of the 3' untranslated regions of about 30 *Plasmodium* genes showed that 77% of them were 500 nt in length (I.C. unpublished observation). Therefore, we analysed in more detail the position in the genome of the 202 tags that were not found within annotated genes. Sixty-
<figure><img src="image_2.png" /><figcaption>**Figure 1**
Flow chart of the rationale used to select Tag sequences from the SAGE library. Shaded boxes correspond to tag sequences retained, clear boxes correspond to tag sequences rejected.</figcaption></figure>

three tag sequences, found within 500 nt downstream of the stop codon of an annotated gene and in the sense orientation, were retained in the analysis as it seemed likely that they came from the 3' end of the adjacent genes. As two genes were identified with two different tag sequences respectively, only one tag was retained for each gene, resulting in a total of 61 (see Figure 1).

In total, 123 identifying tag sequences fulfilled all our criteria (see Table 1 and Additional file 2). BLASTN analysis of the 123 genes identified by these tags showed that 55 aligned with known transcripts from *P. yoelii* or *P. berghei* sporozoites, whereas the remaining 68 did not. Among the latter, proteomic approaches have detected

PB000876.00.0 in sporozoites and a study by Ishino *et al.*[21] showed that PB000892.00.0 (*Pbs36*) is also present. It should be noted that the criteria used have resulted in the elimination of the majority of tag sequences, including several tags corresponding to genes known to be expressed in SG sporozoites (see Additional file 1). For example, two tags aligned with PB001026.00.0, coding for circumsporozoite protein (CS). The most abundant tag, found 77 times, gave 4 hits in the genome and the second tag sequence was in the antisense orientation. Although it is likely that the most abundant tag originated from the CS transcript, it was not taken into account in order to be consistent in the analysis. Another example concerned two tags aligning with PB000892.03.0 (*UIS3*), a gene known to be upregulated in sporozoites [14]. Again the most abundant tag gave multiple hits in the genome and the second was in the antisense orientation.

Taken together this SAGE analysis has identified 66 novel sporozoite expressed genes (indicated as *SIS* on Table 1).

### Validation of SAGE data

SAGE, like microarrays, is designed to give quantitative information on gene expression and the interpretation of the results depends on correct gene identification. This identification is not always straightforward in *P. berghei* since genomic sequence clusters are generally shorter than in *P. falciparum*; this results in the prediction of a large number of truncated genes, often lacking 5' or 3' ends. In addition, there are less EST ressources.

As a first step towards confirming the gene expression in sporozoites and validating tag annotation, we obtained more accurate data on gene structure and organization, by clustering ESTs and comparing with *P. yoelii* and *P. chabaudi* orthologous sequences. In Table 1, 15 genes (indicated by *) have been manually reannotated resulting in a different structure and a longer ORF than that predicted in the databases. These structures have been confirmed for 13 genes by RT-PCR on sporozoite RNA (see Additional file 2). For genes that are split between two genomic clusters, the intervening fragment was cloned and sequenced (see Additional file 2).

To confirm the identification of the 66 novel genes expressed in salivary gland sporozoites, RT-PCR experiments were performed. PCR fragments were obtained for all of them indicating that these genes are truly expressed (Additional file 3).

The number of times a tag sequence is identified in a SAGE library is expected to correlate with the relative abundance of the steady state mRNA [17]. To determine whether the SAGE data correctly reflects transcript abundance
<table><caption>Table 1: List of 123 genes identified by SAGE in P. berghei salivary gland sporozoites</caption><thead><tr><th rowspan="2"></th><th rowspan="2">Tag ID</th><th colspan="3">times found</th><th rowspan="2">GeneDB ID</th><th rowspan="2">Known as</th><th rowspan="2">Sporozoite expression</th><th rowspan="2">Pf orthologue</th><th rowspan="2">Predicted Pb/Pf protein domain and/or function</th></tr><tr><th>d14</th><th>d18</th><th>total</th></tr></thead><tbody><tr><td rowspan="5">group 1</td><td>tag 335</td><td>60</td><td>39</td><td>99</td><td>PB402615.00.0*</td><td></td><td>d</td><td>Pfgenefinder_63r</td><td></td></tr><tr><td>tag 183</td><td>6</td><td>89</td><td>95</td><td>PB100551.00.0</td><td>UIS4</td><td>a, d</td><td>--</td><td>SP, TM</td></tr><tr><td>tag 190</td><td>9</td><td>59</td><td>68</td><td>PB107461.00.0</td><td>UIS7</td><td>d</td><td>--</td><td>TM</td></tr><tr><td>tag 269</td><td>6</td><td>25</td><td>31</td><td>PB000374.03.0</td><td>TRAP</td><td>d, I, P</td><td>PF13_0201</td><td>SP, vWA-like, TSP1</td></tr><tr><td>tag 226</td><td>11</td><td>13</td><td>24</td><td>PB001063.00.0</td><td>S23</td><td>d, I</td><td>PF08_0088</td><td>Pb-reticulocyte binding, SP, TM</td></tr><tr><td rowspan="10">group 2</td><td>tag 435</td><td>10</td><td>7</td><td>17</td><td>PB000252.01.0</td><td>S13, SPECT2</td><td>d, Pb, I, P</td><td>PFD0430c</td><td>SP, MAC/Perforin domain</td></tr><tr><td>tag 168</td><td>4</td><td>12</td><td>16</td><td>PB000317.00.0*</td><td></td><td>d, I, P</td><td>PF08_0008</td><td>SP, TM, GPI</td></tr><tr><td>tag 146</td><td>8</td><td>8</td><td>16</td><td>PB000881.01.0</td><td>S21/PbTRSP</td><td>d, I, P</td><td>PFA0200w</td><td>TM, TSP1</td></tr><tr><td>tag 112</td><td>11</td><td>2</td><td>13</td><td>PB000649.01.0</td><td>ECPI</td><td>Δ, I, P</td><td>PFB0325c</td><td>Papain family cysteine protease</td></tr><tr><td>tag 110</td><td>2</td><td>11</td><td>13</td><td>PB104017.00.0</td><td>S12</td><td>d</td><td>--</td><td>SP</td></tr><tr><td>tag 200</td><td>2</td><td>9</td><td>11</td><td>PB000400.01.0*</td><td>SIS1</td><td>Δ, I</td><td>PFE0395c</td><td>Pf38, Sexual stage antigen, SP, TM, GPI</td></tr><tr><td>tag 75</td><td>3</td><td>8</td><td>11</td><td>PB000484.00.0</td><td>UIS10/PbPL</td><td>a</td><td>PFF1420w</td><td>Phosphatidylcholine-sterol acyltransferase</td></tr><tr><td>tag 50</td><td>4</td><td>6</td><td>10</td><td>PB000891.00.0</td><td>Pbs36p</td><td>d, I</td><td>PFD0215c</td><td>Pf52, Sexual stage antigen, SP</td></tr><tr><td rowspan="24">group 3</td><td>tag 76</td><td>1</td><td>8</td><td>9</td><td>PB104912.00.0*</td><td></td><td>d, I</td><td>PF11_0283</td><td>SP, TM</td></tr><tr><td>tag 264</td><td>3</td><td>5</td><td>8</td><td>AB125370</td><td>SPECT</td><td>d, I, P</td><td>MAL13P1.212</td><td>SP</td></tr><tr><td>tag 278</td><td>4</td><td>4</td><td>8</td><td>PB000757.02.0</td><td>SIS2</td><td>Δ, I, P</td><td>PF14_0425</td><td>Fructose-1, 6-bisphosphate aldolase</td></tr><tr><td>tag 319</td><td>1</td><td>7</td><td>8</td><td>PB108190.00.0*</td><td>UIS1</td><td>a</td><td>--</td><td>Serine/Threonine protein kinases, 9 TM</td></tr><tr><td>tag 136</td><td>2</td><td>5</td><td>7</td><td>PB100963.00.0</td><td></td><td>d</td><td>--</td><td>Guanylyl cyclase, TM</td></tr><tr><td>tag 245</td><td>0</td><td>6</td><td>6</td><td>PB000780.01.0*</td><td></td><td>d, I</td><td>PF11230c</td><td>4 TM, asparagine-rich</td></tr><tr><td>tag 530</td><td>0</td><td>5</td><td>5</td><td>PB000264.00.0*</td><td></td><td>d</td><td>PF14_0044</td><td>SP</td></tr><tr><td>tag 118</td><td>0</td><td>5</td><td>5</td><td>PB000411.03.0*</td><td>SIS3</td><td>Δ</td><td>--</td><td></td></tr><tr><td>tag 294</td><td>1</td><td>4</td><td>5</td><td>PB000464.00.0</td><td>S25</td><td>d, I, P</td><td>PFL0545w</td><td>Kinesin-related protein</td></tr><tr><td>tag 280</td><td>0</td><td>5</td><td>5</td><td>PB000817.02.0</td><td></td><td>d, I, P</td><td>PF08_0054</td><td>Heat shock 70 kDa protein</td></tr><tr><td>tag 268</td><td>1</td><td>4</td><td>5</td><td>PB001070.03.0</td><td></td><td>d</td><td>PFF0685c</td><td>Fe-only hydrogenase</td></tr><tr><td>tag 193</td><td>1</td><td>4</td><td>5</td><td>PB101172.00.0</td><td></td><td>d</td><td>--</td><td></td></tr><tr><td>tag 92</td><td>2</td><td>2</td><td>4</td><td>PB000161.03.0</td><td></td><td>d, I</td><td>PF10955w</td><td>sugar transporter, 6TM</td></tr><tr><td>tag 201</td><td>0</td><td>4</td><td>4</td><td>PB000787.03.0</td><td>SIS4</td><td>Δ, I</td><td>PF10320w</td><td>Arginase family</td></tr><tr><td>tag 101</td><td>0</td><td>4</td><td>4</td><td>PB105090.00.0</td><td>SIS5</td><td>Δ</td><td>--</td><td>TM</td></tr><tr><td>tag 233</td><td>1</td><td>3</td><td>4</td><td>PB106293.00.0*</td><td></td><td>d</td><td>MAL13P1.365</td><td>SNARE</td></tr><tr><td>tag188</td><td>2</td><td>2</td><td>4</td><td>PB107843.00.0*</td><td>SIS6</td><td>Δ, I</td><td>PF10_0290</td><td>TM</td></tr><tr><td>tag 62</td><td>3</td><td>1</td><td>4</td><td>PB108391.00.0</td><td></td><td>d, I, P</td><td>PF08_0072</td><td>TM</td></tr><tr><td>tag 276</td><td>1</td><td>3</td><td>4</td><td>PB000632.00.0</td><td>SIS7</td><td>Δ, I</td><td>PF14_0518</td><td>nifU protein, putative</td></tr><tr><td>tag 354</td><td>3</td><td>0</td><td>3</td><td>PB000175.03.0*</td><td>SIS8</td><td>Δ, I</td><td>PFD0915w</td><td>TM</td></tr><tr><td>tag 171</td><td>1</td><td>2</td><td>3</td><td>PB000416.02.0</td><td></td><td>d, I</td><td>PFD0825c</td><td>RNA-binding protein of pumilio/ mpt5 family</td></tr><tr><td>tag 326</td><td>2</td><td>1</td><td>3</td><td>PB000422.00.0</td><td></td><td>d, I, P</td><td>PF07_0126</td><td></td></tr><tr><td>tag 213</td><td>0</td><td>3</td><td>3</td><td>PB000468.02.0</td><td>SIS9</td><td>Δ, I</td><td>PFL0335c</td><td>Eukaryotic translation initiation factor 5</td></tr><tr><td>tag 197</td><td>0</td><td>3</td><td>3</td><td>PB000626.00.0</td><td></td><td>d, I</td><td>PF14_0224</td><td>PP1-like protein serine/threonine phosphatase</td></tr><tr><td>tag 68</td><td>1</td><td>2</td><td>3</td><td>PB000674.01.0</td><td>SIS10</td><td>Δ, I, P</td><td>PF13_0058</td><td>RNA recognition motif</td></tr><tr><td>tag 170</td><td>0</td><td>3</td><td>3</td><td>PB000876.00.0</td><td></td><td>Pb, I, P</td><td>MAL8P1.78</td><td>Small heat shock protein</td></tr><tr><td>tag 236</td><td>0</td><td>3</td><td>3</td><td>PB000892.00.0</td><td>Pbs 36</td><td>Δ, I</td><td>PFD0210c</td><td>Pbs36, Sexual stage antigen, PS</td></tr><tr><td>tag 20</td><td>0</td><td>3</td><td>3</td><td>PB000911.00.0*</td><td>SIS11</td><td>Δ, I</td><td>PFE0965c</td><td>Vacuolar ATP synthetase, 4 TM</td></tr><tr><td>tag 8</td><td>1</td><td>2</td><td>3</td><td>PB001101.02.0</td><td>SIS12</td><td>ΔΔ</td><td>PF14_0395</td><td></td></tr><tr><td>tag 119</td><td>0</td><td>3</td><td>3</td><td>PB001143.02.0</td><td>SIS13</td><td>Δ, I, P</td><td>PF14_0083</td><td>Ribosomal protein S8e</td></tr><tr><td>tag 122</td><td>1</td><td>2</td><td>3</td><td>PB001247.02.0</td><td>S14</td><td>d, I, P</td><td>PFL0370w</td><td></td></tr><tr><td>tag 127</td><td>0</td><td>3</td><td>3</td><td>PB001432.02.0</td><td>SIS14</td><td>Δ, I</td><td>PFE0340c</td><td>Peptidase S54, rhomboid, 7 TM</td></tr><tr><td>tag 155</td><td>3</td><td>0</td><td>3</td><td>PB001575.02.0</td><td>SIS15</td><td>Δ, I</td><td>PFB0270w</td><td>SP, TM, Fe-S metabolism associated domain</td></tr><tr><td>tag 150</td><td>0</td><td>3</td><td>3</td><td>PB104165.00.0</td><td></td><td>d, I</td><td>PF10675w</td><td></td></tr><tr><td>tag 305</td><td>0</td><td>3</td><td>3</td><td>PB104461.00.0</td><td></td><td>d</td><td>--</td><td>SP, TM</td></tr></tbody></table>
<table><caption>Table 1: List of 123 genes identified by SAGE in P. berghei salivary gland sporozoites (Continued)</caption><thead><tr><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th></tr></thead><tbody><tr><td>tag 350</td><td>0</td><td>3</td><td>3</td><td>PB105607.00.0</td><td>d</td><td>--</td><td>2 TM</td></tr><tr><td>tag 106</td><td>1</td><td>2</td><td>3</td><td>PB105732.00.0</td><td>d</td><td>--</td><td>TM</td></tr><tr><td>tag 196</td><td>0</td><td>3</td><td>3</td><td>PB106033.00.0</td><td>SIS16</td><td>Δ, I</td><td>PF10_0296</td></tr><tr><td>tag 41</td><td>1</td><td>2</td><td>3</td><td>PB107193.00.0*</td><td>SIS17</td><td>ΔΔ</td><td>PFF1195c</td></tr><tr><td>tag 121</td><td>1</td><td>2</td><td>3</td><td>PB108098.00.0</td><td>d</td><td>--</td><td>TM</td></tr><tr><td>tag 157</td><td>2</td><td>1</td><td>3</td><td>PB300784.00.0</td><td>d, I</td><td>MAL8P1.154</td><td></td></tr><tr><td>tag 362</td><td>0</td><td>2</td><td>2</td><td>PY01277**</td><td>SIS18</td><td>ΔΔ</td><td>PFL2080c</td></tr><tr><td>tag 142</td><td>0</td><td>2</td><td>2</td><td>PB000006.00.0</td><td>SIS19</td><td>Δ, I</td><td>PFE0995c</td></tr><tr><td>tag 265</td><td>0</td><td>2</td><td>2</td><td>PB000064.03.0</td><td>SIS20</td><td>ΔΔ</td><td>PFB0510w</td></tr><tr><td>tag 143</td><td>0</td><td>2</td><td>2</td><td>PB000092.02.0</td><td>UIS2</td><td>a, I, P</td><td>PF14_0614</td></tr><tr><td>tag 355</td><td>1</td><td>1</td><td>2</td><td>PB000094.03.0</td><td>SIS21</td><td>Δ, I, P</td><td>MAL12P1.180</td></tr><tr><td>tag 15</td><td>1</td><td>1</td><td>2</td><td>PB000111.01.0</td><td>d, I, P</td><td></td><td>MAL12P1.63</td></tr><tr><td>tag 57</td><td>1</td><td>1</td><td>2</td><td>PB000120.01.0</td><td>SIS22</td><td>Δ, P</td><td>PFC1065w</td></tr><tr><td>tag 26</td><td>1</td><td>1</td><td>2</td><td>PB000133.00.0</td><td>SIS23</td><td>Δ</td><td>--</td></tr><tr><td>tag 59</td><td>0</td><td>2</td><td>2</td><td>PB000171.01.0</td><td>SIS24</td><td>Δ</td><td>--</td></tr><tr><td>tag 339</td><td>0</td><td>2</td><td>2</td><td>PB000178.02.0</td><td>SIS25</td><td>Δ, I</td><td>PFE0810c</td></tr><tr><td>tag 139</td><td>0</td><td>2</td><td>2</td><td>PB000207.01.0</td><td>SIS26</td><td>Δ, I</td><td>PFL1050w</td></tr><tr><td>tag 221</td><td>0</td><td>2</td><td>2</td><td>PB000297.01.0</td><td>SIS27</td><td>Δ, I</td><td>PF11_0251</td></tr><tr><td>tag 12</td><td>0</td><td>2</td><td>2</td><td>PB000307.01.0</td><td>UIS9-7</td><td>a, I</td><td>PFE0935c</td></tr><tr><td>tag 134</td><td>0</td><td>2</td><td>2</td><td>PB000352.00.0</td><td>SIS28</td><td>Δ, I</td><td>PF11_0150</td></tr><tr><td>tag 51</td><td>2</td><td>0</td><td>2</td><td>PB000352.02.0</td><td>d, I, P</td><td>PFL2440w</td><td>DNA repair protein rhp16</td></tr><tr><td>tag 317</td><td>1</td><td>1</td><td>2</td><td>PB000354.03.0</td><td>d, I</td><td>PF10_0283</td><td>TM</td></tr><tr><td>tag 285</td><td>0</td><td>2</td><td>2</td><td>PB000396.01.0</td><td>SIS29</td><td>Δ, I</td><td>PFC0330w</td></tr><tr><td>tag 292</td><td>0</td><td>2</td><td>2</td><td>PB000412.02.0</td><td>SIS30</td><td>Δ, I</td><td>PF10750c</td></tr><tr><td>tag 17</td><td>0</td><td>2</td><td>2</td><td>PB000424.02.0</td><td>d, I</td><td>PF11_0358</td><td>DNA-directed RNA polymerase, beta subunit</td></tr><tr><td>tag 320</td><td>0</td><td>2</td><td>2</td><td>PB000450.03.0</td><td>SIS31</td><td>Δ, I</td><td>MAL13P1.148</td></tr><tr><td>tag 54</td><td>0</td><td>2</td><td>2</td><td>PB000530.00.0</td><td>d, I</td><td>PFE0400w</td><td>myosin</td></tr><tr><td>tag 277</td><td>0</td><td>2</td><td>2</td><td>PB000553.03.0</td><td>d, I</td><td>PF13_0161</td><td>ankyrin</td></tr><tr><td>tag 303</td><td>1</td><td>1</td><td>2</td><td>PB000591.00.0</td><td>SIS32</td><td>ΔΔ</td><td>PFC0125w</td></tr><tr><td>tag 184</td><td>0</td><td>2</td><td>2</td><td>PB000599.03.0</td><td>SIS33</td><td>Δ, I</td><td>PF14_0681</td></tr><tr><td>tag 34</td><td>0</td><td>2</td><td>2</td><td>PB000604.01.0</td><td>d</td><td>PFF0840w</td><td></td></tr><tr><td>tag 87</td><td>0</td><td>2</td><td>2</td><td>PB000610.02.0</td><td>SIS34</td><td>Δ, I</td><td>PFL2390c</td></tr><tr><td>tag 204</td><td>0</td><td>2</td><td>2</td><td>PB000624.01.0</td><td>SIS35</td><td>Δ, I</td><td>PF13_0261</td></tr><tr><td>tag 128</td><td>0</td><td>2</td><td>2</td><td>PB000650.00.0</td><td>SIS36</td><td>Δ, I, P</td><td>PF14_0446</td></tr><tr><td>tag 58</td><td>0</td><td>2</td><td>2</td><td>PB000675.02.0</td><td>SIS37</td><td>ΔΔ</td><td>PFF1450w</td></tr><tr><td>tag 144</td><td>1</td><td>1</td><td>2</td><td>PB000709.00.0</td><td>SIS38</td><td>Δ, I</td><td>PF14_0181</td></tr><tr><td>tag 289</td><td>2</td><td>0</td><td>2</td><td>PB000798.02.0</td><td>MAEBL</td><td>d, I, P</td><td>PF11_0486</td></tr><tr><td>tag 88</td><td>0</td><td>2</td><td>2</td><td>PB000807.01.0</td><td>SIS39</td><td>Δ, I</td><td>PFD0485w</td></tr><tr><td>tag 353</td><td>0</td><td>2</td><td>2</td><td>PB000811.01.0</td><td>d, I</td><td></td><td>MAL7P1.125</td></tr><tr><td>tag 80</td><td>1</td><td>1</td><td>2</td><td>PB000837.02.0</td><td>SIS40</td><td>Δ, I</td><td>PF14_0695</td></tr><tr><td>tag 116</td><td>0</td><td>2</td><td>2</td><td>PB000857.01.0</td><td>SIS41</td><td>Δ, I, P</td><td>PF11_0183</td></tr><tr><td>tag 238</td><td>0</td><td>2</td><td>2</td><td>PB000874.02.0</td><td>SIS42</td><td>Δ, I</td><td>MAL8P1.142</td></tr><tr><td>tag 79</td><td>0</td><td>2</td><td>2</td><td>PB000903.01.0</td><td>SIS43</td><td>Δ, I</td><td>PF11_0197</td></tr><tr><td>tag 227</td><td>0</td><td>2</td><td>2</td><td>PB000904.01.0</td><td>SIS44</td><td>Δ, I</td><td>PF11_0196</td></tr><tr><td>tag 206</td><td>0</td><td>2</td><td>2</td><td>PB000927.03.0</td><td>SIS45</td><td>Δ, I</td><td>MAL12P1.113</td></tr><tr><td>tag 210</td><td>0</td><td>2</td><td>2</td><td>PB000941.01.0</td><td>SIS46</td><td>Δ, I</td><td>PF10_0096</td></tr><tr><td>tag 349</td><td>0</td><td>2</td><td>2</td><td>PB000994.02.0*</td><td>d, I</td><td>PFL1690w</td><td>Octapeptide repeat</td></tr><tr><td>tag 218</td><td>0</td><td>2</td><td>2</td><td>PB001066.02.0</td><td>d, I</td><td>PF10_0140</td><td>SPRY</td></tr><tr><td>tag 228</td><td>0</td><td>2</td><td>2</td><td>PB001080.02.0</td><td>SIS47</td><td>Δ, I</td><td>PFL0870w</td></tr><tr><td>tag 267</td><td>1</td><td>1</td><td>2</td><td>PB001118.02.0</td><td>UIS4-6</td><td>d, I</td><td>PFA0170c</td></tr><tr><td>tag 27</td><td>0</td><td>2</td><td>2</td><td>PB001120.00.0</td><td>SIS48</td><td>Δ, I</td><td>MAL13P1.150</td></tr><tr><td>tag 6</td><td>0</td><td>2</td><td>2</td><td>PB001130.00.0</td><td>SIS49</td><td>Δ, I, P</td><td>MAL8P1.133</td></tr><tr><td>tag 124</td><td>1</td><td>1</td><td>2</td><td>PB001157.02.0</td><td>SIS50</td><td>Δ, I</td><td>MAL8P1.67</td></tr><tr><td>tag 301</td><td>0</td><td>2</td><td>2</td><td>PB001178.02.0</td><td>SIS51</td><td>Δ, I</td><td>PFD0300w</td></tr><tr><td>tag 123</td><td>0</td><td>2</td><td>2</td><td>PB001253.00.0</td><td>SIS52</td><td>Δ, I</td><td>PFD0360w</td></tr><tr><td>tag 331</td><td>0</td><td>2</td><td>2</td><td>PB001305.00.0</td><td>SIS53</td><td>Δ, I, P</td><td>PF11_0168</td></tr></tbody></table>
<table><caption>Table 1: List of 123 genes identified by SAGE in P. berghei salivary gland sporozoites (Continued)</caption><thead><tr><th>tag</th><th>70</th><th>1</th><th>1</th><th>2</th><th>PB001374.02.0</th><th>d, 1, P</th><th>MAL13P1.304</th><th>Malaria antigen, metalloendopeptidase</th></tr></thead><tbody><tr><td>tag 31</td><td>0</td><td>2</td><td>2</td><td>PB001507.02.0</td><td>SIS54</td><td>ΔΔ</td><td>PFB0620w</td><td></td></tr><tr><td>tag 312</td><td>0</td><td>2</td><td>2</td><td>PB001517.02.0</td><td></td><td>d</td><td>PFF0800w</td><td>SP, TM, vWFA, TSP1</td></tr><tr><td>tag 43</td><td>1</td><td>1</td><td>2</td><td>PB100341.00.0</td><td>SIS55</td><td>Δ</td><td>--</td><td></td></tr><tr><td>tag 93</td><td>1</td><td>1</td><td>2</td><td>PB101083.00.0</td><td></td><td>d</td><td>--</td><td></td></tr><tr><td>tag 61</td><td>2</td><td>0</td><td>2</td><td>PB101148.00.0</td><td>SIS56</td><td>Δ</td><td>--</td><td></td></tr><tr><td>tag 114</td><td>2</td><td>0</td><td>2</td><td>PB101408.00.0</td><td>SIS57</td><td>Δ, 1</td><td>PFE1260c</td><td>TM</td></tr><tr><td>tag 56</td><td>2</td><td>0</td><td>2</td><td>PB101596.00.0</td><td>SIS58</td><td>Δ</td><td>--</td><td></td></tr><tr><td>tag 342</td><td>0</td><td>2</td><td>2</td><td>PB102211.00.0</td><td>SIS59</td><td>Δ</td><td>--</td><td></td></tr><tr><td>tag 318</td><td>1</td><td>1</td><td>2</td><td>PB102277.00.0</td><td></td><td>d</td><td>--</td><td></td></tr><tr><td>tag 71</td><td>0</td><td>2</td><td>2</td><td>PB103980.00.0</td><td>SIS60</td><td>Δ</td><td>--</td><td>2 TM</td></tr><tr><td>tag 271</td><td>0</td><td>2</td><td>2</td><td>PB105487.00.0</td><td></td><td>d</td><td>--</td><td>2 TM</td></tr><tr><td>tag 35</td><td>0</td><td>2</td><td>2</td><td>PB106484.00.0</td><td>SIS61</td><td>Δ, 1</td><td>PF14_0449</td><td></td></tr><tr><td>tag 208</td><td>0</td><td>2</td><td>2</td><td>PB108086.00.0</td><td>SIS62</td><td>Δ</td><td>--</td><td></td></tr><tr><td>tag 151</td><td>0</td><td>2</td><td>2</td><td>PB000618.02.0</td><td>SIS63</td><td>Δ</td><td>PFD0725c</td><td>arsenical pump-driving ATPase, putative</td></tr><tr><td>tag 179</td><td>1</td><td>1</td><td>2</td><td>PB300567.00.0</td><td>SIS64</td><td>ΔΔ</td><td>PFL2420w</td><td></td></tr><tr><td>tag 85</td><td>0</td><td>2</td><td>2</td><td>PB300953.00.0</td><td>SIS65</td><td>ΔΔ</td><td>MAL8P1.121</td><td></td></tr><tr><td>tag 33</td><td>1</td><td>1</td><td>2</td><td>PB301273.00.0</td><td></td><td>d, 1</td><td>PFL1695c</td><td></td></tr><tr><td>tag 47</td><td>0</td><td>2</td><td>2</td><td>PB000460.01.0</td><td>SIS66</td><td>Δ, 1</td><td>PFC0750w</td><td>SP</td></tr></tbody></table>

(a) : mRNA detected in P. berghei sporozoites [14].

(d) : mRNA detected in P. yoelii sporozoites [13].

Pb : detected by proteomics in P. berghei sporozoites [11].

P : detected by proteomics in P. falciparum sporozoites [10].

I : detected in P. falciparum sporozoites with microarrays [9].

SP : signal peptide.

TM : transmembrane domain.

** : orthologue of a P. yoelii gene for which a P. berghei gene is not yet annotated.

SIS : new *Sporozoite expressed gene* Identified by SAGE, this work.

Δ : no rodent Plasmodium sporozoite ESTs or proteomics data.

ΔΔ : no sporozoite ESTs or detection by microarrays or proteomics in *Plasmodium* species.

dances, we selected eighteen known or novel genes, predicted to be expressed at high or low levels, and quantified their RNA by qPCR in sporozoites isolated from SG at d18 of infection. The values were normalised to the reference gene PB001026.00.0 (CS) and plotted against the number of times the identifying tag was found in the SAGE library at d18 of infection (Figure 2). A good correlation (R² = 0.8) was found and was even greater (R² = 0.96, not shown), when a second abundant tag sequence for two genes, UIS4 and TRAP, was taken into account. Thus, the number of times a tag was found in our SAGE data correctly reflects the gene expression levels in sporozoites.

## Hierarchical classification of sporozoite expressed genes

In Table 1, the 123 genes unambiguously identified by this SAGE as being expressed in SG sporozoites have been classified into three groups according to the number of times the identifying tag sequence was found in the combined d14 and d18 libraries (with the caveat that for some genes there may be an underestimation due to the criteria used above).

### Group I

The first group of highly expressed genes, defined by tags found more than 20 times in the combined libraries, con-

tains five genes of which four have been described previously: PB000374.03.0, which codes for TRAP, a major sporozoite protein having an essential role in sporozoite motility [22]; PB100551.00.0 and PB107461.00.0, which code for UIS4 and UIS7 respectively, and were identified by SSH as genes whose expression is upregulated in SG sporozoites compared to midgut sporozoites [14]; PB001063.00.0, also known as S23, which codes for a reticulocyte binding protein and was identified by SSH as a gene upregulated in SG sporozoites compared to blood stages [15].

Surprisingly, the gene PB402615.00.0, for which the identifying tag (tag 335) was found the most often, has not been described previously. The gene sequence aligns with numerous ESTs of *P. yoelii* sporozoites and parasites developing in the absence of host cells [13,23]. To further characterize the gene structure, a manual clustering of these ESTs was performed. Several tags were found along this cluster, both in sense and antisense orientations and, as expected, the most 3' tag was also the most abundant (Figure 3A). The other internal tags, one of which is highly represented (tag 189, see additional file 1), may be due to either alternative polyadenylation events or priming at internal polyA stretches during cDNA synthesis, as
<figure><img src="image_3.png" /><figcaption>**Figure 2**
Correlation between the levels of gene expression as determined by qPCR and the number of times the identifying tag sequence was found in the SAGE library. R² is the coefficient of correlation.</figcaption></figure>

hypothesized by others [24]. To rule out the possibility that the cluster covers two different genes, a RT-PCR was performed using primers at the 5' and 3' end of the cluster and this detected a unique transcript of the predicted size (data not shown). Based on our results and annotation, we propose that the ORF predicted for PB402615.00.0 during automatic annotation is incorrect and that this gene codes for a 17 KDa protein, rich in tyrosine and basic amino acids. Furthermore, this new annotation indicates that there are orthologous sequences in *P. yoelii* (PY02432), *P. chabaudi* (PC400629.00.0) and also in *P. falciparum* (PFgenefinder_63r.) (Figure 3B).

Interestingly, PB402615.00.0 is not the only gene characterized by several abundant tags. For instance, two tags are characteristic of PB100551.00.0, encoding UIS4: the most abundant (tag 186, found 95 times) is located at the end of the ORF, while the other (tag 153, found 74 times), is located 650 nt downstream and probably defines a transcript with a longer 3' UTR. Such observations may be useful when defining the gene structure with the complete potential regulatory sequences.

### Group 2
The second group contains eight moderately expressed genes defined by tags found 10 to 20 times. Among these, there is one *UIS* gene (*UIS10*/PB000484.00.0) and three *S* genes *S13/SPECT2* (PB000252.01.0), *S21/PbTRSP* (PB000881.01.0) and *S12* (PB104017.00.0). *UIS10/PbPL* codes for a lecithine-cholesterol acyl-transferase involved in cell traversal [25]. The *S13/SPECT2* protein is characterized by a MAC/Perforin domain and was found to be essential for membrane attack during cell invasion [26]. The protein *S21/PbTRSP*, which contains a thrombospondin type 1 domain, has recently been shown to have a role

in host cell invasion [27]. Finally, *S12* codes for a protein with a predicted signal peptide, but whose function has yet to be characterized.

The gene PB000649.01.0 codes for the cysteine protease ECP1, involved in egress of sporozoites from oocysts [28]. Two genes, PB000891.00.0 (*Pbs36p*) and PB000400.01.0, encode potential GPI- anchored proteins of the Pfs48/45 family, orthologues of *P. falciparum* Pf52 and Pfs38 proteins [29]. Proteins of this family, characterized by a 6-Cys domain, are involved in intra-hepatic development and/ or gametocyte fertilization [21,30,31]. Pfs38 expression was previously observed in *P. falciparum* sporozoites [9] but not in rodent sporozoites. We have named PB000400.01.0, *SIS1*, for Sporozoite expressed gene Identified by SAGE1 (see Table 1). Finally, PB000317.00.0, whose structure has been manually revised, codes for a potential secreted and GPI-anchored protein and therefore constitutes a novel candidate that could be involved in host cell recognition/invasion.

### Group 3
Of the remaining 110 genes, expressed at low levels and defined by tags found less than ten times, 43 aligned with *P. yoelii* or *P. berghei* sporozoite ESTs, one was detected in *P. berghei* sporozoites by proteomics and PB000892.00.0/*Pbs36* has been described in sporozoites [21]. Thus 65 genes in this group are shown for the first time to be expressed in this stage of the *P. berghei* parasite (indicated as *SIS2-66* on Table 1). Orthologous *P. falciparum* genes are predicted for 90 genes in this group, and expression in sporozoites has been detected for 74 by microarrays (of which 20 also by proteomics) and one by proteomics only. The 20 *P. berghei* predicted genes for which there are no obvious *P. falciparum* orthologues could correspond to incomplete genes and require more precise annotation. Alternatively, they may be highly divergent genes or genes whose prediction has been missed during *P. falciparum* genome annotation.

Nine genes in this group have been described previously: AB125370 (*SPECT*) and PB000892.00.0 (encoding Pbs36) are important for liver infectivity [21,32] and PB000798.02.0, encoding MAEBL, has a role in the initial binding of sporozoites to the mosquito SG [33]. Six genes, *UIS1*, *UIS2*, *UIS9-7/UIS29-1*, *UIS4-6*, *S25*, *S14*, were already identified as being expressed in SG sporozoites [14,15].

In addition, two other genes could have a potential role in adhesion/invasion of host cells: PB001517.02.0 codes for a protein with a Thrombospondin 1 domain and a von Willebrand type A domain, and PB001080.02.0, which encodes a protein defined as thrombospondin-related 3 (also found as AF375983 in Genbank). Proteins involved
<figure><img src="image_4.png" /><figcaption>**Figure 3**</figcaption></figure>

New annotation of PB402615.00.0. A) Continuous black line corresponds to the *P. berghei* genomic DNA contig; the position of tags are indicated by small arrows (→: sense tag); (←: antisense tag); thick black arrow: PB402615.00.0 as annotated in PlasmoDB; thin black arrow: EST cluster; thick white arrow: ORF of *P. yoelii* orthologue; thick grey arrow: new annotation of PB402615.00.0 ORF. B) Alignment of protein sequences of PB402615.00.0 orthologues in *P. yoelii*, *P. chabaudi* and *P. falciparum*.

in the molecular motor needed for sporozoite motility are also represented by aldolase (PB000757.02.0), which provides the link between TRAP and myosinA [34,35], and by a kinesin-related protein (PB000464.00.0/S25). While myosinA, which is thought to play a major role in the motility of Apicomplexa zoites, is absent from our description (with a tag found only once), another gene encoding a potential myosin is identified (PB000450.03.0). Two genes encoding proteases of the rhomboid family (PB001432.02.0/SIS14 and PB000352.00.0/SIS28) are also found. Interestingly, these proteins are the orthologues of PfROM4 and PfROM1, respectively, which are able to cleave adhesins, such as TRAP, AMA1, MAEBL and others, that are involved in interactions with host-cell receptors [36].

The expression of only two genes encoding ribosomal proteins has been detected, confirming, as previously reported [13], their under-representation in the sporozoite transcriptome. Interestingly, two genes coding for proteins with RNA-binding domains are identified, one of which has a pumilio domain suggesting a role in the negative regulation of translation [37]. These proteins were

previously described as being upregulated in *P. falciparum* gametocytes and sporozoites [9] and may be involved in the regulation of sporozoite protein expression.

Other genes found in this group code for proteins with diverse biological functions for example energy metabolism, signal transduction and protein modification. Interestingly, a gene encoding a putative sugar transporter (PB000161.03.0), whose *P. falciparum* orthologue is expressed specifically in sporozoites, is present [9]. We previously reported an increase in mRNA levels of an *Anopheles* sugar transporter in infected SG [16]. These two transporters may play a role in providing the sporozoite with sufficient energy resources for its journey from the oocysts to the SG and from the bite site to the liver in the vertebrate host. The presence of three genes coding for proteins involved in iron-sulphur cluster formation and iron homeostasis can also be noted. Iron-sulphur cluster formation is essential for a wide variety of processes, including facilitation of electron transfer in oxidative phosphorylation and enzymatic activities in mitochondria, cytoplasm or nucleus as well as sensing of intracellular iron and/or oxidant levels. Expression of these genes
may be preparing the sporozoite for high mitochondrial activities related to motility and/or for the future iron-rich blood environment.

## Expression profiles of sporozoite expressed genes

To further characterize a number of genes that were not previously studied, namely PB402615.00.0 (group1), PB000317.00.0 and PB000400.01.0 (group2), PB001080.02.0 (Thrombospondin-related 3 protein), PB000757.02.0 (aldolase), PB001432.02.0 (rhomboid protein), PB000787.03.0 (arginase), PB000468.02.0 (eukaryotic translation initiation factor 5) and PB000675.02.0 (sec14-like protein) (group 3), we compared their expression profiles with those of known transcripts such as CS, *UIS1*, *UIS3*, *UIS4*, *UIS7*, *TRAP*, *S23*, *Pbs36*, *Pbs36p*, *MAEBL* and *SPECT2*. The RNA levels were quantified by qPCR in ookinetes, midgut and SG sporozoites and blood stages and the results were normalised with respect to *hsp70* which is expressed at all stages (Figure 4).

All genes showed higher levels of expression in midgut and/or SG sporozoites compared to ookinete or blood stages, with the exception of PB000787.03.0 whose expression is clearly not sporozoite specific. Other genes, for example PB000400.01.0, PB000757.02.0 and PB001432.02.0 are easily detectable in blood stages, suggesting they may also have a role at this stage. Finally, PB001432.02.0, PB000317.00.0, PB001080.02.0 and PB000400.01.0 are expressed, like CS, at relatively high levels in the ookinete stage, although this does not mean that the protein is produced.

<figure><img src="image_5.png" /><figcaption>**Figure 4** Histogram representation of qPCR analysis of gene expression in mixed blood stages, ookinetes, d18 midgut sporozoites and d18 salivary gland sporozoites. X axis: genes tested; Y axis: log scale of mean normalised expression.</figcaption></figure>

As observed in Table 1, the number of times a tag is found in the SG libraries at d18 of infection often differs to that at d14. This is due in part to the higher number of tags sequenced in the d18 library (2 fold) and to the increase (2 fold) in the number of sporozoites inside the glands at d18. However, the tags of some genes, for example *UIS4*, increase up to 15 fold between d14 and d18 whereas they decrease for others. This suggested that there might be variations in gene expression during sporozoite storage in the glands. To further characterize these variations, the relative levels of gene expression for the same panel of genes as above were determined in sporozoites isolated from SG and midguts of *A. gambiae* mosquitoes at d14 or d18 of infection. The values were calculated from the geometric averages and normalised to the geometric mean of PB001026.00.0 (CS), as this gene was determined to be the best reference using GeNORM software.

The results show important increases (4–66 times) in RNA amounts in SG compared to midgut sporozoites for *UIS4*, *UIS1*, *UIS3*, *SPECT2*, *UIS7* and PB402615.00.0 at d18 (see Table 2). These results were expected for the *UIS* genes as they were found in an SSH library between SG and midgut sporozoites [14]. They also suggest that PB402615.00.0 may be important for preparing the sporozoites for infection. On the contrary, there is significantly less *MAEBL* RNA in SG compared to midgut sporozoites in agreement with *MAEBL* having a role in sporozoite adhesion to mosquito salivary glands [33].

In addition, the ratios obtained between d18 and d14 SG sporozoites indicate that there is a significant increase in expression over this period for *UIS4*, *UIS1*, *UIS3* and *UIS7*, whereas there is no substantial change in RNA quantities for the other genes. Similar increases in expression were also seen for *UIS3* and *UIS4* in *A. stephensi* infected salivary glands suggesting that *P. berghei* sporozoites develop similarly in *A. gambiae* and *A. stephensi* (not shown).

Taken together, the qPCR data show that sporozoites are capable of modulating their RNA amounts between the midgut and salivary glands, as well as during their storage within the latter.

## Discussion

Although the transcriptional repertoire of *Plasmodium* sporozoites has been investigated in several laboratories using different techniques, including cDNA libraries, SSH, microarrays and proteomics, only the microarray data has provided information concerning the level of gene expression [9]. We have obtained new data concerning the transcriptional repertoire of *P. berghei* sporozoites using SAGE on *A. gambiae* infected salivary glands. SAGE does not require *a priori* knowledge of the sequence of genes to be
<table><caption>Table 2: Variation of gene expression in *P. berghei* sporozoites developing in *A. gambiae*.</caption><thead><tr><th rowspan="2">Gene ID</th><th rowspan="2">Gene name</th><th colspan="2">between tissues</th><th colspan="2">with time</th></tr><tr><th>ratios at d14 SG/midgut</th><th>ratios at d18 SG/midgut</th><th>ratio SG d18/d14</th><th>ratio midgut d18/d14</th></tr></thead><tbody><tr><td>PB100551.00.0</td><td>UIS4</td><td>37.07*</td><td>66.31**</td><td>4.17*</td><td>2.33*</td></tr><tr><td>PB108190.00.0</td><td>UIS1</td><td>20.72*</td><td>56.52**</td><td>2.15*</td><td>0.79*</td></tr><tr><td>PB000892.03.0</td><td>UIS3</td><td>20.88*</td><td>54.45**</td><td>3.77*</td><td>1.45*</td></tr><tr><td>PB402615.00.0</td><td></td><td>8.22*</td><td>10.67**</td><td>1.57*</td><td>1.21*</td></tr><tr><td>PB000252.01.0</td><td>SPECT2</td><td>7.55*</td><td>7.69**</td><td>1.24*</td><td>1.22*</td></tr><tr><td>PB107461.00.0</td><td>UIS7</td><td>5.81*</td><td>4.59**</td><td>2.55*</td><td>3.24*</td></tr><tr><td>PB000891.00.0</td><td>Pbs36p</td><td>4.67*</td><td>3.79**</td><td>0.77*</td><td>0.95*</td></tr><tr><td>PB000787.03.0</td><td></td><td>1.13*</td><td>2.32**</td><td>1.54*</td><td>0.75*</td></tr><tr><td>PB000675.02.0</td><td></td><td>2.50*</td><td>2.13**</td><td>0.76*</td><td>0.90*</td></tr><tr><td>PB000400.01.0</td><td></td><td>7.52*</td><td>1.97**</td><td>0.49*</td><td>1.88*</td></tr><tr><td>PB000892.00.0</td><td>Pbs36</td><td>1.54*</td><td>1.89**</td><td>1.28*</td><td>1.04*</td></tr><tr><td>PB000374.03.0</td><td>TRAP</td><td>1.57*</td><td>1.88**</td><td>1.22*</td><td>1.01*</td></tr><tr><td>PB000757.02.0</td><td></td><td>1.71*</td><td>1.64**</td><td>0.71*</td><td>0.74*</td></tr><tr><td>PB001432.02.0</td><td></td><td>1.54*</td><td>1.35**</td><td>0.75*</td><td>0.86*</td></tr><tr><td>PB001063.00.0</td><td>S23</td><td>1.76*</td><td>1.00**</td><td>0.80*</td><td>1.41*</td></tr><tr><td>PB000317.00.0</td><td></td><td>1.12*</td><td>0.87**</td><td>0.77*</td><td>0.99*</td></tr><tr><td>PB001080.02.0</td><td></td><td>1.20*</td><td>0.70**</td><td>0.75*</td><td>1.29*</td></tr><tr><td>PB000468.02.0</td><td></td><td>0.52*</td><td>0.44**</td><td>0.90*</td><td>1.05*</td></tr><tr><td>PB000798.02.0</td><td>MAEBL</td><td>0.04*</td><td>0.04**</td><td>0.87*</td><td>0.87*</td></tr></tbody></table>

analysed and provides quantification of gene expression by the number of times a tag sequence is obtained [17]. However, since the sequence and annotation of the *P. berghei* genome is incomplete, the attribution of a tag sequence to a gene was not straightforward.

Several criteria were applied before retaining a tag sequence as a gene identifier: it was found twice or more in the cumulative SAGE libraries, the BLASTN analysis gave a unique hit in the genome, which was in the sense orientation of an annotated gene and was located at the most 3' predicted NlaIII site or within 500 nt of the stop codon of a neighbouring gene. These combined criteria resulted in the unambiguous identification of 123 genes expressed in SG sporozoites, of which 16 were already known to be sporozoite-expressed genes: *TRAP*, several *UIS* and *S* genes, *SPECT*, *Pbs36p*, *Pbs36* and *MAEBL*. It should be noted that this list of 123 *P. berghei* SG sporozoite-expressed genes is not exhaustive.

Of the 2154 unique tag sequences found at least twice in the libraries from infected mosquitoes, only 25% (530) matched perfectly to the present version of the *P. berghei* genome. The 1624 remaining tags may derive from *A. gambiae* cDNAs. Indeed, 587 match on *A. gambiae* ESTs and 943 on the *A. gambiae* genome (see Additional file 1, sheet 2). Alternatively, they may derive from genes not yet sequenced in the *P. berghei* genome. In addition, the

absence of matches of some tags could be due to polymorphisms between the *P. berghei* strain sequenced, ANKA, and the NK65 strain used in our SAGE.

Also, 202 (55%) of the tag sequences that gave a single hit in the *P. berghei* genome were not found within an annotated gene. Upon more detailed analysis of their position in the contigs, 63 were found to be within 500 nt of the stop codon of an annotated gene, in the sense orientation and were considered to derive from transcripts of these genes. However, we were unable to place 139 tag sequences either within an annotated ORF or within 500 nt of a stop codon. These tags should help in refining the present annotation of the *P. berghei* genome and may, in the future, be formally proven as corresponding to sporozoite-expressed genes.

Among the 368 tag sequences that gave a single match in the *P. berghei* genome, 64 were found in the antisense orientation, either within annotated genes (56) or within 500 nt of the stop codon of an annotated gene (8). Antisense RNAs have been described previously in *P. falciparum* and they are suggested to be involved in transcriptional regulation [38-40]. At the present time, there is insufficient *P. berghei* cDNA data and no microarrays using sense and antisense probes, to establish whether or not these tags correspond to antisense RNAs or to transcripts from a gene on the opposite strand.

<footnote>* :p = 0.05 ** :p = 0.01</footnote>
The stringent criteria described above resulted in the elim-
ination of the vast majority of tag sequences, retaining
0.7% (123/17181) of the unique tag sequences found
only in the infected libraries, representing 3% (668/
21721) abundance. Among the tag sequences rejected
there are a number that matched with genes known to be
expressed in sporozoites, such as CS, UIS3, RON4 and
AMA-1. Consequently, other genes expressed in sporo-
zoites may have been eliminated. Furthermore, due to the
SAGE cloning procedure, genes for which the transcripts
have no *NlalIII* site will also be missing. The data supplied
in the supplemental tables should be useful to other
researchers interested in genes expressed in sporozoites.

Amongst the *P. berghei* *SIS* genes identified in this study as
being expressed in sporozoites, eighteen are predicted to
encode proteins with one or more transmembrane
regions and/or a signal peptide sequence, suggesting that
they are membrane associated or secreted proteins. Other
genes of potential interest are PB000464.00.0/S25, a
kinesin-related protein, which could be involved in motil-
ity; PB000416.02.0, a putative RNA-binding protein of
the pumilio/mpt5 family known for their role in repres-
sion of gene expression; PB000903.01.0/SIS43 which
contains an ankyrin repeat suggesting a role in protein-
protein interaction; PB000650.00.0/SIS36, which con-
tains a fasciclin domain, may have a role in cell adhesion.
Finally, PB402615.00.0, for which the identifying tag was
the most abundant in this SAGE analysis, is annotated as
a hypothetical protein and aligns with *P. yoelii* ESTs from
sporozoites as well as parasites developing in the absence
of host cells. Our qPCR data show that this gene is
expressed in sporozoites but not in ookinetes or blood
stages and that there is a substantial increase in the
amount of RNA for this gene between midgut and salivary
gland sporozoites. This differential regulation between
organs and the mRNA abundance in SG suggest a role for
this gene in sporozoite infectivity in the mammalian host.

Several properties (motility, infectivity, etc) differ
between midgut and SG sporozoites, but it is not known
whether these developmental changes are time and/or tis-
sue dependent nor which signalling factors are involved.
Interestingly, the qPCR data presented here indicate that,
at least for *UIS3* and *UIS4*, not only the tissue (midgut ver-
sus salivary gland) but also the time spent in the SG signif-
icantly influences the level of expression of individual
genes. This change in expression of genes that are essential
for development in the liver, is in agreement with the
increase in infectivity of SG sporozoites between d14 and
d18 post blood meal [41].

Since SAGE provides a quantitative read out of gene
expression, and as our qPCR analysis confirmed this, the
123 genes were classified into 3 groups according to the

number of times the tag was found in the libraries. Among
the thirteen genes presenting the most abundant tags
/groups 1 and 2), there are seven, *UIS4*, *TRAP*, *S13/*
*SPECT2*, *S21/Pb TRSP*, *ECP1*, *UIS10/Pb PL*, and *Pbs36p*,
that have been shown via gene knockout experiments to
be essential for the sporozoite or liver stage development
[7,21,22,25-28]. Furthermore, among the genes with less
abundant tags there are three, *SPECT*, *Pbs36* and *MAEBL*,
which have also been shown to have essential roles
[21,32,33]. This indicates that our approach has identified
several genes known to be essential for sporozoites and
points to additional genes that may be required at this
stage. It will be of interest to inactivate the novel sporo-
zoite-expressed genes identified in this paper to define
their function in the parasite. Mutants that are defective in
their development in the mammalian host would be of
particular interest as they could provide new tools to
probe the host immune response to *Plasmodium* infection.

Conclusion

The SAGE described here has lead to the identification of
66 novel genes expressed in *P. berghei* sporozoites. These
novel sporozoite expressed genes, especially those
expressed at high levels in salivary gland sporozoites, are
likely to play a role in *Plasmodium* infectivity in the mam-
malian host.

Methods

Mosquito Infections

A. gambiae (Yaounde strain) and *Anopheles stephensi* (A. stephensi) (Sda500 strain) mosquitoes were reared at the Centre for Production and Infection of Anopheles at the Pasteur Institute under standard conditions [42]. 3-4 day old female mosquitoes were fed on anaesthetised *P. berghei* infected mice. A recombinant *P. berghei* strain expressing GFP in late midgut stages and SG of mosquitoes was used [43].

Construction of SAGE libraries

The construction of the SAGE libraries has been described
elsewhere [16]. Briefly, 2 libraries were constructed from
SG of infected *A. gambiae* mosquitoes: one 14 days (early
stage of infection) and the other 18 days (stage at which
sporozoites are considered to be fully infectious) after
feeding on *P. berghei* infected mice. Two control libraries
were made from SG of mosquitoes 14 and 15-18 days
post-feeding on uninfected mice.

Analysis of Tag Sequences

The analysis of tag sequences has been described in detail
in [16]. Tag sequences found only in libraries from the
infected mosquitoes were considered as potential *P*
*berghei* sequences. BLASTN comparisons were performed
against the library of *P. berghei* contigs and annotated
genes obtained from the ftp site of the Sanger Institute *P*.
berghei genome project (version 12 August 2005). Sequences of genes unambiguously identified by tags were aligned to EST_others (i. e. non-human, non-mouse EST sequences) at the National Center for Biotechnology Information via the *Plasmodium* genome consortium PlasmoDB.

## Real-time Quantitative RT-PCR

Real time qPCR was performed on cDNA preparations using the SYBR green detection system and the ABI Prism 7900 sequence detector (Applied Biosystems) according to manufacturer's instructions. Midgut and SG sporozoites were isolated from *A. gambiae* mosquitoes infected with *P. berghei* at d14 and d18 after the blood meal. Under these conditions, the mean number of sporozoites isolated from midguts was 5500 ± 200 (d14) and 7500 ± 1700 (d18), and the mean number of sporozoites from SG was 1500 ± 20 (d14) and 3900 ± 1200 (d18). SG sporozoites were also isolated from *A. stephensi* mosquitoes infected with *P. berghei* with slightly higher mean numbers of sporozoites per gland: 3300 ± 250 on d14 and 17000 ± 1800 on d18. RNA was extracted with TriReagent, DNAs treated, and cDNAs synthesised with Superscript II reverse transcriptase (Invitrogen) using random primers. Three independent RNA preparations were made for each sample. cDNAs of *P. berghei* ookinetes and asynchronous blood stages were prepared as described in [44]. PCR conditions were 1 cycle of 50°C for 2 min, followed by 95°C for 10 min, then 40 cycles of (95°C 15 sec, 60°C 1 min). qPCR was performed in triplicate. The standard curve was analysed for all primers and gave amplification efficiencies of 90-100%. Data was analysed with SDS 2.1 software and with geNORM Visual Basic application for Microsoft_Excel as described by Vandesompele et al. to determine the best reference gene for the qPCR [45]. This approach relies on the principle that the expression of a perfect reference gene should be identical in all samples, independent of experimental conditions or cell type. For the comparison between midgut and salivary gland samples the Circumsporozoite gene (CS), PB001026.00.0, was determined to be the best reference gene and the expression of genes of interest were subsequently normalized to the geometric mean of this reference. The normalized gene expression was then calculated with the Q-Gene software [46]. For the qPCR analysis between different tissues the gene coding for hsp70 was found to be the best reference. The list of primers used for qPCR analysis is provided in Additional file 4.

## RT-PCR

cDNA samples were prepared as for qPCR except that half of the mRNA was not treated with reverse transcriptase. PCR was performed with 1 unit of Platinum Taq polymerase, in manufacturer's buffer (Invitrogen) and 1.5 mM MgCl₂. PCR conditions were : 1 cycle at 94°C for 2 min,

40 cycles of (94°C for 20 sec, 55°C or 52°C for 35 sec, 68°C for 1 min) and 1 cycle at 68°C for 10 min. Primers used for RT-PCR are shown in Additional file 4.

## Abbreviations

SAGE: Serial Analysis of Gene Expression

SG: salivary glands

cDNA: complementary DNA

ESTs: expressed sequence tags

BLAST: Basic Local Alignment Search Tool

ORF: open reading frame

qPCR: quantitative polymerase chain reaction

KDa: kilodaltons

## Authors' contributions

IR-C was responsible for the conception of the SAGE, performed experiments, analysed data and wrote the manuscript.

TC, BB, SP and EB performed experiments, acquired and interpreted data and provided tables, figures and revised the manuscript.

JB and PC constructed the SAGE libraries and revised the manuscript.

RM and PBrey critically revised the manuscript

PBaldacci analysed the data and wrote the paper.

All authors read and approved the final manuscript.

## Additional material

### Additional file 1

List and annotation of the 530 tag sequences with hits in the *P. berghei* genome (sheet 1). BLASTN analysis against *A. gambiae* genome and EST sequences of the tags found specifically in infected salivary gland libraries that did not align with *P. berghei* genome (sheet 2).

Click here for file

[http://www.biomedcentral.com/content/supplementary/1471-2164-8-466-S1.xls]
### Additional file 2

Compilation of information in databases with BLASTN alignments and annotation of the 123 P. berghei genes with their putative function, predicted molecular weight and domains, expression data and P. falciparum orthologues.

Click here for file

[http://www.biomedcentral.com/content/supplementary/1471-2164-8-466-S2.xls]

### Additional file 3

Confirmation by RT-PCR of the expression of 66 SIS genes in salivary gland sporozoites. The SIS genes tested are indicated above the wells on cDNA samples (+RT) and controls for DNA contamination (-RT). The size marker is 100 nt Smart ladder (Eurogentec). G refers to gel number as indicated in Additional file 4.

Click here for file

[http://www.biomedcentral.com/content/supplementary/1471-2164-8-466-S3.PDF]

### Additional file 4

Primers used for annotation and qPCR. Sheet1: List of primers used for annotation and results of RT-PCR analysis. Some PCR fragments have been sequenced and their sequences submitted to EMBL database (accession numbers are indicated in columns L and M). Column "E" refers to the gel in Additional file 3.

Click here for file

[http://www.biomedcentral.com/content/supplementary/1471-2164-8-466-S4.xls]

## Acknowledgements

We thank the CEPIA for mosquito rearing and Jean Sautereau for infections. The work was supported by grants from the Institut Pasteur (Grand Programme Horizontal Anophèle), CNRS-PAL+, Research Program on Parasitic Diseases Sanofi-Aventis/Ministry for Research, the European Commission (FP6 BioMalPar Network of Excellence) and Howard Hughes Medical Institute. T. C. and B. B. are post-doctoral fellows funded by the Research Program on Parasitic Diseases Sanofi-Aventis/Ministry for Research. E. B. is funded by the Délégation Générale à l'Armement.

## References

1. Baldacci P, Menard R: The elusive malaria sporozoite in the mammalian host. *Mol Microbiol* 2004, **54**:298-306.

2. Frevert U, Engelmann S, Zougbede S, Stange J, Ng B, Matuschewski K, Liebes L, Yee H: *Intravital observation of Plasmodium berghei sporozoite infection of the liver*. *PLoS Biol* 2005, **3**:e192.

3. Amino R, Thiberge S, Martin B, Celli S, Shorte S, Frischknecht F, Menard R: *Quantitative imaging of Plasmodium transmission from mosquito to mammal*. *Nat Med* 2006, **12**:220-224.

4. Tarun AS, Baer K, Dumpit RF, Gray S, Lejarcegui N, Frevert U, Kappe SH: *Quantitative isolation and in vivo imaging of malaria parasite liver stages*. *Int J Parasitol* 2006.

5. Sturm A, Amino R, van de SC, Regen T, Retzlaff S, Rennenberg A, Krueger A, Pollok JM, Menard R, Heussler VT: *Manipulation of host hepatocytes by the malaria parasite for delivery into liver sinusoids*. *Science* 2006, **313**:1287-1290.

6. Nussenzweig R, Vanderberg J, Most H: *Protective immunity produced by the injection of x-irradiated sporozoites of Plasmodium berghei*. **IV. Dose response, specificity and humoral immunity**. *Mil Med* 1969, **134**:1176-1182.

7. Mueller AK, Camargo N, Kaiser K, Andorfer C, Frevert U, Matuschewski K, Kappe SH: *Plasmodium liver stage developmental arrest by depletion of a protein at the parasite-host interface*. *Proc Natl Acad Sci U S A* 2005, **102**:3022-3027.

8. Mueller AK, Labaied M, Kappe SH, Matuschewski K: *Genetically modified Plasmodium parasites as a protective experimental malaria vaccine*. *Nature* 2005, **433**:164-167.

9. Le Roch KG, Zhou Y, Blair PL, Grainger M, Moch JK, Haynes JD, De La Vega P, Holder AA, Batalov S, Carucci DJ, Winzeler EA: *Discovery of gene function by expression profiling of the malaria parasite life cycle*. *Science* 2003, **301**:1503-1508.

10. Florens L, Washburn MP, Raine JD, Anthony RM, Grainger M, Haynes JD, Moch JK, Muster N, Sacci JB, Tabb DL, Witney AA, Wolters D, Wu Y, Gardner MJ, Holder AA, Sinden RE, Yates JR, Carucci DJ: *A proteomic view of the Plasmodium falciparum life cycle*. *Nature* 2002, **419**:520-526.

11. Hall N, Karras M, Raine JD, Carlton JM, Kooij TW, Berriman M, Florens L, Janssen CS, Pain A, Christophides GK, James K, Rutherford K, Harris B, Harris C, Churcher C, Quail MA, Ormond D, Doggett J, Trueman HE, Mendoza J, Bidwell SL, Rajandream MA, Carucci DJ, Yates JR 3rd, Kafatos FC, Janse CJ, Barrell B, Turner CM, Waters AP, Sinden RE: *A comprehensive survey of the Plasmodium life cycle by genomic, transcriptomic, and proteomic analyses*. *Science* 2005, **307**:82-86.

12. Carlton JM, Angiuoli SV, Suh BB, Kooij TW, Pertea M, Silva JC,ERMolaeva MD, Allen JE, Selengut JD, Koo HL, Peterson JD, Pop M, Kosack DS, Shumway MF, Bidwell SL, Shallom SJ, van Aken SE, Riedmuller SB, Feldblyum TV, Cho JK, Quackenbush J, Sedegah M, Shoaibi A, Cummings LM, Florens L, Yates JR, Raine JD, Sinden RE, Harris MA, Cunningham DA, Preiser PR, Bergman LW, Vaidya AB, van Lin LH, Janse CJ, Waters AP, Smith HO, White OR, Salzberg SL, Venter JC, Fraser CM, Hoffman SL, Gardner MJ, Carucci DJ: *Genome sequence and comparative analysis of the model rodent malaria parasite Plasmodium yoelii yoelii*. *Nature* 2002, **419**:512-519.

13. Kappe SH, Gardner MJ, Brown SM, Ross J, Matuschewski K, Ribeiro JM, Adams JH, Quackenbush J, Cho J, Carucci DJ, Hoffman SL, Nussenzweig V: *Exploring the transcriptome of the malaria sporozoite stage*. *Proc Natl Acad Sci U S A* 2001, **98**:9895-9900.

14. Matuschewski K, Ross J, Brown SM, Kaiser K, Nussenzweig V, Kappe SH: *Infectivity-associated changes in the transcriptional repertoire of the malaria parasite sporozoite stage*. *J Biol Chem* 2002, **277**:41948-41953.

15. Kaiser K, Matuschewski K, Camargo N, Ross J, Kappe SH: *Differential transcriptome profiling identifies Plasmodium genes encoding pre-erythrocytic stage-specific proteins*. *Mol Microbiol* 2004, **51**:1221-1232.

16. Rosinski-Chupin I, Briolay J, Brouilly P, Perrot S, Gomez SM, Chem temps T, Roth CW, Keime C, Gandrillon O, Double P, Brey PT: *SAGE analysis of mosquito salivary gland transcriptomes during Plasmodium invasion*. *Cell Microbiol* 2007, **9**:708-724.

17. Velculescu VE, Zhang L, Vogelstein B, Kinzler KW: *Serial analysis of gene expression*. *Science* 1995, **270**:484-487.

18. Knox DP, Skuce PJ: *SAGE and the quantitative analysis of gene expression in parasites*. *Trends Parasitol* 2005, **21**:322-326.

19. Velculescu VE, Vogelstein B, Kinzler KW: *Analysing uncharted transcriptomes with SAGE*. *Trends Genet* 2000, **16**:423-425.

20. Pleasance ED, Marra MA, Jones SJ: *Assessment of SAGE in transcript identification*. *Genome Res* 2003, **13**:1203-1215.

21. Ishino T, Chinzei Y, Yuda M: *Two proteins with 6-cys motifs are required for malarial parasites to commit to infection of the hepatocyte*. *Mol Microbiol* 2005, **58**:1264-1275.

22. Sultan AA, Thathy V, Frevert U, Robson KJ, Crisanti A, Nussenzweig V, Nussenzweig RS, Menard R: *TRAP is necessary for gliding motility and infectivity of plasmodium sporozoites*. *Cell* 1997, **90**:511-522.

23. Wang Q, Brown S, Roos DS, Nussenzweig V, Bhanot P: *Transcriptome of axenic liver stages of Plasmodium yoelii*. *Mol Biochem Parasitol* 2004, **137**:161-168.

24. Munasinghe A, Patankar S, Cook BP, Madden SL, Martin RK, Kyle DE, Shoaibi A, Cummings LM, Wirth DF: *Serial analysis of gene expression (SAGE) in Plasmodium falciparum: application of the technique to A-T rich genomes*. *Mol Biochem Parasitol* 2001, **113**:23-34.

25. Bhanot P, Schauer K, Coppens I, Nussenzweig V: *A surface phospholipase is involved in the migration of plasmodium sporozoites through cells*. *J Biol Chem* 2005, **280**:6752-6760.

26. Ishino T, Chinzei Y, Yuda M: *A Plasmodium sporozoite protein with a membrane attack complex domain is required for breaching the liver sinusoidal cell layer prior to hepatocyte infection*. *Cell Microbiol* 2005, **7**:199-208.
27. Labaied M, Camargo N, Kappe SH: Depletion of the Plasmodium berghei thrombospondin-related sporozoite protein reveals a role in host cell entry by sporozoites. *Mol Biochem Parasitol* 2007, **153**:158-166.

28. Aly AS, Matuschewski K: A malarial cysteine protease is necessary for Plasmodium sporozoite egress from oocysts. *J Exp Med* 2005, **202**:225-230.

29. Thompson J, Janse CJ, Waters AP: Comparative genomics in Plasmodium: a tool for the identification of genes and functional analysis. *Mol Biochem Parasitol* 2001, **118**:147-154.

30. van Dijk MR, Douradinha B, Franke-Fayard B, Heussler V, van Dooren MW, van Schaijk B, van Gemert GJ, Sauerwein RW, Mota MM, Waters AP, Janse CJ: Genetically attenuated, P36p-deficient malarial sporozoites induce protective immunity and apoptosis of infected liver cells. *Proc Natl Acad Sci U S A* 2005, **102**:12194-12199.

31. van Dijk MR, Janse CJ, Thompson J, Waters AP, Braks JA, Dodemont HJ, Stunnenberg HG, van Gemert GJ, Sauerwein RW, Eling W: A central role for P48/45 in malaria parasite male gamete fertility. *Cell* 2001, **104**:153-164.

32. Ishino T, Yano K, Chinzei Y, Yuda M: Cell-passage activity is required for the malarial parasite to cross the liver sinusoidal cell layer. *PLoS Biol* 2004, **2**:

33. Kariu T, Yuda M, Yano K, Chinzei Y: MAEBL is essential for malarial sporozoite infection of the mosquito salivary gland. *J Exp Med* 2002, **195**:1317-1323.

34. Buscaglia CA, Coppens I, Hol WG, Nussenzweig V: Sites of interaction between aldolase and thrombospondin-related anomalous protein in plasmodium. *Mol Biol Cell* 2003, **14**:4947-4957.

35. Jewett TJ, Sibley LD: Aldolase forms a bridge between cell surface adhesins and the actin cytoskeleton in apicomplexan parasites. *Molecular cell* 2003, **11**:885-894.

36. Baker RP, Wijetilaka R, Urban S: Two Plasmodium rhomboid proteases preferentially cleave different adhesins implicated in all invasive stages of malaria. *PLoS Pathog* 2006, **2**:e113.

37. Cui L, Fan Q, Li J: The malaria parasite Plasmodium falciparum encodes members of the Puf RNA-binding protein family with conserved RNA binding activity. *Nucleic Acids Res* 2002, **30**:4607-4617.

38. Patankar S, Munasinghe A, Shoaibi A, Cummings LM, Wirth DF: Serial analysis of gene expression in Plasmodium falciparum reveals the global expression profile of erythrocytic stages and the presence of anti-sense transcripts in the malarial parasite. *Mol Biol Cell* 2001, **12**:3114-3125.

39. Gunasekera AM, Patankar S, Schug J, Eisen G, Kissinger J, Roos D, Wirth DF: Widespread distribution of antisense transcripts in the Plasmodium falciparum genome. *Mol Biochem Parasitol* 2004, **136**:35-42.

40. Ralph SA, Bischoff E, Mattei D, Sismeiro O, Dillies MA, Guigon G, Coppee JY, David PH, Scherf A: Transcriptome analysis of antigenic variation in Plasmodium falciparum–var silencing is not dependent on antisense RNA. *Genome Biol* 2005, **6**:R93.

41. Vanderberg JP: Development of infectivity by the Plasmodium berghei sporozoite. *J Parasitol* 1975, **61**:43-50.

42. Frischknecht F, Martin B, Thiery I, Bourgouin C, Menard R: Using green fluorescent malaria parasites to screen for permissive vector mosquitoes. *Malar J* 2006, **5**:23.

43. Natarajan R, Thathy V, Mota MM, Hafalla JC, Menard R, Vernick KD: Fluorescent Plasmodium berghei sporozoites and pre-erythrocytic stages: a new tool to study mosquito and mammalian host interactions with malaria parasites. *Cell Microbiol* 2001, **3**:371-379.

44. Raibaud A, Brahimi K, Roth CW, Brey PT, Faust DM: Differential gene expression in the ookinete stage of the malaria parasite Plasmodium berghei. *Mol Biochem Parasitol* 2006, **150**:107-113.

45. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. *Genome Biol* 2002, **3**:RESEARCH0034.

46. Muller PY, Janovjak H, Miserez AR, Dobbie Z: Processing of gene expression data generated by quantitative real-time RT-PCR. *BioTechniques* 2002, **32**:1372-1374. 1376, 1378–1379

Publish with **BioMed Central** and every
scientist can read your work free of charge

"BioMed Central will be the most significant development for disseminating the results of biomedical research in our lifetime."

Sir Paul Nurse, Cancer Research UK

Your research papers will be:

• available free of charge to the entire biomedical community

• peer reviewed and published immediately upon acceptance

• cited in PubMed and archived on PubMed Central

• yours — you keep the copyright

Submit your manuscript here:
http://www.biomedcentral.com/info/publishing_adv.asp

<figure><img src="image_6.png" /></figure>

BioMedcentral