Sequence of Plasmodium falciparum chromosomes 1, 3-9 and 13

Plasmodium falciparum Alba6 exhibits DNase activity and participates in stress response

Alba domain proteins, owing to their functional plasticity, play a significant role in organisms. Here, we report an intrinsic DNase activity of PfAlba6 from Plasmodium falciparum, an etiological agent responsible for human malignant malaria. We identified that tyrosine28 plays a critical role in the Mg2+ driven 5'-3' DNase activity of PfAlba6. PfAlba6 cleaves both dsDNA as well as ssDNA. We also characterized PfAlba6-DNA interaction and observed concentration-dependent oligomerization in the presence of DNA, which is evident from size exclusion chromatography and single molecule AFM-imaging. PfAlba6 mRNA expression level is up-regulated several folds following heat stress and treatment with artemisinin, indicating a possible role in stress response. PfAlba6 has no human orthologs and is expressed in all intra-erythrocytic stages; thus, this protein can potentially be a new anti-malarial drug target.

Plasmodium falciparum importation does not sustain malaria transmission in a semi-arid region of Kenya

Human movement impacts the spread and transmission of infectious diseases. Recently, a large reservoir of Plasmodium falciparum malaria was identified in a semi-arid region of northwestern Kenya historically considered unsuitable for malaria transmission. Understanding the sources and patterns of transmission attributable to human movement would aid in designing and targeting interventions to decrease the unexpectedly high malaria burden in the region. Toward this goal, polymorphic parasite genes (ama1, csp) in residents and passengers traveling to Central Turkana were genotyped by amplicon deep sequencing. Genotyping and epidemiological data were combined to assess parasite importation. The contribution of travel to malaria transmission was estimated by modelling case reproductive numbers inclusive and exclusive of travelers. P. falciparum was detected in 6.7% (127/1891) of inbound passengers, including new haplotypes which were later detected in locally-transmitted infections. Case reproductive numbers approximated 1 and did not change when travelers were removed from transmission networks, suggesting that transmission is not fueled by travel to the region but locally endemic. Thus, malaria is not only prevalent in Central Turkana but also sustained by local transmission. As such, interrupting importation is unlikely to be an effective malaria control strategy on its own, but targeting interventions locally has the potential to drive down transmission.

Assessment of genetic diversity of Plasmodium falciparum circumsporozoite protein in Sudan: the RTS,S leading malaria vaccine candidate

Background

The currently used malaria vaccine, RTS,S, is designed based on the Plasmodium falciparum circumsporozoite protein (PfCSP). The pfcsp gene, besides having different polymorphic patterns, can vary between P. falciparum isolates due to geographical origin and host immune response. Such aspects are essential when considering the deployment of the RTS,S vaccine in a certain region. Therefore, this study assessed the genetic diversity of P. falciparum in Sudan based on the pfcsp gene by investigating the diversity at the N-terminal, central repeat, and the C-terminal regions.

Methods

A cross-sectional molecular study was conducted; P. falciparum isolates were collected from different health centres in Khartoum State between January and December 2019. During the study period, a total of 261 febrile patients were recruited. Malaria diagnosis was made by expert microscopists using Giemsa-stained thick and thin blood films. DNA samples were examined by the semi-nested polymerase chain reaction (PCR). Single clonal infection of the confirmed P. falciparum cases, were used to amplify the pfcsp gene. The amplified amplicons of pfcsp have been sequenced using the Sanger dideoxy method. The obtained sequences of pfcsp nucleotide diversity parameters including the numbers of haplotypes (Hap), haplotypes diversity (Hapd), the average number of nucleotide differences between two sequences (p), and the numbers of segregating sites (S) were obtained. The haplotype networks were constructed using the online tcsBU software. Natural selection theory was also tested on pfcsp using Fuand Li’s D, Fuand Li’s F statistics, and Tajima’s D test using DnaSP.

Results

In comparison with the different pfcsp reference strains, the Sudanese isolates showed high similarity with other African isolates. The results of the N-terminal region showed the presence of 2 different haplotypes with a Hapd of 0.425 ± 0.00727. The presence of the unique insertion of NNNGDNGREGKDEDKRDGNN was reported. The KLKQP motif was conserved in all the studied isolates. At the central repeat region, 11 haplotypes were seen with a Hapd of 0.779 ± 0.00097. The analysis of the genetic diversity in the C-terminal region showed the presence of 10 haplotypes with a Hapd of 0.457 ± 0.073. Several non-synonymous amino acids changes were also seen at the Th2R and the Th3R T-cell epitope regions including T317K, E317K, Q318E, K321N, I322K, T322K, R322K, K324Q, I327L, G352N, S354P, R355K, N356D, Q357E, and E361A.

Conclusions

In this study, the results indicated a high conservation at the pfcsp gene. This may further contribute in understanding the genetic polymorphisms of P. falciparum prior to the deployment of the RTS,S vaccine in Sudan.

Exposure to diverse Plasmodium falciparum genotypes shapes the risk of symptomatic malaria in incident and persistent infections: A longitudinal molecular epidemiologic study in Kenya

BACKGROUND

Repeated exposure to malaria infections could protect against symptomatic progression, as people develop adaptive immunity to infections acquired over time.

METHODS

We investigated how new, recurrent, and persistent Plasmodium falciparum infections were associated with the odds of developing symptomatic compared to asymptomatic malaria. Using a 14-month longitudinal cohort in Western Kenya, we used amplicon deep sequencing of two polymorphic genes (pfama1 and pfcsp) to assess overlap of parasite genotypes (represented by haplotypes) acquired within an individual's successive infections. We hypothesized infections with novel haplotypes would increase the odds of symptomatic malaria.

RESULTS

After excluding initial infections, we observed 534 asymptomatic and 88 symptomatic infections across 186 people. We detected 109 pfcsp haplotypes, and each infection was classified as harboring novel, recurrent or persistent haplotypes. Incident infections with only new haplotypes had higher odds of symptomatic malaria when compared to infections with only recurrent haplotypes [odds ratio (OR): 3.24, 95% confidence interval (CI): 1.20 to 8.78], but infections with both new and recurrent haplotypes [OR: 0.64, 95% CI: 0.15 to 2.65] did not. Assessing persistent infections, those with mixed (persistent with new or recurrent) haplotypes [OR: 0.77, 95% CI: 0.21 to 2.75] had no association with symptomatic malaria compared to infections with only persistent haplotypes. Results were similar for pfama1.

CONCLUSIONS

These results confirm that incident infections with only novel haplotypes were associated with increased odds of symptomatic malaria compared to infections with only recurrent haplotypes but this relationship was not seen when haplotypes persisted over time in consecutive infections.

Direct Nanopore Sequencing of mRNA Reveals Landscape of Transcript Isoforms in Apicomplexan Parasites

Alternative splicing is a widespread phenomenon in metazoans by which single genes are able to produce multiple isoforms of the gene product. However, this has been poorly characterized in apicomplexans, a major phylum of some of the most important global parasites. Efforts have been hampered by atypical transcriptomic features, such as the high AU content of Plasmodium RNA, but also the limitations of short-read sequencing in deciphering complex splicing events. In this study, we utilized the long read direct RNA sequencing platform developed by Oxford Nanopore Technologies to survey the alternative splicing landscape of Toxoplasma gondii and Plasmodium falciparum. We find that while native RNA sequencing has a reduced throughput, it allows us to obtain full-length or nearly full-length transcripts with comparable quantification to Illumina sequencing. By comparing these data with available gene models, we find widespread alternative splicing, particularly intron retention, in these parasites. Most of these transcripts contain premature stop codons, suggesting that in these parasites, alternative splicing represents a pathway to transcriptomic diversity, rather than expanding proteomic diversity. Moreover, alternative splicing rates are comparable between parasites, suggesting a shared splicing machinery, despite notable transcriptomic differences between the parasites. This study highlights a strategy in using long-read sequencing to understand splicing events at the whole-transcript level and has implications in the future interpretation of transcriptome sequencing studies.

IMPORTANCE We have used a novel nanopore sequencing technology to directly analyze parasite transcriptomes. The very long reads of this technology reveal the full-length genes of the parasites that cause malaria and toxoplasmosis. Gene transcripts must be processed in a process called splicing before they can be translated to protein. Our analysis reveals that these parasites very frequently only partially process their gene products, in a manner that departs dramatically from their human hosts.

The repertoire of serine rhomboid proteases of piroplasmids of importance to animal and human health

Babesia, Theileria and Cytauxzoon are tick-borne apicomplexan protozoans of the order Piroplasmida, notorious for the diseases they cause in livestock, pets and humans. Host cell invasion is their Achilles heel, allowing for the development of drug or vaccine-based therapies. In other apicomplexans, cleavage of the transmembrane domain of adhesins by the serine rhomboid proteinase ROM4 is required for successful completion of invasion. In this study, we record and classify the rhomboid repertoire encoded in the genomes of 10 piroplasmid species pertaining to the lineages Babesia sensu stricto (s.s., Clade VI), Theileria sensu stricto (Clade IV), Theileria equi (Clade IV), Cytauxzoon felis (Clade IIIb) and Babesia microti (Clade I), as defined by Schnittger et al. (2012). Fifty-six piroplasmid rhomboid-like proteins were assigned by phylogenetic analysis and bidirectional best hit to the ROM4, ROM6, ROM7 or ROM8 groups, and their crucial motifs for conformation and function were identified. Forty-four of these rhomboids had either been incorrectly classified or misannotated. Babesia s.s. encode five or three ROM4 proteinase paralogs, whereas the remaining piroplasmids encode two ROM4 paralogs. All piroplasmids encode a single ROM6, ROM7 and ROM8. Thus, an increased paralog number of ROM4 is the only feature distinguishing Babesia s.s. from other piroplasmid lineages. Piroplasmid ROM6 is related to the mammalian mitochondrial rhomboid and, accordingly, N-terminal mitochondrial targeting signal sequences was found in some cases. ROM6 is the only rhomboid encoded by piroplasmids that is ubiquitous in other organisms. ROM8 represents a pseudoproteinase that is highly conserved between studied piroplasmids, suggesting that it is important in regulatory functions. ROM4, ROM6, ROM7 and ROM8 are exclusively present in Aconoidasida, which comprises piroplasmids and Plasmodium, suggesting a relevant functional role in erythrocyte invasion. The correct classification and designation of piroplasmid rhomboids presented in this study facilitates an informed choice for future in-depth study of their functions.

Genotyping cognate Plasmodium falciparum in humans and mosquitoes to estimate onward transmission of asymptomatic infections

Malaria control may be enhanced by targeting reservoirs of Plasmodium falciparum transmission. One putative reservoir is asymptomatic malaria infections and the scale of their contribution to transmission in natural settings is not known. We assess the contribution of asymptomatic malaria to onward transmission using a 14-month longitudinal cohort of 239 participants in a high transmission site in Western Kenya. We identify P. falciparum in asymptomatically- and symptomatically-infected participants and naturally-fed mosquitoes from their households, genotype all parasites using deep sequencing of the parasite genes pfama1 and pfcsp, and use haplotypes to infer participant-to-mosquito transmission through a probabilistic model. In 1,242 infections (1,039 in people and 203 in mosquitoes), we observe 229 (pfcsp) and 348 (pfama1) unique parasite haplotypes. Using these to link human and mosquito infections, compared with symptomatic infections, asymptomatic infections more than double the odds of transmission to a mosquito among people with both infection types (Odds Ratio: 2.56; 95% Confidence Interval (CI): 1.36-4.81) and among all participants (OR 2.66; 95% CI: 2.05-3.47). Overall, 94.6% (95% CI: 93.1-95.8%) of mosquito infections likely resulted from asymptomatic infections. In high transmission areas, asymptomatic infections are the major contributor to mosquito infections and may be targeted as a component of transmission reduction.

N-Glycosylation in Piroplasmids: Diversity within Simplicity

N-glycosylation has remained mostly unexplored in Piroplasmida, an order of tick-transmitted pathogens of veterinary and medical relevance. Analysis of 11 piroplasmid genomes revealed three distinct scenarios regarding N-glycosylation: Babesia sensu stricto (s.s.) species add one or two N-acetylglucosamine (NAcGlc) molecules to proteins; Theileria equi and Cytauxzoon felis add (NAcGlc)2-mannose, while B. microti and Theileria s.s. synthesize dolichol-P-P-NAcGlc and dolichol-P-P-(NAcGlc)2 without subsequent transfer to proteins. All piroplasmids possess the gene complement needed for the synthesis of the N-glycosylation substrates, dolichol-P and sugar nucleotides. The oligosaccharyl transferase of Babesia species, T. equi and C. felis, is predicted to be composed of only two subunits, STT3 and Ost1. Occurrence of short N-glycans in B. bovis merozoites was experimentally demonstrated by fluorescence microscopy using a NAcGlc-specific lectin. In vitro growth of B. bovis was significantly impaired by tunicamycin, an inhibitor of N-glycosylation, indicating a relevant role for N-glycosylation in this pathogen. Finally, genes coding for N-glycosylation enzymes and substrate biosynthesis are transcribed in B. bovis blood and tick stages, suggesting that this pathway is biologically relevant throughout the parasite life cycle. Elucidation of the role/s exerted by N-glycans will increase our understanding of these successful parasites, for which improved control measures are needed.

DNA analysis reveals non-falciparum malaria in the Democratic Republic of the Congo

BACKGROUND

The World Health Organization (WHO) attributes the entirety of malaria infection and transmission in the Democratic Republic of the Congo (DRC) to Plasmodium falciparum, one of the several species of malaria known to infect humans. Recent studies have put forth some evidence that transmission of Plasmodium vivax may also be occurring in the DRC. As interventions and treatments differ between malaria species, it is crucial to maintain the most accurate understanding of malaria species diversity in each region.

METHODS

Blood samples were taken from aymptomatic children 0-5 years old living in three regions of the DRC in 2014. For this study, samples were taken from a larger pool of samples, collected as part of a population-based survey in three regions. Plasmodium infection was screened for using nested polymerase chain reaction (PCR) assays and species were confirmed by cloning and DNA sequencing.

RESULTS

Of 336 samples screened by PCR, 62.2% (n=209) initially tested positive for P. falciparum and 14.6% (n=49) initially tested positive for P. vivax. Sanger sequencing was performed on PCR-positive Plasmodium samples to confirm identity of Plasmodium species. Sequencing showed Plasmodium malariae in one blood sample and Plasmodium ovale in another sample. Plasmodium vivax was detected in 12/65 cases (18.5%). Overall, 14/65 sequenced cases (21.5%) were infected with a non-falciparum malaria. 330bp 18s P. vivax DNA sequences were obtained.

CONCLUSIONS

This study reveals Plasmodium vivax and other non-falciparum malaria across several regions of the DRC, and enforces the importance of further testing and more precise diagnostics when testing for and treating malaria in the DRC. Here, we find a higher proportion of cases of P. vivax malaria than found in previous studies. This is the most robust DNA sequencing of Plasmodium vivax in the DRC to date.

Characterization of Plasmodium falciparum NEDD8 and identification of cullins as its substrates

A variety of post-translational modifications of Plasmodium falciparum proteins, including phosphorylation and ubiquitination, are shown to have key regulatory roles during parasite development. NEDD8 is a ubiquitin-like modifier of cullin-RING E3 ubiquitin ligases, which regulates diverse cellular processes. Although neddylation is conserved in eukaryotes, it is yet to be characterized in Plasmodium and related apicomplexan parasites. We characterized P. falciparum NEDD8 (PfNEDD8) and identified cullins as its physiological substrates. PfNEDD8 is a 76 amino acid residue protein without the C-terminal tail, indicating that it can be readily conjugated. The wild type and mutant (Gly75Ala/Gly76Ala) PfNEDD8 were expressed in P. falciparum. Western blot of wild type PfNEDD8-expressing parasites indicated multiple high molecular weight conjugates, which were absent in the parasites expressing the mutant, indicating conjugation of NEDD8 through Gly76. Immunoprecipitation followed by mass spectrometry of wild type PfNEDD8-expressing parasites identified two putative cullins. Furthermore, we expressed PfNEDD8 in mutant S. cerevisiae strains that lacked endogenous NEDD8 (rub1Δ) or NEDD8 conjugating E2 enzyme (ubc12Δ). The PfNEDD8 immunoprecipitate also contained S. cerevisiae cullin cdc53, further substantiating cullins as physiological substrates of PfNEDD8. Our findings lay ground for investigation of specific roles and drug target potential of neddylation in malaria parasites.

Characteristic features of the SERA multigene family in the malaria parasite

Serine repeat antigen (SERA) is conserved among species of the genus Plasmodium. Sera genes form a multigene family and are generally tandemly clustered on a single chromosome. Although all Plasmodium species encode multiple sera genes, the number varies between species. Among species, the members share similar sequences and gene organization. SERA possess a central papain-like cysteine protease domain, however, in some members, the active site cysteine residue is substituted with a serine. Recent studies implicate this gene family in a number of aspects in parasite biology and induction of protective immune response. This review summarizes the current understanding on this important gene family in several Plasmodium species. The Plasmodium falciparum (Pf)-sera family, for example, consists of nine gene members. Unlike other multigene families in Plasmodium species, Pf-sera genes do not exhibit antigenic variation. Pf-sera5 nucleotide diversity is also low. Moreover, although Pf-sera5 is highly transcribed during the blood stage of malaria infection, and a large amount is released into the host blood following schizont rupture, in malaria endemic countries the sero-positive rates for Pf-SERA5 are low, likely due to Pf-SERA5 binding of host proteins to avoid immune recognition. As an antigen, the N-terminal 47 kDa domain of Pf-SERA5 is a promising vaccine candidate currently undergoing clinical trials. Pf-SERA5 and Pf-SERA6, as well as P. berghei (Pb)-SERA3, and Pb-SERA5, have been investigated for their roles in parasite egress. Two P. yoelii SERA, which have a serine residue at the protease active center, are implicated in parasite virulence. Overall, these studies provide insight that during the evolution of the Plasmodium parasite, the sera gene family members have increased by gene duplication, and acquired various functions that enable the parasite to survive and successfully maintain infection in the host.

Biosynthesis of heme O in intraerythrocytic stages of Plasmodium falciparum and potential inhibitors of this pathway

A number of antimalarial drugs interfere with the electron transport chain and heme-related reactions; however, the biosynthesis of heme derivatives in Plasmodium parasites has not been fully elucidated. Here, we characterized the steps that lead to the farnesylation of heme. After the identification of a gene encoding heme O synthase, we identified heme O synthesis in blood stage parasites through the incorporation of radioactive precursors. The presence of heme O synthesis in intraerythrocytic stages of Plasmodium falciparum was confirmed by mass spectrometry. Inabenfide and uniconazole-P appeared to interfere in heme synthesis, accordingly, parasite growth was also affected by the addition of these drugs. We conclude that heme O synthesis occurs in blood stage-P. falciparum and this pathway could be a potential target for antimalarial drugs.

Adaptation of Translational Machinery in Malaria Parasites to Accommodate Translation of Poly-Adenosine Stretches Throughout Its Life Cycle

Malaria is caused by unicellular apicomplexan parasites of the genus Plasmodium, which includes the major human parasite Plasmodium falciparum. The complex cycle of the malaria parasite in both mosquito and human hosts has been studied extensively. There is tight control of gene expression in each developmental stage, and at every level of gene synthesis: from RNA transcription, to its subsequent translation, and finally post-translational modifications of the resulting protein. Whole-genome sequencing of P. falciparum has laid the foundation for significant biological advances by revealing surprising genomic information. The P. falciparum genome is extremely AT-rich (∼80%), with a substantial portion of genes encoding intragenic polyadenosine (polyA) tracks being expressed throughout the entire parasite life cycle. In most eukaryotes, intragenic polyA runs act as negative regulators of gene expression. Recent studies have shown that translation of mRNAs containing 12 or more consecutive adenosines results in ribosomal stalling and frameshifting; activating mRNA surveillance mechanisms. In contrast, P. falciparum translational machinery can efficiently and accurately translate polyA tracks without activating mRNA surveillance pathways. This unique feature of P. falciparum raises interesting questions: (1) How is P. falciparum able to efficiently and correctly translate polyA track transcripts, and (2) What are the specifics of the translational machinery and mRNA surveillance mechanisms that separate P. falciparum from other organisms? In this review, we analyze possible evolutionary shifts in P. falciparum protein synthesis machinery that allow efficient translation of an AU rich-transcriptome. We focus on physiological and structural differences of P. falciparum stage specific ribosomes, ribosome-associated proteins, and changes in mRNA surveillance mechanisms throughout the complete parasite life cycle, with an emphasis on the mosquito and liver stages.

Crystal structures of Babesia microti lactate dehydrogenase BmLDH reveal a critical role for Arg99 in catalysis

The tick- and transfusion-transmitted human pathogen Babesia microti infects host erythrocytes to cause the pathologic symptoms associated with human babesiosis, an emerging disease with worldwide distribution and potentially fatal clinical outcome. Drugs currently recommended for the treatment of babesiosis are associated with a high failure rate and significant adverse events, highlighting the urgent need for more-effective and safer babesiosis therapies. Unlike other apicomplexan parasites, B. microti lacks a canonical lactate dehydrogenase (LDH) but instead expresses a unique enzyme, B. microti LDH (BmLDH), acquired through evolution by horizontal transfer from a mammalian host. Here, we report the crystal structures of BmLDH in apo state and ternary complex (enzyme-NADH-oxamate) solved at 2.79 and 1.89 Å. Analysis of these structures reveals that upon binding to the coenzyme and substrate, the active pocket of BmLDH undergoes a major conformational change from an opened and disordered to a closed and stabilized state. Biochemical assays using wild-type and mutant B. microti and human LDHs identified Arg99 as a critical residue for the catalytic activity of BmLDH but not its human counterpart. Interestingly, mutation of Arg99 to Ala had no impact on the overall structure and affinity of BmLDH to NADH but dramatically altered the closure of the enzyme's active pocket. Together, these structural and biochemical data highlight significant differences between B. microti and human LDH enzymes and suggest that BmLDH could be a suitable target for the development of selective antibabesial inhibitors.-Yu, L., Shen, Z., Liu, Q., Zhan, X., Luo, X., An, X., Sun, Y., Li, M., Wang, S., Nie, Z., Ao, Y., Zhao, Y., Peng, G., Ben Mamoun, C., He, L., Zhao, J. Crystal structures of Babesia microti lactate dehydrogenase BmLDH reveal a critical role for Arg99 in catalysis.

Whole genome sequencing and microsatellite analysis of the Plasmodium falciparum E5 NF54 strain show that the var, rifin and stevor gene families follow Mendelian inheritance

Background

Plasmodium falciparum exhibits a high degree of inter-isolate genetic diversity in its variant surface antigen (VSA) families: P. falciparum erythrocyte membrane protein 1, repetitive interspersed family (RIFIN) and subtelomeric variable open reading frame (STEVOR). The role of recombination for the generation of this diversity is a subject of ongoing research. Here the genome of E5, a sibling of the 3D7 genome strain is presented. Short and long read whole genome sequencing (WGS) techniques (Ilumina, Pacific Bioscience) and a set of 84 microsatellites (MS) were employed to characterize the 3D7 and non-3D7 parts of the E5 genome. This is the first time that VSA genes in sibling parasites were analysed with long read sequencing technology.

Results

Of the 5733 E5 genes only 278 genes, mostly var and rifin/stevor genes, had no orthologues in the 3D7 genome. WGS and MS analysis revealed that chromosomal crossovers occurred at a rate of 0–3 per chromosome. var, stevor and rifin genes were inherited within the respective non-3D7 or 3D7 chromosomal context. 54 of the 84 MS PCR fragments correctly identified the respective MS as 3D7- or non-3D7 and this correlated with var and rifin/stevor gene inheritance in the adjacent chromosomal regions. E5 had 61 var and 189 rifin/stevor genes. One large non-chromosomal recombination event resulted in a new var gene on chromosome 14. The remainder of the E5 3D7-type subtelomeric and central regions were identical to 3D7.

Conclusions

The data show that the rifin/stevor and var gene families represent the most diverse compartments of the P. falciparum genome but that the majority of var genes are inherited without alterations within their respective parental chromosomal context. Furthermore, MS genotyping with 54 MS can successfully distinguish between two sibling progeny of a natural P. falciparum cross and thus can be used to investigate identity by descent in field isolates.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2503-2) contains supplementary material, which is available to authorized users.

PhyloChromoMap, a Tool for Mapping Phylogenomic History along Chromosomes, Reveals the Dynamic Nature of Karyotype Evolution in Plasmodium falciparum

The genome of Plasmodium falciparum, the causative agent of malaria in Africa, has been extensively studied since it was first fully sequenced in 2002. However, many open questions remain, including understanding the chromosomal context of molecular evolutionary changes (e.g., relationship between chromosome map and phylogenetic conservation, patterns of gene duplication, and patterns of selection). Here, we present PhyloChromoMap, a method that generates a phylogenomic map of chromosomes from a custom-built bioinformatics pipeline. Using P. falciparum 3D7 as a model, we analyze 2,116 genes with homologs in up to 941 diverse eukaryotic, bacterial and archaeal lineages. We estimate the level of conservation along chromosomes based on conservation across clades, and identify “young” regions (i.e., those with recent or fast evolving genes) that are enriched in subtelomeric regions as compared with internal regions. We also demonstrate that patterns of molecular evolution for paralogous genes differ significantly depending on their location as younger paralogs tend to be found in subtelomeric regions whereas older paralogs are enriched in internal regions. Combining these observations with analyses of synteny, we demonstrate that subtelomeric regions are actively shuffled among chromosome ends, which is consistent with the hypothesis that these regions are prone to ectopic recombination. We also assess patterns of selection by comparing dN/dS ratios of gene family members in subtelomeric versus internal regions, and we include the important antigenic gene family var. These analyses illustrate the highly dynamic nature of the karyotype of P. falciparum, and provide a method for exploring genome dynamics in other lineages.

Evolutionary Dynamics of Cryptophyte Plastid Genomes

Cryptophytes are an ecologically important group of largely photosynthetic unicellular eukaryotes. This lineage is of great interest to evolutionary biologists because their plastids are of red algal secondary endosymbiotic origin and the host cell retains four different genomes (host nuclear, mitochondrial, plastid, and red algal nucleomorph). Here, we report a comparative analysis of plastid genomes from six representative cryptophyte genera. Four newly sequenced cryptophyte plastid genomes of Chroomonas mesostigmatica, Ch. placoidea, Cryptomonas curvata, and Storeatula sp. CCMP1868 share a number of features including synteny and gene content with the previously sequenced genomes of Cryptomonas paramecium, Rhodomonas salina, Teleaulax amphioxeia, and Guillardia theta. Our analysis of these plastid genomes reveals examples of gene loss and intron insertion. In particular, the chlB/chlL/chlN genes, which encode light-independent (dark active) protochlorophyllide oxidoreductase (LIPOR) proteins have undergone recent gene loss and pseudogenization in cryptophytes. Comparison of phylogenetic trees based on plastid and nuclear genome data sets show the introduction, via secondary endosymbiosis, of a red algal derived plastid in a lineage of chlorophyll-c containing algae. This event was followed by additional rounds of eukaryotic endosymbioses that spread the red lineage plastid to diverse groups such as haptophytes and stramenopiles.

Whole-Genome Restriction Mapping by "Subhaploid"-Based RAD Sequencing: An Efficient and Flexible Approach for Physical Mapping and Genome Scaffolding

Assembly of complex genomes using short reads remains a major challenge, which usually yields highly fragmented assemblies. Generation of ultradense linkage maps is promising for anchoring such assemblies, but traditional linkage mapping methods are hindered by the infrequency and unevenness of meiotic recombination that limit attainable map resolution. Here we develop a sequencing-based "in vitro" linkage mapping approach (called RadMap), where chromosome breakage and segregation are realized by generating hundreds of "subhaploid" fosmid/bacterial-artificial-chromosome clone pools, and by restriction site-associated DNA sequencing of these clone pools to produce an ultradense whole-genome restriction map to facilitate genome scaffolding. A bootstrap-based minimum spanning tree algorithm is developed for grouping and ordering of genome-wide markers and is implemented in a user-friendly, integrated software package (AMMO). We perform extensive analyses to validate the power and accuracy of our approach in the model plant Arabidopsis thaliana and human. We also demonstrate the utility of RadMap for enhancing the contiguity of a variety of whole-genome shotgun assemblies generated using either short Illumina reads (300 bp) or long PacBio reads (6-14 kb), with up to 15-fold improvement of N50 (∼816 kb-3.7 Mb) and high scaffolding accuracy (98.1-98.5%). RadMap outperforms BioNano and Hi-C when input assembly is highly fragmented (contig N50 = 54 kb). RadMap can capture wide-range contiguity information and provide an efficient and flexible tool for high-resolution physical mapping and scaffolding of highly fragmented assemblies.

Trans-acting GC-rich non-coding RNA at var expression site modulates gene counting in malaria parasite

Monoallelic expression of the var multigene family enables immune evasion of the malaria parasite Plasmodium falciparum in its human host. At a given time only a single member of the 60-member var gene family is expressed at a discrete perinuclear region called the 'var expression site'. However, the mechanism of var gene counting remains ill-defined. We hypothesize that activation factors associating specifically with the expression site play a key role in this process. Here, we investigate the role of a GC-rich non-coding RNA (ncRNA) gene family composed of 15 highly homologous members. GC-rich genes are positioned adjacent to var genes in chromosome-central gene clusters but are absent near subtelomeric var genes. Fluorescence in situ hybridization demonstrates that GC-rich ncRNA localizes to the perinuclear expression site of central and subtelomeric var genes in trans. Importantly, overexpression of distinct GC-rich ncRNA members disrupts the gene counting process at the single cell level and results in activation of a specific subset of var genes in distinct clones. We identify the first trans-acting factor targeted to the elusive perinuclear var expression site and open up new avenues to investigate ncRNA function in antigenic variation of malaria and other protozoan pathogens.

Exploration of 3-methylisoquinoline-4-carbonitriles as protein kinase A inhibitors of Plasmodium falciparum

A series of isoquinolines have been evaluated in a homology model of Plasmodium falciparum Protein Kinase A (PfPKA) using molecular dynamics. Synthesis of these compounds was then undertaken to investigate their structure-activity relationships. One compound was found to inhibit parasite growth in an in vitro assay and provides a lead to further develop 3-methylisoquinoline-4-carbonitriles as antimalarial compounds. Development of a potent and selective PfPKA inhibitor would provide a useful tool to shed further insight into the mechanisms enabling malaria parasites to establish infection.

Integrated analysis of the Plasmodium species transcriptome

The genome sequence available for different Plasmodium species is a valuable resource for understanding malaria parasite biology. However, comparative genomics on its own cannot fully explain all the species-specific differences which suggests that other genomic aspects such as regulation of gene expression play an important role in defining species-specific characteristics. Here, we developed a comprehensive approach to measure transcriptional changes of the evolutionary conserved syntenic orthologs during the intraerythrocytic developmental cycle across six Plasmodium species. We show significant transcriptional constraint at the mid-developmental stage of Plasmodium species while the earliest stages of parasite development display the greatest transcriptional variation associated with critical functional processes. Modeling of the evolutionary relationship based on changes in transcriptional profile reveal a phylogeny pattern of the Plasmodium species that strictly follows its mammalian hosts. In addition, the work shows that transcriptional conserved orthologs represent potential future targets for anti-malaria intervention as they would be expected to carry out key essential functions within the parasites. This work provides an integrated analysis of orthologous transcriptome, which aims to provide insights into the Plasmodium evolution thereby establishing a framework to explore complex pathways and drug discovery in Plasmodium species with broad host range.

•

Comparison of variations in mRNA abundance across six different Plasmodium species.

•

Transcriptional conservation and divergence of Plasmodium syntenic orthologs.

•

Pattern of Plasmodium transcriptome evolution are established.

•

Transcriptionally conserved orthologs represent attractive intervention targets.

Malaria remains a major public health concern despite global efforts in the fight against this disease. The intraerythrocytic stage of the malaria parasites is currently in the spotlight for anti-malarial intervention and vaccine targets. The primary goal of this study is to generate a comprehensive and directly comparable transcriptome dataset across multiple Plasmodium species originating from different hosts. We establish that specific pathways and intraerythrocytic stages are more transcriptionally diverged than others, reflecting transcriptional evolutionary diversity. We further propose a panel of transcriptionally conserved genes as potential drug targets.

Thyroglobulin From Molecular and Cellular Biology to Clinical Endocrinology

Thyroglobulin (Tg) is a vertebrate secretory protein synthesized in the thyrocyte endoplasmic reticulum (ER), where it acquires N-linked glycosylation and conformational maturation (including formation of many disulfide bonds), leading to homodimerization. Its primary functions include iodide storage and thyroid hormonogenesis. Tg consists largely of repeating domains, and many tyrosyl residues in these domains become iodinated to form monoiodo- and diiodotyrosine, whereas only a small portion of Tg structure is dedicated to hormone formation. Interestingly, evolutionary ancestors, dependent upon thyroid hormone for development, synthesize thyroid hormones without the complete Tg protein architecture. Nevertheless, in all vertebrates, Tg follows a strict pattern of region I, II-III, and the cholinesterase-like (ChEL) domain. In vertebrates, Tg first undergoes intracellular transport through the secretory pathway, which requires the assistance of thyrocyte ER chaperones and oxidoreductases, as well as coordination of distinct regions of Tg, to achieve a native conformation. Curiously, regions II-III and ChEL behave as fully independent folding units that could function as successful secretory proteins by themselves. However, the large Tg region I (bearing the primary T4-forming site) is incompetent by itself for intracellular transport, requiring the downstream regions II-III and ChEL to complete its folding. A combination of nonsense mutations, frameshift mutations, splice site mutations, and missense mutations in Tg occurs spontaneously to cause congenital hypothyroidism and thyroidal ER stress. These Tg mutants are unable to achieve a native conformation within the ER, interfering with the efficiency of Tg maturation and export to the thyroid follicle lumen for iodide storage and hormonogenesis.

Enzymatic Characterization of Recombinant Food Vacuole Plasmepsin 4 from the Rodent Malaria Parasite Plasmodium berghei

The rodent malaria parasite Plasmodium berghei is a practical model organism for experimental studies of human malaria. Plasmepsins are a class of aspartic proteinase isoforms that exert multiple pathological effects in malaria parasites. Plasmepsins residing in the food vacuole (FV) of the parasite hydrolyze hemoglobin in red blood cells. In this study, we cloned PbPM4, the FV plasmepsin gene of P. berghei that encoded an N-terminally truncated pro-segment and the mature enzyme from genomic DNA. We over-expressed this PbPM4 zymogen as inclusion bodies (IB) in Escherichia coli, and purified the protein following in vitro IB refolding. Auto-maturation of the PbPM4 zymogen to mature enzyme was carried out at pH 4.5, 5.0, and 5.5. Interestingly, we found that the PbPM4 zymogen exhibited catalytic activity regardless of the presence of the pro-segment. We determined the optimal catalytic conditions for PbPM4 and studied enzyme kinetics on substrates and inhibitors of aspartic proteinases. Using combinatorial chemistry-based peptide libraries, we studied the active site preferences of PbPM4 at subsites S1, S2, S3, S1’, S2’ and S3’. Based on these results, we designed and synthesized a selective peptidomimetic compound and tested its inhibition of PbPM4, seven FV plasmepsins from human malaria parasites, and human cathepsin D (hcatD). We showed that this compound exhibited a >10-fold selectivity to PbPM4 and human malaria parasite plasmepsin 4 orthologs versus hcatD. Data from this study furthesr our understanding of enzymatic characteristics of the plasmepsin family and provides leads for anti-malarial drug design.

Plasmodium Apicoplast Gln-tRNAGln Biosynthesis Utilizes a Unique GatAB Amidotransferase Essential for Erythrocytic Stage Parasites

The malaria parasite Plasmodium falciparum apicoplast indirect aminoacylation pathway utilizes a non-discriminating glutamyl-tRNA synthetase to synthesize Glu-tRNA(Gln) and a glutaminyl-tRNA amidotransferase to convert Glu-tRNA(Gln) to Gln-tRNA(Gln). Here, we show that Plasmodium falciparum and other apicomplexans possess a unique heterodimeric glutamyl-tRNA amidotransferase consisting of GatA and GatB subunits (GatAB). We localized the P. falciparum GatA and GatB subunits to the apicoplast in blood stage parasites and demonstrated that recombinant GatAB converts Glu-tRNA(Gln) to Gln-tRNA(Gln) in vitro. We demonstrate that the apicoplast GatAB-catalyzed reaction is essential to the parasite blood stages because we could not delete the Plasmodium berghei gene encoding GatA in blood stage parasites in vivo. A phylogenetic analysis placed the split between Plasmodium GatB, archaeal GatE, and bacterial GatB prior to the phylogenetic divide between bacteria and archaea. Moreover, Plasmodium GatA also appears to have emerged prior to the bacterial-archaeal phylogenetic divide. Thus, although GatAB is found in Plasmodium, it emerged prior to the phylogenetic separation of archaea and bacteria.

Identification of Immunodominant B-cell Epitope Regions of Reticulocyte Binding Proteins in Plasmodium vivax by Protein Microarray Based Immunoscreening

Plasmodium falciparum can invade all stages of red blood cells, while Plasmodium vivax can invade only reticulocytes. Although many P. vivax proteins have been discovered, their functions are largely unknown. Among them, P. vivax reticulocyte binding proteins (PvRBP1 and PvRBP2) recognize and bind to reticulocytes. Both proteins possess a C-terminal hydrophobic transmembrane domain, which drives adhesion to reticulocytes. PvRBP1 and PvRBP2 are large (> 326 kDa), which hinders identification of the functional domains. In this study, the complete genome information of the P. vivax RBP family was thoroughly analyzed using a prediction server with bioinformatics data to predict B-cell epitope domains. Eleven pvrbp family genes that included 2 pseudogenes and 9 full or partial length genes were selected and used to express recombinant proteins in a wheat germ cell-free system. The expressed proteins were used to evaluate the humoral immune response with vivax malaria patients and healthy individual serum samples by protein microarray. The recombinant fragments of 9 PvRBP proteins were successfully expressed; the soluble proteins ranged in molecular weight from 16 to 34 kDa. Evaluation of the humoral immune response to each recombinant PvRBP protein indicated a high antigenicity, with 38-88% sensitivity and 100% specificity. Of them, N-terminal parts of PvRBP2c (PVX_090325-1) and PvRBP2 like partial A (PVX_090330-1) elicited high antigenicity. In addition, the PvRBP2-like homologue B (PVX_116930) fragment was newly identified as high antigenicity and may be exploited as a potential antigenic candidate among the PvRBP family. The functional activity of the PvRBP family on merozoite invasion remains unknown.

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.

Representing virus-host interactions and other multi-organism processes in the Gene Ontology

BACKGROUND

The Gene Ontology project is a collaborative effort to provide descriptions of gene products in a consistent and computable language, and in a species-independent manner. The Gene Ontology is designed to be applicable to all organisms but up to now has been largely under-utilized for prokaryotes and viruses, in part because of a lack of appropriate ontology terms.

METHODS

To address this issue, we have developed a set of Gene Ontology classes that are applicable to microbes and their hosts, improving both coverage and quality in this area of the Gene Ontology. Describing microbial and viral gene products brings with it the additional challenge of capturing both the host and the microbe. Recognising this, we have worked closely with annotation groups to test and optimize the GO classes, and we describe here a set of annotation guidelines that allow the controlled description of two interacting organisms.

CONCLUSIONS

Building on the microbial resources already in existence such as ViralZone, UniProtKB keywords and MeGO, this project provides an integrated ontology to describe interactions between microbial species and their hosts, with mappings to the external resources above. Housing this information within the freely-accessible Gene Ontology project allows the classes and annotation structure to be utilized by a large community of biologists and users.

The Plastid Genome of the Cryptomonad Teleaulax amphioxeia

Teleaulax amphioxeia is a photosynthetic unicellular cryptophyte alga that is distributed throughout marine habitats worldwide. This alga is an important plastid donor to the dinoflagellate Dinophysis caudata through the ciliate Mesodinium rubrum in the marine food web. To better understand the genomic characteristics of T. amphioxeia, we have sequenced and analyzed its plastid genome. The plastid genome sequence of T. amphioxeia is similar to that of Rhodomonas salina, and they share significant synteny. This sequence exhibits less similarity to that of Guillardia theta, the representative plastid genome of photosynthetic cryptophytes. The gene content and order of the three photosynthetic cryptomonad plastid genomes studied is highly conserved. The plastid genome of T. amphioxeia is composed of 129,772 bp and includes 143 protein-coding genes, 2 rRNA operons and 30 tRNA sequences. The DNA polymerase III gene (dnaX) was most likely acquired via lateral gene transfer (LGT) from a firmicute bacterium, identical to what occurred in R. salina. On the other hand, the psbN gene was independently encoded by the plastid genome without a reverse transcriptase gene as an intron. To clarify the phylogenetic relationships of the algae with red-algal derived plastids, phylogenetic analyses of 32 taxa were performed, including three previously sequenced cryptophyte plastid genomes containing 93 protein-coding genes. The stramenopiles were found to have branched out from the Chromista taxa (cryptophytes, haptophytes, and stramenopiles), while the cryptophytes and haptophytes were consistently grouped into sister relationships with high resolution.

Polymorphisms in the K13-propeller gene in artemisinin-susceptible Plasmodium falciparum parasites from Bougoula-Hameau and Bandiagara, Mali

Artemisinin-resistant Plasmodium falciparum malaria has been documented in southeast Asia and may already be spreading in that region. Molecular markers are important tools for monitoring the spread of antimalarial drug resistance. Recently, single-nucleotide polymorphisms (SNPs) in the PF3D7_1343700 kelch propeller (K13-propeller) domain were shown to be associated with artemisinin resistance in vivo and in vitro. The prevalence and role of K13-propeller mutations are poorly known in sub-Saharan Africa. K13-propeller mutations were genotyped by direct sequencing of nested polymerase chain reaction (PCR) amplicons from dried blood spots of pre-treatment falciparum malaria infections collected before and after the use of artemisinin-based combination therapy (ACT) as first-line therapy in Mali. Although K13-propeller mutations previously associated with delayed parasite clearance in Cambodia were not identified, 26 K13-propeller mutations were identified in both recent samples and pre-ACT infections. Parasite clearance time was comparable between infections with non-synonymous K13-propeller mutations and infections with the reference allele. These findings suggest that K13-propeller mutations are present in artemisinin-sensitive parasites and that they preceded the wide use of ACTs in Mali.

Novel approaches to identify protective malaria vaccine candidates

Efforts to develop vaccines against malaria have been the focus of substantial research activities for decades. Several categories of candidate vaccines are currently being developed for protection against malaria, based on antigens corresponding to the pre-erythrocytic, blood stage, or sexual stages of the parasite. Long lasting sterile protection from Plasmodium falciparum sporozoite challenge has been observed in human following vaccination with whole parasite formulations, clearly demonstrating that a protective immune response targeting predominantly the pre-erythrocytic stages can develop against malaria. However, most of vaccine candidates currently being investigated, which are mostly subunits vaccines, have not been able to induce substantial (>50%) protection thus far. This is due to the fact that the antigens responsible for protection against the different parasite stages are still yet to be known and relevant correlates of protection have remained elusive. For a vaccine to be developed in a timely manner, novel approaches are required. In this article, we review the novel approaches that have been developed to identify the antigens for the development of an effective malaria vaccine.

Noncoding RNAs as emerging regulators of Plasmodium falciparum virulence gene expression

The eukaryotic unicellular pathogen Plasmodium falciparum tightly regulates gene expression, both during development and in adaptation to dynamic host environments. This regulation is evident in the mutually exclusive expression of members of clonally variant virulence multigene families. While epigenetic regulators have been selectively identified at active or repressed virulence genes, their specific recruitment remains a mystery. In recent years, noncoding RNAs (ncRNAs) have emerged as lynchpins of eukaryotic gene regulation; by binding to epigenetic regulators, they provide target specificity to otherwise non-specific enzyme complexes. Not surprisingly, there is great interest in understanding the role of ncRNA in P. falciparum, in particular, their contribution to the mutually exclusive expression of virulence genes. The current repertoire of P. falciparum ncRNAs includes, but is not limited to, subtelomeric ncRNAs, virulence gene-associated ncRNAs and natural antisense RNA transcripts. Continued improvement in high-throughput sequencing methods is sure to expand this repertoire. Here, we summarize recent advances in P. falciparum ncRNA biology, with an emphasis on ncRNA-mediated epigenetic modes of gene regulation.

Identifying novel cell cycle proteins in Apicomplexa parasites through co-expression decision analysis

Hypothetical proteins comprise roughly half of the predicted gene complement of Toxoplasma gondii and Plasmodium falciparum and represent the largest class of uniquely functioning proteins in these parasites. Following the idea that functional relationships can be informed by the timing of gene expression, we devised a strategy to identify the core set of apicomplexan cell division cycling genes with important roles in parasite division, which includes many uncharacterized proteins. We assembled an expanded list of orthologs from the T. gondii and P. falciparum genome sequences (2781 putative orthologs), compared their mRNA profiles during synchronous replication, and sorted the resulting set of dual cell cycle regulated orthologs (744 total) into protein pairs conserved across many eukaryotic families versus those unique to the Apicomplexa. The analysis identified more than 100 ortholog gene pairs with unknown function in T. gondii and P. falciparum that displayed co-conserved mRNA abundance, dynamics of cyclical expression and similar peak timing that spanned the complete division cycle in each parasite. The unknown cyclical mRNAs encoded a diverse set of proteins with a wide range of mass and showed a remarkable conservation in the internal organization of ordered versus disordered structural domains. A representative sample of cyclical unknown genes (16 total) was epitope tagged in T. gondii tachyzoites yielding the discovery of new protein constituents of the parasite inner membrane complex, key mitotic structures and invasion organelles. These results demonstrate the utility of using gene expression timing and dynamic profile to identify proteins with unique roles in Apicomplexa biology.

Diversity and population structure of Plasmodium falciparum in Thailand based on the spatial and temporal haplotype patterns of the C-terminal 19-kDa domain of merozoite surface protein-1

Background

The 19-kDa C-terminal region of the merozoite surface protein-1 of the human malaria parasite Plasmodium falciparum (PfMSP-1₁₉) constitutes the major component on the surface of merozoites and is considered as one of the leading candidates for asexual blood stage vaccines. Because the protein exhibits a level of sequence variation that may compromise the effectiveness of a vaccine, the global sequence diversity of PfMSP-1₁₉ has been subjected to extensive research, especially in malaria endemic areas. In Thailand, PfMSP-1₁₉ sequences have been derived from a single parasite population in Tak province, located along the Thailand-Myanmar border, since 1995. However, the extent of sequence variation and the spatiotemporal patterns of the MSP-1₁₉ haplotypes along the Thai borders with Laos and Cambodia are unknown.

Methods

Sixty-three isolates of P. falciparum from five geographically isolated populations along the Thai borders with Myanmar, Laos and Cambodia in three transmission seasons between 2002 and 2008 were collected and culture-adapted. The msp-1 gene block 17 was sequenced and analysed for the allelic diversity, frequency and distribution patterns of PfMSP-1₁₉ haplotypes in individual populations. The PfMSP-1₁₉ haplotype patterns were then compared between parasite populations to infer the population structure and genetic differentiation of the malaria parasite.

Results

Five conserved polymorphic positions, which accounted for five distinct haplotypes, of PfMSP-1₁₉ were identified. Differences in the prevalence of PfMSP-1₁₉ haplotypes were detected in different geographical regions, with the highest levels of genetic diversity being found in the Kanchanaburi and Ranong provinces along the Thailand-Myanmar border and Trat province located at the Thailand-Cambodia border. Despite this variability, the distribution patterns of individual PfMSP-1₁₉ haplotypes seemed to be very similar across the country and over the three malarial transmission seasons, suggesting that gene flow may operate between parasite populations circulating in Thailand and the three neighboring countries.

Conclusion

The major MSP-1₁₉ haplotypes of P. falciparum populations in all endemic populations during three transmission seasons in Thailand were identified, providing basic information on the common haplotypes of MSP-1₁₉ that is of use for malaria vaccine development and inferring the population structure of P. falciparum populations in Thailand.

SigHunt: horizontal gene transfer finder optimized for eukaryotic genomes

MOTIVATION

Genomic islands (GIs) are DNA fragments incorporated into a genome through horizontal gene transfer (also called lateral gene transfer), often with functions novel for a given organism. While methods for their detection are well researched in prokaryotes, the complexity of eukaryotic genomes makes direct utilization of these methods unreliable, and so labour-intensive phylogenetic searches are used instead.

RESULTS

We present a surrogate method that investigates nucleotide base composition of the DNA sequence in a eukaryotic genome and identifies putative GIs. We calculate a genomic signature as a vector of tetranucleotide (4-mer) frequencies using a sliding window approach. Extending the neighbourhood of the sliding window, we establish a local kernel density estimate of the 4-mer frequency. We score the number of 4-mer frequencies in the sliding window that deviate from the credibility interval of their local genomic density using a newly developed discrete interval accumulative score (DIAS). To further improve the effectiveness of DIAS, we select informative 4-mers in a range of organisms using the tetranucleotide quality score developed herein. We show that the SigHunt method is computationally efficient and able to detect GIs in eukaryotic genomes that represent non-ameliorated integration. Thus, it is suited to scanning for change in organisms with different DNA composition.

AVAILABILITY AND IMPLEMENTATION

Source code and scripts freely available for download at http://www.iba.muni.cz/index-en.php?pg=research-data-analysis-tools-sighunt are implemented in C and R and are platform-independent.

CONTACT

376090@mail.muni.cz or martinkova@ivb.cz.

Malaria proteomics: insights into the parasite-host interactions in the pathogenic space

Proteomics is improving malaria research by providing global information on relevant protein sets from the parasite and the host in connection with its cellular structures and specific functions. In the last decade, reports have described biologically significant elements in the proteome of Plasmodium, which are selectively targeted and quantified, allowing for sensitive and high-throughput comparisons. The identification of molecules by which the parasite and the host react during the malaria infection is crucial to the understanding of the underlying pathogenic mechanisms. Hence, proteomics is playing a major role by defining the elements within the pathogenic space between both organisms that change across the parasite life cycle in association with the host transformation and response. Proteomics has identified post-translational modifications in the parasite and the host that are discussed in terms of functional interactions in malaria parasitism. Furthermore, the contribution of proteomics to the investigation of immunogens for potential vaccine candidates is summarized. The malaria-specific technological advances in proteomics are particularly suited now for identifying host-parasite interactions that could lead to promising targets for therapy, diagnosis or prevention. In this review, we examine the knowledge gained on the biology, pathogenesis, immunity and diagnosis of Plasmodium infection from recent proteomic studies. This article is part of a Special Issue entitled: Trends in Microbial Proteomics.

A nondiscriminating glutamyl-tRNA synthetase in the plasmodium apicoplast: the first enzyme in an indirect aminoacylation pathway

The malaria parasite Plasmodium falciparum and related organisms possess a relict plastid known as the apicoplast. Apicoplast protein synthesis is a validated drug target in malaria because antibiotics that inhibit translation in prokaryotes also inhibit apicoplast protein synthesis and are sometimes used for malaria prophylaxis or treatment. We identified components of an indirect aminoacylation pathway for Gln-tRNA(Gln) biosynthesis in Plasmodium that we hypothesized would be essential for apicoplast protein synthesis. Here, we report our characterization of the first enzyme in this pathway, the apicoplast glutamyl-tRNA synthetase (GluRS). We expressed the recombinant P. falciparum enzyme in Escherichia coli, showed that it is nondiscriminating because it glutamylates both apicoplast tRNA(Glu) and tRNA(Gln), determined its kinetic parameters, and demonstrated its inhibition by a known bacterial GluRS inhibitor. We also localized the Plasmodium berghei ortholog to the apicoplast in blood stage parasites but could not delete the PbGluRS gene. These data show that Gln-tRNA(Gln) biosynthesis in the Plasmodium apicoplast proceeds via an essential indirect aminoacylation pathway that is reminiscent of bacteria and plastids.

Selective inhibition of an apicoplastic aminoacyl-tRNA synthetase from Plasmodium falciparum

The resistance of malaria parasites to available drugs continues to grow, and this makes the need for new antimalarial therapies pressing. Aminoacyl-tRNA synthetases (ARSs) are essential enzymes and well-established antibacterial targets and so constitute a promising set of targets for the development of new antimalarials. Despite their potential as drug targets, apicoplastic ARSs remain unexplored. We have characterized the lysylation system of Plasmodium falciparum, and designed, synthesized, and tested a set of inhibitors based on the structure of the natural substrate intermediate: lysyl-adenylate. Here we demonstrate that selective inhibition of apicoplastic ARSs is feasible and describe new compounds that that specifically inhibit Plasmodium apicoplastic lysyl-tRNA synthetase and show antimalarial activities in the micromolar range.

Enhanced purification and characterization of the PfeIF4A (PfH45) helicase from Plasmodium falciparum using a codon-optimised clone

With the intention of investigating the DNA strand displacement properties of Plasmodium falciparum helicase PfeIF4A (formerly known as PfH45) a codon-optimized gene for expression in Escherichia coli has been produced. Several histidine-containing proteins with intrinsic helicase activity were captured from the bacterial sonicate by initial Ni(2+)-chromatography. Heparin and size-exclusion steps were subsequently required for unambiguous PfeIF4A purification. This strategy generated an active recombinant protein of significantly improved yield in comparison to previously published studies (~4.2 mg/g wet weight of cells). Helicase unwinding assays confirmed a bipolar activity, but revealed a preference for unwinding a free 3'-end, with a rate of displacement in the 3'-5' direction 2-fold higher than that in the 5'-3' direction. DNA constructs with two, three or four blunt ends were not unwound. Studies confirmed the enzyme to be Mg(2+)-dependent, optimally active at 37°C and had a background ATP turnover rate of 23.16±1.74 pmol/min, which in the presence of single- or double-stranded DNA doubled to 42.92±3.21 pmol/min.

Centromere plasmid: a new genetic tool for the study of Plasmodium falciparum

The introduction of transgenes into Plasmodium falciparum, a highly virulent human malaria parasite, has been conducted either by single crossover recombination or by using episomal plasmids. However, these techniques remain insufficient because of the low transfection efficiency and the low frequency of recombination. To improve the genetic manipulation of P. falciparum, we developed the centromere plasmid as a new genetic tool. First, we attempted to clone all of the predicted centromeres from P. falciparum into E. coli cells but failed because of the high A/T contents of these sequences. To overcome this difficulty, we identified the common sequence features of the centromere of Plasmodium spp. and designed a small centromere that retained those features. The centromere plasmid constructed with the small centromere sequence, pFCEN, segregated into daughter parasites with approximately 99% efficiency, resulting in the stable maintenance of this plasmid in P. falciparum even in the absence of drug selection. This result demonstrated that the small centromere sequence harboured in pFCEN could function as an actual centromere in P. falciparum. In addition, transgenic parasites were more rapidly generated when using pFCEN than when using the control plasmid, which did not contain the centromere sequence. Furthermore, in contrast to the control plasmid, pFCEN did not form concatemers and, thus, was maintained as a single copy over multiple cell divisions. These unique properties of the pFCEN plasmid will solve the current technical limitations of the genetic manipulation of P. falciparum, and thus, this plasmid will become a standard genetic tool for the study of this parasite.

A high-coverage artificial chromosome library for the genome-wide screening of drug-resistance genes in malaria parasites

The global spread of drug-resistant parasites is a serious problem for the treatment of malaria. Although identifying drug-resistance genes is crucial for the efforts against resistant parasites, an effective approach has not yet been developed. Here, we report a robust method for identifying resistance genes from parasites by using a Plasmodium artificial chromosome (PAC). Large genomic DNA fragments (10–50 kb) from the drug-resistant rodent malaria parasite Plasmodium berghei were ligated into the PAC and directly introduced into the drug-sensitive (i.e., wild-type) parasite by electroporation, resulting in a PAC library that encompassed the whole genomic sequence of the parasite. Subsequently, the transformed parasites that acquired resistance were selected by screening with the drug, and the resistance gene in the PAC was successfully identified. Furthermore, the drug-resistance gene was identified from a PAC library that was made from the pyrimethamine-resistant parasite Plasmodium chabaudi, further demonstrating the utility of our method. This method will promote the identification of resistance genes and contribute to the global fight against drug-resistant parasites.

Identification of intracellular and plasma membrane calcium channel homologues in pathogenic parasites

Ca²⁺ channels regulate many crucial processes within cells and their abnormal activity can be damaging to cell survival, suggesting that they might represent attractive therapeutic targets in pathogenic organisms. Parasitic diseases such as malaria, leishmaniasis, trypanosomiasis and schistosomiasis are responsible for millions of deaths each year worldwide. The genomes of many pathogenic parasites have recently been sequenced, opening the way for rational design of targeted therapies. We analyzed genomes of pathogenic protozoan parasites as well as the genome of Schistosoma mansoni, and show the existence within them of genes encoding homologues of mammalian intracellular Ca²⁺ release channels: inositol 1,4,5-trisphosphate receptors (IP₃Rs), ryanodine receptors (RyRs), two-pore Ca²⁺ channels (TPCs) and intracellular transient receptor potential (Trp) channels. The genomes of Trypanosoma, Leishmania and S. mansoni parasites encode IP₃R/RyR and Trp channel homologues, and that of S. mansoni additionally encodes a TPC homologue. In contrast, apicomplexan parasites lack genes encoding IP₃R/RyR homologues and possess only genes encoding TPC and Trp channel homologues (Toxoplasma gondii) or Trp channel homologues alone. The genomes of parasites also encode homologues of mammalian Ca²⁺influx channels, including voltage-gated Ca²⁺ channels and plasma membrane Trp channels. The genome of S. mansoni also encodes Orai Ca²⁺ channel and STIM Ca²⁺ sensor homologues, suggesting that store-operated Ca²⁺ entry may occur in this parasite. Many anti-parasitic agents alter parasite Ca²⁺ homeostasis and some are known modulators of mammalian Ca²⁺ channels, suggesting that parasite Ca²⁺ channel homologues might be the targets of some current anti-parasitic drugs. Differences between human and parasite Ca²⁺ channels suggest that pathogen-specific targeting of these channels may be an attractive therapeutic prospect.

Maturation of thyroglobulin protein region I

In vertebrates, the thyroglobulin (Tg) gene product must be exported to the lumen of thyroid follicles for thyroid hormone synthesis. In toto, Tg is composed of multiple type-1 repeats connected by linker and hinge (altogether considered as "region I," nearly 1,200 residues); regions II-III (~720 residues); and cholinesterase-like (ChEL) domain (~570 residues). Regions II-III and ChEL rapidly acquire competence for secretion, yet regions I-II-III require 20 min to become a partially mature disulfide isomer; stabilization of a fully oxidized form requires ChEL. Transition from partially mature to mature Tg occurs as a discrete "jump" in mobility by nonreducing SDS-PAGE, suggesting formation of at most a few final pairings of Cys residues that may be separated by significant intervening primary sequence. Using two independent approaches, we have investigated which portion of Tg is engaged in this late stage of its maturation. First, we demonstrate that this event is linked to oxidation involving region I. Introduction of the Tg-C1245R mutation in the hinge (identical to that causing human goitrous hypothyroidism) inhibits this maturation, although the Cys-1245 partner remains unidentified. Second, we find that Tg truncated after its fourth type-1 repeat is a fully independent secretory protein. Together, the data indicate that final acquisition of secretory competence includes conformational maturation in the interval between linker and hinge segments of region I.

Single-molecule analysis of genome rearrangements in cancer

Rearrangements of the genome can be detected by microarray methods and massively parallel sequencing, which identify copy-number alterations and breakpoint junctions, but these techniques are poorly suited to reconstructing the long-range organization of rearranged chromosomes, for example, to distinguish between translocations and insertions. The single-DNA-molecule technique HAPPY mapping is a method for mapping normal genomes that should be able to analyse genome rearrangements, i.e. deviations from a known genome map, to assemble rearrangements into a long-range map. We applied HAPPY mapping to cancer cell lines to show that it could identify rearrangement of genomic segments, even in the presence of normal copies of the genome. We could distinguish a simple interstitial deletion from a copy-number loss at an inversion junction, and detect a known translocation. We could determine whether junctions detected by sequencing were on the same chromosome, by measuring their linkage to each other, and hence map the rearrangement. Finally, we mapped an uncharacterized reciprocal translocation in the T-47D breast cancer cell line to about 2 kb and hence cloned the translocation junctions. We conclude that HAPPY mapping is a versatile tool for determining the structure of rearrangements in the human genome.

Iterative Correction of Reference Nucleotides (iCORN) using second generation sequencing technology

MOTIVATION

The accuracy of reference genomes is important for downstream analysis but a low error rate requires expensive manual interrogation of the sequence. Here, we describe a novel algorithm (Iterative Correction of Reference Nucleotides) that iteratively aligns deep coverage of short sequencing reads to correct errors in reference genome sequences and evaluate their accuracy.

RESULTS

Using Plasmodium falciparum (81% A + T content) as an extreme example, we show that the algorithm is highly accurate and corrects over 2000 errors in the reference sequence. We give examples of its application to numerous other eukaryotic and prokaryotic genomes and suggest additional applications.

AVAILABILITY

The software is available at http://icorn.sourceforge.net

The role of platelets in the pathogenesis of cerebral malaria

Malaria is a major cause of morbidity and mortality in the developing world and cerebral malaria is responsible for the majority of malaria-associated deaths. There is a strong association between thrombocytopenia and outcome in malaria, suggesting a role for platelets in the pathogenesis of malaria. This thrombocytopenia is likely due to platelet activation possibly through an interaction between PfEMP1 on plasmodium and CD36 on platelets. Platelet activation by plasmodium has two potential consequences. It can lead to the formation of micro-aggregates of infected red blood cells and platelets which can occlude blood vessels and it also leads to binding to and activation of the endothelium.

High antibody titer against apical membrane antigen-1 is required to protect against malaria in the Aotus model

A Plasmodium falciparum 3D7 strain Apical Membrane Antigen-1 (AMA1) vaccine, formulated with AS02_A adjuvant, slowed parasite growth in a recent Phase 1/2a trial, however sterile protection was not observed. We tested this AS02_A, and a Montanide ISA720 (ISA) formulation of 3D7 AMA1 in Aotus monkeys. The 3D7 parasite does not invade Aotus erythrocytes, hence two heterologous strains, FCH/4 and FVO, were used for challenge, FCH/4 AMA1 being more homologous to 3D7 than FVO AMA1. Following three vaccinations, the monkeys were challenged with 50,000 FCH/4 or 10,000 FVO parasites. Three of the six animals in the AMA+ISA group were protected against FCH/4 challenge. One monkey did not become parasitemic, another showed only a short period of low level parasitemia that self-cured, and a third animal showed a delay before exhibiting its parasitemic phase. This is the first protection shown in primates with a recombinant P. falciparum AMA1 without formulation in Freund's complete adjuvant. No animals in the AMA+AS02_A group were protected, but this group exhibited a trend towards reduced growth rate. A second group of monkeys vaccinated with AMA+ISA vaccine was not protected against FVO challenge, suggesting strain-specificity of AMA1-based protection. Protection against FCH/4 strain correlated with the quantity of induced antibodies, as the protected animals were the only ones to have in vitro parasite growth inhibitory activity of >70% at 1∶10 serum dilution; immuno-fluorescence titers >8,000; ELISA titers against full-length AMA1 >300,000 and ELISA titer against AMA1 domains1+2 >100,000. A negative correlation between log ELISA titer and day 11 cumulative parasitemia (Spearman rank r = −0.780, p value = 0.0001), further confirmed the relationship between antibody titer and protection. High titers of cross-strain inhibitory antibodies against AMA1 are therefore critical to confer solid protection, and the Aotus model can be used to down-select future AMA1 formulations, prior to advanced human trials.

A comprehensive Plasmodium falciparum protein interaction map reveals a distinct architecture of a core interactome

We derive a map of protein interactions in the parasite Plasmodium falciparum from conserved interactions in Saccharomyces cerevisiae, Caenorhabditis elegans, Drosophila melanogaster, and Escherichia coli and pool them with experimental interaction data. The application of a clique-percolation algorithm allows us to find overlapping clusters, strongly correlated with yeast specific conserved protein complexes. Such clusters contain core activities that govern gene expression, largely dominated by components of protein production and degradation processes as well as RNA metabolism. A critical role of protein hubs in the interactome of P. falciparum is supported by their appearance in multiple clusters and the tendencies of their interactions to reach into many distinct protein clusters. Parasite proteins with a human ortholog tend to appear in single complexes. Annotating each protein with the stage where it is maximally expressed we observe a high level of cluster integrity in the ring stage. While we find no signal in the trophozoite phase, expression patterns are reversed in the schizont phase, implying a preponderance of parasite specific functions in this late, invasive schizont stage. As such, the inference of potential protein interactions and their analysis contributes to our understanding of the parasite, indicating basic pathways and processes as unique targets for therapeutic intervention.