Single-Cell RNA Sequencing Reveals Cellular Heterogeneity and Stage Transition under Temperature Stress in Synchronized Plasmodium falciparum Cells

The malaria parasite has a complex life cycle exhibiting phenotypic and morphogenic variations in two different hosts by existing in heterogeneous developmental states. To investigate this cellular heterogeneity of the parasite within the human host, we performed single-cell RNA sequencing of synchronized Plasmodium cells under control and temperature treatment conditions. Using the Malaria Cell Atlas (https://www.sanger.ac.uk/science/tools/mca) as a guide, we identified 9 subtypes of the parasite distributed across known intraerythrocytic stages. Interestingly, temperature treatment results in the upregulation of the AP2-G gene, the master regulator of sexual development in a small subpopulation of the parasites. Moreover, we identified a heterogeneous stress-responsive subpopulation (clusters 5, 6, and 7 [−10% of the total population]) that exhibits upregulation of stress response pathways under normal growth conditions. We also developed an online exploratory tool that will provide new insights into gene function under normal and temperature stress conditions. Thus, our study reveals important insights into cell-to-cell heterogeneity in the parasite population under temperature treatment that will be instrumental toward a mechanistic understanding of cellular adaptation and population dynamics in Plasmodium falciparum. IMPORTANCE The malaria parasite has a complex life cycle exhibiting phenotypic variations in two different hosts accompanied by cell-to-cell variability that is important for stress tolerance, immune evasion, and drug resistance. To investigate cellular heterogeneity determined by gene expression, we performed single-cell RNA sequencing (scRNA-seq) of about 12,000 synchronized Plasmodium cells under physiologically relevant normal (37°C) and temperature stress (40°C) conditions phenocopying the cyclic bouts of fever experienced during malarial infection. In this study, we found that parasites exhibit transcriptional heterogeneity in an otherwise morphologically synchronized culture. Also, a subset of parasites is continually committed to gametocytogenesis and stress-responsive pathways. These observations have important implications for understanding the mechanisms of drug resistance generation and vaccine development against the malaria parasite.</p

Identification of Co-Existing Mutations and Gene Expression Trends Associated With K13-Mediated Artemisinin Resistance in Plasmodium falciparum

Plasmodium falciparum infects millions and kills thousands of people annually the world over. With the emergence of artemisinin and/or multidrug resistant strains of the pathogen, it has become even more challenging to control and eliminate the disease. Multiomics studies of the parasite have started to provide a glimpse into the confounding genetics and mechanisms of artemisinin resistance and identified mutations in Kelch13 (K13) as a molecular marker of resistance. Over the years, thousands of genomes and transcriptomes of artemisinin-resistant/sensitive isolates have been documented, supplementing the search for new genes/pathways to target artemisinin-resistant isolates. This meta-analysis seeks to recap the genetic landscape and the transcriptional deregulation that demarcate artemisinin resistance in the field. To explore the genetic territory of artemisinin resistance, we use genomic single-nucleotide polymorphism (SNP) datasets from 2,517 isolates from 15 countries from the MalariaGEN Network (The Pf3K project, pilot data release 4, 2015) to dissect the prevalence, geographical distribution, and co-existing patterns of genetic markers associated with/enabling artemisinin resistance. We have identified several mutations which co-exist with the established markers of artemisinin resistance. Interestingly, K13-resistant parasites harbor α-ß hydrolase and putative HECT domain–containing protein genes with the maximum number of SNPs. We have also explored the multiple, publicly available transcriptomic datasets to identify genes from key biological pathways whose consistent deregulation may be contributing to the biology of resistant parasites. Surprisingly, glycolytic and pentose phosphate pathways were consistently downregulated in artemisinin-resistant parasites. Thus, this meta-analysis highlights the genetic and transcriptomic features of resistant parasites to propel further exploratory studies in the community to tackle artemisinin resistance.</p

Dynamic association of the H3K64 trimethylation mark with genes encoding exported proteins in Plasmodium falciparum

Epigenetic modifications have emerged as critical regulators of virulence genes and stage-specific gene expression in Plasmodium falciparum. However, the specific roles of histone core epigenetic modifications in regulating the stage-specific gene expression are not well understood. In this study, we report an unconventional trimethylation at lysine 64 on histone 3 (H3K64me3) and characterize its functional relevance in P. falciparum. We show that PfSET4 and PfSET5 proteins of P. falciparum methylate H3K64 and that they prefer the nucleosome as a substrate over free histone 3 proteins. Structural analysis of PfSET5 revealed that it interacts with the nucleosome as a dimer. The H3K64me3 mark is dynamic, being enriched in the ring and trophozoite stages and drastically reduced in the schizont stages. Stage-specific global chromatin immunoprecipitation -sequencing analysis of the H3K64me3 mark revealed the selective enrichment of this methyl mark on the genes of exported family proteins in the ring and trophozoite stages and a significant reduction of the same in the schizont stages. Collectively, our data identify a novel epigenetic mark that is associated with the subset of genes encoding for exported proteins, which may regulate their expression in different stages of P. falciparum.</p

LHX2 Interacts with the NuRD Complex and Regulates Cortical Neuron Subtype Determinants Fezf2 and Sox11

In the developing cerebral cortex, sequential transcriptional programs take neuroepithelial cells from proliferating progenitors to differentiated neurons with unique molecular identities. The regulatory changes that occur in the chromatin of the progenitors are not well understood. During deep layer neurogenesis, we show that transcription factor LHX2 binds to distal regulatory elements of Fezf2 and Sox11, critical determinants of neuron subtype identity in the mouse neocortex. We demonstrate that LHX2 binds to the nucleosome remodeling and histone deacetylase histone remodeling complex subunits LSD1, HDAC2, and RBBP4, which are proximal regulators of the epigenetic state of chromatin. When LHX2 is absent, active histone marks at the Fezf2 and Sox11 loci are increased. Loss of LHX2 produces an increase, and overexpression of LHX2 causes a decrease, in layer 5 Fezf2 and CTIP2-expressing neurons. Our results provide mechanistic insight into how LHX2 acts as a necessary and sufficient regulator of genes that control cortical neuronal subtype identity

Regulation of vascular signalling by nuclear Sprouty2 in fetal lung epithelial cells:Implications for co-ordinated airway and vascular branching in lung development

This work was supported by a Wellcome Trust project grant 088032/Z/08/Z to SCL.Sprouty2 (Spry2) acts as a central regulator of tubular growth and branch patterning in the developing mammalian lung by controlling both magnitude and duration of growth factor signalling. To determine if this protein coordinates airway and vascular growth factor signalling, we tested the hypothesis that Spry2 links the primary cue for airway outgrowth, fibroblast growth factor-10 (FGF-10), to genomic events underpinning the expression and release of vascular endothelial growth factor-A (VEGF-A). Using primary fetal distal lung epithelial cells (FDLE) from rat, and immortalised human bronchial epithelial cells (16HBE14o-), we identified a nuclear sub-population of Spry2 which interacted with regions of the rat and human VEGF-A promoter spanning the hypoxia response element (HRE) and adjacent 3′ sites. In FDLE cultured at the PO2 of the fetal lung, FGF-10 relieved the Spry2 interaction at the HRE region by promoting clearance of a 39 kDa form and this was accompanied by histone-3 S10K14 phosphoacetylation, promoter de-methylation, hypoxia inducible factor-1α activation and VEGF-A expression. This repressive characteristic of nuclear Spry2 was relieved in 16HBE14o- by shRNA knockdown, and stable expression of mutants (C218A; C221A) that do not interact with the VEGF-A promoter HRE region. We conclude that nuclear Spry2 acts as a molecular link which co-ordinates airway and vascular growth of the cardiopulmonary system. This identifies Spry2 as a contributing determinant of design optimality in the mammalian lung.Publisher PDFPeer reviewe

The double PHD finger domain of MOZ/MYST3 induces α-helical structure of the histone H3 tail to facilitate acetylation and methylation sampling and modification

Histone tail modifications control many nuclear processes by dictating the dynamic exchange of regulatory proteins on chromatin. Here we report novel insights into histone H3 tail structure in complex with the double PHD finger (DPF) of the lysine acetyltransferase MOZ/MYST3/KAT6A. In addition to sampling H3 and H4 modification status, we show that the DPF cooperates with the MYST domain to promote H3K9 and H3K14 acetylation, although not if H3K4 is trimethylated. Four crystal structures of an extended DPF alone and in complex with unmodified or acetylated forms of the H3 tail reveal the molecular basis of crosstalk between H3K4me3 and H3K14ac. We show for the first time that MOZ DPF induces α-helical conformation of H3K4-T11, revealing a unique mode of H3 recognition. The helical structure facilitates sampling of H3K4 methylation status, and proffers H3K9 and other residues for modification. Additionally, we show that a conserved double glycine hinge flanking the H3 tail helix is required for a conformational change enabling docking of H3K14ac with the DPF. In summary, our data provide the first observations of extensive helical structure in a histone tail, revealing the inherent ability of the H3 tail to adopt alternate conformations in complex with chromatin regulators

PBRM1 Regulates Stress Response in Epithelial Cells

Polybromo1 (PBRM1) is a chromatin remodeler subunit highly mutated in cancer, particularly clear cell renal carcinoma. PBRM1 is a member of the SWI/SNF subcomplex, PBAF (PBRM1-Brg1/Brm-associated factors), and is characterized by six tandem bromodomains. Here we establish a role for PBRM1 in epithelial cell maintenance through the expression of genes involved in cell adhesion, metabolism, stress response, and apoptosis. In support of a general role for PBRM1 in stress response and apoptosis, we observe that loss of PBRM1 results in an increase in reactive oxygen species generation and a decrease in cellular viability under stress conditions. We find that loss of PBRM1 promotes cell growth under favorable conditions but is required for cell survival under conditions of cellular stress

Analyses of co-operative transitions in Plasmodium falciparum β-ketoacyl acyl carrier protein reductase upon co-factor and acyl carrier protein binding

The type II fatty acid synthase pathway of Plasmodium falciparum is a validated unique target for developing novel antimalarials because of its intrinsic differences from the type I pathway operating in humans. β -Ketoacyl-acyl carrier protein reductase is the only enzyme of this pathway that has no isoforms and thus selective inhibitors can be developed for this player of the pathway. We report here intensive studies on the direct interactions of Plasmodiumβ -ketoacyl-acyl carrier protein reductase with its cofactor, NADPH, acyl carrier protein, acetoacetyl-coenzyme A and other ligands in solution, by monitoring the intrinsic fluorescence (γ<SUB>max</SUB> 334 nm) of the protein as a result of its lone tryptophan, as well as the fluorescence of NADPH (γ<SUB>max</SUB> 450 nm) upon binding to the enzyme. Binding of the reduced cofactor makes the enzyme catalytically efficient, as it increases the binding affinity of the substrate, acetoacetyl-coenzyme A, by 16-fold. The binding affinity of acyl carrier protein to the enzyme also increases by approximately threefold upon NADPH binding. Plasmodiumβ -ketoacyl-acyl carrier protein reductase exhibits negative, homotropic co-operative binding for NADPH, which is enhanced in the presence of acyl carrier protein. Acyl carrier protein increases the accessibility of NADPH to β-ketoacyl-acyl carrier protein reductase, as evident from the increase in the accessibility of the tryptophan of β -ketoacyl-acyl carrier protein reductase to acrylamide, from 81 to 98%. In the presence of NADP<SUP>+</SUP>, the reaction proceeds in the reverse direction (K<SUB>a</SUB> = 23.17 μ m<SUP>-1</SUP>). These findings provide impetus for exploring the influence of ligands on the structure-activity relationship of Plasmodiumβ-ketoacyl-acyl carrier protein reductase