#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Genome-Wide Identification of the Target Genes of AP2-O, a AP2-Family Transcription Factor


Although malarial parasites have a complex life cycle, they harbor only 30 transcription factors in their genome. The majority of these transcription factors belong to a single family referred to as the AP2 family. Our previous study suggested that stage-specific AP2 family transcription factors have critical roles in maintaining the Plasmodium parasite life cycle. However, it remains fairly elusive as to how these transcription factors regulate each stage. AP2-O is an AP2 family transcription factor that is expressed during the mosquito midgut-invading stage, the ookinete, and is essential for normal development of this stage. In the present study, we identified the entire set of AP2-O target genes to elucidate how this AP2 family transcription factor contributes to the formation of this stage. Our results showed that AP2-O directly regulates 10% of the parasite genome and is involved in the whole process of mosquito midgut-invasion by ookinetes. The global and comprehensive regulation by the AP2 family transcription factor that we revealed provides a model for transcriptional regulation of this parasite and may explain how malarial parasites regulate their complex life cycle using a small number of sequence-specific transcription factors.


Published in the journal: . PLoS Pathog 11(5): e32767. doi:10.1371/journal.ppat.1004905
Category: Research Article
doi: https://doi.org/10.1371/journal.ppat.1004905

Summary

Although malarial parasites have a complex life cycle, they harbor only 30 transcription factors in their genome. The majority of these transcription factors belong to a single family referred to as the AP2 family. Our previous study suggested that stage-specific AP2 family transcription factors have critical roles in maintaining the Plasmodium parasite life cycle. However, it remains fairly elusive as to how these transcription factors regulate each stage. AP2-O is an AP2 family transcription factor that is expressed during the mosquito midgut-invading stage, the ookinete, and is essential for normal development of this stage. In the present study, we identified the entire set of AP2-O target genes to elucidate how this AP2 family transcription factor contributes to the formation of this stage. Our results showed that AP2-O directly regulates 10% of the parasite genome and is involved in the whole process of mosquito midgut-invasion by ookinetes. The global and comprehensive regulation by the AP2 family transcription factor that we revealed provides a model for transcriptional regulation of this parasite and may explain how malarial parasites regulate their complex life cycle using a small number of sequence-specific transcription factors.

Introduction

Malarial parasites require two host animals during their life cycle and undergo multiple developmental changes in each host. According to these changes in the life cycle, parasites remarkably alter their repertoire of gene expression [1]. However, the corresponding regulatory mechanisms of gene expression remain poorly understood. In contrast to the lifecycle, malarial parasites have only a small set of sequence-specific transcription factors in their genome. The majority of the transcription factors belong to a single transcription factor family known as the Apetala2 (AP2) family, and 26–27 genes in this family have been detected in the genome [2,3]. The total number of sequence-specific transcription factors is exceptionally small compared with that in other eukaryotic organisms [46], suggesting that malaria parasites have a unique gene regulation system. Previous studies by us and other groups suggest that AP2-family transcription factors are involved in stage-specific gene regulation and are essential for normal development of the stages in which they are expressed [711]. However, only partial information has been obtained about their target genes; thus, it remains elusive how these transcription factors contribute to the development of each stage.

Ookinetes are motile forms of malarial parasites that are generated in the midgut of a mosquito after ingestion of an infected blood meal. The ookinetes promptly invade midgut epithelial cells and arrive at the basal lamina. There, they transform into oocysts, in which sporozoites, the liver-invading form, develop. We reported previously that the AP2-O AP2-family transcription factor is expressed in developing ookinetes of Plasmodium berghei [7]. Targeting experiments demonstrated that the disruptants display abnormal morphologies and completely lose infectivity to mosquitoes. We explored AP2-O targets by microarray analysis and identified 19 genes as targets. They included genes that encode microneme proteins and major surface proteins [7]. However, these genes do not explain the major abnormal morphogenesis phenotype of AP2-O disruptants, suggesting that most targets of this transcription factor remain to be identified.

The aim of this study was to investigate the basic features of transcriptional regulation in malaria parasites through elucidating the role of P. berghei AP2-O in this motile stage; determining the types of target genes controlled by the transcription factor and the extent to which they are responsible for gene regulation in this stage. We performed chromatin immunoprecipitation-high-throughput sequencing (ChIP-seq) and determined the whole range of P. berghei AP2-O target genes [12,13]. The results revealed that AP2-O regulates hundreds of genes directly and oversees the transcriptional regulation of this stage as a master regulator. Based on this result, we discuss the possibility that this centralized gene regulatory system represents a basic feature of transcriptional regulation in malarial parasites and could explain the paucity of transcription factors in this parasite.

Results

AP2-O has more than 1,000 binding sites on the genome

ChIP-seq was performed with transgenic P. berghei that expressed green fluorescent protein (GFP)-tagged AP2-O, using ChIP conditions established previously for ChIP-quantitative polymerase chain reaction analysis in ookinetes [7]. Two independent ChIP-seq analyses were performed, and the data were compared to confirm experimental reproducibility. The comparison showed that the results were quite reproducible (Fig 1A). Approximately 90% of the significantly-enriched peaks of the second experiment were found within 90-bp from a significantly-enriched peak in the first experiment (Fig 1B), suggesting that these peaks are common to both. The numbers of target genes predicted in the two experiments were 541 and 573, and 465 genes were common to both (S1 and S2 Tables in S1 File). In a subsequent analysis, we primarily used the results of the second experiment, which was performed more recently with a next-generation sequencing platform that is higher in both read length and read number.

Fig. 1. AP2-O targets over 500 genes in the genome.
AP2-O targets over 500 genes in the genome.
A. ChIP-seq experiments were performed independently using two different sequence platforms: an Illumina Genome analyzer and an ABI SOLiD 5500 system. This figure shows AP2-O peaks in each experiment within a 200-kb region (350–550-kb) of the fifth chromosome. The Integrative Genomic Viewer [56] was used for generating this image from bedgraph files (S ChIP-seq bedgraph). Peak call was performed with the MACS2 program [57], and 1,540 and 1,111 peaks were identified (FDR < 0.01, fold enrichments > 5). B. Distance of each peak in experiment 2 (1,111 peaks) to the nearest peak in experiment 1 (1,570 peaks) was calculated. Numbers of peaks that have a matching peak within the selected distance were plotted. The graph indicates that nearly 90% of peaks in experiment 2 have counterparts in those in experiment 1. C. Sequences enriched around the summits of AP2-O peaks. Logos were generated using WebLogo 3.3 (http://weblogo.threeplusone.com/create.cgi) [58]. D. Distances between the predicted summits of the AP2-O peaks and the motif sequences. The horizontal axis indicates the distance from the summit to the nearest motif sequence (the bin size is 20-bp). The vertical axis indicates the number of peaks in each region (total number of peaks = 959). E. The distribution of AP2-O peaks in the upstream regions of the target genes. The horizontal axis indicates the distance between the summits of the AP2-O peaks and the first methionine codon. The data were obtained from 53 putative target genes identified using microarray analysis. F. The P. berghei genes were divided into six groups containing 1,000 genes each (except for the 6th group, which contained 900 genes) according to RPKM values estimated by RNA-seq. The number of target genes in each group is shown as a histogram. The horizontal axis indicates the groups ordered according to their expression levels. G. A pie chart showing the functional categories of target genes (282 genes in total). Hypothetical protein genes are not shown. The number of members in each group is shown on the chart.

Analysis of the second set of ChIP-seq data with the MACS2 peak-calling program identified 1,111 peaks [false discovery rate (FDR) < 0.01 and fold enrichment over input > 5] [14]. Regions located within 100 base pairs (bp) of the predicted summit of each peak were sorted from the genome to identify the AP2-O-binding sequences. DNA sequences appearing frequently in these regions were analyzed using Fisher’s exact test (S3 and S4 Tables in S1 File). When six-base sequences were ranked according to their p-values, 17 of the top 20 sequences overlapped, thereby yielding the [TC][AG]GC[TC][AG] binding motif. This was the same motif predicted in our previous study [7] (Fig 1C), and 86% of the peaks had this motif within the peak region. Fig 1D summarizes the distances between the summits of these peaks and the binding motif. These results confirm that the motif predicted in our previous study is used for AP2-O binding in vivo. We further investigated the remaining 14% of the peaks that did not contain this motif. Statistical analysis using these peaks showed that another motif sequence, TG[ATC]ACA, was highly concentrated around the summit (S5 and S6 Tables in S1 File). The two nucleotides on either end of this motif, viz. TG and CA, were the same as those in the primary motif, suggesting that they are minor variants of the binding motif. After including those sequences, 99% of the peaks (1,097 of 1,111) had either motif within their regions. The distances of these motifs from the predicted summits averaged 26.6-bp.

Target gene prediction revealed over 500 putative target genes

Next, the positional relationship of the binding sites with the target was investigated to predict the target genes based on the binding sites. Our previous study suggested that AP2-O binding sequences are usually located in the 1-kb upstream regions of the target genes, on the basis of data of 19 genes that were mainly (15 of 19 genes) identified by microarray analysis between wild-type parasites and AP2-O disruptants [7]. However, the number of genes seemed too small to define the relationship (the array was designed on our expressed sequence tag data). Accordingly, to identify more candidates as targets, we performed another microarray analysis using an array designed on the P. berghei genome sequence [11]. In that analysis, we identified 63 genes as candidate AP2-O targets (a twofold decrease in the AP2-O disruptants was observed). The 63 genes contained all 15 genes identified in the previous study. The ChIP-seq data showed that 53 of these genes had AP2-O peaks in the upstream intergenic region (Table 1). Among these genes, approximately 90% of the AP2-O peaks (49/53) existed within a 1,200-bp region upstream of the start codon and most frequently at 400–600-bp (Fig 1E). Based on this result, we concluded that it was appropriate to define the predicted target genes as those with AP2-O binding sites within the 1,200-bp region upstream of the start codon. However, this result indicated that target genes would be missed in such a prediction when they harbored AP2-O binding sites further upstream of this 1,200-bp region. Thus, in the following study, the targets were further manually investigated using the ChIP-seq and high-throughput cDNA sequencing (RNA-seq) data (See also Table 1, Table 2, and S7 Table in S1 File).

Tab. 1. List of genes whose expression decreased in AP2-O(–) ookinetes.
List of genes whose expression decreased in <i>AP2-O</i>(–) ookinetes.
Microarray analysis was performed between wild-type and AP2-O(–) parasites. Genes whose expression was decreased over two-fold are listed. Asterisks indicate genes that have peaks over 1200-bp upstream, in the intron, or in the exon of adjacent genes. Target genes identified in experiment 2 are indicated by “Yes” in the third column.

Tab. 2. Gene expression analysis of P. berghei ookinetes by RNA-seq.
Gene expression analysis of <i>P</i>. <i>berghei</i> ookinetes by RNA-seq.
RNA-seq analysis was performed in P. berghei ookinetes cultured for 24 h, and genes were ordered according to reads per kilobase of coding sequence per million reads (RPKM) value. Asterisks indicate genes that have peaks over 1200-bp upstream or in the exon of adjacent genes.

According to the prediction rule defined above, 541 genes (approximately 10% of all P. berghei genes) were predicted to be AP2-O target genes (S1 Table in S1 File). These predicted target genes (later target genes) included 18 of the 19 target genes identified in our previous study [7], and almost all genes that have been reported to be ookinete-specific genes or to be important in midgut-infection of the parasite, as descried later.

Of the remaining 10 genes that did not show AP2-O peaks in the upstream intergenic region, one gene showed a peak in the first intron and four genes, which all harbored a short upstream intergenic region, showed a peak in the exon of the adjacent gene (S1 Fig). All these peaks accompanied the transcripts downstream. These results suggest that exons and introns are occasionally used as promoters in Plasmodium. However, AP2-O summits on the exons and introns were not used for predicting target genes in subsequent analyses because whether such cryptic transcripts would be finally translated into proteins with biological functions was unclear.

AP2-O target genes in the ookinete transcriptome and proteome

An RNA-seq analysis was performed to investigate the expression levels of AP2-O target genes in the ookinete transcriptome. As summarized in a histogram (Fig 1F), targets were biased towards highly expressed genes. Among the 50 genes showing the highest number of reads per kilobase of coding sequence per million reads (RPKM), 38 genes were AP2-O targets (Table 2), suggesting that AP2-O contributes to a stage-specific gene expression pattern in ookinetes. It is notable that the majority of the AP2-O target genes within these 50 genes were not predicted to be targets by microarray analysis (Table 2). This was probably because a high background of transcripts carried from the female gametocytes (e.g., p28 transcripts) made it difficult to predict targets using expression differences [7,15]. This clearly shows that ChIP-seq has an advantage in target prediction.

Proteomics data are now available for the ookinete stage [1,1618]. These data provide a useful resource for exploring genes involved in midgut invasion [19]. We compared the list of AP2-O target genes with the data of microneme proteome which contains 330 proteins (S8 Table in S1 File) [16]. Microneme proteome contained six known microneme proteins, two microneme proteins newly identified in this proteome study, and six putative secreted ookinete proteins (PSOP). ChIP-seq identified 11 of these 14 proteins, and only three proteins were missed in the list of AP2-O targets. On the other hand, the list of AP2-O targets contains two genes that encode the membrane-attack ookinete protein (MAOP) [or perforin-like protein 3 (PPLP3)] and PPLP5. These proteins are not present in the microneme proteome (S7 and S8 Tables in S1 File). These results indicate that the target gene data are as comprehensive as the proteome data. The proteome data contain a number of proteins specific to other stages, including eight rhoptry proteins, five merozoite surface proteins, and eight proteins exported to the erythrocyte (21 proteins in total). These contaminants were probably included because it is difficult to completely remove blood stage parasites from ookinete samples by using differences in density. This means that ookinete proteome data inevitably contain contaminants from other stages. In clear contrast, the list of AP2-O targets contained none of these proteins (S8 Table in S1 File). This lack of contamination from other stages may result from the fact that ChIP-seq analysis is based solely on information about the binding of a stage-specific transcription factor.

Overview of target gene functions

A total of 541 target genes were identified in the analyses described above. In addition, as described later, one novel gene was identified as a target gene in the genomic region that had been thought to be intergenic. All 542 genes were classified into several major groups, and this classification was based on their functional annotations in PlasmoDB (http://PlasmoDB.org), their corresponding structures (such as an N-terminal signal sequence), as well as similarities to functionally annotated genes in other apicomplexan parasites (Fig 1G and S7 Table in S1 File). This classification showed that the targets included several genes related to the morphogenesis of, and midgut invasion by ookinetes, together with genes involved in general functions such as translation, transcription, and DNA replication. In the following sections, we describe the target genes belonging to the categories related to midgut invasion and further explore novel members of these categories among the target genes using the list of identified targets as a resource for gene exploration.

Exploration for possible missing genes in the P. berghei genome

The target prediction described above was based on the genome annotation in PlasmoDB. However, this annotation is still underway, and some small-sized genes may have been missed. Therefore, we used our data to look for target genes still missing in the P. berghei genome. We selected AP2-O peaks present in the intergenic regions that had no corresponding targets and then manually looked for open reading frames (ORFs) in the nearby regions using the BLAST program and the RNA-seq data. Using this screening technique, an ORF was identified in the region flanked by the 3′-portions of two genes (PBANKA_141090 and PBANKA_141100, respectively) (Fig 2A). It encoded a small putative protein comprising 60 amino acids that is conserved in Plasmodium spp. (Fig 2B).

Fig. 2. Identification of a missing gene in the P. berghei genome.
Identification of a missing gene in the <i>P</i>. <i>berghei</i> genome.
A. Mapped view of the ChIP-seq and RNA-seq reads of the identified gene in ookinetes. The RNA-seq reads were located in the intergenic region between the 3′-portions of neighboring genes (PBANKA 141090 and PBANKA 141100). The identified ORF is indicated by a red rectangle. B. The amino acid sequences encoded by the corresponding ORFs were aligned in different Plasmodium species. The alignment was performed using ClustalW2.1 (http://www.ebi.ac.uk/Tools/msa/clustalw2). Pb, P. berghei; Pv. P. vivax; Pf, P. falciparum. C. A fluorescence image of P. berghei ookinetes expressing the GFP-tagged protein. Ookietes were cultured in vitro and observed by fluorescence microscopy at 24 h after fertilization. Apical ends of ookinetes are indicated by arrowheads. Left, merged image of GFP and nuclear staining with Hoechst 33342. Right, optical transmission image. Scale bar, 10 μm. D. Immunoelectron microscopy image of an ookinete expressing the putative protein. Immunoelectron microscopy was performed with anti-GFP antibodies. Left, sagittal section. Colloidal gold particles (15 nm) are located on the apical polar ring, a low-electron density area between the collar (edges of the color are indicated by closed arrowheads) and the striated structure of microtubules (indicated by open arrowheads. Right, cross-section image. Colloidal gold particles (15 nm) are located on the apical polar ring, a low-electron density area between the collar (high-electron density area), and microtubules (cross-sections of microtubules are indicated by arrowheads). The colloidal gold particles are also located on the fibrous tissue observed in the lower-left part of the section, which could be the microtubules adhering to the apical ring and the cytoskeletal fibers surrounding them. Scale bars, 0.5 μm.

To examine whether this hypothetical ORF was translated into a protein, it was expressed as a GFP-tagged protein under control of the original promoter, using the centromere plasmid pCen-GFP [8]. Expression of the tagged protein was found at the apical tip of the ookinetes (Fig 2C). Immunoelectron microscopy showed that the GFP-tagged protein was localized to the apical polar ring, corresponding to the low electron density structure beneath the collar (Fig 2D, left). Transverse section images demonstrated that this protein was also distributed at the boundaries between subpellicular microtubules and the apical ring, suggesting that the protein plays a role connecting the above two structures (Fig 2D, right). This is the first Plasmodium protein reported to localize to the apical polar ring, and we designated this protein apical ring associated protein 1 (ARA1). No ARA1 homologs have been found in organisms other than Plasmodium parasites.

AP2-O induces genes involved in ookinete morphogenesis and gliding motility

The pellicle is a thin structure constituted from the plasma membrane and the closely apposed inner membrane complex (IMC). The IMC is a complex of the membranous structure (IMC membrane) and the subpellicular network (SPN) of cytoskeletal proteins, such as the alveolins/IMC1 proteins. The pellicle underlies the entire parasite plasma membrane, except for the apical portion, and is supported by a row of subpellicular microtubules that originate from the apical ring. The pellicle space contains an actin-myosin motor that generates the driving force for gliding motility. Therefore, this structure is essential for the motile stages of these parasites. Pellicular/IMC proteins have been studied intensively in Toxoplasma gondii tachyzoites [20,21], and several pellicle-associated proteins have been reported. Most of these genes have orthologs in the Plasmodium genome. However, none of these genes have been identified as AP2-O target genes [7]. In the present study, 22 genes encoding pellicullar/subpellicular proteins were identified as AP2-O target genes (S7 Table in S1 File). The constellation of these target gene products in the pellicle and subpellicular structure is illustrated in Fig 3A.

Fig. 3. Genes encoding pellicular proteins are AP2-O targets.
Genes encoding pellicular proteins are AP2-O targets.
A. The pellicle structure and its components are illustrated. Targets are highlighted in red. In addition to the proteins mentioned in the text, targets contained several genes encoding putative pellicular/IMC proteins, viz. apicortin, tubulin-tyrosine ligase (TTL), and small heat shock-related protein 20 (HSP20). Apicortin is involved in the stabilization of microtubules [59]. TTL, an enzyme adding a tyrosine to the carboxyl end of tubulin, marks the plus ends of growing microtubules and regulates microtubule growth at this site [60]. HSP20 is involved in ookinete motility [51,61]. MyoA, myosin A; ADF1, actin-depolymerizing factor 1; SPM, subpellicular microtubule. B. IMC1i was expressed as a GFP-tagged protein using a P. berghei centromere plasmid. Ookinetes were cultured in vitro and observed by fluorescence microscopy at 24 h after fertilization. The apical end of the ookinete is indicated by an arrowhead. Left, merged image of GFP and nuclear staining with Hoechst 33342. Right, optical transmission image. Scale bar, 10 μm. C. Giemsa-stained image of wild-type (left) and IMC1i(-) (right) ookinetes, 24 h after fertilization. Scale bar, 10 μm. D. Fluorescence microscopy image of ookinetes expressing GFP-tagged PUA26. Left, merged image of GFP and nuclear staining with Hoechst 33342. Right, optical transmission image. Scale bar, 10 μm. E. Immunoelectron microscopy image. Cross-sections of two ookinetes expressing GFP-tagged PUA26 are shown. Immunoelectron microscopy was performed with anti-GFP antibodies. Colloidal gold particles (15 nm) are mainly localized at the IMC, which is the layer of high-electron density beneath the plasma membrane. Subpellicular microtubules are indicated by arrowheads. Scale bar, 1 μm. F. Giemsa-stained image of wild-type (left) and PUA26(-) (right) ookinetes at 24 h after fertilization. Scale bar, 10 μm. G. Height—width ratios were compared between wild-type and PUA26(-) ookinetes. Ookinetes were cultured for 24 h and stained with Giemsa. Micrographs were obtained with an Olympus BX60 fluorescence microscope. Height—width ratios of ookinetes were measured with the AquaCosmos software (Hamamatsu Photonic System). In total, 100 ookinetes were analyzed in each parasite population. Bars represent mean ± SE. H. The numbers of parasites associated with the midgut were compared between wild-type and PUA26(-) parasites at 24 h after an infective blood meal by mosquitoes. Data are the mean ± SE of three independent experiments using 20 mosquitoes each. Only parasites on a single side of the midgut were counted.

The IMC1/alveolin-like genes constitute the largest group of paralogous genes in this category. Eight IMC1 genes (IMC1a—IMC1h) have been identified in the P. berghei genome [22], and there are five additional members of the alveolin family (tentatively named IMC1i—IMC1m) in it (S2 Fig and S7 Table in S1 File). Of these 13 Plasmodium IMC1-like genes, eight were identified as AP2-O targets (Fig 3A and S2 Fig). This result indicates that the SPN of this motile stage is composed of a stage-specific cocktail of IMC1 proteins. We generated parasites expressing IMC1i as a C-terminal GFP fusion protein to examine if newly identified IMC1-like genes actually encode IMC proteins (Fig 3B and S3 Fig). Signals were observed along the surface, but not at the apical end of the mature ookinete. This distribution pattern is characteristic of SPN proteins. We further disrupted the IMC1i gene and investigated the resulting phenotype (S3 Fig). The mutants developed into ookinetes with normal conversion rates (S9 Table in S2 File); however, nearly all of them (99.8%, n = 520) displayed abnormal morphologies (Fig 3C). To exclude the possibility that unfertilized gametocytes were mistakenly counted as immature ookinete, immunofluorescent staining with antibody against circumsporozoite- and TRAP-related protein (CTRP) was performed. Number of ookinetes per microscopic field observed by the immunofluorescent staining [14.7 ± 12.5 (mean ± SE; n = 10)] was essentially same as that observed by Giemsa staining [15.0 ± 11.9 (mean ± SE; n = 10)], which confirmed that round ookinetes-like parasites observed by Giemsa staining were genuine ookinetes. The morphologies were similar to those of ookinetes depleted of IMC1b, another ookinete protein in this family [23], and to those of AP2-O-depleted ookinetes [7]. They nearly lost their ability to infect mosquitoes (S10 Table in S2 File). These results suggest that the reduced production of these putative cytoskeletal proteins may have caused the abnormal morphology seen in AP2-O-depleted ookinetes.

Along with these IMC1 genes, the list of targets included orthologs of the T. gondii pellicullar/IMC proteins: P. berghei PhIL1 (photosensitized INA-labeled protein 1), IMC sub-compartment protein 1 (ISP1), ISP3, subpellicular microtubule-associated protein 1 (SPM1), and SPM2 [20,2426] (Fig 3A). PhIL1 and ISP1 are localized mainly with IMC and the apical cap of T. gondii tachyzoites, which is an apical structure linked to the IMC membrane that covers the parasite’s apical protrusion. We investigated localization of P. berghei PhIL1 by generating parasites expressing a GFP-tagged protein. The tagged protein spread along the cell surface of mature ookinetes, but still localized predominantly to the apical protrusion (S4 Fig). This distribution pattern was similar to the pattern reported for T. gondii PhIL1 [20,21], suggesting that ookinetes have the apical structure corresponding to the apical cap of T. gondii tachyzoites.

The pellicle also contains the glideosome, a complex of motor proteins (Fig 3A). The glideosome is linked to the cytoplasmic domains of adhesins integrated into the plasma membrane, facilitating the gliding motility of apicomplexan parasites. According to T. gondii studies, the glideosome is composed of actin filaments, myosin A, the myosin A tail domain-interacting protein, aldolase, glideosome-associated proteins (GAPs) [27], and GAPs with multiple membrane spans [28]. All of the corresponding genes have orthologs in Plasmodium parasites, and the majority of them, excluding actin, myosin A, and aldolase, are AP2-O targets (The GAP 50 gene shows a ChIP peak in the 1.4-kbp upstream). Taken together, these results suggest that AP2-O is involved in morphogenesis and motility of ookinetes.

Exploring novel genes encoding ookinete pellicle proteins

The pellicle proteins described above are linage-specific; they have homologs solely in apicomplexan parasites or alveolates (alveolates are higher-order groups of apicomplexan parasites). Therefore, we explored genes encoding novel pellicle proteins among the target genes that have not been functionally annotated, but which are uniquely conserved in organisms of this lineage. The list of target genes contained 36 genes satisfying these criteria. These genes were tentatively designated PUA (protein unique to apicomplexan parasites) 1–36 (S7 Table in S1 File). We expressed three small PUA genes (PUA17, PUA19, and PUA26) as GFP-tagged proteins using pCen-GFP and examined their subcellular localization (Fig 3D and S5 Fig). Among these proteins, GFP-tagged PUA26 displayed a distribution pattern characteristic of SPN proteins (localization at the parasite surface, except at the apical end, as shown in Fig 3D). Therefore, we further investigated its localization and function using immunoelectron microscopy and gene targeting (S3 Fig), respectively. By immunoelectron microscopy gold particles were localized predominantly along the parasite surface, except at the apical region (S6 Fig). Under higher magnification, it was evident that the particles were localized not to the plasma membrane, but on the electron-dense structure beneath it, indicating that PUA26 is located on IMC (Fig 3E). Targeting this gene did not affect ookinete conversion rates (S9 Table in S2 File); however, the morphologies of generated ookinetes appeared somewhat laterally longer than those of the wild-type parasites (Fig 3F). Calculation of the height—width ratios of these ookinetes confirmed that they had abnormal morphologies (Fig 3G). The number of oocysts observed in the midgut was also significantly reduced in these parasites (Table 3). We prepared GFP-expressing parasites from these disruptants and investigated the number of ookinetes that successfully reached the midgut lamina at 24 h after an infective blood meal by mosquitoes. Ookinetes observed in the midgut were significantly decreased in the disruptants (Fig 3H), while the majority of them (84.5 ± 5.4%) began to transform into spherical early oocysts, as observed in wild-type parasites (75.5 ± 5.4%). These results suggest that the disruption impaired ookinete cytoskeletal structure and reduced ookinete motility, as reported for other IMC proteins [29]. PUA26 encodes a 96-amino-acid protein that possesses no known functional motifs. Orthologs of this gene exist in T. gondii (TGME49_014220), Neospora caninum(XP_003884749), and Eimeria spp (CDJ60835). We designated this protein as IMC-associated apicomplexan protein (IAAP).

Tab. 3. Oocyst formation in mutant parasites.
Oocyst formation in mutant parasites.
Two independently prepared mutants were used for assessing the phenotype in each gene. The number of oocysts was counted 14 days after an infective blood meal by mosquitoes. The value is an average from 25 mosquitoes. SE, standard error.

Putative microneme proteins among the targets

Secreted proteins constituted the largest group among the targets; in total, 43 genes encoded a putative secreted protein (S7 Table in S1 File). The majority of them could be microneme proteins, as ookinetes have no secretory organelles for intracellular parasitic infection. The list of targets included all seven proteins reported to be localized to the ookinete microneme: chitinase, CTRP, secreted ookinete adhesive protein, von Willebrand factor A-domain-related protein, cell traversal protein for ookinetes and sporozoites, GPI-anchored micronemal antigen (GAMA), and MAOP (PPLP3), and the seven putative microneme proteins detected during the ookinete proteomic studies: PPLP4, PPLP5, PSOP1, PSOP2, PSOP6, PSOP7, and PSOP12 [19,3037] (Fig 4A and S7 Table in S1 File). These results indicate that microneme proteins are induced as a set by AP2-O in developing ookinetes. Another 25 unannotated genes in this category were tentatively designated here as POM (putative ookinete microneme proteins) 1–25 (S7 Table in S1 File).

Fig. 4. Target genes required for midgut invasion and oocyst formation.
Target genes required for midgut invasion and oocyst formation.
A. Overview of the target genes involved in midgut invasion and oocyst formation. B. The numbers of parasites associated with the midgut were compared between wild-type and PPLP4(-) parasites at 24 h after an infective blood meal by mosquitoes, as in Fig 3G. Data are the mean ± SE of three independent experiments using 20 mosquitoes each. C. Ratios of okinetes to total parasites associated with midguts at 24 h after an infective blood meal by mosquitoes. Only ookinetes of fully elongated shape were judged as ookinetes. Data are the mean ± SE of three independent experiments using 20 mosquitoes each. D. Oocysts were counted 2 and 14 dpi. Three independent experiments were performed with each clone. Data are represented as mean ± SE. E. Dot plots of diameters of oocysts at 14 dpi. In total, 200 oocysts were used for measurements in each experiment. Bars represent the mean ± SE of diameters. Mosquitoes were fed on mice infected with POS8(-) parasites constitutively expressing GFP. F. The oocysts in the midgut were counted at different time points (2, 10, and 14 dpi), using the same GFP-expressing parasites as in H. Three independent experiments were performed with each clone. Data are represented as mean ± SE. G. Dot plots of diameters of oocysts at 14 dpi. In total, 200 oocysts were used for measurements in each experiment. Bars represent the mean ± SE of diameters. Mosquitoes were fed on mice infected with CYC3(-) parasites constitutively expressing GFP. H. Fluorescence microscopy image of a mosquito midgut infected with wild-type (left) and CYC3(-) (right) parasites expressing GFP at 14 dpi. Scale bar, 300 μm.

In addition to these genes, the list of predicted targets included two genes encoding secreted proteins of the osmiophilic body, which is the secretory organelle of gametocytes. We expressed one of the genes, viz. gamete egress and sporozoite traversal protein (GEST), as a GFP-tagged protein and demonstrated that it is expressed in the ookinete as a microneme protein (S7 Fig). As the osmiophilic body is involved in egress of the gametocyte from the host erythrocyte, its expression during motile stages suggests that parasites employ common mechanisms for egression and invasion into host cells. This similarity between the two organelles was also suggested from the finding that perforin-like proteins are used for egression by gametocytes [38] and for invasion by ookinetes and sporozoites, respectively [35,39].

We generated disruptants in five target genes belonging to this category whose roles in ookinetes remain to be elucidated: CS domain protein, PPLP4, POM2, POM7, and POM16. They were all successfully disrupted, but a clear reduction in oocysts was observed only with PPLP4 (Table 3 and S11 Table in S2 File). PPLP4 encodes one of the five PPLP genes (PPLP1PPLP5) identified in the P. berghei genome that contains a membrane attack complex/perforin (MACPF) domain. Targeting PPLP4 did not affect ookinete conversion rates (S9 Table in S2 File) and morphologies of ookinetes but resulted in the complete loss of infectivity to mosquitoes (Table 3). We generated parasites that constitutively expressed GFP from these disruptant populations and assessed midgut invasion by ookinetes in vivo at 24 h after an infective blood meal by mosquitoes (Fig 4B and 4C). A small number of ookinetes were observed in the midgut (Fig 4B). However, almost all these ookinetes still displayed elongated shapes (99.5 ± 1.4%), whereas in wild-type parasites, the majority of ookinetes associated with the midgut had already started to transform into early oocysts (Fig 4C). This means that the ookinetes associated with the midguts observed in the parasites may not have arrived at the basal lamina but may have rather attached to the apical side of the epithelial cells, as we reported in PPLP3/MAOP disruptants [35]. Three PPLPs, PPLP3PPLP5, were detected in ookinete micronemes during a proteomic study [19], and PPLP3 and PPLP5 are essential for ookinete infection of the midgut [35,40]. The present results demonstrated that all three ookinete PPLP genes were AP2-O targets and are essential for infection of the midgut. This is in clear contrast with the redundancy of other target genes in this category, raising the possibility that these PPLPs assemble into a protein complex and form MAC on the plasma membrane of a target cell, as do human complement proteins C6–C9, which all contain a MACPF domain [41].

In addition to secretory proteins, the list of targets contained genes involved in vesicular protein transport (Fig 4A and S7 Table in S1 File). Considering that the predominant secretory proteins in this stage are all microneme proteins, these genes could be involved in transport and secretion of microneme proteins. In particular, several proteins, possibly constituting the soluble N-ethylmaleimide-sensitive factor activating protein receptor (SNARE) machinery, were included in the list of targets. These included vesicle-associated SNARE protein, target-associated SNARE protein, syntaxin, syntaxin-binding protein, mammalian uncoordinated protein 18-related protein, and a C2-domain protein with a transmembrane region (S7 Table in S1 File). The functions of these genes during this stage remain elusive, but our data could serve as a resource for future investigations into the mechanism of ookinete microneme secretion.

Putative plasma membrane-associated proteins on the list of targets

We classified genes encoding proteins with an N-terminal signal sequence and a structure for anchoring to the plasma membrane, such as the membrane-integrated domain and the GPI modification site, as putative plasma membrane-associated proteins. This group could include proteins that are targeted to the parasite surface after being secreted from micronemes, such as CTRP and GAMA, as well as proteins directly targeted to the cell surface. Genes in this group are important for ookinete biology because their products could be involved in parasite interactions with the midgut epithelium, and are still largely unknown. The list of targets contained 14 of these genes, including two genes that have been reported to be expressed during this stage: P25 and P28 (Fig 4A and S7 Table in S1 File, respectively). The MSP10 gene may have been incorrectly predicted as a target, mostly because the AP2-O binding site is very close to the start codon (< 50-bp) and probably exists in the 5′-UTR. This finding suggests that the region close to the start codon must be excluded from target prediction. Ten other unannotated genes in this category were tentatively designated as putative ookinete surface-associated protein (POS) 1–10 (S7 Table in S1 File).

We performed targeting experiments with three genes of unknown function in this category (POS7–9) (S3 Fig), and we investigated their involvement in midgut infection. During the targeting experiments, a significant decrease in the number of oocysts was observed only for the POS8 disruptant (Table 3 and S12 Table in S2 File). We further analyzed the phenotype of this disruptant using independently prepared disruptant populations. In these disruptants, the number of oocysts decreased more than 20-fold (Table 3). The size of parasite oocysts was clearly smaller than that of the wild-type parasite oocysts [as observed by phase contrast microscopy at 14 days post-infection (dpi)]. Therefore, to identify the step in which they decreased in number, we generated parasites that constitutively expressed GFP from these disruptant populations and counted the oocysts shortly after ookinete invasion of the midgut at 2 dpi (Fig 4D). The number of oocysts was already approximately 10-fold smaller than that of wild-type parasites. We further measured oocyst diameters by epifluorescence microscopy at 14 dpi. As shown in Fig 4E, the average oocyst diameter was approximately 60% of that of the wild-type oocysts, suggesting that oocyst development was impaired by this disruption. This phenotype indicates that POS8 plays a critical role in ookinete-mediated midgut invasion as well as in subsequent oocyst development. We speculate that this gene product might participate in oocyst development by establishing the foothold necessary for subsequent oocyst development. To determine the localization of this protein at the surface of the ookinetes and to explain the phenotype of the disruptants by in vivo observation, we transfected wild-type parasites with a centromere plasmid containing a construct for expressing POS8 as a GFP-fused protein. However, we could not detect GFP signals in the ookinetes of these parasites. At present, we speculate that the fusion of GFP to the putative membrane-associated region of the protein may affect its targeting to the plasma membrane. Further study is required to explain why the abnormal phenotypes were observed over the two stages.

The glutathione redox system is activated by AP2-O for survival in mosquitoes

Reactive oxygen species (ROS) are produced in mosquitoes after a blood meal and in response to penetration of the midgut epithelium by ookinetes. Epithelial cells play a critical role eliminating parasites from the midgut [42,43]. The parasite reduction-oxidation (redox) system is essential for parasites to combat these oxidative stressors and is therefore critically important for their midgut survival. Parasites have two redox systems, a thioredoxin system and a glutathione system [44]. Of these two systems, the glutathione system plays a central role in midgut parasite infection and oocyst development [45,46]. This system is maintained by three genes that produce gamma-glutamylcysteine synthetase (gammaGCS), glutathione synthetase, and glutathione reductase. Glutathione is synthesized de novo from glutamine and cysteine by gammaGCS and glutathione synthetase. Glutathione reductase reduces oxidized glutathione to glutathione. A ChIP-seq analysis showed that all of these genes are AP2-O targets (Fig 4A and S7 Table in S1 File), indicating that these genes are “programmed” to be induced during this stage prior to midgut invasion, thereby protecting ookinetes from the mosquito immune system during midgut penetration. This is in clear contrast to the mechanisms for sensing and responding to oxidative stressors observed in other eukaryotes. The glutathione system is also important for oocyst development. Considering that transcripts of these genes are abundant in the ookinete stage, it is possible that transcripts of these genes in ookinetes are prepared in part to protect early oocysts from ROS produced in the midgut.

AP2-O induces the genes required for oocyst formation

After arrival at the basal side of the midgut, ookinetes lose their elongated shape and transform into spherical oocysts. The list of AP2-O targets contained the Cap380 gene, which encodes the oocyst major capsule protein [47]. RNA-seq analysis showed that the transcripts of this gene are abundant in ookinetes (S13 Table in S2 File), and DNA microarray analysis showed that its expression decreased more than 10-fold in AP2-O(-) parasites (Table 1). This finding suggests that the AP2-O target genes could include genes for oocyst development and that expression of a portion of target genes is post-transcriptionally regulated in ookinetes for the subsequent oocyst stage as reported in other stages, such as the female gametocyte stage or the sporozoite stage [15,48]. The list of targets contained a number of gene products with RNA-interacting motifs, such as the zinc finger motif and RNA recognition motifs (Fig 1G and S7 Table in S1 File). They could be involved in post-transcriptional regulation of such genes and contribute to the progression of ookinetes into oocysts.

We found another target gene with an important role in oocyst development. This gene (PBANKA_123320), which encodes a cyclin-like protein, was tentatively named CYC3, to be consistent with the annotation of its ortholog in P. falciparum (PlasmoDB). Among the five genes annotated as cyclins in the P. berghei genome, CYC3 is most abundantly expressed in ookinetes (S13 Table in S2 File). Parasites in which this gene was disrupted (S3 Fig) showed a normal ookinete conversion rate (S9 Table in S2 File) and formed morphologically normal ookinetes. However, the numbers of oocysts and oocyst sporozoites at 14 dpi were several-fold lower than those present in the wild-type (Table 3). Moreover, oocysts observed with phase-contrast microscopy were smaller than corresponding wild-type oocysts. Following this, we prepared disruptants constitutively expressing GFP from the same set of parasites and conducted a time-course study of the number of oocysts. The oocyst number of the disruptants was approximately 75% of that of wild-type parasites at 2 dpi and decreased to 20%–30% of that of the wild-type parasites at 14 dpi (Fig 4F). A high proportion of small-size oocysts, which were too small to be detected by phase-contrast microscopy, were observed by epifluorescence microscopy at 14 dpi (Fig 4G and 4H). The average diameter of the oocysts at this time point was approximately 65% of that of wild-type oocysts (Fig 4G). By nuclear staining with cell permeable Hoechst 33342 at 4 dpi, nuclei of most disruptants were not detected (27 of 30 oocysts), suggesting that nuclear division scarcely proceeded in them. To exclude the possibility that uncompleted meiosis during zygote to ookinetes transition affected the subsequent oocyst development [49,50], we measured the content of nuclear DNA in mature ookinetes by staining with Hoechst 33342. No differences in DNA contents were observed between wild-type ookinetes and disruptants (S8 Fig), indicating that meiosis was normally completed in the disruptants. Collectively, these results strongly suggest that CYC3 controls cell cycle progression in early oocysts.

Discussion

We explored the AP2-O targets with ChIP-seq analyses, identified approximately 1,100 binding sites on the genome, and predicted over 500 genes as targets. The list of targets included a set of genes necessary for pellicle formation, which finally explained why ookinetes with disrupted AP2-O display abnormal morphologies. The targets also included a series of genes necessary for midgut infection by the parasite. The targets comprised genes encoding proteins for gliding motility, microneme proteins, and surface proteins; genes necessary for the redox system; and genes for oocyst development. Further, although not mentioned specifically, the list of targets included genes encoding protein kinases (Fig 1G and S7 Table in S1 File), three of which, viz. CDPK1, CDPK3, and PBANKA_146050, are essential for midgut infection of ookinetes [5153]. These results show that a single transcription factor is involved in all processes of ookinete-mediated midgut invasion and that transcriptional regulation during this stage is centralized to this transcription factor.

This transcriptional regulation in ookinetes is in clear contrast to that reported in model eukaryotes, in which a network system composed of a number of transcription factors regulates gene expression. A transcription factor usually regulates functionally related genes in budding yeast [6], and a bundle of these modules constitutes a transcriptional regulatory system [54]. Therefore, different biological pathways can be activated independently, enabling the organism to create different gene expression patterns according to environmental conditions. In contrast, transcriptional regulation in ookinetes definitely lacks the flexibility necessary for adapting to environmental change. In this stage, parasites may not be able to control a group of genes separately from other target genes. Hundreds of target genes should be induced in a set according to the program defined in advance for this stage. This observation suggests that the parasites can survive only under limited environmental conditions.

Obviously, stable and predictable environments are prerequisite for this inflexible system, and their parasitic lifestyle seems to ensure them such environments. It would be reasonable to speculate that parasitism relieved them from the necessity to respond to environmental change and allowed them to reduce the number of regulatory genes including transcription factors. If the ookinete regulatory system is an evolutionary consequence of the parasitic lifestyle, it would not be surprising if these parasites adopt similar regulation systems in other lifecycle stages.

Major families of transcription factors different from the AP2-family have not been demonstrated, supporting the assumption that the AP2-family is the sole sequence-specific transcription factor family of this parasite and suggests that the total number of transcription factors in these parasites is about 30, which is an order smaller than those in ordinary eukaryotes, because transcription factors usually constitute 5–10% of the eukaryote genome and these parasites possess >5000 genes in their genome. This paucity of transcription factors and their complicated lifecycle appear to oppose each other and suggest that these parasites have a unique gene regulatory system. In this study, we provide the first information regarding the parasite transcriptional regulatory system and revealed that a single master transcription factor directly regulates the broad range of targets in ookinetes. This finding suggests one simple gene regulatory model in these parasites as follows: Each lifecycle stage possesses a stage-specific master transcription factor. That factor directly induces a number of target genes during the stage, creating a stage-specific gene expression repertoire. This simple model could explain the paucity of transcription factors in these parasites [2,55]. To examine if this model could be extended to other stages, we are now performing ChIP-seq analyses of AP2 transcription factors in several lifecycle stages including the proliferation and motility stages.

AP2-O was first observed in the nucleus of the developing zygote/ookinete approximately 8 h after fertilization [7]. This observation suggests that AP2-O would not be involved in transcriptional regulation in the early development of ookinetes. This is also supported by the fact that parasites disrupted with AP2-O can develop into retort forms [7]. Therefore, it is possible that other sequence-specific transcription factors are expressed during this period and participate in transcriptional regulation of this stage. However, at present, the most likely regulation involved in this early development would be translational regulation. It has been reported that large amounts of transcripts are stored in female gametocytes for development in this period and that the development of parasites lacking RNA helicase DOZI, the major component of this translation regulation system, is completely halted during the early phase [15]. The expression of AP2-O is also regulated by this regulation system, and as demonstrated by cross-fertilization experiments, its transcripts are mainly derived from the female gametocyte [7]. Therefore, it is possible that gene expression in this period largely depends on transcripts prepared in female gametocytes. On the other hand, RNA-seq analysis in the present study showed that most highly expressed genes in mature ookinetes are targets of AP2-O, strongly suggesting that the development in the later stage would depend on transcriptional regulation by AP2-O. Therefore, it seems that corroboratory regulation by these two gene regulation systems contributes to the transition from gametocytes to ookinetes in the mosquito midgut. Elucidation of gene regulation in this transition stage is a next important theme because it may deepen our understanding of the parasite lifecycle that is constituted from serial stage conversions in host animals.

We demonstrated that target information obtained by ChIP-seq provided an overview of the molecular events proceeding in ookinetes. Based on this target information, we identified novel genes important for midgut infection by this parasite. Importantly, two independent experiments demonstrated consistency of peaks and calculated targets, demonstrating the robustness of this method. This result suggests that analyzing the entire set of target genes of a transcription factor (or targetome analysis) is a powerful way to study the biology of this parasite. Obviously, success in this attempt with ookinetes occurred largely because of the unique transcriptional regulatory features described above. In fact, if ookinete gene expression was regulated by a network of sequence-specific transcription factors, only a partial view of this stage would have been possible. Our results suggest that this method could be a potent omics tool for studying the biology of malarial parasites.

In conclusion, we report the first application of ChIP-seq to genome-wide identification of transcription factor targets in a malaria parasite. We determined the entire set of target genes of a malaria parasite transcription factor and elucidated how the AP2 family transcription factor contributes to formation of the motile stage. In addition, we revealed a unique gene regulatory system that is employed in ookinetes and possibly in other stages of the parasitic lifecycle.

Materials and Methods

Ethics statement

This study was carried out in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. The protocol was approved by the Committee on the Ethics of Animal Experiments of the Mie University (Permit Number: 23–29). All efforts were made to minimize animal suffering during the course of these studies.

Parasite preparations

The ANKA strain of P. berghei was maintained in female BALB/c mice (6–10 weeks old; Japan SLC, Inc., Hamamatsu, Japan). Ookinete culturing was performed as previously described [7]. To examine the number of oocysts, infected mice were subjected to Anopheles stephensi mosquitoes. Fully engorged mosquitoes were selected and maintained at 20°C. The number of oocysts and oocyst sporozoites was evaluated 14 days after an infective blood meal.

ChIP-seq

For the ookinete culture, six mice were pre-treated with phenylhydrazine and then infected with P. berghei expressing GFP-fused AP2-O. Sulfadiazine was added to their drinking water (10 mg/L) in order to deplete asexual blood-stage parasites. When the exflagellation of the male gametes reached approximately 300 per 105 red blood cells, the blood was harvested and cultured in an ookinete culture medium for 16 h, and then fixed with 1% paraformaldehyde. Erythrocytes in the fixed culture were removed by lysis in 0.84% NH4Cl, and the remaining ookinetes were subjected to ChIP. ChIP was performed using the ChIP Assay Kit (Millipore) according to the manufacturer’s protocol. Briefly, samples in the lysis buffer were sonicated with a Bioruptor (Tosho Denki, Yokohama, Japan) until chromatin DNA was fragmented to 300–500-bp in size for sequencing with an Illumina Genome Analyzer and to 150-bp for sequencing with a SOLiD 5500 system (Life Technologies). Immunoprecipitation (IP) was performed with anti-GFP antibodies, and the harvested DNA fragments were subjected to sequencing. Input DNAs were obtained from the chromatin without IP. Anti-GFP antibodies used for ChIP were purchased from Clontech and Abcam.

Analysis of ChIP-seq data

In experiment 1, sequence data obtained with an Illumina Genome Analyzer were mapped onto the P. berghei genome sequence (PlasmoDB, version 12.0) using the Bowtie program under conditions allowing one mismatch within 35-bp. In experiment 2, sequence data obtained with a SOLiD 5500 system were mapped onto the P. berghei genome sequence with a lifescope program equipped with the sequencing system in the default conditions and then filtered under more stringent conditions allowing no mismatches within 60-bp using an in-house program. The mapping data were analyzed with the MACS2 peak-calling algorithm using approximately 0.22 × 107 reads for IP and 0.77 × 107 reads for input control in experiment 1 or 0.22 × 107 reads for IP and 1.12 × 107 reads for input in experiment 2. Conditions for peak calling included FDR < 0.01 and fold enrichments over input control > 5 in both analyses. Genes were identified as AP2-O targets when their 1.2-kbp upstream regions contained the nearest binding motifs from the predicted summits of the ChIP-seq peaks. When the upstream region was less than 1.2-kbp, the entire intergenic region was used for target prediction. The ChIP-seq data have been deposited to Gene Expression Omnibus (GEO) with the accession no. GSE58584. Gene ID and functional annotation were attributed to each gene according to those in PlasmoDB ver.12.0.

To identify the binding sequences of AP2-O, six-base sequences concentrated around the predicted summits were investigated. Fisher’s exact test was performed between the 200-bp regions that have summits in the center and 200-bp regions excised from the genome excluding the former regions, to cover the entire genome sequence. Six-base sequences were ordered according to the calculated p-values, and common sequence motifs were searched among the sequences with the least p-values.

RNA-seq

Ookinetes were obtained from asexual parasite-depleted infected mouse blood that had been cultured for 24 h in an ookinete culture medium. Total RNA was extracted using the RNeasy mini kit (Qiagen, Hilden, Germany) according to the manufacturer’s protocols. Poly (A)+ RNA was purified using the Oligotex-dT30 mRNA purification kit (Takara Bio, Japan). The harvested RNA was used for sequencing with the SOLiD sequencing system. cDNA libraries for sequencing were prepared according to manufacturer protocols. Reads were mapped on the P. berghei genome, allowing no mismatches within 60-bp. The total number of reads mapped onto the genome were 15,740,028. The RNA-seq data have been deposited to GEO with the accession no. GSE58584.

DNA microarray

RNA samples used for the DNA microarray analysis were identical as those used in the previous study (seven biologically independent AP2-O(-) and five biologically independent wild-type ookinete samples) [7]. The DNA microarray experiments were performed with the one color method using a custom chip designed on a Agilent platform, as previously described [11]. The data were analyzed by using the GeneSpring software (Agilent Technologies, Santa Clara, CA), and genes whose expression was reduced more than two-fold relative to the wild-type were selected. Genes located at subtelomeric regions were excluded from the analysis. The microarray data have been deposited to GEO with the accession no. GSE58584.

Gene-targeting experiments

Gene-targeting experiments were carried out essentially using the same procedure as described previously [7]. The genotypes of cloned parasites were checked by PCR. Primers used for the preparation of targeting constructs and for the purpose of genotyping are listed in S14 Table in S2 File. The infectivity of parasites was estimated by examining the oocyst number. When the oocyst number was normal in one disruptant clone, the phenotype was not investigated further. When this number decreased, relative to that of the corresponding wild-type parasites, mutant parasites were obtained by another transfection experiment, and the phenotype and the genotype were subsequently determined.

Analysis of protein targeting using P. berghei centromere plasmid

Protein targeting was investigated by expressing each gene as a GFP-tagged protein, using the centromere plasmid pCen-GFP. Target genes with the upstream regulatory region (1.1–1.5-kbp upstream of the first methionine codon) were amplified by PCR using genomic DNA as the template, and then subcloned into the region immediately upstream of the GFP gene of the pCen vector. Merozoites were transfected with these constructs as reported previously [7]. The primers used for the preparation of the pCen constructs are listed in S14 Table in S2 File. The micrographs were obtained with an Olympus BX60 fluorescence microscope (Olympus, Tokyo, Japan) with a DXM1200C digital color camera (Nikon Corporation, Tokyo, Japan).

Preparation of the centromere plasmid vector with a sulfadiazine-resistant selectable marker

The sulfadiazine-resistant P. falciparum DHPS (Dihydropteroate synthase) gene (Lys-460 to Glu) was generated by PCR using the primers listed in S14 Table in S2 File, using P. falciparum 3D7 genomic DNA as a template. The P. berghei centromere plasmid pCen-GFP-mDHPS was prepared from the P. berghei heat shock protein promoter, the GFP gene, the P. berghei heat shock protein 70 termination sequence, the P. berghei elongation factor 1 alpha promoter, the Sulfadiazine-resistant P. falciparum DHPS gene, the P. berghei DHFR-ts (Dihydrofolate reductase-thymidylate synthase) transcription termination sequence, and the plasmid vector pBluescript II SK+ (Agilent Technologies). Pyrimethamine-resistant mutant parasites were transfected with this centromere plasmid, and transfectants were selected in mice by adding sulfadiazine to the drinking water (10 mg/mL). Micrographs were obtained with an Olympus BX60 fluorescence microscope. The oocyst diameter was measured with the AquaCosmos software (Hamamatsu Photonic System).

Immunoelectron microscopy

Immunoelectron microscopy was performed as described previously [37]. Briefly, ookinetes were fixed in 0.1 M phosphate buffer (pH 7.4) containing 1% paraformaldehyde and 0.1% glutaraldehyde (TAAB), after a 24-hour culture. They were dehydrated in ethanol and embedded in LR Gold resin (London Resin Company, UK). Ultrathin sections were blocked for 30 min in PBS containing 0.01% Tween 20 and 5% non-fat dry milk, incubated with anti-GFP antibodies, and then with goat anti-rabbit IgG conjugated to gold particles of 15 nm diameter (Amersham Pharmacia Biotech) diluted in a blocking buffer. Finally, the sections were fixed with 2.5% glutaraldehyde for 10 min and stained with 2% uranyl acetate and Reynold’s lead citrate. The anti-GFP antibodies used were identical to those used in the ChIP assay.

Supporting Information

Attachment 1

Attachment 2

Attachment 3

Attachment 4

Attachment 5

Attachment 6

Attachment 7

Attachment 8

Attachment 9

Attachment 10

Attachment 11


Zdroje

1. Hall N, Karras M, Raine JD, Carlton JM, Kooij TW, et al. (2005) A comprehensive survey of the Plasmodium life cycle by genomic, transcriptomic, and proteomic analyses. Science 307: 82–86. 15637271

2. Balaji S, Babu MM, Iyer LM, Aravind L (2005) Discovery of the principal specific transcription factors of Apicomplexa and their implication for the evolution of the AP2-integrase DNA binding domains. Nucleic Acids Res 33: 3994–4006. 16040597

3. De Silva EK, Gehrke AR, Olszewski K, Leon I, Chahal JS, et al. (2008) Specific DNA-binding by apicomplexan AP2 transcription factors. Proc Natl Acad Sci U S A 105: 8393–8398. doi: 10.1073/pnas.0801993105 18541913

4. Reece-Hoyes JS, Deplancke B, Shingles J, Grove CA, Hope IA, et al. (2005) A compendium of Caenorhabditis elegans regulatory transcription factors: a resource for mapping transcription regulatory networks. Genome Biol 6: R110. 16420670

5. Adryan B, Teichmann SA (2006) FlyTF: a systematic review of site-specific transcription factors in the fruit fly Drosophila melanogaster. Bioinformatics 22: 1532–1533. 16613907

6. Hughes TR, de Boer CG (2013) Mapping yeast transcriptional networks. Genetics 195: 9–36. doi: 10.1534/genetics.113.153262 24018767

7. Yuda M, Iwanaga S, Shigenobu S, Mair GR, Janse CJ, et al. (2009) Identification of a transcription factor in the mosquito-invasive stage of malaria parasites. Mol Microbiol 71: 1402–1414. doi: 10.1111/j.1365-2958.2009.06609.x 19220746

8. Yuda M, Iwanaga S, Shigenobu S, Kato T, Kaneko I (2010) Transcription factor AP2-Sp and its target genes in malarial sporozoites. Mol Microbiol 75: 854–863. doi: 10.1111/j.1365-2958.2009.07005.x 20025671

9. Kafsack BF, Rovira-Graells N, Clark TG, Bancells C, Crowley VM, et al. (2014) A transcriptional switch underlies commitment to sexual development in malaria parasites. Nature 507: 248–252. doi: 10.1038/nature12920 24572369

10. Sinha A, Hughes KR, Modrzynska KK, Otto TD, Pfander C, et al. (2014) A cascade of DNA-binding proteins for sexual commitment and development in Plasmodium. Nature 507: 253–257. doi: 10.1038/nature12970 24572359

11. Iwanaga S, Kaneko I, Kato T, Yuda M (2012) Identification of an AP2-family protein that is critical for malaria liver stage development. PLoS One 7: e47557. doi: 10.1371/journal.pone.0047557 23144823

12. Cheng C, Min R, Gerstein M (2011) TIP: a probabilistic method for identifying transcription factor target genes from ChIP-seq binding profiles. Bioinformatics 27: 3221–3227. doi: 10.1093/bioinformatics/btr552 22039215

13. Wu S, Wang J, Zhao W, Pounds S, Cheng C (2010) ChIP-PaM: an algorithm to identify protein-DNA interaction using ChIP-Seq data. Theor Biol Med Model 7: 18. doi: 10.1186/1742-4682-7-18 20525272

14. Feng J, Liu T, Qin B, Zhang Y, Liu XS (2012) Identifying ChIP-seq enrichment using MACS. Nat Protocols 7: 1728–1740. doi: 10.1038/nprot.2012.101 22936215

15. Mair GR, Braks JA, Garver LS, Wiegant JC, Hall N, et al. (2006) Regulation of sexual development of Plasmodium by translational repression. Science 313: 667–669. 16888139

16. Lal K, Prieto JH, Bromley E, Sanderson SJ, Yates JR 3rd, et al. (2009) Characterisation of Plasmodium invasive organelles; an ookinete microneme proteome. Proteomics 9: 1142–1151. doi: 10.1002/pmic.200800404 19206106

17. Patra KP, Johnson JR, Cantin GT, Yates JR 3rd, Vinetz JM (2008) Proteomic analysis of zygote and ookinete stages of the avian malaria parasite Plasmodium gallinaceum delineates the homologous proteomes of the lethal human malaria parasite Plasmodium falciparum. Proteomics 8: 2492–2499. doi: 10.1002/pmic.200700727 18563747

18. Hall N, Carlton J (2005) Comparative genomics of malaria parasites. Curr Opin Genet Dev 15: 609–613. 16182520

19. Ecker A, Bushell ES, Tewari R, Sinden RE (2008) Reverse genetics screen identifies six proteins important for malaria development in the mosquito. Mol Microbiol 70: 209–220. doi: 10.1111/j.1365-2958.2008.06407.x 18761621

20. Gilk SD, Raviv Y, Hu K, Murray JM, Beckers CJ, et al. (2006) Identification of PhIL1, a novel cytoskeletal protein of the Toxoplasma gondii pellicle, through photosensitized labeling with 5-[125I]iodonaphthalene-1-azide. Eukaryot Cell 5: 1622–1634. 17030994

21. Barkhuff WD, Gilk SD, Whitmarsh R, Tilley LD, Hunter C, et al. (2011) Targeted disruption of TgPhIL1 in Toxoplasma gondii results in altered parasite morphology and fitness. PLoS One 6: e23977. doi: 10.1371/journal.pone.0023977 21901148

22. Khater EI, Sinden RE, Dessens JT (2004) A malaria membrane skeletal protein is essential for normal morphogenesis, motility, and infectivity of sporozoites. J Cell Biol 167: 425–432. 15533999

23. Tremp AZ, Khater EI, Dessens JT (2008) IMC1b is a putative membrane skeleton protein involved in cell shape, mechanical strength, motility, and infectivity of malaria ookinetes. J Biol Chem 283: 27604–27611. doi: 10.1074/jbc.M801302200 18650444

24. Tran JQ, Li C, Chyan A, Chung L, Morrissette NS (2012) SPM1 stabilizes subpellicular microtubules in Toxoplasma gondii. Eukaryot Cell 11: 206–216. doi: 10.1128/EC.05161-11 22021240

25. Beck JR, Rodriguez-Fernandez IA, Cruz de Leon J, Huynh MH, Carruthers VB, et al. (2010) A novel family of Toxoplasma IMC proteins displays a hierarchical organization and functions in coordinating parasite division. PLoS Pathog 6: e1001094. doi: 10.1371/journal.ppat.1001094 20844581

26. Poulin B, Patzewitz EM, Brady D, Silvie O, Wright MH, et al. (2013) Unique apicomplexan IMC sub-compartment proteins are early markers for apical polarity in the malaria parasite. Biol Open 2: 1160–1170. doi: 10.1242/bio.20136163 24244852

27. Frenal K, Polonais V, Marq JB, Stratmann R, Limenitakis J, et al. (2010) Functional dissection of the apicomplexan glideosome molecular architecture. Cell Host Microbe 8: 343–357. doi: 10.1016/j.chom.2010.09.002 20951968

28. Bullen HE, Tonkin CJ, O'Donnell RA, Tham WH, Papenfuss AT, et al. (2009) A novel family of Apicomplexan glideosome-associated proteins with an inner membrane-anchoring role. J Biol Chem 284: 25353–25363. doi: 10.1074/jbc.M109.036772 19561073

29. Tremp AZ, Dessens JT (2011) Malaria IMC1 membrane skeleton proteins operate autonomously and participate in motility independently of cell shape. J Biol Chem 286: 5383–5391. doi: 10.1074/jbc.M110.187195 21098480

30. Kaiser K, Camargo N, Coppens I, Morrisey JM, Vaidya AB, et al. (2004) A member of a conserved Plasmodium protein family with membrane-attack complex/perforin (MACPF)-like domains localizes to the micronemes of sporozoites. Mol Biochem Parasitol 133: 15–26. 14668008

31. Yuda M, Sawai T, Chinzei Y (1999) Structure and expression of an adhesive protein-like molecule of mosquito invasive-stage malarial parasite. J Exp Med 189: 1947–1952. 10377190

32. Dessens JT, Siden-Kiamos I, Mendoza J, Mahairaki V, Khater E, et al. (2003) SOAP, a novel malaria ookinete protein involved in mosquito midgut invasion and oocyst development. Mol Microbiol 49: 319–329. 12828632

33. Yuda M, Yano K, Tsuboi T, Torii M, Chinzei Y (2001) von Willebrand Factor A domain-related protein, a novel microneme protein of the malaria ookinete highly conserved throughout Plasmodium parasites. Mol Biochem Parasitol 116: 65–72. 11463467

34. Hinds L, Green JL, Knuepfer E, Grainger M, Holder AA (2009) Novel putative glycosylphosphatidylinositol-anchored micronemal antigen of Plasmodium falciparum that binds to erythrocytes. Eukaryot Cell 8: 1869–1879. doi: 10.1128/EC.00218-09 19820120

35. Kadota K, Ishino T, Matsuyama T, Chinzei Y, Yuda M (2004) Essential role of membrane-attack protein in malarial transmission to mosquito host. Proc Natl Acad Sci U S A 101: 16310–16315. 15520375

36. Huber M, Cabib E, Miller LH (1991) Malaria parasite chitinase and penetration of the mosquito peritrophic membrane. Proc Natl Acad Sci U S A 88: 2807–2810. 2011589

37. Kariu T, Ishino T, Yano K, Chinzei Y, Yuda M (2006) CelTOS, a novel malarial protein that mediates transmission to mosquito and vertebrate hosts. Mol Microbiol 59: 1369–1379. 16468982

38. Deligianni E, Morgan RN, Bertuccini L, Wirth CC, de Monerri NC, et al. (2013) A perforin-like protein mediates disruption of the erythrocyte membrane during egress of Plasmodium berghei male gametocytes. Cell Microbiol 15: 1438–1455. doi: 10.1111/cmi.12131 23461714

39. Ishino T, Chinzei Y, Yuda M (2005) 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 7: 199–208. 15659064

40. Ecker A, Pinto SB, Baker KW, Kafatos FC, Sinden RE (2007) Plasmodium berghei: Plasmodium perforin-like protein 5 is required for mosquito midgut invasion in Anopheles stephensi. Exp Parasitol 116: 504–508. 17367780

41. Kondos SC, Hatfaludi T, Voskoboinik I, Trapani JA, Law RH, et al. (2010) The structure and function of mammalian membrane-attack complex/perforin-like proteins. Tissue Antigens 76: 341–351. doi: 10.1111/j.1399-0039.2010.01566.x 20860583

42. Marois E (2011) The multifaceted mosquito anti-Plasmodium response. Curr Opin Microbiol 14: 429–435. doi: 10.1016/j.mib.2011.07.016 21802348

43. Molina-Cruz A, DeJong RJ, Charles B, Gupta L, Kumar S, et al. (2008) Reactive oxygen species modulate Anopheles gambiae immunity against bacteria and Plasmodium. J Biol Chem 283: 3217–3223. 18065421

44. Jortzik E, Becker K (2012) Thioredoxin and glutathione systems in Plasmodium falciparum. Int J Med Microbiol 302: 187–194. doi: 10.1016/j.ijmm.2012.07.007 22939033

45. Vega-Rodriguez J, Franke-Fayard B, Dinglasan RR, Janse CJ, Pastrana-Mena R, et al. (2009) The glutathione biosynthetic pathway of Plasmodium is essential for mosquito transmission. PLoS Pathog 5: e1000302. doi: 10.1371/journal.ppat.1000302 19229315

46. Pastrana-Mena R, Dinglasan RR, Franke-Fayard B, Vega-Rodriguez J, Fuentes-Caraballo M, et al. (2010) Glutathione reductase-null malaria parasites have normal blood stage growth but arrest during development in the mosquito. J Biol Chem 285: 27045–27056. doi: 10.1074/jbc.M110.122275 20573956

47. Srinivasan P, Fujioka H, Jacobs-Lorena M (2008) PbCap380, a novel oocyst capsule protein, is essential for malaria parasite survival in the mosquito. Cell Microbiol 10: 1304–1312. doi: 10.1111/j.1462-5822.2008.01127.x 18248630

48. Zhang M, Fennell C, Ranford-Cartwright L, Sakthivel R, Gueirard P, et al. (2010) The Plasmodium eukaryotic initiation factor-2alpha kinase IK2 controls the latency of sporozoites in the mosquito salivary glands. J Exp Med 207: 1465–1474. doi: 10.1084/jem.20091975 20584882

49. Ning J, Otto TD, Pfander C, Schwach F, Brochet M, et al. (2013) Comparative genomics in Chlamydomonas and Plasmodium identifies an ancient nuclear envelope protein family essential for sexual reproduction in protists, fungi, plants, and vertebrates. Genes Dev 27: 1198–1215. doi: 10.1101/gad.212746.112 23699412

50. Bushell ES, Ecker A, Schlegelmilch T, Goulding D, Dougan G, et al. (2009) Paternal effect of the nuclear formin-like protein MISFIT on Plasmodium development in the mosquito vector. PLoS Pathog 5: e1000539. doi: 10.1371/journal.ppat.1000539 19662167

51. Ishino T, Orito Y, Chinzei Y, Yuda M (2006) A calcium-dependent protein kinase regulates Plasmodium ookinete access to the midgut epithelial cell. Mol Microbiol 59: 1175–1184. 16430692

52. Tewari R, Straschil U, Bateman A, Bohme U, Cherevach I, et al. (2010) The systematic functional analysis of Plasmodium protein kinases identifies essential regulators of mosquito transmission. Cell Host Microbe 8: 377–387. doi: 10.1016/j.chom.2010.09.006 20951971

53. Sebastian S, Brochet M, Collins MO, Schwach F, Jones ML, et al. (2012) A Plasmodium calcium-dependent protein kinase controls zygote development and transmission by translationally activating repressed mRNAs. Cell Host Microbe 12: 9–19. doi: 10.1016/j.chom.2012.05.014 22817984

54. Segal E, Shapira M, Regev A, Pe'er D, Botstein D, et al. (2003) Module networks: identifying regulatory modules and their condition-specific regulators from gene expression data. Nat Genet 34: 166–176. 12740579

55. Iyer LM, Anantharaman V, Wolf MY, Aravind L (2008) Comparative genomics of transcription factors and chromatin proteins in parasitic protists and other eukaryotes. Int J Parasitol 38: 1–31. 17949725

56. Thorvaldsdottir H, Robinson JT, Mesirov JP (2013) Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform 14: 178–192. doi: 10.1093/bib/bbs017 22517427

57. Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, et al. (2008) Model-based analysis of ChIP-Seq (MACS). Genome Biol 9: R137. doi: 10.1186/gb-2008-9-9-r137 18798982

58. Crooks GE, Hon G, Chandonia JM, Brenner SE (2004) WebLogo: a sequence logo generator. Genome Res 14: 1188–1190. 15173120

59. Orosz F (2009) Apicortin, a unique protein, with a putative cytoskeletal role, shared only by apicomplexan parasites and the placozoan Trichoplax adhaerens. Infect Genet Evol 9: 1275–1286. doi: 10.1016/j.meegid.2009.09.001 19778640

60. Peris L, Thery M, Faure J, Saoudi Y, Lafanechere L, et al. (2006) Tubulin tyrosination is a major factor affecting the recruitment of CAP-Gly proteins at microtubule plus ends. J Cell Biol 174: 839–849. 16954346

61. Montagna GN, Buscaglia CA, Munter S, Goosmann C, Frischknecht F, et al. (2012) Critical role for heat shock protein 20 (HSP20) in migration of malarial sporozoites. J Biol Chem 287: 2410–2422. doi: 10.1074/jbc.M111.302109 22139844

Štítky
Hygiena a epidemiologie Infekční lékařství Laboratoř

Článek vyšel v časopise

PLOS Pathogens


2015 Číslo 5
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy

Zvyšte si kvalifikaci online z pohodlí domova

Důležitost adherence při depresivním onemocnění
nový kurz
Autoři: MUDr. Eliška Bartečková, Ph.D.

Koncepce osteologické péče pro gynekology a praktické lékaře
Autoři: MUDr. František Šenk

Sekvenční léčba schizofrenie
Autoři: MUDr. Jana Hořínková, Ph.D.

Hypertenze a hypercholesterolémie – synergický efekt léčby
Autoři: prof. MUDr. Hana Rosolová, DrSc.

Multidisciplinární zkušenosti u pacientů s diabetem
Autoři: Prof. MUDr. Martin Haluzík, DrSc., prof. MUDr. Vojtěch Melenovský, CSc., prof. MUDr. Vladimír Tesař, DrSc.

Všechny kurzy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#