- Open Access
Alternative splicing is required for stage differentiation in malaria parasites
Genome Biology volume 20, Article number: 151 (2019)
In multicellular organisms, alternative splicing is central to tissue differentiation and identity. Unicellular protists lack multicellular tissue but differentiate into variable cell types during their life cycles. The role of alternative splicing in transitions between cell types and establishing cellular identity is currently unknown in any unicellular organism.
To test whether alternative splicing in unicellular protists plays a role in cellular differentiation, we conduct RNA-seq to compare splicing in female and male sexual stages to asexual intraerythrocytic stages in the rodent malaria parasite Plasmodium berghei. We find extensive changes in alternative splicing between stages and a role for alternative splicing in sexual differentiation. Previously, general gametocyte differentiation was shown to be modulated by specific transcription factors. Here, we show that alternative splicing establishes a subsequent layer of regulation, controlling genes relating to consequent sex-specific differentiation of gametocytes.
We demonstrate that alternative splicing is reprogrammed during cellular differentiation of a unicellular protist. Disruption of an alternative splicing factor, PbSR-MG, perturbs sex-specific alternative splicing and decreases the ability of the parasites to differentiate into male gametes and oocysts, thereby reducing transmission between vertebrate and insect hosts. Our results reveal alternative splicing as an integral, stage-specific phenomenon in these protists and as a regulator of cellular differentiation that arose early in eukaryotic evolution.
Alternative splicing is a phenomenon where different mature transcripts are created from a single pre-mRNA species. In some multicellular eukaryotes, alternative splicing is essential for cellular differentiation to form specific tissues. Alternative splicing generates transcript variants needed to drive cells down a developmental pathway and to maintain the tissue after differentiation has occurred [1, 2]. Indeed, most of the alternative splicing that affects 90% of human genes is tissue-specific [3, 4]. Many unicellular protists undergo stage-specific cellular differentiation, but reported levels of alternative splicing are much lower than in metazoans. For example, the obligate intracellular parasitic eukaryotes of the phylum Apicomplexa [5, 6] have reported splicing rates ranging from 4.5% in Plasmodium falciparum  to over 20% in Toxoplasma gondii . The lower levels of alternative splicing observed may reflect real biological differences, but could also reflect the fact that only a single life stage was used in earlier studies, masking higher levels of alternative splicing that occur during life stage progression.
The transcriptome and proteome of P. falciparum are not congruent. Abundance and timing of mRNA appearance correlate poorly with the expression of corresponding proteins [9–11]. This suggests that post-transcriptional regulation processes, perhaps including alternative splicing, are central to apicomplexan biology. Hence, we tested the hypothesis that alternative splicing in unicellular protists plays an important role in life stage differentiation, a phenomenon analogous to tissue-specific splicing in metazoans.
Plasmodium species, which cause malaria, are ideal model organisms in which to address this question. These parasites have a complicated life cycle, spanning multiple, morphologically distinct stages in their mammalian and insect hosts. The symptomatic stage in the vertebrate host is the asexual blood cycle, where parasites proliferate rapidly in erythrocytes. Some blood-stage parasites irreversibly convert into sexually dimorphic gametocytes . After ingestion by mosquitoes, female gametocytes activate and egress from their host cell, while male gametocytes produce eight motile gametes (sperm) in a process known as exflagellation . Gametes fuse to form a zygote, which metamorphoses into an ookinete, and subsequently establishes an oocyst in the midgut wall . Sporozoites develop within the oocyst; these migrate to the mosquito salivary glands and are capable of infection of a new vertebrate host  via invasion and development of parasites within liver cells in mammalian hosts . The mouse malaria model, P. berghei, provides laboratory access to all these life stages, facilitating analysis of the regulation of stage differentiation. In addition, many introns are conserved between Plasmodium species, making this an ideal model organism to provide insights into human parasites.
Transcriptional reprogramming plays a significant role in the transition from asexual blood stage to gametocytes, with thousands of genes differentially expressed between the stages . 13.1% of all differentially expressed genes in gametocytes were differentially expressed simultaneously in both female and male gametocytes compared to their asexual developmental precursors; this excludes genes that differ when comparing female and male gametocytes directly and, hence, indicates the “common” genes differentially expressed in initial sexual commitment (Fig. 1a, data derived from ). Initial differentiation of asexual parasites into gametocytes is triggered by the expression of transcription factors AP2-G and AP2-G2, which regulate whole-gene expression . The importance of these transcription factors in initiating stage transition is emphasized by the increase in asexual growth seen after the ablation of AP2-G . This increase in asexual growth results from a loss of non-replicative gametocytes from the population .
Results and discussion
We used the well-established differential splicing tool, DEXSeq, to examine alternative splicing between these stages . We detected changes in alternative splicing in thousands of genes in both gametocyte sexes compared to their asexual precursors (see Additional file 1: Table S1 and Additional file 2: able S2 for the full list). We detected changes in alternative splicing common to both sexes in 11.4% of genes, again excluding changes between both sexes when compared pairwise (Fig. 1b). Splicing alterations included diverse categories of alternative splicing, including exon skipping, alternative 5′ and 3′ splice sites, and intron retention. In many cases, the splice changes are predicted to result in truncated or aberrant proteins. These changes in splicing between stages and sexes are viewable for mapped reads preloaded for individual genes using the GBrowse genome browser at www.plasmodb.org . Almost 90% of the changes in alternative splicing events were unique to either male or female gametocytes, suggesting an important role for alternative splicing in sexual determination after the initial commitment to gametocyte development. Indeed, while general gametocyte differentiation is clearly modulated by specific transcriptional factors , our analysis indicates that alternative splicing may establish a subsequent layer of regulation, controlling genes relating to consequent sex-specific differentiation of gametocytes.
We attempted to validate changes in alternative splicing by analysis of 15 genes using qRT-PCR. Five samples failed due to poor, non-specific amplification, but ten reactions worked well. In four samples, we predicted higher abundance in a specific sub-exonic bin in female compared to asexual parasites. We predicted lower abundance in one sample, and no change in five control genes. qRT-PCR analysis verified our predictions in all ten cases (Fig. 1c, Additional file 3: Table S3, Additional file 4: Figure S1).
Analysis of sex-specific splicing changes in individual genes reveals changes in several genes that are required for appropriate differentiation (see Additional file 1: Table S1, Additional file 2: Table S2). In female gametocytes, changes in alternative splicing were observed in a major epigenetic regulator histone acetyltransferase HAT1 (PBANKA_0718400) , such that retention of the final intron is decreased compared to asexual stages (Fig. 1b). Histone acetyltransferases have been implicated in cell cycle regulation , which is consistent with a role in differentiation into female gametocytes.
Differentiation to male gametocytes is accompanied by alternative splicing changes in the male gamete fusion factor HAP2, where we detected relatively more transcripts with a longer first intron in male gametocytes compared to blood stages (PBANKA_1212600, Fig. 1b). In P. berghei, HAP2 is essential for fusion of gamete membranes and, hence, fertilisation . HAP2 is upregulated 20-fold in male gametocytes , and thus, alternative splicing appears to regulate the function of this gene on an additional level.
To explicitly test whether alternative splicing was required to regulate differentiation of gametes, we genetically ablated a putative gametocyte-specific alternative splicing factor, using the PlasmoGEM system, replacing the entire coding sequence of the gene [23, 24]. Serine-arginine–rich (SR) proteins are major regulators of alternative splicing, with regular expression levels necessary for maintaining normal alternative splicing and viability [8, 25]. One SR protein, PBANKA_0921600, is specifically upregulated in male gametocytes of P. berghei and downregulated in female gametocytes , and phylogenetic analysis identifies this gene as an orthologue of a known splicing factor in T. gondii . We named PBANKA_0921600 PbSR-MG (P.berghei SR protein male gametocytes). Localizing an integrated, hemagglutinin-tagged (HA-tagged) PbSR-MG confirmed its expected nuclear localization (Fig. 2a).
To examine the role of PbSR-MG in regulating splicing, we ablated the gene, obtaining monoclonal parasites (Fig. 2b). We were then able to compare changes in alternative splicing between clonal mutant lines (PbSR-MG KO) and parental lines. RNA-seq analysis revealed differences in alternative splicing to be heavily skewed by sex. Ten percent of genes had modified alternative splicing in male gametocytes, while only 3.9% of female-specific and 0.90% of asexual-stage genes were affected (Fig. 2c, see Additional file 5: Table S4, Additional file 6: Table S5, Additional file 7: Table S6 for full list). Loss of PbSR-MG had minimal impact on transcription, altering the overall expression of only 0.5% of genes (see Additional file 8: Table S7 for full list).
We could now compare all male-specific changes in alternative splicing (Fig. 1b, 264+974=1238 genes total) to male genes specifically regulated by PbSR-MG (Fig. 2c, black 10% region, 314 genes total). There were 247 genes that had both male-specific changes in alternative splicing and were male genes regulated by PbSR-MG. Hence, 79% (247/314) of the genes regulated by PbSR-MG were originally identified as having alternative splicing changes specific to male gametocytes. From the opposing perspective, only 20% (247/1238) of all male-specific alternative splicing changes are regulated by PbSR-MG. This contrasts with female gametocytes, where only 5.4% (73/1340) of the female-specific alternative splicing changes are affected by the deletion of PbSR-MG. Presumably, there are other unidentified factors regulating female-specific alternative splicing.
In the analyses presented above, we investigated the changes in alternative splicing between pairs of samples. To assess the total levels of alternative splicing in individual stages, we employed our custom-built JunctionJuror computational tool . This method estimates the abundance of alternative splicing in single samples, although it explicitly underestimates the total amount, because it excludes instances of intron retention. In wild-type parasites, we detected alternative splicing in 800–1100 genes, after the analysis of each sample independently (Fig. 2d, see Additional file 9: Tables S8 for full list). When PbSR-MG was ablated, the number of alternatively spliced genes detected was relatively unchanged for asexual parasites, modestly increased for female gametocytes, and most increased for male gametocytes (Fig. 2d). Many genes were alternatively spliced in both wild-type and mutant parasites (Fig. 2d, striped regions), suggesting that much alternative splicing is regulated by factors additional to PbSR-MG in all three stages.
Ablation of PbSR-MG is associated with increased alternative splicing in males (Fig. 2d, male, WT vs. KO). There are several interpretations consistent with this. PbSR-MG may enhance spliceosomal specificity at alternative splicing sites. When PbSR-MG is absent, the spliceosome acts imprecisely, resulting in aberrant alternative splicing and the production of additional isoforms. Alternatively, PbSR-MG may act as a general suppressor of alternative splicing but this would imply that suppression of alternative splicing is important for male gametocytes. We observed no clear decrease in male alternative splicing compared to the other stages (Fig. 2d, WT, male vs. asexual), so this seems a less compelling mechanism, although it is possible that other splicing factors also regulate alternative splicing in male gametocytes.
Having established changes in alternative splicing as a major phenomemon underlying stage differentiation via the screens above, we wanted to quantify the types of alternative splicing involved. There are presently no integrated computational tools for this kind of analysis, so we applied two different approaches. Unlike our analyses above (Fig. 2c, d), which are agnostic to canonical gene models, these new analyses rely on detecting the differences between observed reads and expected gene models. There are likely to be some differences in genes identified because of these different approaches, but the gene model-based methods do allow us to categorize which types of alternative splicing events occur. We quantified intron retention utilizing featureCounts, which counts the reads that map to defined regions of a gene , and quantified alternative 5′ and 3′ splice site usage and exon skipping using Junction annotation, which allows us to identify junctions that match or conflict with the canonical gene model . We attempted to utilize consistent thresholds with these two methods, so that results between the two software packages were comparable.
In single stages, intron retention accounted for over 80% of alternatively spliced genes (Fig. 3, colored bar graphs; Additional file 10: Tables S9). Similarly, when looking at alternative splicing that changed between stages, there was also a strong over-representation of intron retention (Fig. 3, gray bar graphs; Additional file 11: Tables S10). In all cases, alternative 5′ and 3′ splice site usage occured in far fewer genes, and exon skipping contributed a very minor proportion. This distribution of alternative splice types is consistent with the results in plants and fungi [28, 29], where intron retention is the dominant form of alternative splicing, and differs from animals, where intron retention is relatively rare , although these analyses did not correlate the potential effect of relative intron size with proportion of intron retention, which might be a confounding factor.
These data were consistent with the broader screen for changes in alternative splicing presented earlier (Additional file 1: Table S1, Additional file 2: Table S2). In the transition from asexual to females, 87% (761/875) of the genes identified in our pipeline were also identified in the screen above. In males, we observed 83% (578/699) overlap.
We also compared wild-type and PbSR-MG KO parasites to analyze the changes in different categories of alternative splicing. We identified only a single gene in both asexual parasites (5′ splice site change) and female gametocytes (exon skipping), but observed 60 genes in male gametocytes, with 78% of these identified as undergoing intron retention. These proportions are consistent with the DEXSeq analysis (Fig. 2c), although numbers are smaller and likely reflect the ability of DEXSeq to detect alternative splicing changes that do not correspond to canonical gene models. The gene model-based methods additionally incorporate a fold change threshold for splicing changes, and this also reduces the number of splice changes that are flagged. Hence, the DEXSeq analyses are more likely to capture lower fold change variations, such as those observed in the PbSR-MG KO female parasites.
To determine the biological consequences of changes to alternative splicing caused by the loss of PbSR-MG, we compared the viability of wild-type and PbSR-MG KO parasites through the entire life cycle (Fig. 4). Deletion of PbSR-MG has no detectable effect in asexual blood stages (Fig. 4a, ) nor in the proportion of female to male gametocytes produced (Fig. 4b). The proportion of female gametocytes that can be activated is also unchanged (Fig. 4c), but the exflagellation of male gametocytes is severely impacted in PbSR-MG KO parasites (Fig. 4d). When wild-type parasites exflagellate, flagella from the nascent gametes protrude conspicuously from the remnant gametocyte body (Fig. 5a). This occurred in up to 36% of wild-type males after in vitro induction (Fig. 4d). In comparison, we observed almost no exflagellation in PbSR-MG KO parasites (Fig. 4d). Flagella stained with an anti-tubulin antibody remain curled within the gametocyte body (Fig. 5a). Further, on the rare occasion when exflagellation of PbSR-MG KO male gametocytes could be observed, the motility of gametes was markedly impaired, almost to the point of immobility. A similar defect in exflagellation was reported in an SR protein kinase and is consistent with the requirement for phosphorylation to activate PbSR-MG proteins .
Ablation of PbSR-MG also resulted in fewer oocysts per mosquito host (Fig. 4e). To test if this difference solely reflected the defect in exflagellation, we repeated the mosquito infection using equal numbers of in vitro generated ookinetes. Despite exposing mosquitoes to equal infective load, we still observed a significant difference in oocyst numbers in PbSR-MG KO compared to wild-type parental parasites (Additional file 4: Figure S2).
No further defects were observed during sporozoite development or in the ability to infect a naive mouse (Fig. 4f, g). Despite the defect in exflagellation, parasites were still able to progress past this stage and complete their life cycle. The PbSR-MG exflagellation phenotype was consistent across clones from independent transfections into three different parasite lines and could be fully complemented by the introduction of full-length PbSR-MG on an artificial chromosome (Fig. 5b, Additional file 4: Figure S3). The oocyst growth defect could also be similarly complemented (Additional file 4: Figure S3).
PbSR-MG clearly plays an important role in stage differentiation, with the development of functional male gametes being significantly impaired by the loss of this splicing factor. There is also a second, less severe effect on the transition from ookinete to oocyst (Fig. 4e). This finding clearly links the changes in alternative splicing resulting from the loss of PbSR-MG with stage-specific development of male gametocytes.
Initial differentiation of asexual parasites into gametocytes is triggered by expression of transcription factors such as AP2-G and AP2-G2, which regulate whole-gene expression . PbSR-MG appears to be responsible for subsequent differentiation into specific sexes, via regulation of alternative splicing. This multi-stage regulation is further supported by asexual blood stage growth being unaffected by ablation of PbSR-MG, unlike AP2-G (, Fig. 4a). This provides additional evidence that PbSR-MG has no role in general gametocytogenesis, instead functioning downstream of this initial commitment.
The central role of post-transcriptional regulation in stage progression has previously been observed in female gametocytes, which regulate gene expression through protein complexes (DOZI and CITH) that stabilize transcripts and temporarily silence their translation . We demonstrate that post-transcriptional regulation by alternative splicing is key to the sexual differentiation of female and male gametes and that a specific alternative splicing factor is essential for the production of viable male gametocytes. Thus, analogous to the regulation of tissue specificity in metazoans, alternative splicing is integral to stage differentiation in these unicellular protists. This reveals alternative splicing as a common mechanism for regulating cellular differentiation in evolutionarily diverse eukaryotes.
Generation of knockout mutants
Vectors were transfected into parasites and transfectant parasites selected as previously described . Male Swiss Webster mice aged 4 to 6 weeks old were used in all experiments. Three independent transfections into three different parental parasite lines were performed. The first was into P. berghei ANKA wild-type parasites. The second was into P. berghei ANKA 820cl1m1cl1 parasites, which produce RFP and GFP in female and male gametocytes, respectively (a kind gift from Andy Waters, University of Glasgow, Glasgow, Scotland) . The final transfection was into P. berghei ANKA parasites constitutively expressing tdTomato .
After recovery of pyrimethamine-resistant parasites, correct integration of the cassette into the PbSR-MG locus was verified by PCR with the antisense primer “PbSR-MG ASP” and either sense primer “PbGEM GW2+ SP” or “PbSR-MG WT SP” (primer sequences in Additional file 12: Table S11). We generated monoclonal lines via limiting dilution in ten mice per mutant (see Fig. 2b and Additional file 4: Figure S4 for PCR results).
Female and male gametocyte mutants based on the 820cl1m1cl1 line , as well as asexual erythrocytic stages, were purified based on differential density and on fluorescent sorting and sequenced as previously reported .
Generation of epitope-tagged mutant
The hemagglutinin epitope-tagged line was created by modification of a pL0006-based plasmid . This plasmid was designed for single-site recombination into P. berghei and includes three hemagglutinin (HA) tags, followed by the 3′-UTR of KAHRP, followed by a DHFR selectable marker. We cloned the 3′ end of PbSR-MG from genomic DNA with the primers “PbSR-MG 3′-HindIII-ApaI SP” and “PbSR-MG 3′-AvrII-BglII ASP” (primer sequences in Additional file 12: Table S11) and ligated this into the vector backbone.
The integrity of the plasmid was confirmed by sequencing. The final vector was linearized with the restriction enzyme BsrGI-HF, and P. berghei ANKA wild-type parasites transfected as above.
Parasites were screened with the sense primer “PbSR-MG test SP” and either antisense primer “PbSR-MG 3′ UTR test ASP” or “pPb3HA 3’ UTR ASP” (primer sequences in Additional file 12: Table S11). The strongest signal in transfected parasites was for the mutant locus (see Additional file 4: Figure S4).
The epitope was visualized in blood-stage parasites via western blot. Samples were run on CF11 columns to deplete leukocytes, before being resuspended in 0.15% w/v saponin in Dulbecco’s phosphate-buffered saline and incubated for 10 min. Samples were washed four times, before resuspension in sample buffer and running on a NuPage SDS-PAGE gel (Invitrogen, Australia). Blots were probed with rat anti-HA (1/200; Roche, Australia, 18674) and anti-rat horseradish peroxidase (1/1000; Pierce, Australia, 31470). The apparent mass of 98 kDa was close to the predicted mass of 89.0 kDa.
Immunofluorescence assays of HA-tagged parasites were performed on asexual blood stages, and assays were also performed on exflagellating male gametocytes, comparing knockout and parental 820cl1m1cl1 parasites, as previously described [36, 37], with the following differences.
For blood-stage assays, samples were fixed in 4% paraformaldehyde and 0.00625% glutaraldehyde. For exflagellating gametocytes, 4% paraformaldehyde and microtubule stabilization buffer  were used. Parasites were permeabilied for 10 min in 0.1% or 0.5% Triton X-100. Blood-stage parasites were incubated in rat anti-HA (1/250; Roche, Australia, 18674) and Alexa Fluor Ⓡ goat anti-rat-488 (1/5000; Life Technologies, Australia, A-11006). Exflagellating gametocytes were incubated in chicken anti-GFP (1/1000; Abcam, USA, ab101863) and mouse anti-alpha-tubulin (1/300; Sigma-Aldrich, Australia, T9026), followed by Alexa Fluor Ⓡ goat anti-chicken-488 (1/2000; Life Technologies, Australia, A-11039) and Alexa Fluor Ⓡ goat anti-mouse-546 (1/1000; Life Technologies, Australia, A-11003). Both stages were incubated in 200 µg/ml Hoechst 33342 (Life Technologies, Australia).
Images were acquired with a Leica SP2 confocal microscope, linearly adjusted for brightness and contrast, and merged using the bundled Leica software.
Wild-type female, male, and asexual transcriptomes were previously published in . RNA-seq mapping and differential gene expression analyses were performed as previously described . Data verification and quality control were performed as previously described  and included generation of multidimensional scaling plots and heatmaps (Additional file 4: Figure S5–S6). Alternative splicing analyses were performed as previously described , using DEXSeq for changes in splicing between the two samples , and JunctionJuror for the presence of alternative splicing in single samples . Full commands and version numbers are presented in Additional file 13: Text 1.
Categorization of alternative splicing type (Fig. 3) was achieved using two tools. Intron retention was quantified with featureCounts  and alternative 5′/ 3′ splicing and exon skipping using Junction annotation, in the RSeQC package . For both approaches, we set a minimum threshold of three reads per replicate for the alternative junction to be accepted for further analysis. The intron read counts were normalized for intron length. We filtered for alternative splicing events where the alternative junction reads were more than 10% of the total splicing observed for each junction.
To analyze splicing changes between stages, we compared the fractions of alternative splicing reads for each junction, and obtained corrected p values using Benjamini-Hochberg false discovery correction, filtering at 0.25 FDR. We again filtered for junctions with a minimum coverage of three reads and changes that exceeded 10% of the total transcripts. We also filtered out differential alternative splicing events where the magnitude of the relative change was less than twofold in abundance. Because differential expression cannot be distinguished from differential splicing when the gene has negligible expression in one condition, we also excluded junctions with an average of fewer than 10 canonical junction per replicate for any condition.
Prior to reverse transcription, RNA samples were treated with DNase I (Qiagen, Australia). This comprized one asexual wild-type sample and one female wild-type sample. cDNA was synthesized using the SensiFASTTM cDNA Synthesis Kit (Bioline, Australia) according to the manufacturer’s instructions. Quantitative PCR was performed with the SensiFASTTM SYBR Ⓡ Hi-ROX Kit on an Applied Biosystems Quant Studio 6 Flex (Life Technologies, Australia). This was repeated three times for each sample. The ΔΔCt for the average of three experiments was calculated with ΔΔCt=ΔCt(asexual)−ΔCt(female), where ΔCt compared two amplicons from the same gene, with the numerator as the amplicon associated with the exonic region with changed expression (the “_T” primer), and the denominator representing the amplicon for an unchanged region (the “_C” primer). Statistical analysis was performed with Welch’s one-sample t test. Primer sequences are listed in Additional file 12: Table S11.
Blood stage growth assays were performed by competition of tdTomato parasites with parasites expressing GFP constitutively under the control of the HSP70 promoter (a kind gift from Robert Ménard; promoter described in ). Naive mice were infected with 5×104 parasites of either PbSR-MG KO or parental tdTomato strains, plus 5×104 parasites of the GFP strain, each combination with four replicates. Parasites were counted by FACS 3–6 days post-infection. For each mouse and each day, the parasitemia of tdTomato strains was normalized to that of the GFP strain. We presented this visually by dividing this normalized growth at 6 days by that at 3 days (Fig. 4). For statistical tests, normalized PbSR-MG KO tdTomato parasites were divided by normalized parental tdTomato parasites for each day. To see if this ratio changed as parasitemia increased, we assayed by simple linear regression, with no change in ratio over time detected (p value = 0.886).
Gametocytes ratios were determined by comparison of female and male counts of 820cl1m1cl1-based parasites, as partitioned by FACS sorting. See  for a representative FACS plot.
Activated females were quantified as previously described ; the numbers reported are the sum of activated females and ookinetes, since the latter are activated females that have been subsequently fertilized. This was performed in tdTomato-based parasites.
The exflagellation rates presented were determined from 820cl1m1cl1-based parasites. In vitro induction of exflagellation was performed as previously described , but using 2.5 µl blood in 50 µl exflagellation media. After 13 min, parasites were fixed as described above for immunofluorescence assays. Males were detected by the presence of the male-specific GFP marker, and exflagellating males were identified by the protrusion of tubulin (Fig. 5a). Exflagellation defects were confirmed in all three knockout lines.
Mosquitoes were infected as previously described . Given the low number of exflagellation events in PbSR-MG KO mutants, we were sometimes forced to infect mosquitoes with fewer than the normal minimum limit of exflagellating gametocytes .
Oocyst and sporozoite numbers were quantified by analysis of tdTomato parental and PbSR-MG KO parasites in paired experiments, as previously described , but without staining. Time to patency in mice reinfected with these parasites was determined as previously described . To control for lack of growth in upstream mosquito parasite stages, mice that exhibited no infection after several weeks were not included in the analysis.
Normalized oocyst growth data (Additional file 4: Figure S2) were obtained by in vitro fertilization of blood-stage parasites (using either 820cl1m1cl1 or tdTomato lines) as previously reported . Parasites were diluted to 460 or 785 ookinetes per microliter blood and fed to mosquitoes via an artificial membrane. Due to the low numbers of ookinetes acquired in the PbSR-MG KO strain, the exact number of normalized ookinetes varied between experiments. Hence, KO and parental strains were experimentally paired and subsequently analyzed with a paired t test.
Complementation was achieved by growing 820cl1m1cl1-based PbSR-MG KO parasites in mice treated with 1 mg/ml 5-fluorocytosine for 4 days, to select for parasites that had lost the DHFR-yFCU cassette due to recombination (Additional file 4: Figure S3A). Parasites without the hDHFR-yFCU cassette were cloned by limiting the dilution, and the recovery of a clone lacking hDHFR-yFCU (PbSR-MG-KO ΔDH) was confirmed by PCR (Additional file 4: Figure S3B). The PlasmoGEM construct containing the full-length PbSR-MG gene on artificial chromosome PbAC0281c08 was linearized with NotI, transfected into PbSR-MG-KO ΔDH parasites, and transfectants selected by treatment with pyrimethamine. The genotype of the complemented parasite line was confirmed by PCR (Additional file 4: Figure S3B). Exflagellation counts were obtained as previously described . Oocyst counts were obtained as above.
For all graphs in Figs. 4 and 5b and Additional file 4: Figure S3, boxplots and medians are based on log-transformed data. Whiskers are restricted to 1.5 times the interquartile range, in log-space. With the exception of the blood stage competition assay, which is described above, statistical significance was inferred using Welch’s two-sample t test on log-transformed data. Counts that were zero (e.g., exflagellation of PbSR-MG KO parasites) were converted to the limit of detection (i.e., 0.5 parasites per sample) for the purposes of these calculations. p values of < 0.05 were considered statistically significant for all tests. All statistical tests were two-sided.
Paronetto MP, Passacantilli I, Sette C. Alternative splicing and cell survival: from tissue homeostasis to disease. Cell Death Differ. 2016; 23:1919–29.
Baralle FE, Giudice J. Alternative splicing as a regulator of development and tissue identity. Nat Rev Mol Cell Biol. 2017; 18:437–51.
Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, Mayr C, Kingsmore SF, Schroth GP, Burge CB. Alternative isoform regulation in human tissue transcriptomes. Nature. 2008; 456(7221):470–6.
Pan Q, Shai O, Lee LJ, Frey BJ, Blencowe BJ. Deep surveying of alternative splicing complexity in the human transcriptome by high-throughput sequencing. Nat Genet. 2008; 40:1413–5.
Cavalier-Smith T. Kingdom protozoa and its 18 phyla. Microbiol Rev. 1993; 57:953–94.
Adl SM, Simpson AGB, Lane CE, Lukeš J, Bass D, Bowser SS, Brown MW, Burki F, Dunthorn M, Hampl V, Heiss A, Hoppenrath M, Lara E, Le Gall L, Lynn DH, McManus H, Mitchell EAD, Mozley-Stanridge SE, Parfrey LW, Pawlowski J, Rueckert S, Shadwick L, Shadwick L, Schoch CL, Smirnov A, Spiegel FW. The revised classification of eukaryotes. J Eukaryot Microbiol. 2012; 59:429–93.
Sorber K, Dimon MT, DeRisi JL. RNA-Seq analysis of splicing in Plasmodium falciparum uncovers new splice junctions, alternative splicing and splicing of antisense transcripts. Nucleic Acids Res. 2011; 39:3820–35.
Yeoh LM, Goodman CD, Hall NE, van Dooren GG, McFadden GI, Ralph SA. A serine-arginine–rich (SR) splicing factor modulates alternative splicing of over a thousand genes in Toxoplasma gondii. Nucleic Acids Res. 2015; 43:4661–75.
Le Roch KG, Johnson JR, Florens L, Zhou Y, Santrosyan A, Grainger M, Yan SF, Williamson KC, Holder AA, Carucci DJ, Yates JR, Winzeler EA. Global analysis of transcript and protein levels across the Plasmodium falciparum life cycle. Genome Res. 2004; 14:2308–18.
Foth BJ, Zhang N, Mok S, Preiser PR, Bozdech Z. Quantitative protein expression profiling reveals extensive post-transcriptional regulation and post-translational modifications in schizont-stage malaria parasites. Genome Biol. 2008; 9:177.
Foth BJ, Zhang N, Chaal BK, Sze SK, Preiser PR, Bozdech Z. Quantitative time-course profiling of parasite and host cell proteins in the human malaria parasite Plasmodium falciparum. Mol Cell Proteomics. 2011; 10:110–006411.
Janse CJ, Waters AP. Sexual development of parasites In: Waters AP, Janse CJ, editors. Malaria parasites: genomes and molecular biology. Norfolk: Caister Academic Press: 2004. p. 445–74.
Sinden RE, Alavi YIH, Butcher GA, Dessens JT, Raine JD, Trueman HE. Ookinete cell biology In: Waters AP, Janse CJ, editors. Malaria parasites: genomes and molecular biology. Norfolk: Caister Academic Press: 2004. p. 445–74.
Vanderberg J, Rdodin J, Yoeli M. Electron microscopic and histochemical studies of sporozoite formation in Plasmodium berghei. J Protozool. 1967; 14(1):82–103.
Shortt HE, Fairley NH, Covell G, Shute PG, Garnham PCC. The pre-erythrocytic stage of Plasmodium falciparum. Trans R Soc Trop Med Hyg. 1951; 44:405–19.
Yeoh LM, Goodman CD, Mollard V, McFadden GI, Ralph SA. Comparative transcriptomics of female and male gametocytes in Plasmodium berghei and the evolution of sex in alveolates. BioMed Cent Genom. 2017; 18:734.
Sinha A, Hughes KR, Modrzynska KK, Otto TD, Pfander C, Dickens NJ, Religa AA, Bushell E, Graham AL, Cameron R, Kafsack BFC, Williams AE, Llinás M, Berriman M, Billker O, Waters AP. A cascade of DNA-binding proteins for sexual commitment and development in Plasmodium. Nature. 2014; 507:253–7.
Bushell E, Gomes AR, Sanderson T, Anar B, Girling G, Herd C, Metcalf T, Modrzynska K, Schwach F, Martin RE, Mather MW, McFadden GI, Parts L, Rutledge GG, Vaidya AB, Wengelnik K, Rayner JC, Billker O. Functional profiling of a Plasmodium genome reveals an abundance of essential genes. Cell. 2017; 170:260–2728.
Anders S, Reyes A, Huber W. Detecting differential usage of exons from RNA-seq data. Genome Res. 2012; 22(10):2008–17.
Aurrecoechea C, Brestelli J, Brunk BP, Dommer J, Fischer S, Gajria B, Gao X, Gingle A, Grant G, Harb OS, Heiges M, Innamorato F, Iodice J, Kissinger JC, Kraemer E, Li W, Miller JA, Nayak V, Pennington C, Pinney DF, Roos DS, Ross C, Stoeckert Jr CJ, Treatman C, Wang H. PlasmoDB: a functional genomic database for malaria parasites. Nucleic Acids Res. 2009; 37(Database issue):539–43.
Miao J, Fan Q, Cui L, Li X, Wang H, Ning G, Reese JC, Cui L. The MYST family histone acetyltransferase regulates gene expression and cell cycle in malaria parasite Plasmodium falciparum. Mol Microbiol. 2010; 78:883–902.
Liu Y, Tewari R, Ning J, Blagborough AM, Garbom S, Pei J, Grishin NV, Steele RE, Sinden RE, Snell WJ, Billker O. The conserved plant sterility gene HAP2 functions after attachment of fusogenic membranes in Chlamydomonas and Plasmodium gametes. Genes Dev. 2008; 22:1051–68.
Schwach F, Bushell E, Gomes AR, Anar B, Girling G, Herd C, Rayner JC, Billker O. PlasmoGEM, a database supporting a community resource for large-scale experimental genetics in malaria parasites. Nucleic Acids Res. 2015; 43:1176–82.
Gomes AR, Bushell E, Schwach F, Girling G, Anar B, Quail MA, Herd C, Pfander C, Modrzynska K, Rayner JC, Billker O. A genome-scale vector resource enables high-throughput reverse genetic screening in a malaria parasite. Cell Host Microbe. 2015; 17:404–13.
Karni R, de Stanchina E, Lowe SW, Sinha R, Mu D, Krainer AR. The gene encoding the splicing factor SF2/ASF is a proto-oncogene. Nat Struct Mol Biol. 2007; 14(3):185–93.
Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014; 30:923–30.
Wang L, Wang S, Li W. RSeQC: quality control of RNA-seq experiments. Bioinformatics. 2012; 28:2184–5.
Kim E, Magen A, Ast G. Different levels of alternative splicing among eukaryotes. Nucleic Acids Res. 2007; 35:125–31.
Sieber P, Voigt K, Kämmer P, Brunke S, Schuster S, Linde J. Comparative study on alternative splicing in human fungal pathogens suggests its involvement during host invasion. Front Microbiol. 2018; 9:2313.
Tewari R, Straschil U, Bateman A, Böhme U, Cherevach I, Gong P, Pain A, Billker O. The systematic functional analysis of Plasmodium protein kinases identifies essential regulators of mosquito transmission. Cell Host Microbe. 2010; 8:377–87.
Mair GR, Lasonder E, Garver LS, Franke-Fayard BMD, Carret CK, Wiegant JCAG, Dirks RW, Dimopoulos G, Janse CJ, Waters AP. Universal features of post-transcriptional gene regulation are critical for Plasmodium zygote development. PLoS Pathog. 2010; 6:1000767.
Janse CJ, Ramesar J, Waters AP. High-efficiency transfection and drug selection of genetically transformed blood stages of the rodent malaria parasite Plasmodium berghei. Nat Protoc. 2006; 1:346–56.
Ponzi M, Sidén-Kiamos I, Bertuccini L, Currà C, Kroeze H, Camarda G, Pace T, Franke-Fayard B, Laurentino EC, Louis C, Waters AP, Janse CJ, Alano P. Egress of Plasmodium berghei gametes from their host erythrocyte is mediated by the MDV-1/PEG3 protein. Cell Microbiol. 2009; 11:1272–88.
Rathnapala UL, Goodman CD, McFadden GI. A novel genetic technique in Plasmodium berghei allows liver stage analysis of genes required for mosquito stage development and demonstrates that de novo heme synthesis is essential for liver stage development in the malaria parasite. PLoS Pathog. 2017; 13:1006396.
Lavazec C, Moreira CK, Mair GR, Waters AP, Janse CJ, Templeton TJ. Analysis of mutant Plasmodium berghei parasites lacking expression of multiple PbCCp genes. Mol Biochem Parasitol. 2009; 163:1–7.
Tonkin CJ, van Dooren GG, Spurck TP, Struck NS, Good RT, Handman E, Cowman AF, McFadden GI. Localization of organellar proteins in Plasmodium falciparum using a novel set of transfection vectors and a new immunofluorescence fixation method. Mol Biochem Parasitol. 2004; 137:13–21.
Hliscs M, Millet C, Dixon MW, Siden-Kiamos I, McMillan P, Tilley L. Organization and function of an actin cytoskeleton in Plasmodium falciparum gametocytes. Cell Microbiol. 2015; 17:207–25.
Hliscs M, Nahar C, Frischknecht F, Matuschewski K. Expression profiling of Plasmodium berghei HSP70 genes for generation of bright red fluorescent parasites. PLoS ONE. 2013; 8:72771.
Sturm A, Mollard V, Cozijnsen A, Goodman CD, McFadden GI. Mitochondrial ATP synthase is dispensable in blood-stage Plasmodium berghei rodent malaria but essential in the mosquito phase. Proc Natl Acad Sci U S A. 2015; 112:10216–23.
Shears MJ, MacRae JI, Mollard V, Goodman CD, Sturm A, Orchard LM, Llinás M, McConville MJ, Botté CY, McFadden GI. Characterization of the Plasmodium falciparum and P. berghei glycerol 3-phosphate acyltransferase involved in FASII fatty acid utilization in the malaria parasite apicoplast. Cell Microbiol. 2017; 19:e12633.
Goodman CD, Siregar JE, Mollard V, Vega-Rodriguez J, Syafruddin D, Matsuoka H, Matsuzaki M, Toyama T, Sturm A, Cozijnsen A, Jacobs-Lorena M, Kita K, Marzuki S, McFadden GI. Parasites resistant to the antimalarial atovaquone fail to transmit by mosquitoes. Science. 2016; 352:349–53.
Yeoh LM, Goodman CD, Mollard V, McHugh E, Lee V V, Sturm A, Cozijnsen A, McFadden GI, Ralph SA. Plasmodium berghei transcriptome for PbSR-MG knock-out parasites: female gametocytes, male gametocytes, and asexual erythrocytic stages. Datasets. Gene Expression Omnibus. 2017. https://www.ncbi.nlm.nih.gov/bioproject/391033. Accessed 16 April 2019.
Yeoh LM, Goodman CD, Mollard V, McHugh E, Lee V V, Sturm A, Cozijnsen A, McFadden GI, Ralph SA. Plasmodium berghei transcriptome for female gametocytes, male gametocytes, and asexual erythrocytic stages. Datasets. Gene Expression Omnibus. 2017. https://www.ncbi.nlm.nih.gov/bioproject/374918. Accessed 16 April 2019.
This research was supported by Melbourne Bioinformatics at the University of Melbourne, grant number [UOM0023]. The authors would like to thank Andy Waters (University of Glasgow) for the kind gift of the 820cl1m1cl1 parasite line and Robert Ménard (Institut Pasteur) for the kind gift of the constitutively expressing HSP70-GFP line. We wish to express our gratitude to Hayley Buchanan and Marion Hliscs (The University of Melbourne) for the technical assistance with parasite culturing and transfection. We thank Malcolm McConville (The University of Melbourne) for helpful discussions.
LMY was supported by an Australian Postgraduate Award. GIM gratefully acknowledges a program grant from the National Health and Medical Research Council (Australia), a Discovery Project from the Australian Research Council, and an Australian Laureate Fellowship from the Australian Research Council. SAR is supported by an Australian National Health and Medical Research Council grant (628704), an Australian Research Council Discovery Project grant (DP160100389), and an NHMRC RD Wright Biomedical fellowship (APP1062504). The funding bodies had no role in the design of the study, nor the collection, analysis, and interpretation of data, nor the writing of the manuscript.
Ethics approval and consent to participate
Experiments were performed in conformity with the Australian Prevention of Cruelty to Animals Act 1986 and the Prevention of Cruelty to Animals Regulations 2008, and reviewed and permitted by the Melbourne University Animal Ethics Committee (AEEC ethics ID: 1413078).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. A tab-delimited text file containing a list of genes where alternative splicing has changed in female gametocytes compared to asexual erythrocytic stages (in wild-type parasites). (TSV 527 kb)
Table S2. A tab-delimited text file containing a list of genes where alternative splicing has changed in male gametocytes compared to asexual erythrocytic stages (in wild-type parasites). (TSV 487 kb)
Table S3. A tab-delimited text file containing qRT-PCR results that verify changes in alternative splicing of ten genes identified by RNA-seq data in Additional file 1: Table S1. (TSV 2 kb)
Figure S1. A PDF image showing changes in alternative splicing between samples validated by qRT-PCR. Figure S2. A PDF image showing growth of PbSR-MG KO compared to parental parasites after normalization of ookinetes in all samples. Figure S3. A PDF image showing complementation of PbSR-MG KO parasites with full-length PbSR-MG. Figure S4. A PDF image showing PCR screens of monoclonal parasites for PbSR-MG KO integration into P. berghei ANKA 820cl1m1cl1 and P. berghei ANKA tdTomato; PCR analysis of polyclonal epitope-tagged PbSR-MG 3′-HA. Figure S5. A PDF image showing a multidimensional scaling (MDS) plot based on gene-level differential expression. Figure S6. A PDF image showing a heatmap of gene-level differential expression. (PDF 1172 kb)
Table S4. A tab-delimited text file containing a list of genes where alternative splicing has changed in PbSR-MG KO parasites compared to parental parasites, for asexual erythrocytic stages. (TSV 6 kb)
Table S5. A tab-delimited text file containing a list of genes where alternative splicing has changed in PbSR-MG KO parasites compared to parental parasites, for female gametocytes. (TSV 20 kb)
Table S6. A tab-delimited text file containing a list of genes where alternative splicing has changed in PbSR-MG KO parasites compared to parental parasites, for male gametocytes. (TSV 86 kb)
Table S7. A tab-delimited text file containing a list of genes with differential expression of full genes in PbSR-MG KO parasites compared to parental parasites, for all stages. (TSV 5 kb)
Tables S8. A zip-compressed file of tab-delimited text files containing six lists of genes that are constitutively alternatively spliced (excluding intron retention), for all samples. (ZIP 15 kb)
Tables S9. An Excel spreadsheet quantifying categories of alternative splicing in single stages. (XLSX 6008 kb)
Tables S10. An Excel spreadsheet quantifying categories of alternative splicing between stages. (XLSX 986 kb)
Table S11. A tab-delimited text file containing a list of primers used in this study. (TSV 3 kb)
Text 1. A text file containing detailed bioinformatics methods with full commands and version numbers. (TXT 4 kb)
About this article
Cite this article
Yeoh, L.M., Goodman, C.D., Mollard, V. et al. Alternative splicing is required for stage differentiation in malaria parasites. Genome Biol 20, 151 (2019) doi:10.1186/s13059-019-1756-6
- Plasmodium berghei
- Alternative splicing
- SR proteins
- Next-generation sequencing