Expression profiling of the schizont and trophozoite stages of Plasmodium falciparum with a long-oligonucleotide microarray

Size-dependent enhancement of gene expression by Plasmodium 5'UTR introns

BACKGROUND

Eukaryotic genes contain introns that are removed by the spliceosomal machinery during mRNA maturation. Introns impose a huge energetic burden on a cell; therefore, they must play an essential role in maintaining genome stability and/or regulating gene expression. Many genes (> 50%) in Plasmodium parasites contain predicted introns, including introns in 5' and 3' untranslated regions (UTR). However, the roles of UTR introns in the gene expression of malaria parasites remain unknown.

METHODS

In this study, an episomal dual-luciferase assay was developed to evaluate gene expression driven by promoters with or without a 5'UTR intron from four Plasmodium yoelii genes. To investigate the effect of the 5'UTR intron on endogenous gene expression, the pytctp gene was tagged with 3xHA at the N-terminal of the coding region, and parasites with or without the 5'UTR intron were generated using the CRISPR/Cas9 system.

RESULTS

We showed that promoters with 5'UTR introns had higher activities in driving gene expression than those without 5'UTR introns. The results were confirmed in recombinant parasites expressing an HA-tagged gene (pytctp) driven by promoter with or without 5'UTR intron. The enhancement of gene expression was intron size dependent, but not the DNA sequence, e.g. the longer the intron, the higher levels of expression. Similar results were observed when a promoter from one strain of P. yoelii was introduced into different parasite strains. Finally, the 5'UTR introns were alternatively spliced in different parasite development stages, suggesting an active mechanism employed by the parasites to regulate gene expression in various developmental stages.

CONCLUSIONS

Plasmodium 5'UTR introns enhance gene expression in a size-dependent manner; the presence of alternatively spliced mRNAs in different parasite developmental stages suggests that alternative slicing of 5'UTR introns is one of the key mechanisms in regulating parasite gene expression and differentiation.

Plasmodium falciparum topoisomerases: Emerging targets for anti-malarial therapy

Different metabolic pathways like DNA replication, transcription, and recombination generate topological constrains in the genome. These topological constraints are resolved by essential molecular machines known as topoisomerases. To bring changes in DNA topology, the topoisomerases create a single or double-stranded nick in the template DNA, hold the nicked ends to let the tangled DNA pass through, and finally re-ligate the breaks. The DNA nicking and re-ligation activities as well as ATPase activities (when present) in topoisomerases are subjected to inhibition by several anticancer and antibacterial drugs, thus establishing these enzymes as successful targets in anticancer and antibacterial therapies. The anti-topoisomerase drugs interfere with the functioning of these enzymes and result in the accumulation of DNA tangles or lethal genomic breaks, thereby promoting host cell (or organism) death. The potential of topoisomerases in the human malarial parasite, Plasmodium falciparum in antimalarial drug development has received little attention so far. Interestingly, the parasite genome encodes orthologs of topoisomerases found in eukaryotes, prokaryotes, and archaea, thus, providing an enormous opportunity for investigating these enzymes for antimalarial therapeutics. This review focuses on the features of Plasmodium falciparum topoisomerases (PfTopos) with respect to their closer counterparts in other organisms. We will discuss overall advances and basic challenges with topoisomerase research in Plasmodium falciparum and our attempts to understand the interaction of PfTopos with classical and new-generation topoisomerase inhibitors using in silico molecular docking approach. The recent episodes of parasite resistance against artemisinin, the only effective antimalarial drug at present, further highlight the significance of investigating new drug targets including topoisomerases in antimalarial therapeutics.

PfMORC protein regulates chromatin accessibility and transcriptional repression in the human malaria parasite, P. falciparum

The environmental challenges the human malaria parasite, Plasmodium falciparum, faces during its progression into its various lifecycle stages warrant the use of effective and highly regulated access to chromatin for transcriptional regulation. Microrchidia (MORC) proteins have been implicated in DNA compaction and gene silencing across plant and animal kingdoms. Accumulating evidence has shed light into the role MORC protein plays as a transcriptional switch in apicomplexan parasites. In this study, using CRISPR/Cas9 genome editing tool along with complementary molecular and genomics approaches, we demonstrate that PfMORC not only modulates chromatin structure and heterochromatin formation throughout the parasite erythrocytic cycle, but is also essential to the parasite survival. Chromatin immunoprecipitation followed by deep sequencing (ChIP-seq) experiments suggest that PfMORC binds to not only sub-telomeric regions and genes involved in antigenic variation but is also most likely a key modulator of stage transition. Protein knockdown experiments followed by chromatin conformation capture (Hi-C) studies indicate that downregulation of PfMORC induces the collapse of the parasite heterochromatin structure leading to its death. All together these findings confirm that PfMORC plays a crucial role in chromatin structure and gene regulation, validating this factor as a strong candidate for novel antimalarial strategies.

Novel insights into the role of long non-coding RNA in the human malaria parasite, Plasmodium falciparum

The complex life cycle of Plasmodium falciparum requires coordinated gene expression regulation to allow host cell invasion, transmission, and immune evasion. Increasing evidence now suggests a major role for epigenetic mechanisms in gene expression in the parasite. In eukaryotes, many lncRNAs have been identified to be pivotal regulators of genome structure and gene expression. To investigate the regulatory roles of lncRNAs in P. falciparum we explore the intergenic lncRNA distribution in nuclear and cytoplasmic subcellular locations. Using nascent RNA expression profiles, we identify a total of 1768 lncRNAs, of which 718 (~41%) are novels in P. falciparum. The subcellular localization and stage-specific expression of several putative lncRNAs are validated using RNA-FISH. Additionally, the genome-wide occupancy of several candidate nuclear lncRNAs is explored using ChIRP. The results reveal that lncRNA occupancy sites are focal and sequence-specific with a particular enrichment for several parasite-specific gene families, including those involved in pathogenesis and sexual differentiation. Genomic and phenotypic analysis of one specific lncRNA demonstrate its importance in sexual differentiation and reproduction. Our findings bring a new level of insight into the role of lncRNAs in pathogenicity, gene regulation and sexual differentiation, opening new avenues for targeted therapeutic strategies against the deadly malaria parasite.

Isoliensinine from Cissampelos pariera rhizomes exhibits potential gametocytocidal and anti-malarial activities against Plasmodium falciparum clinical isolates

BACKGROUND

The unmet demand for effective malaria transmission-blocking agents targeting the transmissible stages of Plasmodium necessitates intensive discovery efforts. In this study, a bioactive bisbenzylisoquinoline (BBIQ), isoliensinine, from Cissampelos pariera (Menispermaceae) rhizomes was identified and characterized for its anti-malarial activity.

METHODS

Malaria SYBR Green I fluorescence assay was performed to evaluate the in vitro antimalarial activity against D6, Dd2, and F32-ART5 clones, and immediate ex vivo (IEV) susceptibility for 10 freshly collected P. falciparum isolates. To determine the speed- and stage-of-action of isoliensinine, an IC₅₀ speed assay and morphological analyses were performed using synchronized Dd2 asexuals. Gametocytocidal activity against two culture-adapted gametocyte-producing clinical isolates was determined using microscopy readouts, with possible molecular targets and their binding affinities deduced in silico.

RESULTS

Isoliensinine displayed a potent in vitro gametocytocidal activity at mean IC₅₀^gam values ranging between 0.41 and 0.69 µM for Plasmodium falciparum clinical isolates. The BBIQ compound also inhibited asexual replication at mean IC₅₀^Asexual of 2.17 µM, 2.22 µM, and 2.39 µM for D6, Dd2 and F32-ART5 respectively, targeting the late-trophozoite to schizont transition. Further characterization demonstrated a considerable immediate ex vivo potency against human clinical isolates at a geometric mean IC₅₀^IEV = 1.433 µM (95% CI 0.917-2.242). In silico analyses postulated a probable anti-malarial mechanism of action by high binding affinities for four mitotic division protein kinases; Pfnek1, Pfmap2, Pfclk1, and Pfclk4. Additionally, isoliensinine was predicted to possess an optimal pharmacokinetics profile and drug-likeness properties.

CONCLUSION

These findings highlight considerable grounds for further exploration of isoliensinine as an amenable scaffold for malaria transmission-blocking chemistry and target validation.

Exploring potential of Plasmodium RUVBL proteins as anti-malarial drug target

Although malaria related cases and deaths have consistently declined over time, growing resistance to existing anti-malarial drugs in Plasmodium remains a matter of extreme concern. Since we rely so heavily on use of chemotherapy for malaria treatment and knowing that all the available anti-malarial drug will become virtually useless in the near future, we have to increase our understanding of basic biology of the parasite as well as characterize new molecular targets that can be exploited for anti-malarial therapy. In the present study, PfRUVBLs (AAA family member proteins) were evaluated for their potential as novel anti-malarial drug target candidates, using computational approaches. Virtual High-throughput screening of various pharmacophore libraries obtained from three different databases (which included, Asinex, ZINC15 & PubChem) followed by extra precision docking, resulted in identification of relevant hit compounds that showed binding affinity with the active region of PfRUVBL1 protein. Based on molecular docking data, MD simulations, and protein-ligand interaction studies, combined with toxicity assessment & ADME profiling data, at least three best hits were eventually identified that could be novel potent inhibitors of PfRUVBL1 protein and can be further tested for anti-malarial activity using in vitro protocols. Communicated by Ramaswamy H. Sarma.

A rising tide of parasite transcriptomics propels pathogen biology

Twenty years ago, the first transcriptome of the intraerythrocytic developmental cycle of the malaria parasite Plasmodium falciparum was published in PLOS Biology. Since then, transcriptomics studies have transformed the study of parasite biology.

Selective expression of variant surface antigens enables Plasmodium falciparum to evade immune clearance in vivo

Plasmodium falciparum has developed extensive mechanisms to evade host immune clearance. Currently, most of our understanding is based on in vitro studies of individual parasite variant surface antigens and how this relates to the processes in vivo is not well-understood. Here, we have used a humanized mouse model to identify parasite factors important for in vivo growth. We show that upregulation of the specific PfEMP1, VAR2CSA, provides the parasite with protection from macrophage phagocytosis and clearance in the humanized mice. Furthermore, parasites adapted to thrive in the humanized mice show reduced NK cell-mediated killing through interaction with the immune inhibitory receptor, LILRB1. Taken together, these findings reveal new insights into the molecular and cellular mechanisms that the parasite utilizes to coordinate immune escape in vivo. Identification and targeting of these specific parasite variant surface antigens crucial for immune evasion provides a unique approach for therapy.

Stochastic expression of invasion genes in Plasmodium falciparum schizonts

Genetically identical cells are known to exhibit differential phenotypes in the same environmental conditions. These phenotypic variants are linked to transcriptional stochasticity and have been shown to contribute towards adaptive flexibility of a wide range of unicellular organisms. Here, we investigate transcriptional heterogeneity and stochastic gene expression in Plasmodium falciparum by performing the quasilinear multiple annealing and looping based amplification cycles (MALBAC) based amplification and single cell RNA sequencing of blood stage schizonts. Our data reveals significant transcriptional variations in the schizont stage with a distinct group of highly variable invasion gene transcripts being identified. Moreover, the data reflects several diversification processes including putative developmental “checkpoint”; transcriptomically distinct parasite sub-populations and transcriptional switches in variable gene families (var, rifin, phist). Most of these features of transcriptional variability are preserved in isogenic parasite cell populations (albeit with a lesser amplitude) suggesting a role of epigenetic factors in cell-to-cell transcriptional variations in human malaria parasites. Lastly, we apply quantitative RT-PCR and RNA-FISH approach and confirm stochastic expression of key invasion genes, such as, msp1, msp3, msp7, eba181 and ama1 which represent prime candidates for invasion-blocking vaccines.

Genetically identical cells can be phenotypically diverse to allow adaptive flexibility in a given environment. This phenotypic diversity is driven by epigenetic and transcriptional variability. Here, Tripathi et al. perform scRNA-seq of isogenic and non-isogenic Plasmodium falciparum schizont populations to explore transcriptional heterogeneity and stochastic gene expression during the course of development.

Contribution of Transcriptome to Elucidate the Biology of Plasmodium spp

In the present review, we discuss some of the new technologies that have been applied to elucidate how Plasmodium spp escape from the immune system and subvert the host physiology to orchestrate the regulation of its biological pathways. Our manuscript describes how techniques such as microarray approaches, RNA-Seq and single-cell RNA sequencing have contributed to the discovery of transcripts and changed the concept of gene expression regulation in closely related malaria parasite species. Moreover, the text highlights the contributions of high-throughput RNA sequencing for the current knowledge of malaria parasite biology, physiology, vaccine target and the revelation of new players in parasite signaling.

Enhanced detection of minimal residual disease by targeted sequencing of phased variants in circulating tumor DNA

Circulating tumor-derived DNA (ctDNA) is an emerging biomarker for many cancers, but the limited sensitivity of current detection methods reduces its utility for diagnosing minimal residual disease. Here we describe phased variant enrichment and detection sequencing (PhasED-seq), a method that uses multiple somatic mutations in individual DNA fragments to improve the sensitivity of ctDNA detection. Leveraging whole-genome sequences from 2,538 tumors, we identify phased variants and their associations with mutational signatures. We show that even without molecular barcodes, the limits of detection of PhasED-seq outperform prior methods, including duplex barcoding, allowing ctDNA detection in the ppm range in participant samples. We profiled 678 specimens from 213 participants with B cell lymphomas, including serial cell-free DNA samples before and during therapy for diffuse large B cell lymphoma. In participants with undetectable ctDNA after two cycles of therapy using a next-generation sequencing-based approach termed cancer personalized profiling by deep sequencing, an additional 25% have ctDNA detectable by PhasED-seq and have worse outcomes. Finally, we demonstrate the application of PhasED-seq to solid tumors.

Plasmodium falciparum Multidrug Resistance Proteins (pfMRPs)

The capacity of the lethal Plasmodium falciparum parasite to develop resistance against anti-malarial drugs represents a central challenge in the global control and elimination of malaria. Historically, the action of drug transporters is known to play a pivotal role in the capacity of the parasite to evade drug action. MRPs (Multidrug Resistance Protein) are known in many phylogenetically diverse groups to be related to drug resistance by being able to handle a large range of substrates, including important endogenous substances as glutathione and its conjugates. P. falciparum MRPs are associated with in vivo and in vitro altered drug response, and might be important factors for the development of multi-drug resistance phenotypes, a latent possibility in the present, and future, combination therapy environment. Information on P. falciparum MRPs is scattered in the literature, with no specialized review available. We herein address this issue by reviewing the present state of knowledge.

Heat Shock Proteins as Targets for Novel Antimalarial Drug Discovery

Plasmodium falciparum, the parasitic agent that is responsible for a severe and dangerous form of human malaria, has a history of long years of cohabitation with human beings with attendant negative consequences. While there have been some gains in the fight against malaria through the application of various control measures and the use of chemotherapeutic agents, and despite the global decline in malaria cases and associated deaths, the continual search for new and effective therapeutic agents is key to achieving sustainable development goals. An important parasite survival strategy, which is also of serious concern to the scientific community, is the rate at which the parasites continually develop resistance to drugs. Among the key players in the parasite's ability to develop resistance, maintain cellular integrity, and survives within an unusual environment of the red blood cells are the molecular chaperones of the heat shock proteins (HSP) family. HSPs constitute a novel avenue for antimalarial drug discovery and by exploring their ubiquitous nature and multifunctional activities, they may be suitable targets for the discovery of multi-targets antimalarial drugs, needed to fight incessant drug resistance. In this chapter, features of selected families of plasmodial HSPs that can be exploited in drug discovery are presented. Also, known applications of HSPs in small molecule screening, their potential usefulness in high throughput drug screening, as well as possible challenges are highlighted.

Ring Distortion of Vincamine Leads to the Identification of Re-Engineered Antiplasmodial Agents

There is a significant need for new agents to combat malaria, which resulted in ∼409,000 deaths globally in 2019. We utilized a ring distortion strategy to create complex and diverse compounds from vincamine with the goal of discovering molecules with re-engineered biological activities. We found compound 8 (V3b) to target chloroquine-resistant Plasmodium falciparum Dd2 parasites (EC₅₀ = 1.81 ± 0.09 μM against Dd2 parasites; EC₅₀ > 40 μM against HepG2 cells) and established structure-activity relationships for 25 related analogues. New analogue 30 (V3ss, Dd2, EC₅₀ = 0.25 ± 0.004 μM; HepG2, EC₅₀ > 25 μM) was found to demonstrate the most potent activity, which prevents exit on the parasite from the schizont stage of intraerythrocytic development and requires >24 h to kill P. falciparum Dd2 cells. These findings demonstrate the potential that vincamine ring distortion has toward the discovery of novel antimalarial agents and other therapies significant to human health.

Identification of Macrophage Polarization-Related Genes as Biomarkers of Chronic Obstructive Pulmonary Disease Based on Bioinformatics Analyses

OBJECTIVES

Chronic obstructive pulmonary disease (COPD) is characterized by lung inflammation and remodeling. Macrophage polarization is associated with inflammation and tissue remodeling, as well as immunity. Therefore, this study attempts to investigate the diagnostic value and regulatory mechanism of macrophage polarization-related genes for COPD by bioinformatics analysis and to provide a new theoretical basis for experimental research.

METHODS

The raw gene expression profile dataset (GSE124180) was collected from the Gene Expression Omnibus (GEO) database. Next, a weighted gene coexpression network analysis (WGCNA) was conducted to screen macrophage polarization-related genes. The differentially expressed genes (DEGs) between the COPD and normal samples were generated using DESeq2 v3.11 and overlapped with the macrophage polarization-related genes. Moreover, functional annotations of overlapped genes were conducted by Database for Annotation, Visualization and Integrated Discovery (DAVID) Bioinformatics Resource. The immune-related genes were selected, and their correlation with the differential immune cells was analyzed by Pearson. Finally, receiver operating characteristic (ROC) curves were used to verify the diagnostic value of genes.

RESULTS

A total of 4922 coexpressed genes related to macrophage polarization were overlapped with the 203 DEGs between the COPD and normal samples, obtaining 25 genes related to COPD and macrophage polarization. GEM, S100B, and GZMA of them participated in the immune response, which were considered the candidate biomarkers. GEM and S100B were significantly correlated with marker genes of B cells which had a significant difference between the COPD and normal samples. Moreover, GEM was highly associated with the genes in the PI3K/Akt/GSK3β signaling pathway, regulation of actin cytoskeleton, and calcium signaling pathway based on a Pearson correlation analysis of the candidate genes and the genes in the B cell receptor signaling pathway. PPI network analysis also indicated that GEM might participate in the regulation of the PI3K/Akt/GSK3β signaling pathway. The ROC curve showed that GEM possessed an excellent accuracy in distinguishing COPD from normal samples.

CONCLUSIONS

The data provide a transcriptome-based evidence that GEM is related to COPD and macrophage polarization likely contributes to COPD diagnosis. At the same time, it is hoped that in-depth functional mining can provide new ideas for exploring the COPD pathogenesis.

Comprehensive Analysis of Transcript and Protein Relative Abundance During Blood Stages of Plasmodium falciparum Infection

Plasmodium falciparum is the main causative agent of human malaria. During the intraerythrocytic development cycle, the P. falciparum morphology changes dramatically from circulating young rings to sequestered mature trophozoites and schizonts. Sequestered forms contribute to the pathophysiology of severe malaria as the infected erythrocytes obstruct the microvascular flow in deep organs and induce local inflammation. However, the sequestration mechanism limits the access to the corresponding parasitic form in the clinical samples from patients infected with P. falciparum. To complement this deficiency, we aimed to evaluate the relevance of mRNA study as a proxy of protein expression in sequestered parasites. To do so, we conducted a proteotranscriptomic analysis using five independent P. falciparum laboratory strain samples. RNA sequencing was performed, and the mRNA expression level was assessed on circulating ring-stage parasites. The level of protein expression were measured by LC-MS/MS on the corresponding sequestered mature forms after 18-24 h of maturation. Overall, our results showed a strong transcriptome/transcriptome and a very strong proteome/proteome correlation between samples. Moreover, positive correlations of mRNA and protein expression levels were found between ring-stage transcriptomes and mature form proteomes. However, twice more transcripts were identified at the ring stage than proteins at the mature trophozoite stage. A high level of transcript expression did not guarantee the detection of the corresponding protein. Finally, we pointed out discrepancies at the individual gene level. Taken together, our results show that transcript and protein expressions are overall correlated. However, mRNA abundance is not a perfect proxy of protein expression at the individual level. Importantly, our study shows limitations of the "blind" use of RNA-seq and the importance of multiomics approaches for P. falciparum blood stage study in clinical samples.

An integrated platform for genome engineering and gene expression perturbation in Plasmodium falciparum

Establishing robust genome engineering methods in the malarial parasite, Plasmodium falciparum, has the potential to substantially improve the efficiency with which we gain understanding of this pathogen's biology to propel treatment and elimination efforts. Methods for manipulating gene expression and engineering the P. falciparum genome have been validated. However, a significant barrier to fully leveraging these advances is the difficulty associated with assembling the extremely high AT content DNA constructs required for modifying the P. falciparum genome. These are frequently unstable in commonly-used circular plasmids. We address this bottleneck by devising a DNA assembly framework leveraging the improved reliability with which large AT-rich regions can be efficiently manipulated in linear plasmids. This framework integrates several key functional genetics outcomes via CRISPR/Cas9 and other methods from a common, validated framework. Overall, this molecular toolkit enables P. falciparum genetics broadly and facilitates deeper interrogation of parasite genes involved in diverse biological processes.

The Mediator Subunit, Med23 Is Required for Embryonic Survival and Regulation of Canonical WNT Signaling During Cranial Ganglia Development

Development of the vertebrate head is a complex and dynamic process, which requires integration of all three germ layers and their derivatives. Of special importance are ectoderm-derived cells that form the cranial placodes, which then differentiate into the cranial ganglia and sensory organs. Critical to a fully functioning head, defects in cranial placode and sensory organ development can result in congenital craniofacial anomalies. In a forward genetic screen aimed at identifying novel regulators of craniofacial development, we discovered an embryonically lethal mouse mutant, snouty, which exhibits malformation of the facial prominences, cranial nerves and vasculature. The snouty mutation was mapped to a single nucleotide change in a ubiquitously expressed gene, Med23, which encodes a subunit of the global transcription co-factor complex, Mediator. Phenotypic analyses revealed that the craniofacial anomalies, particularly of the cranial ganglia, were caused by a failure in the proper specification of cranial placode neuronal precursors. Molecular analyses determined that defects in cranial placode neuronal differentiation in Med23 ^sn/sn mutants were associated with elevated WNT/β-catenin signaling, which can be partially rescued through combined Lrp6 and Wise loss-of-function. Our work therefore reveals a surprisingly tissue specific role for the ubiquitously expressed mediator complex protein Med23 in placode differentiation during cranial ganglia development. This highlights the importance of coupling general transcription to the regulation of WNT signaling during embryogenesis.

N6-Adenosine methylation on mRNA is recognized by YTH2 domain protein of human malaria parasite Plasmodium falciparum

Background

Plasmodium falciparum exhibits high translational plasticity during its development in RBCs, yet the regulation at the post-transcriptional level is not well understood. The N⁶-methyl adenosine (m6A) is an important epigenetic modification primarily present on mRNA that controls the levels of transcripts and efficiency of translation in eukaryotes. Recently, the dynamics of m6A on mRNAs at all three developmental stages of P. falciparum in RBCs have been profiled; however, the proteins that regulate the m6A containing mRNAs in the parasites are unknown.

Results

Using sequence analysis, we computationally identified that the P. falciparum genome encodes two putative YTH (YT521-B Homology) domain-containing proteins, which could potentially bind to m6A containing mRNA. We developed a modified methylated RNA immunoprecipitation (MeRIP) assay using PfYTH2 and find that it binds selectively to m6A containing transcripts. The PfYTH2 has a conserved aromatic amino acid cage that forms the methyl-binding pocket. Through site-directed mutagenesis experiments and molecular dynamics simulations, we show that F98 residue is important for m6A binding on mRNA. Fluorescence depolarization assay confirmed that PfYTH2 binds to methylated RNA oligos with high affinity. Further, MeRIP sequencing data revealed that PfYTH2 has more permissive sequence specificity on target m6A containing mRNA than other known eukaryotic YTH proteins. Taken together, here we identify and characterize PfYTH2 as the major protein that could regulate m6A containing transcripts in P. falciparum.

Conclusion

Plasmodium spp. lost the canonical m6A-specific demethylases in their genomes, however, the YTH domain-containing proteins seem to be retained. This study presents a possibility that the YTH proteins are involved in post-transcriptional control in P. falciparum, and might orchestrate the translation of mRNA in various developmental stages of P. falciparum. This is perhaps the first characterization of the methyl-reading function of YTH protein in any parasites.

Rapid activation of distinct members of multigene families in Plasmodium spp

The genomes of Plasmodium spp. encode a number of different multigene families that are thought to play a critical role for survival. However, with the exception of the P. falciparum var genes, very little is known about the biological roles of any of the other multigene families. Using the recently developed Selection Linked Integration method, we have been able to activate the expression of a single member of a multigene family of our choice in Plasmodium spp. from its endogenous promoter. We demonstrate the usefulness of this approach by activating the expression of a unique var, rifin and stevor in P. falciparum as well as yir in P. yoelii. Characterization of the selected parasites reveals differences between the different families in terms of mutual exclusive control, co-regulation, and host adaptation. Our results further support the application of the approach for the study of multigene families in Plasmodium and other organisms.

Omelianczyk, Loh et al. activate the expression of a single member of a multigene family in Plasmodium spp. from its endogenous promoter, identifying differences between the different families. This study supports the application of the Selection Linked Integration method for studying multigene families in Plasmodium.

Re-Engineering of Yohimbine's Biological Activity through Ring Distortion: Identification and Structure-Activity Relationships of a New Class of Antiplasmodial Agents

Select natural products are ideal starting points for ring distortion, or the dramatic altering of inherently complex molecules through short synthetic pathways, to generate an array of novel compounds with diverse skeletal architectures. A major goal of our ring distortion approach is to re-engineer the biological activity of indole alkaloids to identify new compounds with diverse biological activities in areas of significance to human health and medicine. In this study, we re-engineered the biological activity of the indole alkaloid yohimbine through ring rearrangement and ring cleavage synthesis pathways to discover new series of antiplasmodial agents. One new compound, Y7j, was found to demonstrate good potency against chloroquine-resistant Plasmodium falciparum Dd2 cells (EC₅₀ = 0.33 μM) without eliciting cytotoxicity against HepG2 cells (EC₅₀ > 40 μM). Y7j demonstrated stage-specific action against parasites at the late ring/trophozoite stage. A series of analogues was synthesized to gain structure-activity relationship insights, and we learned that both benzyl groups of Y7j are required for activity and fine-tuning of antiplasmodial activities could be accomplished by changing substitution patterns on the benzyl moieties. This study demonstrates the potential for ring distortion to drive new discoveries and change paradigms in chemical biology and drug discovery.

Nucleic Acid Platform Technologies

Every microarray experiment is based on a common format. First, a large number of nucleotide "spots" are arrayed onto a substrate, typically a glass slide, a silicon chip, or microbeads. Second, a complex population of nucleic acids (isolated from cells, selected from in vitro-synthesized libraries, or obtained from another source) is labeled, typically with fluorescent dyes. Third, the labeled nucleic acids are allowed to hybridize to their complementary spot(s) on the microarray. Fourth, the hybridized microarray is washed, allowing the amount of hybridized label to then be quantified. Analysis of the raw data generates a readout of the levels of each species of RNA in the original complex population. This introduction includes several examples of microarray applications and provides a discussion of the basic steps of most microarray experiments.

Partners in Mischief: Functional Networks of Heat Shock Proteins of Plasmodium falciparum and Their Influence on Parasite Virulence

The survival of the human malaria parasite Plasmodium falciparum under the physiologically distinct environments associated with their development in the cold-blooded invertebrate mosquito vectors and the warm-blooded vertebrate human host requires a genome that caters to adaptability. To this end, a robust stress response system coupled to an efficient protein quality control system are essential features of the parasite. Heat shock proteins constitute the main molecular chaperone system of the cell, accounting for approximately two percent of the malaria genome. Some heat shock proteins of parasites constitute a large part (5%) of the 'exportome' (parasite proteins that are exported to the infected host erythrocyte) that modify the host cell, promoting its cyto-adherence. In light of their importance in protein folding and refolding, and thus the survival of the parasite, heat shock proteins of P. falciparum have been a major subject of study. Emerging evidence points to their role not only being cyto-protection of the parasite, as they are also implicated in regulating parasite virulence. In undertaking their roles, heat shock proteins operate in networks that involve not only partners of parasite origin, but also potentially functionally associate with human proteins to facilitate parasite survival and pathogenicity. This review seeks to highlight these interplays and their roles in parasite pathogenicity. We further discuss the prospects of targeting the parasite heat shock protein network towards the developments of alternative antimalarial chemotherapies.

Parasite-Related Genetic and Epigenetic Aspects and Host Factors Influencing Plasmodium falciparum Invasion of Erythrocytes

Malaria, a disease caused by Plasmodium parasites, is widespread throughout tropical and sub-tropical regions worldwide; it mostly affects children and pregnant woman. Eradication has stalled despite effective prevention measures and medication being available for this disease; this has mainly been due to the parasite's resistance to medical treatment and the mosquito vector's resistance to insecticides. Tackling such resistance involves using renewed approaches and techniques for accruing a deep understanding of the parasite's biology, and developing new drugs and vaccines. Studying the parasite's invasion of erythrocytes should shed light on its ability to switch between invasion phenotypes related to the expression of gene sets encoding proteins acting as ligands during target cell invasion, thereby conferring mechanisms for evading a particular host's immune response and adapting to changes in target cell surface receptors. This review considers some factors influencing the expression of such phenotypes, such as Plasmodium's genetic, transcriptional and epigenetic characteristics, and explores some host-related aspects which could affect parasite phenotypes, aiming at integrating knowledge regarding this topic and the possible relationship between the parasite's biology and host factors playing a role in erythrocyte invasion.

Candidate Biomarkers and Molecular Mechanism Investigation for Glioblastoma Multiforme Utilizing WGCNA

To reveal the potential molecular mechanism of glioblastoma multiforme (GBM) and provide the candidate biomarkers for GBM gene therapy. Microarray dataset GSE50161 was obtained from GEO database. The differentially expressed genes (DEGs) were identified between GBM samples and control samples, followed by the module partition analysis based on WGCNA. Then, the pathway and functional enrichment analyses of DEGs were performed. The hub genes were further investigated, followed by the survival analysis and data validation. A total of 1913 DEGs were investigated between two groups, followed by analysis of 5 modules using WGCNA. These DEGs were mainly enriched in functions like inflammatory response. The hub genes including upregulated N-Myc and STAT Interactor (NMI), Capping Actin Protein-Gelsolin Like (CAPG), and Proteasome Subunit Beta 8 (PSMB8) were revealed as potential liquid biopsy molecules for GBM diagnose. Moreover, Nucleolar and Spindle Associated Protein 1 (NUSAP1) and G Protein-Coupled Receptor 65 (GPR65) were outstanding genes in survival analysis. Our results suggested that CPNE6, HAPLN2, CMTM3, NMI, CAPG, and PSMB8 might be used as potential molecules for liquid biopsy of GBM. NUSAP1 and GPR65 might be novel prognostic targets for GBM gene therapy. Furthermore, the upregulated NMI might play an important role in GBM progression via inflammatory response.

Gene copy number variation in natural populations of Plasmodium falciparum in Eastern Africa

BACKGROUND

Gene copy number variants (CNVs), which consist of deletions and amplifications of single or sets of contiguous genes, contribute to the great diversity in the Plasmodium falciparum genome. In vitro studies in the laboratory have revealed their important role in parasite fitness phenotypes such as red cell invasion, transmissibility and cytoadherence. Studies of natural parasite populations indicate that CNVs are also common in the field and thus may facilitate adaptation of the parasite to its local environment.

RESULTS

In a survey of 183 fresh field isolates from three populations in Eastern Africa with different malaria transmission intensities, we identified 94 CNV loci using microarrays. All CNVs had low population frequencies (minor allele frequency < 5%) but each parasite isolate carried an average of 8 CNVs. Nine CNVs showed high levels of population differentiation (F_ST > 0.3) and nine exhibited significant clines in population frequency across a gradient in transmission intensity. The clearest example of this was a large deletion on chromosome 9 previously reported only in laboratory-adapted isolates. This deletion was present in 33% of isolates from a population with low and highly seasonal malaria transmission, and in < 9% of isolates from populations with higher transmission. Subsets of CNVs were strongly correlated in their population frequencies, implying co-selection.

CONCLUSIONS

These results support the hypothesis that CNVs are the target of selection in natural populations of P. falciparum. Their environment-specific patterns observed here imply an important role for them in conferring adaptability to the parasite thus enabling it to persist in its highly diverse ecological environment.

Plasmodium falciparum strains spontaneously switch invasion phenotype in suspension culture

The extensive redundancy in the use of invasion ligands by Plasmodium falciparum, and its unique ability to switch between invasion pathways have hampered vaccine development. P. falciparum strains Dd2 and W2mef have been shown to change from sialic acid (SA)-dependent to SA-independent phenotypes when selected on neuraminidase-treated erythrocytes. Following an observation of increasing ability of Dd2 to invade neuraminidase-treated cells when cultured for several weeks, we systematically investigated this phenomenon by comparing invasion phenotypes of Dd2, W2mef and 3D7 strains of P. falciparum that were cultured with gentle shaking (Suspended) or under static (Static) conditions. While Static Dd2 and W2mef remained SA-dependent for the entire duration of the investigation, Suspended parasites spontaneously and progressively switched to SA-independent phenotype from week 2 onwards. Furthermore, returning Suspended cultures to Static conditions led to a gradual reversal to SA-dependent phenotype. The switch to SA-independent phenotype was accompanied by upregulation of the key invasion ligand, reticulocyte-binding homologue 4 (RH4), and the increased invasion was inhibited by antibodies to the RH4 receptor, CR1. Our data demonstrates a novel mechanism for inducing the switching of invasion pathways in P. falciparum parasites and may provide clues for understanding the mechanisms involved.

Single-cell RNA-seq reveals hidden transcriptional variation in malaria parasites

Single-cell RNA-sequencing is revolutionising our understanding of seemingly homogeneous cell populations but has not yet been widely applied to single-celled organisms. Transcriptional variation in unicellular malaria parasites from the Plasmodium genus is associated with critical phenotypes including red blood cell invasion and immune evasion, yet transcriptional variation at an individual parasite level has not been examined in depth. Here, we describe the adaptation of a single-cell RNA-sequencing (scRNA-seq) protocol to deconvolute transcriptional variation for more than 500 individual parasites of both rodent and human malaria comprising asexual and sexual life-cycle stages. We uncover previously hidden discrete transcriptional signatures during the pathogenic part of the life cycle, suggesting that expression over development is not as continuous as commonly thought. In transmission stages, we find novel, sex-specific roles for differential expression of contingency gene families that are usually associated with immune evasion and pathogenesis.

Oxidative stress and protein damage responses mediate artemisinin resistance in malaria parasites

Due to their remarkable parasitocidal activity, artemisinins represent the key components of first-line therapies against Plasmodium falciparum malaria. However, the decline in efficacy of artemisinin-based drugs jeopardizes global efforts to control and ultimately eradicate the disease. To better understand the resistance phenotype, artemisinin-resistant parasite lines were derived from two clones of the 3D7 strain of P. falciparum using a selection regimen that mimics how parasites interact with the drug within patients. This long term in vitro selection induced profound stage-specific resistance to artemisinin and its relative compounds. Chemosensitivity and transcriptional profiling of artemisinin-resistant parasites indicate that enhanced adaptive responses against oxidative stress and protein damage are associated with decreased artemisinin susceptibility. This corroborates our previous findings implicating these cellular functions in artemisinin resistance in natural infections. Genomic characterization of the two derived parasite lines revealed a spectrum of sequence and copy number polymorphisms that could play a role in regulating artemisinin response, but did not include mutations in pfk13, the main marker of artemisinin resistance in Southeast Asia. Taken together, here we present a functional in vitro model of artemisinin resistance that is underlined by a new set of genetic polymorphisms as potential genetic markers.

Malaria Pathogenesis

In the mosquito-human life cycle, the six species of malaria parasites infecting humans (Plasmodium falciparum, Plasmodium vivax, Plasmodium ovale wallickeri, Plasmodium ovale curtisi, Plasmodium malariae, and Plasmodium knowlesi) undergo 10 or more morphological states, replicate from single to 10,000+ cells, and vary in total population from one to many more than 10⁶ organisms. In the human host, only a small number of these morphological stages lead to clinical disease and the vast majority of all malaria-infected patients in the world produce few (if any) symptoms in the human. Human clinical disease (e.g., fever, anemia, coma) is the result of the parasite preprogrammed biology in concert with the human pathophysiological response. Caveats and corollaries that add variation to this host-parasite interaction include parasite genetic diversity of key proteins, coinfections, comorbidities, delays in treatment, human polymorphisms, and environmental determinants.

Comparison of mitochondrial gene expression and polysome loading in different tobacco tissues

BACKGROUND

To investigate translational regulation of gene expression in plant mitochondria, a mitochondrial polysome isolation protocol was established for tobacco to investigate polysomal mRNA loading as a proxy for translational activity. Furthermore, we developed an oligonucleotide based microarray platform to determine the level of Nicotiana tabacum and Arabidopsis thaliana mitochondrial mRNA.

RESULTS

Microarray analysis of free and polysomal mRNAs was used to characterize differences in the levels of free transcripts and ribosome-bound mRNAs in various organs of tobacco plants. We have observed higher mitochondrial transcript levels in young leaves, flowers and floral buds as compared to fully expanded leaves and roots. A similar pattern of abundance was observed for ribosome-bound mitochondrial mRNAs in these tissues. However, the accumulation of the mitochondrial protein COX2 was found to be inversely related to that of its ribosome-bound mRNA.

CONCLUSIONS

Our results indicate that the association of mitochondrial mRNAs to ribosomes is largely determined by the total transcript level of a gene. However, at least for Cox2, we demonstrated that the level of ribosome-bound mRNA is not reflected by the amount of COX2 protein.

Conserved Binding Regions Provide the Clue for Peptide-Based Vaccine Development: A Chemical Perspective

Synthetic peptides have become invaluable biomedical research and medicinal chemistry tools for studying functional roles, i.e., binding or proteolytic activity, naturally-occurring regions’ immunogenicity in proteins and developing therapeutic agents and vaccines. Synthetic peptides can mimic protein sites; their structure and function can be easily modulated by specific amino acid replacement. They have major advantages, i.e., they are cheap, easily-produced and chemically stable, lack infectious and secondary adverse reactions and can induce immune responses via T- and B-cell epitopes. Our group has previously shown that using synthetic peptides and adopting a functional approach has led to identifying Plasmodium falciparumconserved regions binding to host cells. Conserved high activity binding peptides’ (cHABPs) physicochemical, structural and immunological characteristics have been taken into account for properly modifying and converting them into highly immunogenic, protection-inducing peptides (mHABPs) in the experimental Aotus monkey model. This article describes stereo–electron and topochemical characteristics regarding major histocompatibility complex (MHC)-mHABP-T-cell receptor (TCR) complex formation. Some mHABPs in this complex inducing long-lasting, protective immunity have been named immune protection-inducing protein structures (IMPIPS), forming the subunit components in chemically synthesized vaccines. This manuscript summarizes this particular field and adds our recent findings concerning intramolecular interactions (H-bonds or π-interactions) enabling proper IMPIPS structure as well as the peripheral flanking residues (PFR) to stabilize the MHCII-IMPIPS-TCR interaction, aimed at inducing long-lasting, protective immunological memory.

New Insights into the Genetic Regulation of Plasmodium Falciparum Obtained by Bayesian Modeling

The most fatal and prevalent form of malaria is caused by the bloodborne pathogen Plasmodium falciparum (henceforth P.f). Annually, approximately three million people died of malaria. Despite P.f devastivating effect globally, the vast majority of its proteins have not been characterized experimentally. In this work, we provide computational insight that explore the modalities of the regulation for some important group of genes of P.f, namely components of the glycolytic pathway, and those involved in apicoplast metabolism. Glycolysis is a crucial pathway in the maintenance of the parasite while the recently discovered apicoplast contains a range of metabolic pathways and housekeeping processes that differ radically to those of the host, which makes it ideal for drug therapy.

We have been able to validate some of our findings from available literature and therefore provide a basis to give theoretical insight for some genes regulations, which has not been characterized experimentally.

Simultaneous genome-wide gene expression and transcript isoform profiling in the human malaria parasite

Gene expression DNA microarrays have been vital for characterizing whole-genome transcriptional profiles. Nevertheless, their effectiveness relies heavily on the accuracy of genome sequences, the annotation of gene structures, and the sequence-dependent performance of individual probes. Currently available gene expression arrays for the malaria parasite Plasmodium falciparum rely on an average of 2 probes per gene, usually positioned near the 3′ end of genes; consequently, existing designs are prone to measurement bias and cannot capture complexities such as the occurrence of transcript isoforms arising from alternative splicing or alternative start/ stop sites. Here, we describe two novel gene expression arrays with exon-focused probes designed with an average of 12 and 20 probes spanning each gene. This high probe density minimizes signal noise inherent in probe-to-probe sequence-dependent hybridization intensity. We demonstrate that these exon arrays accurately profile genome-wide expression, comparing favorably to currently available arrays and RNA-seq profiling, and can detect alternatively spliced transcript isoforms as well as non-coding RNAs (ncRNAs). Of the 964 candidate alternate splicing events from published RNA-seq studies, 162 are confirmed using the exon array. Furthermore, the exon array predicted 330 previously unidentified alternate splicing events. Gene expression microarrays continue to offer a cost-effective alternative to RNA-seq for the simultaneous monitoring of gene expression and alternative splicing events. Microarrays may even be preferred in some cases due to their affordability and the rapid turn-around of results when hundreds of samples are required for fine-scale systems biology investigations, including the monitoring of the networks of gene co-expression in the emergence of drug resistance.

Tube-Wise Diagnostic Microarray for the Multiplex Characterization of the Complex Plant Pathogen Ralstonia solanacearum

Ralstonia solanacearum is a well-known agricultural and ecological threat worldwide. The complexity of the R. solanacearum species complex (Rssc) represents a challenge for the accurate characterization of epidemiological strains by official services and research laboratories. The majority of protocols only focus on a narrow range of strains; however, this species complex includes strains that represent major constraints and are under strict regulation. The main drawback associated with the current methods of detecting and characterizing Rssc strains is their reliance on combining different protocols to properly characterize the strains at the ecotype level, which require time and money. Therefore, we used microarray technology (ArrayTube) to develop a standard protocol, which characterizes 17 major groups of interest in the Rssc, in a single multiplex reaction. These 17 majors groups are linked with a phylogenetic assignation (phylotypes, sequevars), but also with an ecotype assignation associated with a range of hosts (e.g., brown rot, Moko). Probes were designed with a 50-mer length constraint and thoroughly evaluated for any flaws or secondary structures. The strains are characterized based on a DNA extraction from pure culture. Validation data showed strong intra-repeatability, inter-repeatability, and reproducibility as well as good specificity. A hierarchical analysis of the probe groups is suitable for an accurate characterization. Compared with single marker detection tests, the method described in this paper addresses efficiently the issue of combining several tests by testing a large number of phylogenetic markers in a single reaction assay. This custom microarray (RsscAT) represents a significant improvement in the epidemiological monitoring of Rssc strains worldwide, and it has the potential to provide insights for phylogenetic incongruence of Rssc strains based on the host of isolation and may be used to indicate potentially emergent strains.

Most of the tight positional conservation of transcription factor binding sites near the transcription start site reflects their co-localization within regulatory modules

Background

Transcription factors (TFs) form complexes that bind regulatory modules (RMs) within DNA, to control specific sets of genes. Some transcription factor binding sites (TFBSs) near the transcription start site (TSS) display tight positional preferences relative to the TSS. Furthermore, near the TSS, RMs can co-localize TFBSs with each other and the TSS. The proportion of TFBS positional preferences due to TFBS co-localization within RMs is unknown, however. ChIP experiments confirm co-localization of some TFBSs genome-wide, including near the TSS, but they typically examine only a few TFs at a time, using non-physiological conditions that can vary from lab to lab. In contrast, sequence analysis can examine many TFs uniformly and methodically, broadly surveying the co-localization of TFBSs with tight positional preferences relative to the TSS.

Results

Our statistics found 43 significant sets of human motifs in the JASPAR TF Database with positional preferences relative to the TSS, with 38 preferences tight (±5 bp). Each set of motifs corresponded to a gene group of 135 to 3304 genes, with 42/43 (98%) gene groups independently validated by DAVID, a gene ontology database, with FDR < 0.05. Motifs corresponding to two TFBSs in a RM should co-occur more than by chance alone, enriching the intersection of the gene groups corresponding to the two TFs. Thus, a gene-group intersection systematically enriched beyond chance alone provides evidence that the two TFs participate in an RM. Of the 903 = 43*42/2 intersections of the 43 significant gene groups, we found 768/903 (85%) pairs of gene groups with significantly enriched intersections, with 564/768 (73%) intersections independently validated by DAVID with FDR < 0.05. A user-friendly web site at http://go.usa.gov/3kjsH permits biologists to explore the interaction network of our TFBSs to identify candidate subunit RMs.

Conclusions

Gene duplication and convergent evolution within a genome provide obvious biological mechanisms for replicating an RM near the TSS that binds a particular TF subunit. Of all intersections of our 43 significant gene groups, 85% were significantly enriched, with 73% of the significant enrichments independently validated by gene ontology. The co-localization of TFBSs within RMs therefore likely explains much of the tight TFBS positional preferences near the TSS.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1354-5) contains supplementary material, which is available to authorized users.

Experimental conditions improving in-solution target enrichment for ancient DNA

High-throughput sequencing has dramatically fostered ancient DNA research in recent years. Shotgun sequencing, however, does not necessarily appear as the best-suited approach due to the extensive contamination of samples with exogenous environmental microbial DNA. DNA capture-enrichment methods represent cost-effective alternatives that increase the sequencing focus on the endogenous fraction, whether it is from mitochondrial or nuclear genomes, or parts thereof. Here, we explored experimental parameters that could impact the efficacy of MYbaits in-solution capture assays of ~5000 nuclear loci or the whole genome. We found that varying quantities of the starting probes had only moderate effect on capture outcomes. Starting DNA, probe tiling, the hybridization temperature and the proportion of endogenous DNA all affected the assay, however. Additionally, probe features such as their GC content, number of CpG dinucleotides, sequence complexity and entropy and self-annealing properties need to be carefully addressed during the design stage of the capture assay. The experimental conditions and probe molecular features identified in this study will improve the recovery of genetic information extracted from degraded and ancient remains.

A Plasmodium yoelii Mei2-Like RNA Binding Protein Is Essential for Completion of Liver Stage Schizogony

Plasmodium parasites employ posttranscriptional regulatory mechanisms as their life cycle transitions between host cell invasion and replication within both the mosquito vector and mammalian host. RNA binding proteins (RBPs) provide one mechanism for modulation of RNA function. To explore the role of Plasmodium RBPs during parasite replication, we searched for RBPs that might play a role during liver stage development, the parasite stage that exhibits the most extensive growth and replication. We identified a parasite ortholog of the Mei2 (Meiosis inhibited 2) RBP that is conserved among Plasmodium species (PlasMei2) and exclusively transcribed in liver stage parasites. Epitope-tagged Plasmodium yoelii PlasMei2 was expressed only during liver stage schizogony and showed an apparent granular cytoplasmic location. Knockout of PlasMei2 (plasmei2(-)) in P. yoelii only affected late liver stage development. The P. yoelii plasmei2(-) liver stage size increased progressively until late in development, similar to wild-type parasite development. However, P. yoelii plasmei2(-) liver stage schizonts exhibited an abnormal DNA segregation phenotype and failed to form exoerythrocytic merozoites. Consequently the cellular integrity of P. yoelii plasmei2(-) liver stages became increasingly compromised late in development and the majority of P. yoelii plasmei2(-) underwent cell death by the time wild-type liver stages mature and release merozoites. This resulted in a complete block of P. yoelii plasmei2(-) transition from liver stage to blood stage infection in mice. Our results show for the first time the importance of a Plasmodium RBP in the coordinated progression of late liver stage schizogony and maturation of new invasive forms.

Landscape and Dynamics of Transcription Initiation in the Malaria Parasite Plasmodium falciparum

A comprehensive map of transcription start sites (TSSs) across the highly AT-rich genome of P. falciparum would aid progress toward deciphering the molecular mechanisms that underlie the timely regulation of gene expression in this malaria parasite. Using high-throughput sequencing technologies, we generated a comprehensive atlas of transcription initiation events at single-nucleotide resolution during the parasite intra-erythrocytic developmental cycle. This detailed analysis of TSS usage enabled us to define architectural features of plasmodial promoters. We demonstrate that TSS selection and strength are constrained by local nucleotide composition. Furthermore, we provide evidence for coordinate and stage-specific TSS usage from distinct sites within the same transcription unit, thereby producing transcript isoforms, a subset of which are developmentally regulated. This work offers a framework for further investigations into the interactions between genomic sequences and regulatory factors governing the complex transcriptional program of this major human pathogen.

Probe Design Strategies for Oligonucleotide Microarrays

Oligonucleotide microarrays have been widely used for gene detection and/or quantification of gene expression in various samples ranging from a single organism to a complex microbial assemblage. The success of a microarray experiment, however, strongly relies on the quality of designed probes. Consequently, probe design is of critical importance and therefore multiple parameters should be considered for each probe in order to ensure high specificity, sensitivity, and uniformity as well as potentially quantitative power. Moreover, to assess the complete gene repertoire of complex biological samples such as those studied in the field of microbial ecology, exploratory probe design strategies must be also implemented to target not-yet-described sequences. To design such probes, two algorithms, KASpOD and HiSpOD, have been developed and they are available via two user-friendly web services. Here, we describe the use of this software necessary for the design of highly effective probes especially in the context of microbial oligonucleotide microarrays by taking into account all the crucial parameters.

Genome-wide transcriptome profiling reveals functional networks involving the Plasmodium falciparum drug resistance transporters PfCRT and PfMDR1

BACKGROUND

The acquisition of multidrug resistance by Plasmodium falciparum underscores the need to understand the underlying molecular mechanisms so as to counter their impact on malaria control. For the many antimalarials whose mode of action relates to inhibition of heme detoxification inside infected erythrocytes, the digestive vacuole transporters PfCRT and PfMDR1 constitute primary resistance determinants.

RESULTS

Using gene expression microarrays over the course of the parasite intra-erythrocytic developmental cycle, we compared the transcriptomic profiles between P. falciparum strains displaying mutant or wild-type pfcrt or varying in pfcrt or pfmdr1 expression levels. To account for differences in the time of sampling, we developed a computational method termed Hypergeometric Analysis of Time Series, which combines Fast Fourier Transform with a modified Gene Set Enrichment Analysis. Our analysis revealed coordinated changes in genes involved in protein catabolism, translation initiation and DNA/RNA metabolism. We also observed differential expression of genes with a role in transport or coding for components of the digestive vacuole. Interestingly, a global comparison of all profiled transcriptomes uncovered a tight correlation between the transcript levels of pfcrt and pfmdr1, extending to dozens of other genes, suggesting an intricate regulatory balance in order to maintain optimal physiological processes.

CONCLUSIONS

This study provides insight into the mechanisms by which P. falciparum adjusts to the acquisition of mutations or gene amplification in key transporter loci that mediate drug resistance. Our results implicate several biological pathways that may be differentially regulated to compensate for impaired transporter function and alterations in parasite vacuole physiology.

The nucleosome landscape of Plasmodium falciparum reveals chromatin architecture and dynamics of regulatory sequences

In eukaryotes, the chromatin architecture has a pivotal role in regulating all DNA-associated processes and it is central to the control of gene expression. For Plasmodium falciparum, a causative agent of human malaria, the nucleosome positioning profile of regulatory regions deserves particular attention because of their extreme AT-content. With the aid of a highly controlled MNase-seq procedure we reveal how positioning of nucleosomes provides a structural and regulatory framework to the transcriptional unit by demarcating landmark sites (transcription/translation start and end sites). In addition, our analysis provides strong indications for the function of positioned nucleosomes in splice site recognition. Transcription start sites (TSSs) are bordered by a small nucleosome-depleted region, but lack the stereotypic downstream nucleosome arrays, highlighting a key difference in chromatin organization compared to model organisms. Furthermore, we observe transcription-coupled eviction of nucleosomes on strong TSSs during intraerythrocytic development and demonstrate that nucleosome positioning and dynamics can be predictive for the functionality of regulatory DNA elements. Collectively, the strong nucleosome positioning over splice sites and surrounding putative transcription factor binding sites highlights the regulatory capacity of the nucleosome landscape in this deadly human pathogen.

The PfAlba1 RNA-binding protein is an important regulator of translational timing in Plasmodium falciparum blood stages

Background

Transcriptome-wide ribosome occupancy studies have suggested that during the intra-erythrocytic lifecycle of Plasmodium falciparum, select mRNAs are post-transcriptionally regulated. A subset of these encodes parasite virulence factors required for invading host erythrocytes, and are currently being developed as vaccine candidates. However, the molecular mechanisms that govern post-transcriptional regulation are currently unknown.

Results

We explore the previously identified DNA/RNA-binding protein PfAlba1, which localizes to multiple foci in the cytoplasm of P. falciparum trophozoites. We establish that PfAlba1 is essential for asexual proliferation, and subsequently investigate parasites overexpressing epitope-tagged PfAlba1 to identify its RNA targets and effects on mRNA homeostasis and translational regulation. Using deep sequencing of affinity-purified PfAlba1-associated RNAs, we identify 1193 transcripts that directly bind to PfAlba1 in trophozoites. For 105 such transcripts, 43 % of which are uncharacterized and 13 % of which encode erythrocyte invasion components, the steady state levels significantly change at this stage, evidencing a role for PfAlba1 in maintaining mRNA homeostasis. Additionally, we discover that binding of PfAlba1 to four erythrocyte invasion mRNAs, Rap1, RhopH3, CDPK1, and AMA1, is linked to translation repression in trophozoites whereas release of these mRNAs from a PfAlba1 complex in mature stages correlates with protein synthesis.

Conclusions

We show that PfAlba1 binds to a sub-population of asexual stage mRNAs and fine-tunes the timing of translation. This mode of post-transcriptional regulation may be especially important for P. falciparum erythrocyte invasion components that have to be assembled into apical secretory organelles in a highly time-dependent manner towards the end of the parasite’s asexual lifecycle.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0771-5) contains supplementary material, which is available to authorized users.

Plasmodium falciparum coronin organizes arrays of parallel actin filaments potentially guiding directional motility in invasive malaria parasites

Background

Gliding motility in Plasmodium parasites, the aetiological agents of malaria disease, is mediated by an actomyosin motor anchored in the outer pellicle of the motile cell. Effective motility is dependent on a parasite myosin motor and turnover of dynamic parasite actin filaments. To date, however, the basis for directional motility is not known. Whilst myosin is very likely orientated as a result of its anchorage within the parasite, how actin filaments are orientated to facilitate directional force generation remains unexplained. In addition, recent evidence has questioned the linkage between actin filaments and secreted surface antigens leaving the way by which motor force is transmitted to the extracellular milieu unknown. Malaria parasites possess a markedly reduced repertoire of actin regulators, among which few are predicted to interact with filamentous (F)-actin directly. One of these, PF3D7_1251200, shows strong homology to the coronin family of actin-filament binding proteins, herein referred to as PfCoronin.

Methods

Here the N terminal beta propeller domain of PfCoronin (PfCor-N) was expressed to assess its ability to bind and bundle pre-formed actin filaments by sedimentation assay, total internal reflection fluorescence (TIRF) microscopy and confocal imaging as well as to explore its ability to bind phospholipids. In parallel a tagged PfCoronin line in Plasmodium falciparum was generated to determine the cellular localization of the protein during asexual parasite development and blood-stage merozoite invasion.

Results

A combination of biochemical approaches demonstrated that the N-terminal beta-propeller domain of PfCoronin is capable of binding F-actin and facilitating formation of parallel filament bundles. In parasites, PfCoronin is expressed late in the asexual lifecycle and localizes to the pellicle region of invasive merozoites before and during erythrocyte entry. PfCoronin also associates strongly with membranes within the cell, likely mediated by interactions with phosphatidylinositol-4,5-bisphosphate (PI(4,5)P₂) at the plasma membrane.

Conclusions

These data suggest PfCoronin may fulfil a key role as the critical determinant of actin filament organization in the Plasmodium cell. This raises the possibility that macro-molecular organization of actin mediates directional motility in gliding parasites.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-015-0801-5) contains supplementary material, which is available to authorized users.

The effects of cell asynchrony on gene expression levels: analysis and application to Plasmodium falciparum

To investigate the intraerythrocytic developmental cycle of Plasmodium falciparum, time-series gene expression data is commonly measured of infected red blood cells. However, the observed data are usually blurred due to cell asynchrony during experiments. In this paper, the effect of cell asynchrony is investigated by conducting numerical experiments. The simulation results suggest that cell asynchrony has varying effects on different intrinsic expression patterns. Specifically, the intrinsic patterns with high expression around the late life stage are more likely to be affected by cell asynchrony. It is also investigated how the effect of cell asynchrony depends on the experimental conditions. From this analysis, the burst rate r% in infection period and the standard deviation σ of growth rate are identified to have a strong impact on the blurring due to cell asynchrony. Consequently, it is important to measure these two parameters during biological experiments in order to deblur time-series gene expression data.

Genome-Wide Gene Expression in relation to Age in Large Laboratory Cohorts of Drosophila melanogaster

Aging is a complex process characterized by a steady decline in an organism's ability to perform life-sustaining tasks. In the present study, two cages of approximately 12,000 mated Drosophila melanogaster females were used as a source of RNA from individuals sampled frequently as a function of age. A linear model for microarray data method was used for the microarray analysis to adjust for the box effect; it identified 1,581 candidate aging genes. Cluster analyses using a self-organizing map algorithm on the 1,581 significant genes identified gene expression patterns across different ages. Genes involved in immune system function and regulation, chorion assembly and function, and metabolism were all significantly differentially expressed as a function of age. The temporal pattern of data indicated that gene expression related to aging is affected relatively early in life span. In addition, the temporal variance in gene expression in immune function genes was compared to a random set of genes. There was an increase in the variance of gene expression within each cohort, which was not observed in the set of random genes. This observation is compatible with the hypothesis that D. melanogaster immune function genes lose control of gene expression as flies age.

Mitochondrial ATP synthase is dispensable in blood-stage Plasmodium berghei rodent malaria but essential in the mosquito phase

Mitochondrial ATP synthase is driven by chemiosmotic oxidation of pyruvate derived from glycolysis. Blood-stage malaria parasites eschew chemiosmosis, instead relying almost solely on glycolysis for their ATP generation, which begs the question of whether mitochondrial ATP synthase is necessary during the blood stage of the parasite life cycle. We knocked out the mitochondrial ATP synthase β subunit gene in the rodent malaria parasite, Plasmodium berghei, ablating the protein that converts ADP to ATP. Disruption of the β subunit gene of the ATP synthase only marginally reduced asexual blood-stage parasite growth but completely blocked mouse-to-mouse transmission via Anopheles stephensi mosquitoes. Parasites lacking the β subunit gene of the ATP synthase generated viable gametes that fuse and form ookinetes but cannot progress beyond this stage. Ookinetes lacking the β subunit gene of the ATP synthase had normal motility but were not viable in the mosquito midgut and never made oocysts or sporozoites, thereby abrogating transmission to naive mice via mosquito bite. We crossed the self-infertile ATP synthase β subunit knockout parasites with a male-deficient, self-infertile strain of P. berghei, which restored fertility and production of oocysts and sporozoites, which demonstrates that mitochondrial ATP synthase is essential for ongoing viability through the female, mitochondrion-carrying line of sexual reproduction in P. berghei malaria. Perturbation of ATP synthase completely blocks transmission to the mosquito vector and could potentially be targeted for disease control.