Abstract
Malaria parasite genes exhibit variation in both sequence and expression level. There is much information on sequence polymorphism, but less resolution on natural variation in transcriptomes of parasites at specific developmental stages. This is largely because it is challenging to obtain highly replicated sampling of transcriptomes to overcome potentially confounding technical and biological variation. We address the issue in the major human parasite Plasmodium falciparum by obtaining RNA-seq profiles of multiple independent replicate preparations of mature schizont-stage parasites from a panel of clinical isolates recently established in culture and from long-term laboratory-adapted clones. With a goal of robustly identifying variably expressed genes, we show that increasing the numbers of biological sample replicates greatly improves the discovery rate. Generally, six independent replicates of each parasite culture is recommendable as being significantly to lower numbers, although for highly expressed genes variable expression can be detected when fewer replicates are available. A broad comparison identifies genes differing in relative expression between cultured clinical isolates and laboratory-adapted clones. Genes more highly expressed in the laboratory-adapted clones include an AP2 transcription factor gene Pf3D7_0420300 and putative methyl transferase genes. The variable expression of several known merozoite invasion ligands is confirmed, and previously uncharacterised genes are shown to be differentially expressed among clinical isolates. New RT-qPCR assays validate the variation in transcript levels of these genes, and allow quantitation of expression to be extended to a wider panel of clinical isolate samples. These variably expressed genes are new candidates for investigation as potential determinants of alternative parasite developmental pathways or targets of immunity.
Author summary Understanding parasite diversity and adaptation may require characterisation of gene expression variation, and is vital if chemotherapeutic or vaccine development is to consider new candidate targets, but it is technically challenging to generate precise data on clinical isolates. Here, we analyse the transcriptomes of mature Plasmodium falciparum schizonts using RNA-sequencing, using large numbers of biological replicate samples to minimise the impact of inter-replicate variation on observed patterns of differential expression. This identifies genes that are differentially expressed in long term laboratory-adapted parasites and recently cultured clinical isolates, as well as among different clinical isolates. In additional samples of schizonts grown in the first cycle ex vivo prior to any erythrocyte invasion, expression levels of a selected panel of these genes vary among isolates, but mean levels are similar to those in the continuously cultured clinical isolates, indicating that the latter are useful for experimental studies requiring biological replication.