Apicoplast ribosomal protein S10-V127M enhances artemisinin resistance of a Kelch13 transgenic Plasmodium falciparum.

BACKGROUND: The resistance of Plasmodium falciparum to artemisinin-based (ART) drugs, the front-line drug family used in artemisinin-based combination therapy (ACT) for treatment of malaria, is of great concern. Mutations in the kelch13 (k13) gene (for example, those resulting in the Cys580Tyr [C580Y] variant) were identified as genetic markers for ART-resistant parasites, which suggests they are associated with resistance mechanisms. However, not all resistant parasites contain a k13 mutation, and clearly greater understanding of resistance mechanisms is required. A genome-wide association study (GWAS) found single nucleotide polymorphisms associated with ART-resistance in fd (ferredoxin), arps10 (apicoplast ribosomal protein S10), mdr2 (multidrug resistance protein 2), and crt (chloroquine resistance transporter), in addition to k13 gene mutations, suggesting that these alleles contribute to the resistance phenotype. The importance of the FD and ARPS10 variants in ART resistance was then studied since both proteins likely function in the apicoplast, which is a location distinct from that of K13. METHODS: The reported mutations were introduced, together with a mutation to produce the k13-C580Y variant into the ART-sensitive 3D7 parasite line and the effect on ART-susceptibility using the 0-3 h ring survival assay (RSA0-3 h) was investigated. RESULTS AND CONCLUSION: Introducing both fd-D193Y and arps10-V127M into a k13-C580Y-containing parasite, but not a wild-type k13 parasite, increased survival of the parasite in the RSA0-3 h. The results suggest epistasis of arps10 and k13, with arps10-V127M a modifier of ART susceptibility in different k13 allele backgrounds.

PTEX is an essential nexus for protein export in malaria parasites.

During the blood stages of malaria, several hundred parasite-encoded proteins are exported beyond the double-membrane barrier that separates the parasite from the host cell cytosol. These proteins have a variety of roles that are essential to virulence or parasite growth. There is keen interest in understanding how proteins are exported and whether common machineries are involved in trafficking the different classes of exported proteins. One potential trafficking machine is a protein complex known as the Plasmodium translocon of exported proteins (PTEX). Although PTEX has been linked to the export of one class of exported proteins, there has been no direct evidence for its role and scope in protein translocation. Here we show, through the generation of two parasite lines defective for essential PTEX components (HSP101 or PTEX150), and analysis of a line lacking the non-essential component TRX2 (ref. 12), greatly reduced trafficking of all classes of exported proteins beyond the double membrane barrier enveloping the parasite. This includes proteins containing the PEXEL motif (RxLxE/Q/D) and PEXEL-negative exported proteins (PNEPs). Moreover, the export of proteins destined for expression on the infected erythrocyte surface, including the major virulence factor PfEMP1 in Plasmodium falciparum, was significantly reduced in PTEX knockdown parasites. PTEX function was also essential for blood-stage growth, because even a modest knockdown of PTEX components had a strong effect on the parasite's capacity to complete the erythrocytic cycle both in vitro and in vivo. Hence, as the only known nexus for protein export in Plasmodium parasites, and an essential enzymic machine, PTEX is a prime drug target.

Inhibition of Plasmepsin V activity demonstrates its essential role in protein export, PfEMP1 display, and survival of malaria parasites.

The malaria parasite Plasmodium falciparum exports several hundred proteins into the infected erythrocyte that are involved in cellular remodeling and severe virulence. The export mechanism involves the Plasmodium export element (PEXEL), which is a cleavage site for the parasite protease, Plasmepsin V (PMV). The PMV gene is refractory to deletion, suggesting it is essential, but definitive proof is lacking. Here, we generated a PEXEL-mimetic inhibitor that potently blocks the activity of PMV isolated from P. falciparum and Plasmodium vivax. Assessment of PMV activity in P. falciparum revealed PEXEL cleavage occurs cotranslationaly, similar to signal peptidase. Treatment of P. falciparum-infected erythrocytes with the inhibitor caused dose-dependent inhibition of PEXEL processing as well as protein export, including impaired display of the major virulence adhesin, PfEMP1, on the erythrocyte surface, and cytoadherence. The inhibitor killed parasites at the trophozoite stage and knockdown of PMV enhanced sensitivity to the inhibitor, while overexpression of PMV increased resistance. This provides the first direct evidence that PMV activity is essential for protein export in Plasmodium spp. and for parasite survival in human erythrocytes and validates PMV as an antimalarial drug target.

Fast processing of microscopic images using object-based extended depth of field.

BACKGROUND: Microscopic analysis requires that foreground objects of interest, e.g. cells, are in focus. In a typical microscopic specimen, the foreground objects may lie on different depths of field necessitating capture of multiple images taken at different focal planes. The extended depth of field (EDoF) technique is a computational method for merging images from different depths of field into a composite image with all foreground objects in focus. Composite images generated by EDoF can be applied in automated image processing and pattern recognition systems. However, current algorithms for EDoF are computationally intensive and impractical, especially for applications such as medical diagnosis where rapid sample turnaround is important. Since foreground objects typically constitute a minor part of an image, the EDoF technique could be made to work much faster if only foreground regions are processed to make the composite image. We propose a novel algorithm called object-based extended depths of field (OEDoF) to address this issue. METHODS: The OEDoF algorithm consists of four major modules: 1) color conversion, 2) object region identification, 3) good contrast pixel identification and 4) detail merging. First, the algorithm employs color conversion to enhance contrast followed by identification of foreground pixels. A composite image is constructed using only these foreground pixels, which dramatically reduces the computational time. RESULTS: We used 250 images obtained from 45 specimens of confirmed malaria infections to test our proposed algorithm. The resulting composite images with all in-focus objects were produced using the proposed OEDoF algorithm. We measured the performance of OEDoF in terms of image clarity (quality) and processing time. The features of interest selected by the OEDoF algorithm are comparable in quality with equivalent regions in images processed by the state-of-the-art complex wavelet EDoF algorithm; however, OEDoF required four times less processing time. CONCLUSIONS: This work presents a modification of the extended depth of field approach for efficiently enhancing microscopic images. This selective object processing scheme used in OEDoF can significantly reduce the overall processing time while maintaining the clarity of important image features. The empirical results from parasite-infected red cell images revealed that our proposed method efficiently and effectively produced in-focus composite images. With the speed improvement of OEDoF, this proposed algorithm is suitable for processing large numbers of microscope images, e.g., as required for medical diagnosis.

Automatic DNA Diagnosis for 1D Gel Electrophoresis Images using Bio-image Processing Technique.

BACKGROUND: DNA gel electrophoresis is a molecular biology technique for separating different sizes of DNA fragments. Applications of DNA gel electrophoresis include DNA fingerprinting (genetic diagnosis), size estimation of DNA, and DNA separation for Southern blotting. Accurate interpretation of DNA banding patterns from electrophoretic images can be laborious and error prone when a large number of bands are interrogated manually. Although many bio-imaging techniques have been proposed, none of them can fully automate the typing of DNA owing to the complexities of migration patterns typically obtained. RESULTS: We developed an image-processing tool that automatically calls genotypes from DNA gel electrophoresis images. The image processing workflow comprises three main steps: 1) lane segmentation, 2) extraction of DNA bands and 3) band genotyping classification. The tool was originally intended to facilitate large-scale genotyping analysis of sugarcane cultivars. We tested the proposed tool on 10 gel images (433 cultivars) obtained from polyacrylamide gel electrophoresis (PAGE) of PCR amplicons for detecting intron length polymorphisms (ILP) on one locus of the sugarcanes. These gel images demonstrated many challenges in automated lane/band segmentation in image processing including lane distortion, band deformity, high degree of noise in the background, and bands that are very close together (doublets). Using the proposed bio-imaging workflow, lanes and DNA bands contained within are properly segmented, even for adjacent bands with aberrant migration that cannot be separated by conventional techniques. The software, called GELect, automatically performs genotype calling on each lane by comparing with an all-banding reference, which was created by clustering the existing bands into the non-redundant set of reference bands. The automated genotype calling results were verified by independent manual typing by molecular biologists. CONCLUSIONS: This work presents an automated genotyping tool from DNA gel electrophoresis images, called GELect, which was written in Java and made available through the imageJ framework. With a novel automated image processing workflow, the tool can accurately segment lanes from a gel matrix, intelligently extract distorted and even doublet bands that are difficult to identify by existing image processing tools. Consequently, genotyping from DNA gel electrophoresis can be performed automatically allowing users to efficiently conduct large scale DNA fingerprinting via DNA gel electrophoresis. The software is freely available from http://www.biotec.or.th/gi/tools/gelect.

Plasmodium parasites mount an arrest response to dihydroartemisinin, as revealed by whole transcriptome shotgun sequencing (RNA-seq) and microarray study.

BACKGROUND: Control of malaria is threatened by emerging parasite resistance to artemisinin and derivative drug (ART) therapies. The molecular detail of how Plasmodium malaria parasites respond to ART and how this could contribute to resistance are not well understood. To address this question, we performed a transcriptomic study of dihydroartemisinin (DHA) response in P. falciparum K1 strain and in P. berghei ANKA strain using microarray and RNA-seq technology. RESULTS: Microarray data from DHA-treated P. falciparum trophozoite stage parasites revealed a response pattern that is overall less trophozoite-like and more like the other stages of asexual development. A meta-analysis of these data with previously published data from other ART treatments revealed a set of common differentially expressed genes. Notably, ribosomal protein genes are down-regulated in response to ART. A similar pattern of trophozoite transcriptomic change was observed from RNA-seq data. RNA-seq data from DHA-treated P. falciparum rings reveal a more muted response, although there is considerable overlap of differentially expressed genes with DHA-treated trophozoites. No genes are differentially expressed in DHA-treated P. falciparum schizonts. The transcriptional response of P. berghei to DHA treatment in vivo in infected mice is similar to the P. falciparum in vitro culture ring and trophozoite responses, in which ribosomal protein genes are notably down-regulated. CONCLUSIONS: Ring and trophozoite stage Plasmodium respond to ART by arresting metabolic processes such as protein synthesis and glycolysis. This response can be protective in rings, as shown by the phenomenon of dormancy. In contrast, this response is not as protective in trophozoites owing to their commitment to a highly active and vulnerable metabolic state. The lower metabolic demands of schizonts could explain why they are less sensitive and unresponsive to ART. The ART response pattern is revealed clearly from RNA-seq data, suggesting that this technology is of great utility for studying drug response in Plasmodium.

Antimalarial target vulnerability of the putative Plasmodium falciparum methionine synthase.

Background: Plasmodium falciparum possesses a cobalamin-dependent methionine synthase (MS). MS is putatively encoded by the PF3D7_1233700 gene, which is orthologous and syntenic in Plasmodium. However, its vulnerability as an antimalarial target has not been assessed. Methods: We edited the PF3D7_1233700 and PF3D7_0417200 (dihydrofolate reductase-thymidylate synthase, DHFR-TS) genes and obtained transgenic P. falciparum parasites expressing epitope-tagged target proteins under the control of the glmS ribozyme. Conditional loss-of-function mutants were obtained by treating transgenic parasites with glucosamine. Results: DHFR-TS, but not MS mutants showed a significant proliferation defect over 96 h, suggesting that P. falciparum MS is not a vulnerable antimalarial target.

Antimalarial mechanism of action of the natural product 9-methoxystrobilurin G.

The natural product 9-methoxystrobilurin G (9MG) from Favolaschia spp basidiomycetes is a potent and selective antimalarial. The mechanism of action of 9MG is unknown. We induced 9MG resistance in Plasmodium falciparum 3D7 and Dd2 strains and identified mutations associated with resistance by genome sequencing. All 9MG-resistant clones possessed missense mutations in the cytochrome b (CYTB) gene, a key component of mitochondrial complex III. The mutations map to the quinol oxidation site of CYTB, which is also the target of antimalarials such as atovaquone. In a complementary approach to identify protein targets of 9MG, a photoactivatable derivative of 9MG was synthesized and applied in chemoproteomic-based target profiling. Three components of mitochondrial complex III (QCR7, QCR9, and COX15) were specifically enriched consistent with 9MG targeting CYTB and complex III function in P. falciparum. Inhibition of complex III activity by 9MG was confirmed by ubiquinone cytochrome c reductase assay using P. falciparum extract. The findings from this study may be useful for developing novel antimalarials targeting CYTB.

Susceptibility of Southeast Asian Plasmodium falciparum isolates to P218.

A major threat to the goal of eliminating malaria, particularly in Southeast Asia, is the spread of Plasmodium falciparum resistant to artemisinin-based combination therapies. P218 is a drug candidate designed to combat antifolate-sensitive and -resistant parasites. However, there is no evidence that P218 is effective against artemisinin-resistant P. falciparum. This report investigated the susceptibilities of 10 parasite isolates from Southeast Asia to P218 and other antimalarial drugs. All isolates with different levels of artemisinin resistance were genetically distinct from one another, although common haplotypes associated with antimalarial resistance were identified. All isolates were highly resistant to pyrimethamine, and none of them were significantly less sensitive to P218 than the pyrimethamine-resistant laboratory strain V1/S. Significant differences in sensitivity to other types of antimalarials (mefloquine, atovaquone and chloroquine) compared with V1/S were found for some isolates, although the differences were not clinically relevant. P218 is thus efficacious against multi-drug (including artemisinin-resistant P. falciparum.

Antagonistic antimalarial properties of a methoxyamino chalcone derivative and 3-hydroxypyridinones in combination with dihydroartemisinin against Plasmodium falciparum.

Background: The spread of artemisinin (ART)-resistant Plasmodium falciparum threatens the control of malaria. Mutations in the propeller domains of P. falciparum Kelch13 (k13) are strongly associated with ART resistance. Ferredoxin (Fd), a component of the ferredoxin/NADP+ reductase (Fd/FNR) redox system, is essential for isoprenoid precursor synthesis in the plasmodial apicoplast, which is important for K13-dependent hemoglobin trafficking and ART activation. Therefore, Fd is an antimalarial drug target and fd mutations may modulate ART sensitivity. We hypothesized that loss of Fd/FNR function enhances the effect of k13 mutation on ART resistance. Methods: In this study, methoxyamino chalcone (C3), an antimalarial compound that has been reported to inhibit the interaction of recombinant Fd and FNR proteins, was used as a chemical inhibitor of the Fd/FNR redox system. We investigated the inhibitory effects of dihydroartemisinin (DHA), C3, and iron chelators including deferiprone (DFP), 1-(N-acetyl-6-aminohexyl)-3-hydroxy-2-methylpyridin-4-one (CM1) and deferiprone-resveratrol hybrid (DFP-RVT) against wild-type (WT), k13 mutant, fd mutant, and k13 fd double mutant P. falciparum parasites. Furthermore, we investigated the pharmacological interaction of C3 with DHA, in which the iron chelators were used as reference ART antagonists. Results: C3 showed antimalarial potency similar to that of the iron chelators. As expected, combining DHA with C3 or iron chelators exhibited a moderately antagonistic effect. No differences were observed among the mutant parasites with respect to their sensitivity to C3, iron chelators, or the interactions of these compounds with DHA. Discussion: The data suggest that inhibitors of the Fd/FNR redox system should be avoided as ART partner drugs in ART combination therapy for treating malaria.

iLOCi: a SNP interaction prioritization technique for detecting epistasis in genome-wide association studies.

BACKGROUND: Genome-wide association studies (GWAS) do not provide a full account of the heritability of genetic diseases since gene-gene interactions, also known as epistasis are not considered in single locus GWAS. To address this problem, a considerable number of methods have been developed for identifying disease-associated gene-gene interactions. However, these methods typically fail to identify interacting markers explaining more of the disease heritability over single locus GWAS, since many of the interactions significant for disease are obscured by uninformative marker interactions e.g., linkage disequilibrium (LD). RESULTS: In this study, we present a novel SNP interaction prioritization algorithm, named iLOCi (Interacting Loci). This algorithm accounts for marker dependencies separately in case and control groups. Disease-associated interactions are then prioritized according to a novel ranking score calculated from the difference in marker dependencies for every possible pair between case and control groups. The analysis of a typical GWAS dataset can be completed in less than a day on a standard workstation with parallel processing capability. The proposed framework was validated using simulated data and applied to real GWAS datasets using the Wellcome Trust Case Control Consortium (WTCCC) data. The results from simulated data showed the ability of iLOCi to identify various types of gene-gene interactions, especially for high-order interaction. From the WTCCC data, we found that among the top ranked interacting SNP pairs, several mapped to genes previously known to be associated with disease, and interestingly, other previously unreported genes with biologically related roles. CONCLUSION: iLOCi is a powerful tool for uncovering true disease interacting markers and thus can provide a more complete understanding of the genetic basis underlying complex disease. The program is available for download at http://www4a.biotec.or.th/GI/tools/iloci.

Stepwise acquisition of pyrimethamine resistance in the malaria parasite.

The spread of high-level pyrimethamine resistance in Africa threatens to curtail the therapeutic lifetime of antifolate antimalarials. We studied the possible evolutionary pathways in the evolution of pyrimethamine resistance using an approach in which all possible mutational intermediates were created by site-directed mutagenesis and assayed for their level of drug resistance. The coding sequence for dihydrofolate reductase (DHFR) from the malaria parasite Plasmodium falciparum was mutagenized, and tests were carried out in Escherichia coli under conditions in which the endogenous bacterial enzyme was selectively inhibited. We studied 4 key amino acid replacements implicated in pyrimethamine resistance: N51I, C59R, S108N, and I164L. Using empirical estimates of the mutational spectrum in P. falciparum and probabilities of fixation based on the relative levels of resistance, we found that the predicted favored pathways of drug resistance are consistent with those reported in previous kinetic studies, as well as DHFR polymorphisms observed in natural populations. We found that 3 pathways account for nearly 90% of the simulated realizations of the evolution of pyrimethamine resistance. The most frequent pathway (S108N and then C59R, N51I, and I164L) accounts for more than half of the simulated realizations. Our results also suggest an explanation for why I164L is detected in Southeast Asia and South America, but not at significant frequencies in Africa.

Transcriptomic complexity of the human malaria parasite Plasmodium falciparum revealed by long-read sequencing.

The Plasmodium falciparum human malaria parasite genome is incompletely annotated and does not accurately represent the transcriptomic diversity of this species. To address this need, we performed long-read transcriptomic sequencing. 5' capped mRNA was enriched from samples of total and nuclear-fractionated RNA from intra-erythrocytic stages and converted to cDNA library. The cDNA libraries were sequenced on PacBio and Nanopore long-read platforms. 12,495 novel isoforms were annotated from the data. Alternative 5' and 3' ends represent the majority of isoform events among the novel isoforms, with retained introns being the next most common event. The majority of alternative 5' ends correspond to genomic regions with features similar to those of the reference transcript 5' ends. However, a minority of alternative 5' ends showed markedly different features, including locations within protein-coding regions. Alternative 3' ends showed similar features to the reference transcript 3' ends, notably adenine-rich termination signals. Distinguishing features of retained introns could not be observed, except for a tendency towards shorter length and greater GC content compared with spliced introns. Expression of antisense and retained intron isoforms was detected at different intra-erythrocytic stages, suggesting developmental regulation of these isoform events. To gain insights into the possible functions of the novel isoforms, their protein-coding potential was assessed. Variants of P. falciparum proteins and novel proteins encoded by alternative open reading frames suggest that P. falciparum has a greater proteomic repertoire than the current annotation. We provide a catalog of annotated transcripts and encoded alternative proteins to support further studies on gene and protein regulation of this pathogen.

Study of large and highly stratified population datasets by combining iterative pruning principal component analysis and structure.

BACKGROUND: The ever increasing sizes of population genetic datasets pose great challenges for population structure analysis. The Tracy-Widom (TW) statistical test is widely used for detecting structure. However, it has not been adequately investigated whether the TW statistic is susceptible to type I error, especially in large, complex datasets. Non-parametric, Principal Component Analysis (PCA) based methods for resolving structure have been developed which rely on the TW test. Although PCA-based methods can resolve structure, they cannot infer ancestry. Model-based methods are still needed for ancestry analysis, but they are not suitable for large datasets. We propose a new structure analysis framework for large datasets. This includes a new heuristic for detecting structure and incorporation of the structure patterns inferred by a PCA method to complement STRUCTURE analysis. RESULTS: A new heuristic called EigenDev for detecting population structure is presented. When tested on simulated data, this heuristic is robust to sample size. In contrast, the TW statistic was found to be susceptible to type I error, especially for large population samples. EigenDev is thus better-suited for analysis of large datasets containing many individuals, in which spurious patterns are likely to exist and could be incorrectly interpreted as population stratification. EigenDev was applied to the iterative pruning PCA (ipPCA) method, which resolves the underlying subpopulations. This subpopulation information was used to supervise STRUCTURE analysis to infer patterns of ancestry at an unprecedented level of resolution. To validate the new approach, a bovine and a large human genetic dataset (3945 individuals) were analyzed. We found new ancestry patterns consistent with the subpopulations resolved by ipPCA. CONCLUSIONS: The EigenDev heuristic is robust to sampling and is thus superior for detecting structure in large datasets. The application of EigenDev to the ipPCA algorithm improves the estimation of the number of subpopulations and the individual assignment accuracy, especially for very large and complex datasets. Furthermore, we have demonstrated that the structure resolved by this approach complements parametric analysis, allowing a much more comprehensive account of population structure. The new version of the ipPCA software with EigenDev incorporated can be downloaded from http://www4a.biotec.or.th/GI/tools/ippca.

Global gene expression profiling of Plasmodium falciparum in response to the anti-malarial drug pyronaridine.

BACKGROUND: Pyronaridine (PN) and chloroquine (CQ) are structurally related anti-malarial drugs with primarily the same mode of action. However, PN is effective against several multidrug-resistant lines of Plasmodium falciparum, including CQ resistant lines, suggestive of important operational differences between the two drugs. METHODS: Synchronized trophozoite stage cultures of P. falciparum strain K1 (CQ resistant) were exposed to 50% inhibitory concentrations (IC50) of PN and CQ, and parasites were harvested from culture after 4 and 24 hours exposure. Global transcriptional changes effected by drug treatment were investigated using DNA microarrays. RESULTS: After a 4 h drug exposure, PN induced a greater degree of transcriptional perturbation (61 differentially expressed features) than CQ (10 features). More genes were found to respond to 24 h treatments with both drugs, and 461 features were found to be significantly responsive to one or both drugs across all treatment conditions. Filtering was employed to remove features unrelated to primary drug action, specifically features representing genes developmentally regulated, secondary stress/death related processes and sexual stage development. The only significant gene ontologies represented among the 46 remaining features after filtering relate to host exported proteins from multi-gene families. CONCLUSIONS: The malaria parasite's molecular responses to PN and CQ treatment are similar in terms of the genes and pathways affected. However, PN appears to exert a more rapid response than CQ. The faster action of PN may explain why PN is more efficacious than CQ, particularly against CQ resistant isolates. In agreement with several other microarray studies of drug action on the parasite, it is not possible, however, to discern mechanism of drug action from the drug-responsive genes.

Formation of catalytically active cross-species heterodimers of thymidylate synthase from Plasmodium falciparum and Plasmodium vivax.

Thymidylate synthase (TS) of Plasmodium dihydrofolate reductase-thymidylate synthase (DHFR-TS) functions as a homodimeric enzyme with two active sites located near the subunit interface. The dimerization is essential for catalysis, since the active site of each subunit contains amino acid residues contributed from the other TS domain. In P. falciparum DHFR-TS, it has been shown that the active sites require Cys-490 from one domain and Arg-470 donated from the other domain. Mutants of these two series can complement one another giving rise to active enzyme. Here, the potential to form cross-species heterodimers between P. falciparum and P. vivax TS has been explored. Formation of cross-species heterodimer was tested by co-transformation of TS-inactive Cys-490 mutants of P. falciparum or P. vivax with corresponding TS-inactive Arg-486 mutants of P. vivax or P. falciparum into thymidine-requiring Escherichia coli. Active heterodimers were detected by subunit complementation and 6-[(3)H]-FdUMP binding assays. All combinations of the mutants tested, except for (Pf)R470A+(Pv)C506Y, were able to form catalytically active cross-species heterodimers. The single active site formed by (Pf)R470D+(Pv)C506Y and (Pv)R486D+(Pf)C490A pairs of cross-species heterodimers has k(cat) and K(m) values similar to those of intra-species heterodimers of P. falciparum and P. vivax. This is the first report to demonstrate that the TS subunit interface between Plasmodium species is sufficiently conserved to allow formation of fully active cross-species heterodimer.

Identifying transcript 5' capped ends in Plasmodium falciparum.

BACKGROUND: The genome of the human malaria parasite Plasmodium falciparum is poorly annotated, in particular, the 5' capped ends of its mRNA transcripts. New approaches are needed to fully catalog P. falciparum transcripts for understanding gene function and regulation in this organism. METHODS: We developed a transcriptomic method based on next-generation sequencing of complementary DNA (cDNA) enriched for full-length fragments using eIF4E, a 5' cap-binding protein, and an unenriched control. DNA sequencing adapter was added after enrichment of full-length cDNA using two different ligation protocols. From the mapped sequence reads, enrichment scores were calculated for all transcribed nucleotides and used to calculate P-values of 5' capped nucleotide enrichment. Sensitivity and accuracy were increased by combining P-values from replicate experiments. Data were obtained for P. falciparum ring, trophozoite and schizont stages of intra-erythrocytic development. RESULTS: 5' capped nucleotide signals were mapped to 17,961 non-overlapping P. falciparum genomic intervals. Analysis of the dominant 5' capped nucleotide in these genomic intervals revealed the presence of two groups with distinctive epigenetic features and sequence patterns. A total of 4,512 transcripts were annotated as 5' capped based on the correspondence of 5' end with 5' capped nucleotide annotated from full-length cDNA data. DISCUSSION: The presence of two groups of 5' capped nucleotides suggests that alternative mechanisms may exist for producing 5' capped transcript ends in P. falciparum. The 5' capped transcripts that are antisense, outside of, or partially overlapping coding regions may be important regulators of gene function in P. falciparum.

PharmVIP: A Web-Based Tool for Pharmacogenomic Variant Analysis and Interpretation.

The increasing availability of next generation sequencing (NGS) for personal genomics could promote pharmacogenomics (PGx) discovery and application. However, current tools for analysis and interpretation of pharmacogenomic variants from NGS data are inadequate, as none offer comprehensive analytic functions in a simple, web-based platform. In addition, no tools exist to analyze human leukocyte antigen (HLA) genes for determining potential risks of immune-mediated adverse drug reaction (IM-ADR). We describe PharmVIP, a web-based PGx tool, for one-stop comprehensive analysis and interpretation of genome-wide variants obtained from NGS platforms. PharmVIP comprises three main interpretation modules covering analyses of pharmacogenes involved in pharmacokinetics, pharmacodynamics and IM-ADR. The Guideline module provides Clinical Pharmacogenetics Implementation Consortium (CPIC) drug guideline recommendations based on the translation of genotypic data in genes having guidelines. The HLA module reports HLA genotypes, potential adverse drug reactions, and the relevant drug guidelines. The Pharmacogenes module is employed for prioritizing variants according to variant effect on gene function. Detailed, customizable reports are provided as exportable files and as an interactive web version. PharmVIP is a new integrated NGS workflow for the PGx community to facilitate discovery and clinical application.

5' and 3' end modifications of spliceosomal RNAs in Plasmodium falciparum.

5' caps provide recognition sequences for the nuclear import of snRNAs. The 5' and 3' ends of snRNAs were studied in Plasmodium falciparum with a modified adapter ligation method, which showed that 5' ends of U1, U2, U4, U5 and U6 snRNAs are capped. In P. falciparum, the 3' ends of U1, U2, U4 and U5 snRNAs have free hydroxyl groups whereas U6 snRNA has a blocked 3' end. An immunoprecipitation assay for trimethyl guanosine caps shows that the cap structures of parasite U1-U5 snRNAs are hypermethylated while U6 snRNA may be gamma-mono-methylated. Bioinformatics analysis of proteins involved in hypermethylation and trafficking of snRNAs indicates that the methyltransferase TGS1 is present in the P. falciparum genome. PfTGS1 is larger than its orthologs and may have transmembrane domains in the C-terminus. Surprisingly, the snRNA trafficking protein Snurportin is absent from the P. falciparum genome suggesting that reminiscent of yeast, parasite snRNAs may be retained in the nucleus.

Iterative pruning PCA improves resolution of highly structured populations.

BACKGROUND: Non-random patterns of genetic variation exist among individuals in a population owing to a variety of evolutionary factors. Therefore, populations are structured into genetically distinct subpopulations. As genotypic datasets become ever larger, it is increasingly difficult to correctly estimate the number of subpopulations and assign individuals to them. The computationally efficient non-parametric, chiefly Principal Components Analysis (PCA)-based methods are thus becoming increasingly relied upon for population structure analysis. Current PCA-based methods can accurately detect structure; however, the accuracy in resolving subpopulations and assigning individuals to them is wanting. When subpopulations are closely related to one another, they overlap in PCA space and appear as a conglomerate. This problem is exacerbated when some subpopulations in the dataset are genetically far removed from others. We propose a novel PCA-based framework which addresses this shortcoming. RESULTS: A novel population structure analysis algorithm called iterative pruning PCA (ipPCA) was developed which assigns individuals to subpopulations and infers the total number of subpopulations present. Genotypic data from simulated and real population datasets with different degrees of structure were analyzed. For datasets with simple structures, the subpopulation assignments of individuals made by ipPCA were largely consistent with the STRUCTURE, BAPS and AWclust algorithms. On the other hand, highly structured populations containing many closely related subpopulations could be accurately resolved only by ipPCA, and not by other methods. CONCLUSION: The algorithm is computationally efficient and not constrained by the dataset complexity. This systematic subpopulation assignment approach removes the need for prior population labels, which could be advantageous when cryptic stratification is encountered in datasets containing individuals otherwise assumed to belong to a homogenous population.