Malaria Parasite Genes Research Articles

Suppression by RNA Polymerase I Inhibitors Varies Greatly Between Distinct RNA Polymerase I Transcribed Genes in Malaria Parasites.

The transcription of ribosomal RNA (rRNA) by RNA Polymerase I (Pol I) is the rate-limiting step in ribosome biogenesis and a major determinant of cellular growth rates. Unlike other eukaryotes, which express identical rRNA from large tandem arrays of dozens to hundreds of identical rRNA genes in every cell, the genome of the human malaria parasite Plasmodium falciparum contains only a handful single-copy 47S rRNA loci that differ substantially from one another in length, sequence, and expression in different cell types. We found that the growth of the malaria parasite was acutely sensitive to the Pol I inhibitors 9-hydroxyellipticine and BMH-21 and demonstrated that they greatly reduce the transcription of 47S rRNAs as well as the transcription of other non-coding RNA genes. This makes P. falciparum only the second known organism where RNA Polymerase I transcribes genes other than the 47S rRNAs. We found that the various types of Pol I-transcribed genes differed by more than two orders of magnitude in their susceptibility to these inhibitors and explored the implications of these findings for the regulation of rRNA in P. falciparum.

Suppression by RNA Polymerase I Inhibitors Varies Greatly Between Distinct RNA Polymerase I Transcribed Genes in Malaria Parasites.

Transcription of ribosomal RNA (rRNA) by RNA Polymerase I (Pol I) is the rate-limiting step in ribosome biogenesis and a major determinant of cellular growth rates. Unlike virtually every other eukaryote, which express identical rRNA from large tandem arrays of dozens to hundreds of identical rRNA genes in every cell, the genome of the human malaria parasite Plasmodium falciparum contains only a handful single-copy 47S rRNA loci that differ substantially from one another in length, sequence and expression in different cell-types. We found that growth of malaria parasite was acutely sensitive to the Pol I inhibitors 9-hydroxyellipticine and BMH-21 and demonstrate that they greatly reduce the transcription of 47S rRNAs as well as transcription of other non-coding RNA genes. Surprisingly, we found that the various types of Pol I-transcribed genes differed by more than two orders of magnitude in their susceptibility to these inhibitors and explore the implications of these findings for regulation of rRNA in P. falciparum.

Novel Ion Channel Genes in Malaria Parasites.

Ion channels serve many cellular functions including ion homeostasis, volume regulation, signaling, nutrient acquisition, and developmental progression. Although the complex life cycles of malaria parasites necessitate ion and solute flux across membranes, the whole-genome sequencing of the human pathogen Plasmodium falciparum revealed remarkably few orthologs of known ion channel genes. Contrasting with this, biochemical studies have implicated the channel-mediated flux of ions and nutritive solutes across several membranes in infected erythrocytes. Here, I review advances in the cellular and molecular biology of ion channels in malaria parasites. These studies have implicated novel parasite genes in the formation of at least two ion channels, with additional ion channels likely present in various membranes and parasite stages. Computational approaches that rely on homology to known channel genes from higher organisms will not be very helpful in identifying the molecular determinants of these activities. Given their unusual properties, novel molecular and structural features, and essential roles in pathogen survival and development, parasite channels should be promising targets for therapy development.

Consistent Community Detection in Inter-Layer Dependent Multi-Layer Networks

Community detection in multi-layer networks, which aims at finding groups of nodes with similar connective patterns among all layers, has attracted tremendous interests in multi-layer network analysis. Most existing methods are extended from those for single-layer networks, which assume that different layers are independent. In this article, we propose a novel community detection method in multi-layer networks with inter-layer dependence, which integrates the stochastic block model (SBM) and the Ising model. The community structure is modeled by the SBM model and the inter-layer dependence is incorporated via the Ising model. An efficient alternative updating algorithm is developed to tackle the resultant optimization task. Moreover, the asymptotic consistencies of the proposed method in terms of both parameter estimation and community detection are established, which are supported by extensive simulated examples and a real example on a multi-layer malaria parasite gene network. Supplementary materials for this article are available online.

Unusually Divergent Ubiquitin Genes and Proteins in Plasmodium Species.

Ubiquitin is an extraordinarily highly conserved 76 amino acid protein encoded by three different types of gene, where the primary translation products are fusions either of ubiquitin with one of two ribosomal proteins (RPs) or of multiple ubiquitin monomers from head to tail. Here, we investigate the evolution of ubiquitin genes in mammalian malaria parasites (Plasmodium species). The ubiquitin encoded by the RPS27a fusion gene is highly divergent, as previously found in a variety of protists. However, we also find that two other forms of divergent ubiquitin sequence, each previously thought to be extremely rare, have arisen recently during the divergence of Plasmodium subgenera. On two occasions, in two distinct lineages, the ubiquitin encoded by the RPL40 fusion gene has rapidly diverged. In addition, in one of these lineages, the polyubiquitin genes have undergone a single codon insertion, previously considered a unique feature of Rhizaria. There has been disagreement whether the multiple ubiquitin coding repeats within a genome exhibit concerted evolution or undergo a birth-and-death process; the Plasmodium ubiquitin genes show clear signs of concerted evolution, including the spread of this codon insertion to multiple repeats within the polyubiquitin gene.

Molecular Profiles of Multiple Antimalarial Drug Resistance Markers in Plasmodium falciparum and Plasmodium vivax in the Mandalay Region, Myanmar.

Emergence and spreading of antimalarial drug resistant malaria parasites are great hurdles to combating malaria. Although approaches to investigate antimalarial drug resistance status in Myanmar malaria parasites have been made, more expanded studies are necessary to understand the nationwide aspect of antimalarial drug resistance. In the present study, molecular epidemiological analysis for antimalarial drug resistance genes in Plasmodium falciparum and P. vivax from the Mandalay region of Myanmar was performed. Blood samples were collected from patients infected with P. falciparum and P. vivax in four townships around the Mandalay region, Myanmar in 2015. Partial regions flanking major mutations in 11 antimalarial drug resistance genes, including seven genes (pfdhfr, pfdhps, pfmdr-1, pfcrt, pfk13, pfubp-1, and pfcytb) of P. falciparum and four genes (pvdhfr, pvdhps, pvmdr-1, and pvk12) of P. vivax were amplified, sequenced, and overall mutation patterns in these genes were analyzed. Substantial levels of mutations conferring antimalarial drug resistance were detected in both P. falciparum and P. vivax isolated in Mandalay region of Myanmar. Mutations associated with sulfadoxine-pyrimethamine resistance were found in pfdhfr, pfdhps, pvdhfr, and pvdhps of Myanmar P. falciparum and P. vivax with very high frequencies up to 90%. High or moderate levels of mutations were detected in genes such as pfmdr-1, pfcrt, and pvmdr-1 associated with chloroquine resistance. Meanwhile, low frequency mutations or none were found in pfk13, pfubp-1, pfcytb, and pvk12 of the parasites. Overall molecular profiles for antimalarial drug resistance genes in malaria parasites in the Mandalay region suggest that parasite populations in the region have substantial levels of mutations conferring antimalarial drug resistance. Continuous monitoring of mutations linked with antimalarial drug resistance is necessary to provide useful information for policymakers to plan for proper antimalarial drug regimens to control and eliminate malaria in the country.

Community detection on mixture multilayer networks via regularized tensor decomposition

We study the problem of community detection in multilayer networks, where pairs of nodes can be related in multiple modalities. We introduce a general framework, that is, mixture multilayer stochastic block model (MMSBM), which includes many earlier models as special cases. We propose a tensor-based algorithm (TWIST) to reveal both global/local memberships of nodes, and memberships of layers. We show that the TWIST procedure can accurately detect the communities with small misclassification error as the number of nodes and/or number of layers increases. Numerical studies confirm our theoretical findings. To our best knowledge, this is the first systematic study on the mixture multilayer networks using tensor decomposition. The method is applied to two real datasets: worldwide trading networks and malaria parasite genes networks, yielding new and interesting findings.

Mutations in Pfcrt and Pfmdr1 genes of Plasmodium falciparum isolates from two sites in Northcentral and Southwest Nigeria.

The ability of malaria parasites to develop resistance to antimalarial drugs has made it necessary to continuously survey malaria parasite populations for resistance markers. Mutations in specific malaria parasite genes confer resistance to antimalarial drugs. The study compared mutations in Pfcrt and Pfmdr1 genes of P. falciparum from two ecologically different areas of Nigeria. Plasmodium falciparum dried blood spots collected from New Bussa (Northcentral Nigeria) and Ijede (Southwest Nigeria) were analysed by PCR-RFLP for Pfcrt, K76T, Pfmdr1, N86Y and Y184F mutations. Pfmdr1 copy number was determined by quantitative-PCR. A total of 145 blood spots [Ijede=55; New Bussa=90 blood spots] were analysed, but Pfcrt gene was successfully amplified in 144 samples while Pfmdr1 was amplified in 132 samples. Overall, prevalence of mutant forms of Pfcrt 76T,Pfmdr1 86Y and 184F were 74.3% (95% CI: 66.4-81.2%), 18.2% (95% CI: 12.0-25.8%) and 35.6% (95% CI: 27.5-44.4%). The frequency of Pfcrt 76T was similar in both study sites [Ijede: 81.8% (95%CI: 69.1-90.9%); New Bussa: 69.7% (95%CI: 59.0-79.0), p=0.105]. However, the frequencies of Pfmdr1 86Y and 184F were significantly higher in Ijede (28.3% and 62.3%) than in New Bussa (11.4% and 17.7%), respectively (P<0.05). Eight parasite genotypes based on three codons of the two genes were identified. The most frequent genotype was TNY 53(40.5%) while the least was KYF 1 (0.8%). The most frequent genotype in Ijede and New Bussa were TNF 18(34.0%) and TNY 40 (51.3%) respectively. The frequency of wild strain KNF in Ijede and New Bussa were 3 (5.7%) and 18 (23.1%), respectively. The distribution of the genotypes differed significantly by location. The genotypes with more than two or more mutations were more in Ijede 32 (60.4%) than in New Bussa 16 (20.5%) (p<0.001). Amplification of Pfmdr1 copy number was not observed in the two study sites. The prevalence of Pfcrt 76T was similar in both locations while Pfmdr1 86Y and 184F differed in both locations. Single nucleotide polymorphisms in the three codons assessed were more in Ijede than in New Bussa.

Differential diagnosis of Plasmodium falciparum and Plasmodium vivax in mixed infection by colorimetric nanogold probes

Malaria is an infectious disease reported mostly in the tropical region. The most severe human malaria is Plasmodium falciparum since it can cause cerebral malaria. Therefore, the presence of P. falciparum either in single or mixed infection needs accurate diagnosis. In some mixed infections, the presence of P. falciparum may be cryptic which cannot be detected by microscopic examination. The molecular diagnosis is required in these cases. Many methods based on amplification of malaria parasite genes have been developed but most of them need sophisticated instruments. Here, we created a colorimetric method using probe immobilized gold nanoparticles (AuNPs) to detect the malaria parasite gene. Color changes rely on salt-induced aggregation of AuNPs in the presence or absence of DNA hybridization. Color changes could be observed either by a naked eye or UV–vis spectrophotometer. By this approach, single infection by the most common malaria parasite, P. falciparum or P. vivax could be differentially identified. Mixed infection of these two malaria species could also be clearly diagnosed including cases of cryptic P. falciparum. The novel nanogold based molecular malaria diagnosis is sensitive, specific, rapid and cheap ($0.94). The prepared nanogold malaria probes are stable for up to 3 months indicating their filed application in remote areas.

Transcriptional Analysis of Tightly Synchronized Plasmodium falciparum Intraerythrocytic Stages by RT-qPCR.

In Plasmodium falciparum, the parasite responsible for the most severe forms of human malaria, many fundamental processes are controlled at the transcriptional level. Studies on diverse aspects of basic parasite biology as well as molecular epidemiology studies often rely on the ability to accurately measure transcript levels, but this is complicated by the cyclic expression patterns of the majority of malaria parasite genes. Here, we provide a complete workflow to measure transcript levels in P. falciparum intraerythrocytic blood stages, overcoming the confounding factors that are commonly encountered. The method described covers all the steps from synchronization of parasite cultures to reverse transcriptase quantitative PCR (RT-qPCR) analysis.

Effect of Dihydroartemisinin on Plasmodium NADH-Dependent Glutamate Synthase: The Implication in Malaria Management.

Artemisinin and its analogues (ARTs) are currently the most effective anti-malarial drugs, but the precise mechanism of action is still highly controversial. Effects of ARTs on Plasmodium genes expression are studied in our Lab. The overexpression of an interesting amidotransferase, NADH-dependent glutamate synthase (NADH-GltS) was found in treated by dihydroartemisinin (DHA). The increased expression occurred not only from global transcriptomics analysis on the human malaria parasite Plasmodium falciparum (P.falciparum) 3D7 and gene expression screening on all of iron-sulphur cluster proteins from . 3D7 invitro but also from Plasmodium berghei (P.berghei) ANKA in mice. Influence of DHA on NADH-GltS was specifically at trophozoite stage of P.falciparum and in a dose-dependent manner below the effective doses. L-glutamine (Gln) and L-glutamate (Glu) are the substrate and product of NADH-GltS respectively. Azaserine (Aza) is specific inhibitor for NADH-GltS. Experimental data showed that Glu levels were significantly decreasing with DHA dose increasing but NADH-GltS enzyme activities were still remained at higher levels in parasites, and appropriate amount of exogenous Glu could significantly reduce anti-malarial action of DHA but excessive amount lost the above effect. Aza alone could inhibit proliferation of P. falciparum and had an additive effect in combination with DHA. Those results could suggest that: Glutamate depletion is one of the anti-malarial actions of DHA; overexpression of NADH-GltS would be a feedback pattern of parasite itself due to glutamate depletion, but not a direct action of DHA; the "feedback pattern" is one of protective strategies of Plasmodium to interfere with the anti-malarial actions of DHA; and specific inhibitor for NADH-GltS as a new type of anti-malarial agents or new partner in ACT might provide a potential.

Multi-population genomic analysis of malaria parasites indicates local selection and differentiation at the gdv1 locus regulating sexual development

Parasites infect hosts in widely varying environments, encountering diverse challenges for adaptation. To identify malaria parasite genes under locally divergent selection across a large endemic region with a wide spectrum of transmission intensity, genome sequences were obtained from 284 clinical Plasmodium falciparum infections from four newly sampled locations in Senegal, The Gambia, Mali and Guinea. Combining these with previous data from seven other sites in West Africa enabled a multi-population analysis to identify discrete loci under varying local selection. A genome-wide scan showed the most exceptional geographical divergence to be at the early gametocyte gene locus gdv1 which is essential for parasite sexual development and transmission. We identified a major structural dimorphism with alternative 1.5 kb and 1.0 kb sequence deletions at different positions of the 3′-intergenic region, in tight linkage disequilibrium with the most highly differentiated single nucleotide polymorphism, one of the alleles being very frequent in Senegal and The Gambia but rare in the other locations. Long non-coding RNA transcripts were previously shown to include the entire antisense of the gdv1 coding sequence and the portion of the intergenic region with allelic deletions, suggesting adaptive regulation of parasite sexual development and transmission in response to local conditions.

De novo synthesis of thiamine (vitamin B1) is the ancestral state in Plasmodium parasites - evidence from avian haemosporidians.

Parasites often have reduced genomes as their own genes become redundant when utilizing their host as a source of metabolites, thus losing their own de novo production of metabolites. Primate malaria parasites can synthesize vitamin B1 (thiamine) de novo but rodent malaria and other genome-sequenced apicomplexans cannot, as the three essential genes responsible for this pathway are absent in their genomes. The unique presence of functional thiamine synthesis genes in primate malaria parasites and their sequence similarities to bacterial orthologues, have led to speculations that this pathway was horizontally acquired from bacteria. Here we show that the genes essential for the de novo synthesis of thiamine are found also in avian Plasmodium species. Importantly, they are also present in species phylogenetically basal to all mammalian and avian Plasmodium parasites, i.e. Haemoproteus. Furthermore, we found that these genes are expressed during the blood stage of the avian malaria infection, indicating that this metabolic pathway is actively transcribed. We conclude that the ability to synthesize thiamine is widespread among haemosporidians, with a recent loss in the rodent malaria species.

PhenoPlasm: a database of disruption phenotypes for malaria parasite genes.

Two decades after the first Plasmodium transfection, attempts have been made to disrupt more than 3,151 genes in malaria parasites, across five Plasmodium species. While results from rodent malaria transfections have been curated and systematised, empowering large-scale analysis, phenotypic data from human malaria parasite transfections currently exists as individual reports scattered across a the literature. To facilitate systematic analysis of published experimental genetic data across Plasmodium species, we have built PhenoPlasm ( http://www.phenoplasm.org), a database of phenotypes generated by transfection experiments in all Plasmodium parasites. The site provides a simple interface linking citation-backed Plasmodium reverse-genetic phenotypes to gene IDs. The database has been populated with phenotypic data on 367 P. falciparum genes, curated from 176 individual publications, as well as existing data on rodent Plasmodium species from RMgmDB and PlasmoGEM. This is the first time that all available data on P. falciparum transfection experiments has been brought together in a single place. These data are presented using ortholog mapping to allow a researcher interested in a gene in one species to see results across other Plasmodium species. The collaborative nature of the database enables any researcher to add new phenotypes as they are discovered. As an example of database utility, we use the currently available datasets to identify RAP (RNA-binding domain abundant in Apicomplexa)-domain containing proteins as crucial to parasite survival.

PhenoPlasm: a database of disruption phenotypes for malaria parasite genes

Two decades after the first Plasmodium transfection, attempts have been made to disrupt more than 3,151 genes in malaria parasites, across five Plasmodium species. While results from rodent malaria transfections have been curated and systematised, empowering large-scale analysis, phenotypic data from human malaria parasite transfections currently exists as individual reports scattered across a the literature. To facilitate systematic analysis of published experimental genetic data across Plasmodium species, we have built PhenoPlasm (http://www.phenoplasm.org), a database of phenotypes generated by transfection experiments in all Plasmodium parasites. The site provides a simple interface linking citation-backed Plasmodium reverse-genetic phenotypes to gene IDs. The database has been populated with phenotypic data on 367 P. falciparum genes, curated from 176 individual publications, as well as existing data on rodent Plasmodium species from RMgmDB and PlasmoGEM. This is the first time that all available data on P. falciparum transfection experiments has been brought together in a single place. These data are presented using ortholog mapping to allow a researcher interested in a gene in one species to see results across other Plasmodium species. The collaborative nature of the database enables any researcher to add new phenotypes as they are discovered. As an example of database utility, we use the currently available datasets to identify RAP (RNA-binding domain abundant in Apicomplexa)-domain containing proteins as crucial to parasite survival.

Recombination events among virulence genes in malaria parasites are associated with G-quadruplex-forming DNA motifs.

BackgroundMalaria parasites of the genus Plasmodium possess large hyper-variable families of antigen-encoding genes. These are often variantly-expressed and are major virulence factors for immune evasion and the maintenance of chronic infections. Recombination and diversification of these gene families occurs readily, and may be promoted by G-quadruplex (G4) DNA motifs within and close to the variant genes. G4s have been shown to cause replication fork stalling, DNA breakage and recombination in model systems, but these motifs remain largely unstudied in Plasmodium.ResultsWe examined the nature and distribution of putative G4-forming sequences in multiple Plasmodium genomes, finding that their co-distribution with variant gene families is conserved across different Plasmodium species that have different types of variant gene families. In P. falciparum, where a large set of recombination events that occurred over time in cultured parasites has been mapped, we found a strong spatial association between these recombination events and putative G4-forming sequences. Finally, we searched Plasmodium genomes for the three classes of helicase that can unwind G4s: Plasmodium spp. have no identifiable homologue of the highly efficient G4 helicase PIF1, but they do encode two putative RecQ helicases and one homologue of the RAD3-family helicase FANCJ.ConclusionsOur analyses, conducted at the whole-genome level in multiple species of Plasmodium, support the concept that G4s are likely to be involved in recombination and diversification of antigen-encoding gene families in this important protozoan pathogen.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-3183-3) contains supplementary material, which is available to authorized users.

Dual regulatory effects of non-coding GC-rich elements on the expression of virulence genes in malaria parasites

As the primary virulence factor of falciparum malaria, var genes harboring mutually exclusive expression pattern lead to antigenic variation and immune evasion of this pathogen in human host. Although various mechanisms contribute to silence of var genes, little is known of transcriptional activation pathways of a single var gene and maintenance of its active state with other silent var loci. Here, we report a monoallelic expression pattern of the non-coding GC-elements flanking chromosomal internal var genes, and transcript from the active one was required for activation of the var gene in the same array. Meanwhile, GFP reporter assays revealed a repressive effect on the adjacent gene induced by DNA motifs of the insulator-like GC-element, which was linked to heterochromatin subnuclear localization. Taken together, these data for the first time provide experimental evidence of the dual cis- and trans-acting regulatory functions of the GC-elements in both silence and activation of var genes, which would advance our understanding of the complex regulatory network of the virulence gene family in P. falciparum.

Antimalarial Drug Resistance: In the Past, Current Status and Future Perspectives

The aim of this study was to review the antimalarial drug treatment and its resistance during the course of therapy. Malaria affects the populations of tropical and subtropical areas world-wide, as well as an increasing number of travellers to these areas. Although malaria is found in over 100 countries, the major burden of disease is carried by the nations of Africa, where over 90% of all deaths from falciparum malaria are recorded and where the high levels of morbidity and transmission place considerable strains on public health services and economic infrastructure. In the absence of effective vaccines, management of the disease has depended largely upon chemotherapy and chemoprophylaxis. The number of available and effective antimalarial drugs is quickly dwindling due to different cause. This is mainly because a number of drug resistance-associated mutations in malaria parasite genes, such as crt, mdr1, dhfr/dhps and others, have led to widespread resistance to all known classes of antimalarial compounds. Unfortunately, malaria parasites have started to exhibit some level of resistance in Southeast Asia even to the most recently introduced class of drugs, artemisinins. Molecular evolutionary and population genetic approaches will greatly facilitate our understanding of the evolution and spread of parasite drug resistance and will contribute to developing strategies for better control of malaria. While there is much need, the antimalarial drug development pipeline remains woefully thin, with little chemical diversity and there is currently no alternative to the precious artemisinins. In general, four basic methods have been routinely used to study or measure antimalarial drug resistance: in vivo, in vitro, animal model studies and molecular characterization. This review endeavors to present the background facts and treatment and more on the current status of antimalarial drug resistance, what their causes, mechanisms, spread and management are and future perspective with antimalarial treatment. Furthermore, findings from this review will be useful to clinical researchers, health planners, policy makers and stakeholders and others concerned for malaria patients.

Plasmodium falciparum gene expression measured directly from tissue during human infection.

BackgroundDuring the latter half of the natural 48-h intraerythrocytic life cycle of human Plasmodium falciparum infection, parasites sequester deep in endothelium of tissues, away from the spleen and inaccessible to peripheral blood. These late-stage parasites may cause tissue damage and likely contribute to clinical disease, and a more complete understanding of their biology is needed. Because these life cycle stages are not easily sampled due to deep tissue sequestration, measuring in vivo gene expression of parasites in the trophozoite and schizont stages has been a challenge.MethodsWe developed a custom nCounter® gene expression platform and used this platform to measure malaria parasite gene expression profiles in vitro and in vivo. We also used imputation to generate global transcriptional profiles and assessed differential gene expression between parasites growing in vitro and those recovered from malaria-infected patient tissues collected at autopsy.ResultsWe demonstrate, for the first time, global transcriptional expression profiles from in vivo malaria parasites sequestered in human tissues. We found that parasite physiology can be correlated with in vitro data from an existing life cycle data set, and that parasites in sequestered tissues show an expected schizont-like transcriptional profile, which is conserved across tissues from the same patient. Imputation based on 60 landmark genes generated global transcriptional profiles that were highly correlated with genome-wide expression patterns from the same samples measured by microarray. Finally, differential expression revealed a limited set of in vivo upregulated transcripts, which may indicate unique parasite genes involved in human clinical infections.ConclusionsOur study highlights the utility of a custom nCounter® P. falciparum probe set, validation of imputation within Plasmodium species, and documentation of in vivo schizont-stage expression patterns from human tissues.Electronic supplementary materialThe online version of this article (doi:10.1186/s13073-014-0110-6) contains supplementary material, which is available to authorized users.

Transcriptome sequencing and analysis of Plasmodium gallinaceum reveals polymorphisms and selection on the apical membrane antigen-1.

BackgroundPlasmodium erythrocyte invasion genes play a key role in malaria parasite transmission, host-specificity and immuno-evasion. However, the evolution of the genes responsible remains understudied. Investigating these genes in avian malaria parasites, where diversity is particularly high, offers new insights into the processes that confer malaria pathogenesis. These parasites can pose a significant threat to birds and since birds play crucial ecological roles they serve as important models for disease dynamics. Comprehensive knowledge of the genetic factors involved in avian malaria parasite invasion is lacking and has been hampered by difficulties in obtaining nuclear data from avian malaria parasites. Thus the first Illumina-based de novo transcriptome sequencing and analysis of the chicken parasite Plasmodium gallinaceum was performed to assess the evolution of essential Plasmodium genes.MethodsWhite leghorn chickens were inoculated intravenously with erythrocytes containing P. gallinaceum. cDNA libraries were prepared from RNA extracts collected from infected chick blood and sequencing was run on the HiSeq2000 platform. Orthologues identified by transcriptome sequencing were characterized using phylogenetic, ab initio protein modelling and comparative and population-based methods.ResultsAnalysis of the transcriptome identified several orthologues required for intra-erythrocytic survival and erythrocyte invasion, including the rhoptry neck protein 2 (RON2) and the apical membrane antigen-1 (AMA-1). Ama-1 of avian malaria parasites exhibits high levels of genetic diversity and evolves under positive diversifying selection, ostensibly due to protective host immune responses.ConclusionErythrocyte invasion by Plasmodium parasites require AMA-1 and RON2 interactions. AMA-1 and RON2 of P. gallinaceum are evolutionarily and structurally conserved, suggesting that these proteins may play essential roles for avian malaria parasites to invade host erythrocytes. In addition, host-driven selection presumably results in the high levels of genetic variation found in ama-1 of avian Plasmodium species. These findings have implications for investigating avian malaria epidemiology and population dynamics. Moreover, this work highlights the P. gallinaceum transcriptome as an important public resource for investigating the diversity and evolution of essential Plasmodium genes.Electronic supplementary materialThe online version of this article (doi:10.1186/1475-2875-13-382) contains supplementary material, which is available to authorized users.