Both host and pathogen factors contribute to disease outcome in Plasmodium falciparum infection. The feasibility of studying the P. falciparum in vivo transcriptome to understand parasite transcriptional response while it resides in the human host is presented.

One therapeutic approach to Duchenne Muscular Dystrophy (DMD) recently entering clinical trials aims to convert DMD phenotypes to that of a milder disease variant, Becker Muscular Dystrophy (BMD), by employing antisense oligonucleotides (AONs) targeting splice sites, to induce exon skipping and restore partial dystrophin function. In order to search for small molecule and genetic modulators of AON-dependent and independent exon skipping, we screened approximately 10,000 known small molecule drugs, >17,000 cDNA clones, and >2,000 kinase- targeted siRNAs against a 5.6 kb luciferase minigene construct, encompassing exon 71 to exon 73 of human dystrophin. As a result, we identified several enhancers of exon skipping, acting on both the reporter construct as well as endogenous dystrophin in mdx cells. Multiple mechanisms of action were identified, including histone deacetylase inhibition, tubulin modulation and pre-mRNA processing. Among others, the nucleolar protein NOL8 and staufen RNA binding protein homolog 2 (Stau2) were found to induce endogenous exon skipping in mdx cells in an AON-dependent fashion. An unexpected but recurrent theme observed in our screening efforts was the apparent link between the inhibition of cell cycle progression and the induction of exon skipping.

Influenza A virus is an RNA virus that encodes up to 11 proteins and this small coding capacity demands that the virus use the host cellular machinery for many aspects of its life cycle. Knowledge of these host cell requirements not only informs us of the molecular pathways exploited by the virus but also provides further targets that could be pursued for antiviral drug development. Here we use an integrative systems approach, based on genome-wide RNA interference screening, to identify 295 cellular cofactors required for early-stage influenza virus replication. Within this group, those involved in kinase-regulated signalling, ubiquitination and phosphatase activity are the most highly enriched, and 181 factors assemble into a highly significant host-pathogen interaction network. Moreover, 219 of the 295 factors were confirmed to be required for efficient wild-type influenza virus growth, and further analysis of a subset of genes showed 23 factors necessary for viral entry, including members of the vacuolar ATPase (vATPase) and COPI-protein families, fibroblast growth factor receptor (FGFR) proteins, and glycogen synthase kinase 3 (GSK3)-beta. Furthermore, 10 proteins were confirmed to be involved in post-entry steps of influenza virus replication. These include nuclear import components, proteases, and the calcium/calmodulin-dependent protein kinase (CaM kinase) IIbeta (CAMK2B). Notably, growth of swine-origin H1N1 influenza virus is also dependent on the identified host factors, and we show that small molecule inhibitors of several factors, including vATPase and CAMK2B, antagonize influenza virus replication.

The development of the human malaria parasite Plasmodium falciparum is controlled by coordinated changes in gene expression throughout its complex life cycle, but the corresponding regulatory mechanisms are incompletely understood. To study the relationship between genome architecture and gene regulation in Plasmodium, we assayed the genome architecture of P. falciparum at three time points during its erythrocytic (asexual) cycle. Using chromosome conformation capture coupled with next-generation sequencing technology (Hi-C), we obtained high-resolution chromosomal contact maps, which we then used to construct a consensus three-dimensional genome structure for each time point. We observed strong clustering of centromeres, telomeres, ribosomal DNA, and virulence genes, resulting in a complex architecture that cannot be explained by a simple volume exclusion model. Internal virulence gene clusters exhibit domain-like structures in contact maps, suggesting that they play an important role in the genome architecture. Midway during the erythrocytic cycle, at the highly transcriptionally active trophozoite stage, the genome adopts a more open chromatin structure with increased chromosomal intermingling. In addition, we observed reduced expression of genes located in spatial proximity to the repressive subtelomeric center, and colocalization of distinct groups of parasite-specific genes with coordinated expression profiles. Overall, our results are indicative of a strong association between the P. falciparum spatial genome organization and gene expression. Understanding the molecular processes involved in genome conformation dynamics could contribute to the discovery of novel antimalarial strategies.

Nucleosomes are the basic elements of chromatin structure. They control the packaging of DNA and play a critical role in gene regulation by allowing physical access to transcription factors. The advent of second-generation sequencing has enabled landmark genome-wide studies of nucleosome positions for several model organisms. Current methods to determine nucleosome positioning first compute an occupancy coverage profile by mapping nucleosome-enriched sequenced reads to a reference genome; then, nucleosomes are placed according to the peaks of the coverage profile. These methods are quite accurate on placing isolated nucleosomes, but they do not properly handle more complex configurations. Also, they can only provide the positions of nucleosomes and their occupancy level, whereas it is very beneficial to supply molecular biologists additional information about nucleosomes like the probability of placement, the size of DNA fragments enriched for nucleosomes and/or whether nucleosomes are well positioned or 'fuzzy' in the sequenced cell sample.

Cytosine DNA methylation is an epigenetic mark in most eukaryotic cells that regulates numerous processes, including gene expression and stress responses. We performed a genome-wide analysis of DNA methylation in the human malaria parasite Plasmodium falciparum. We mapped the positions of methylated cytosines and identified a single functional DNA methyltransferase (Plasmodium falciparum DNA methyltransferase; PfDNMT) that may mediate these genomic modifications. These analyses revealed that the malaria genome is asymmetrically methylated and shares common features with undifferentiated plant and mammalian cells. Notably, core promoters are hypomethylated, and transcript levels correlate with intraexonic methylation. Additionally, there are sharp methylation transitions at nucleosome and exon-intron boundaries. These data suggest that DNA methylation could regulate virulence gene expression and transcription elongation. Furthermore, the broad range of action of DNA methylation and the uniqueness of PfDNMT suggest that the methylation pathway is a potential target for antimalarial strategies.

Autophagy is a catabolic pathway typically induced by nutrient starvation to recycle amino acids, but can also function in removing damaged organelles. In addition, this pathway plays a key role in eukaryotic development. To date, not much is known about the role of autophagy in apicomplexan parasites and more specifically in the human malaria parasite Plasmodium falciparum. Comparative genomic analysis has uncovered some, but not all, orthologs of autophagy-related (ATG) genes in the malaria parasite genome. Here, using a genome-wide in silico analysis, we confirmed that ATG genes whose products are required for vesicle expansion and completion are present, while genes involved in induction of autophagy and cargo packaging are mostly absent. We subsequently focused on the molecular and cellular function of P. falciparum ATG8 (PfATG8), an autophagosome membrane marker and key component of the autophagy pathway, throughout the parasite asexual and sexual erythrocytic stages. In this context, we showed that PfATG8 has a distinct and atypical role in parasite development. PfATG8 localized in the apicoplast and in vesicles throughout the cytosol during parasite development. Immunofluorescence assays of PfATG8 in apicoplast-minus parasites suggest that PfATG8 is involved in apicoplast biogenesis. Furthermore, treatment of parasite cultures with bafilomycin A 1 and chloroquine, both lysosomotropic agents that inhibit autophagosome and lysosome fusion, resulted in dramatic morphological changes of the apicoplast, and parasite death. Furthermore, deep proteomic analysis of components associated with PfATG8 indicated that it may possibly be involved in ribophagy and piecemeal microautophagy of the nucleus. Collectively, our data revealed the importance and specificity of the autophagy pathway in the malaria parasite and offer potential novel therapeutic strategies.

Several systems-level datasets designed to dissect host-pathogen interactions during influenza A infection have been reported. However, apparent discordance among these data has hampered their full utility toward advancing mechanistic and therapeutic knowledge. To collectively reconcile these datasets, we performed a meta-analysis of data from eight published RNAi screens and integrated these data with three protein interaction datasets, including one generated within the context of this study. Further integration of these data with global virus-host interaction analyses revealed a functionally validated biochemical landscape of the influenza-host interface, which can be queried through a simplified and customizable web portal (http://www.metascape.org/IAV). Follow-up studies revealed that the putative ubiquitin ligase UBR4 associates with the viral M2 protein and promotes apical transport of viral proteins. Taken together, the integrative analysis of influenza OMICs datasets illuminates a viral-host network of high-confidence human proteins that are essential for influenza A virus replication.

Microbial resistance to chemotherapy has caused countless deaths where malaria is endemic. Chemotherapy may fail either due to pre-existing resistance or evolution of drug-resistant parasites. Here we use a diverse set of antimalarial compounds to investigate the acquisition of drug resistance and the degree of cross-resistance against common resistance alleles. We assess cross-resistance using a set of 15 parasite lines carrying resistance-conferring alleles in pfatp4, cytochrome bc1, pfcarl, pfdhod, pfcrt, pfmdr, pfdhfr, cytoplasmic prolyl t-RNA synthetase or hsp90. Subsequently, we assess whether resistant parasites can be obtained after several rounds of drug selection. Twenty-three of the 48 in vitro selections result in resistant parasites, with time to resistance onset ranging from 15 to 300 days. Our data indicate that pre-existing resistance may not be a major hurdle for novel-target antimalarial candidates, and focusing our attention on fast-killing compounds may result in a slower onset of clinical resistance.

Despite extensive efforts to discover transcription factors and their binding sites in the human malaria parasite Plasmodium falciparum, only a few transcription factor binding motifs have been experimentally validated to date. As a consequence, gene regulation in P. falciparum is still poorly understood. There is now evidence that the chromatin architecture plays an important role in transcriptional control in malaria.

As genome-wide experiments and annotations become more prevalent, researchers increasingly require tools to help interpret data at this scale. Many functional genomics experiments involve partitioning the genome into labeled segments, such that segments sharing the same label exhibit one or more biochemical or functional traits. For example, a collection of ChlP-seq experiments yields a compendium of peaks, each labeled with one or more associated DNA-binding proteins. Similarly, manually or automatically generated annotations of functional genomic elements, including cis-regulatory modules and protein-coding or RNA genes, can also be summarized as genomic segmentations.

Discovering novel genes involved in immune evasion and drug resistance in the human malaria parasite, Plasmodium falciparum, is of critical importance to global health. Such knowledge may assist in the development of new effective vaccines and in the appropriate use of antimalarial drugs. By performing a full-genome scan of allelic variability in 14 field and laboratory strains of P. falciparum, we comprehensively identified approximately 500 genes evolving at higher than neutral rates. The majority of the most variable genes have paralogs within the P. falciparum genome and may be subject to a different evolutionary clock than those without. The group of 211 variable genes without paralogs contains most known immunogens and a few drug targets, consistent with the idea that the human immune system and drug use is driving parasite evolution. We also reveal gene-amplification events including one surrounding pfmdr1, the P. falciparum multidrug-resistance gene, and a previously uncharacterized amplification centered around the P. falciparum GTP cyclohydrolase gene, the first enzyme in the folate biosynthesis pathway. Although GTP cyclohydrolase is not the known target of any current drugs, downstream members of the pathway are targeted by several widely used antimalarials. We speculate that an amplification of the GTP cyclohydrolase enzyme in the folate biosynthesis pathway may increase flux through this pathway and facilitate parasite resistance to antifolate drugs.

Despite essential roles played by long noncoding RNAs (lncRNAs) in development and disease, methods to determine lncRNA cis-elements are lacking. Here, we developed a screening method named "Tiling CRISPR" to identify lncRNA functional domains. Using this approach, we identified Xist A-Repeats as the silencing domain, an observation in agreement with published work, suggesting Tiling CRISPR feasibility. Mechanistic analysis suggested a novel function for Xist A-repeats in promoting Xist transcription. Overall, our method allows mapping of lncRNA functional domains in an unbiased and potentially high-throughput manner to facilitate the understanding of lncRNA functions.

Eukaryotic flagella are conserved microtubule-based organelles that drive cell motility. Plasmodium, the causative agent of malaria, has a single flagellate stage: the male gamete in the mosquito. Three rounds of endomitotic division in male gametocyte together with an unusual mode of flagellum assembly rapidly produce eight motile gametes. These processes are tightly coordinated, but their regulation is poorly understood. To understand this important developmental stage, we studied the function and location of the microtubule-based motor kinesin-8B, using gene-targeting, electron microscopy, and live cell imaging. Deletion of the kinesin-8B gene showed no effect on mitosis but disrupted 9+2 axoneme assembly and flagellum formation during male gamete development and also completely ablated parasite transmission. Live cell imaging showed that kinesin-8B-GFP did not co-localise with kinetochores in the nucleus but instead revealed a dynamic, cytoplasmic localisation with the basal bodies and the assembling axoneme during flagellum formation. We, thus, uncovered an unexpected role for kinesin-8B in parasite flagellum formation that is vital for the parasite life cycle.

Empirical evidence suggests that the malaria parasite Plasmodium falciparum employs a broad range of mechanisms to regulate gene transcription throughout the organism's complex life cycle. To better understand this regulatory machinery, we assembled a rich collection of genomic and epigenomic data sets, including information about transcription factor (TF) binding motifs, patterns of covalent histone modifications, nucleosome occupancy, GC content, and global 3D genome architecture. We used these data to train machine learning models to discriminate between high-expression and low-expression genes, focusing on three distinct stages of the red blood cell phase of the Plasmodium life cycle. Our results highlight the importance of histone modifications and 3D chromatin architecture in Plasmodium transcriptional regulation and suggest that AP2 transcription factors may play a limited regulatory role, perhaps operating in conjunction with epigenetic factors.

Respiratory syncytial virus (RSV) infection can cause mucus overproduction and bronchiolitis in infants leading to severe disease and hospitalization. As a therapeutic strategy, immune modulatory agents may help prevent RSV-driven immune responses that cause severe airway disease. We developed a high throughput screen to identify compounds that reduced RSV-driven mucin 5AC (Muc5AC) expression and identified dexamethasone. Despite leading to a pronounced reduction in RSV-driven Muc5AC, dexamethasone increased RSV infection in vitro and delayed viral clearance in mice. This correlated with reduced expression of a subset of immune response genes and reduced lymphocyte infiltration in vivo. Interestingly, dexamethasone increased RSV infection levels without altering antiviral interferon signaling. In summary, the immunosuppressive activities of dexamethasone had favorable inhibitory effects on RSV-driven mucus production yet prevented immune defense activities that limit RSV infection in vitro and in vivo. These findings offer an explanation for the lack of efficacy of glucocorticoids in RSV-infected patients.

Proteins interacting with DNA are fundamental for mediating processes such as gene expression, DNA replication and maintenance of genome integrity. Accumulating evidence suggests that the chromatin of apicomplexan parasites, such as Plasmodium falciparum, is highly organized, and this structure provides an epigenetic mechanism for transcriptional regulation. To investigate how parasite chromatin structure is being regulated, we undertook comparative genomics analysis using 12 distinct eukaryotic genomes. We identified conserved and parasite-specific chromatin-associated domains (CADs) and proteins (CAPs). We then used the chromatin enrichment for proteomics (ChEP) approach to experimentally capture CAPs in P. falciparum. A topological scoring analysis of the proteomics dataset revealed stage-specific enrichments of CADs and CAPs. Finally, we characterized, two candidate CAPs: a conserved homologue of the structural maintenance of chromosome 3 protein and a homologue of the crowded-like nuclei protein, a plant-like protein functionally analogous to animal nuclear lamina proteins. Collectively, our results provide a comprehensive overview of CAPs in apicomplexans, and contribute to our understanding of the complex molecular components regulating chromatin structure and genome architecture in these deadly parasites.

Eukaryotic cell proliferation requires chromosome replication and precise segregation to ensure daughter cells have identical genomic copies. Species of the genus Plasmodium, the causative agents of malaria, display remarkable aspects of nuclear division throughout their life cycle to meet some peculiar and unique challenges to DNA replication and chromosome segregation. The parasite undergoes atypical endomitosis and endoreduplication with an intact nuclear membrane and intranuclear mitotic spindle. To understand these diverse modes of Plasmodium cell division, we have studied the behaviour and composition of the outer kinetochore NDC80 complex, a key part of the mitotic apparatus that attaches the centromere of chromosomes to microtubules of the mitotic spindle. Using NDC80-GFP live-cell imaging in Plasmodium berghei, we observe dynamic spatiotemporal changes during proliferation, including highly unusual kinetochore arrangements during sexual stages. We identify a very divergent candidate for the SPC24 subunit of the NDC80 complex, previously thought to be missing in Plasmodium, which completes a canonical, albeit unusual, NDC80 complex structure. Altogether, our studies reveal the kinetochore to be an ideal tool to investigate the non-canonical modes of chromosome segregation and cell division in Plasmodium.

Burkholderia mallei, the causative agent of glanders, is a gram-negative intracellular bacterium. Depending on different routes of infection, the disease is manifested by pneumonia, septicemia, and chronic infections of the skin. B. mallei poses a serious biological threat due to its ability to infect via aerosol route, resistance to multiple antibiotics and to date there are no US Food and Drug Administration (FDA) approved vaccines available. Induction of innate immunity, inflammatory cytokines and chemokines following B. mallei infection, have been observed in in vitro and small rodent models; however, a global characterization of host responses has never been systematically investigated using a non-human primate (NHP) model. Here, using a liquid chromatography-tandem mass spectrometry (LC-MS/MS) approach, we identified alterations in expression levels of host proteins in peripheral blood mononuclear cells (PBMCs) originating from naïve rhesus macaques (Macaca mulatta), African green monkeys (Chlorocebus sabaeus), and cynomolgus macaques (Macaca fascicularis) exposed to aerosolized B. mallei. Gene ontology (GO) analysis identified several statistically significant overrepresented biological annotations including complement and coagulation cascade, nucleoside metabolic process, vesicle-mediated transport, intracellular signal transduction and cytoskeletal protein binding. By integrating an LC-MS/MS derived proteomics dataset with a previously published B. mallei host-pathogen interaction dataset, a statistically significant predictive protein-protein interaction (PPI) network was constructed. Pharmacological perturbation of one component of the PPI network, specifically ezrin, reduced B. mallei mediated interleukin-1β (IL-1β). On the contrary, the expression of IL-1β receptor antagonist (IL-1Ra) was upregulated upon pretreatment with the ezrin inhibitor. Taken together, inflammasome activation as demonstrated by IL-1β production and the homeostasis of inflammatory response is critical during the pathogenesis of glanders. Furthermore, the topology of the network reflects the underlying molecular mechanism of B. mallei infections in the NHP model.

The sexual stages of malarial parasites are essential for the mosquito transmission of the disease and therefore are the focus of transmission-blocking drug and vaccine development. In order to better understand genes important to the sexual development process, the transcriptomes of high-purity stage I-V Plasmodium falciparum gametocytes were comprehensively profiled using a full-genome high-density oligonucleotide microarray. The interpretation of this transcriptional data was aided by applying a novel knowledge-based data-mining algorithm termed ontology-based pattern identification (OPI) using current information regarding known sexual stage genes as a guide. This analysis resulted in the identification of a sexual development cluster containing 246 genes, of which approximately 75% were hypothetical, exhibiting highly-correlated, gametocyte-specific expression patterns. Inspection of the upstream promoter regions of these 246 genes revealed putative cis-regulatory elements for sexual development transcriptional control mechanisms. Furthermore, OPI analysis was extended using current annotations provided by the Gene Ontology Consortium to identify 380 statistically significant clusters containing genes with expression patterns characteristic of various biological processes, cellular components, and molecular functions. Collectively, these results, available as part of a web-accessible OPI database (http://carrier.gnf.org/publications/Gametocyte), shed light on the components of molecular mechanisms underlying parasite sexual development and other areas of malarial parasite biology.