Epigenetically regulated RNA-binding proteins signify malaria hypnozoite dormancy.

Dormancy enables relapsing malaria parasites, such as Plasmodium vivax and cynomolgi, to survive unfavorable conditions. It is enabled by hypnozoites, parasites remaining quiescent inside hepatocytes before reactivating and establishing blood-stage infection. We integrate omics approaches to explore gene-regulatory mechanisms underlying hypnozoite dormancy. Genome-wide profiling of activating and repressing histone marks identifies a few genes that get silenced by heterochromatin during hepatic infection of relapsing parasites. By combining single-cell transcriptomics, chromatin accessibility profiling, and fluorescent in situ RNA hybridization, we show that these genes are expressed in hypnozoites and that their silencing precedes parasite development. Intriguingly, these hypnozoite-specific genes mainly encode proteins with RNA-binding domains. We hence hypothesize that these likely repressive RNA-binding proteins keep hypnozoites in a developmentally competent but dormant state and that heterochromatin-mediated silencing of the corresponding genes aids reactivation. Exploring the regulation and exact function of these proteins hence could provide clues for targeted reactivation and killing of these latent pathogens.

A genome-wide map of DNA replication at single-molecule resolution in the malaria parasite Plasmodium falciparum.

The malaria parasite Plasmodium falciparum replicates via schizogony: an unusual type of cell cycle involving asynchronous replication of multiple nuclei within the same cytoplasm. Here, we present the first comprehensive study of DNA replication origin specification and activation during Plasmodium schizogony. Potential replication origins were abundant, with ORC1-binding sites detected every ∼800 bp. In this extremely A/T-biased genome, the sites were biased towards areas of higher G/C content, and contained no specific sequence motif. Origin activation was then measured at single-molecule resolution using newly developed DNAscent technology: a powerful method of detecting replication fork movement via base analogues in DNA sequenced on the Oxford Nanopore platform. Unusually, origins were preferentially activated in areas of low transcriptional activity, and replication forks also moved fastest through lowly transcribed genes. This contrasts with the way that origin activation is organised in other systems, such as human cells, and suggests that P. falciparum has evolved its S-phase specifically to minimise conflicts between transcription and origin firing. This may be particularly important to maximise the efficiency and accuracy of schizogony, with its multiple rounds of DNA replication and its absence of canonical cell-cycle checkpoints.

The troubled puberty of malaria parasites.

Sexual differentiation of malaria parasites is essential for transmission, yet the underlying mechanisms are poorly understood. Russell et al. elegantly combined a loss-of-function screen with single-cell RNA-sequencing to identify key factors in this process. Gomes et al. further characterized one of them, MD1, as a regulator contributing to male fate determination.

Gene-by-gene screen of the unknown proteins encoded on Plasmodium falciparum chromosome 3.

Taxon-specific proteins are key determinants defining the biology of all organisms and represent prime drug targets in pathogens. However, lacking comparability with proteins in other lineages makes them particularly difficult to study. In malaria parasites, this is exacerbated by technical limitations. Here, we analyzed the cellular location, essentiality, function, and, in selected cases, interactome of all unknown non-secretory proteins encoded on an entire P. falciparum chromosome. The nucleus was the most common localization, indicating that it is a hotspot of parasite-specific biology. More in-depth functional studies with four proteins revealed essential roles in DNA replication and mitosis. The mitosis proteins defined a possible orphan complex and a highly diverged complex needed for spindle-kinetochore connection. Structure-function comparisons indicated that the taxon-specific proteins evolved by different mechanisms. This work demonstrates the feasibility of gene-by-gene screens to elucidate the biology of malaria parasites and reveal critical parasite-specific processes of interest as drug targets.

Human branching cholangiocyte organoids recapitulate functional bile duct formation.

Human cholangiocyte organoids show great promise for regenerative therapies and in vitro modeling of bile duct development and diseases. However, the cystic organoids lack the branching morphology of intrahepatic bile ducts (IHBDs). Here, we report establishing human branching cholangiocyte organoid (BRCO) cultures. BRCOs self-organize into complex tubular structures resembling the IHBD architecture. Single-cell transcriptomics and functional analysis showed high similarity to primary cholangiocytes, and importantly, the branching growth mimics aspects of tubular development and is dependent on JAG1/NOTCH2 signaling. When applied to cholangiocarcinoma tumor organoids, the morphology changes to an in vitro morphology like primary tumors. Moreover, these branching cholangiocarcinoma organoids (BRCCAOs) better match the transcriptomic profile of primary tumors and showed increased chemoresistance to gemcitabine and cisplatin. In conclusion, BRCOs recapitulate a complex process of branching morphogenesis in vitro. This provides an improved model to study tubular formation, bile duct functionality, and associated biliary diseases.

Cholangiocyte organoids from human bile retain a local phenotype and can repopulate bile ducts in vitro.

The well-established 3D organoid culture method enabled efficient expansion of cholangiocyte-like cells from intrahepatic (IHBD) and extrahepatic bile duct (EHBD) tissue biopsies. The extensive expansion capacity of these organoids enables various applications, from cholangiocyte disease modelling to bile duct tissue engineering. Recent research demonstrated the feasibility of culturing cholangiocyte organoids from bile, which was minimal-invasive collected via endoscopic retrograde pancreaticography (ERCP). However, a detailed analysis of these bile cholangiocyte organoids (BCOs) and the cellular region of origin was not yet demonstrated. In this study, we characterize BCOs and mirror them to the already established organoids initiated from IHBD- and EHBD-tissue. We demonstrate successful organoid-initiation from extrahepatic bile collected from gallbladder after resection and by ERCP or percutaneous transhepatic cholangiopathy from a variety of patients. BCOs initiated from these three sources of bile all show features similar to in vivo cholangiocytes. The regional-specific characteristics of the BCOs are reflected by the exclusive expression of regional common bile duct genes (HOXB2 and HOXB3) by ERCP-derived BCOs and gallbladder-derived BCOs expressing gallbladder-specific genes. Moreover, BCOs have limited hepatocyte-fate differentiation potential compared to intrahepatic cholangiocyte organoids. These results indicate that organoid-initiating cells in bile are likely of local (extrahepatic) origin and are not of intrahepatic origin. Regarding the functionality of organoid initiating cells in bile, we demonstrate that BCOs efficiently repopulate decellularized EHBD scaffolds and restore the monolayer of cholangiocyte-like cells in vitro. Bile samples obtained through minimally invasive procedures provide a safe and effective alternative source of cholangiocyte organoids. The shedding of (organoid-initiating) cholangiocytes in bile provides a convenient source of organoids for regenerative medicine.

The ApiAP2 factor PfAP2-HC is an integral component of heterochromatin in the malaria parasite Plasmodium falciparum.

Malaria parasites undergo a complex life cycle in the human host and the mosquito vector. The ApiAP2 family of DNA-binding proteins plays a dominant role in parasite development and life cycle progression. Most ApiAP2 factors studied to date act as transcription factors regulating stage-specific gene expression. Here, we characterized an ApiAP2 factor in Plasmodium falciparum that we termed PfAP2-HC. We demonstrate that PfAP2-HC specifically binds to heterochromatin throughout the genome. Intriguingly, PfAP2-HC does not bind DNA in vivo and recruitment of PfAP2-HC to heterochromatin is independent of its DNA-binding domain but strictly dependent on heterochromatin protein 1. Furthermore, our results suggest that PfAP2-HC functions neither in the regulation of gene expression nor in heterochromatin formation or maintenance. In summary, our findings reveal PfAP2-HC as a core component of heterochromatin in malaria parasites and identify unexpected properties and substantial functional divergence among the members of the ApiAP2 family of regulatory proteins.

An epigenetic map of malaria parasite development from host to vector.

The malaria parasite replicates asexually in the red blood cells of its vertebrate host employing epigenetic mechanisms to regulate gene expression in response to changes in its environment. We used chromatin immunoprecipitation followed by sequencing in conjunction with RNA sequencing to create an epigenomic and transcriptomic map of the developmental transition from asexual blood stages to male and female gametocytes and to ookinetes in the rodent malaria parasite Plasmodium berghei. Across the developmental stages examined, heterochromatin protein 1 associates with variantly expressed gene families localised at subtelomeric regions and variant gene expression based on heterochromatic silencing is observed only in some genes. Conversely, the euchromatin mark histone 3 lysine 9 acetylation (H3K9ac) is abundant in non-heterochromatic regions across all developmental stages. H3K9ac presents a distinct pattern of enrichment around the start codon of ribosomal protein genes in all stages but male gametocytes. Additionally, H3K9ac occupancy positively correlates with transcript abundance in all stages but female gametocytes suggesting that transcription in this stage is independent of H3K9ac levels. This finding together with known mRNA repression in female gametocytes suggests a multilayered mechanism operating in female gametocytes in preparation for fertilization and zygote development, coinciding with parasite transition from host to vector.

A Kelch13-defined endocytosis pathway mediates artemisinin resistance in malaria parasites.

Artemisinin and its derivatives (ARTs) are the frontline drugs against malaria, but resistance is jeopardizing their effectiveness. ART resistance is mediated by mutations in the parasite's Kelch13 protein, but Kelch13 function and its role in resistance remain unclear. In this study, we identified proteins located at a Kelch13-defined compartment. Inactivation of eight of these proteins, including Kelch13, rendered parasites resistant to ART, revealing a pathway critical for resistance. Functional analysis showed that these proteins are required for endocytosis of hemoglobin from the host cell. Parasites with inactivated Kelch13 or a resistance-conferring Kelch13 mutation displayed reduced hemoglobin endocytosis. ARTs are activated by degradation products of hemoglobin. Hence, reduced activity of Kelch13 and its interactors diminishes hemoglobin endocytosis and thereby ART activation, resulting in parasite resistance.

Gene knockdown in malaria parasites via non-canonical RNAi.

The lack of endogenous RNAi machinery in the malaria parasite Plasmodium hampers gene annotation and hence antimalarial drug and vaccine development. Here, we engineered rodent Plasmodium berghei to express a minimal, non-canonical RNAi machinery that solely requires Argonaute 2 (Ago2) and a modified short hairpin RNA, so-called AgoshRNA. Using this strategy, we achieved robust and specific gene knockdown throughout the entire parasite life cycle. We also successfully silenced the endogenous gene perforin-like protein 2, phenocopying a full gene knockout. Transcriptionally restricting Ago2 expression to the liver stage further enabled us to perform a stage-specific gene knockout. The RNAi-competent Plasmodium lines reported here will be a valuable resource for loss-of-function phenotyping of the many uncharacterized genes of Plasmodium in low or high throughput, without the need to engineer the target gene locus. Thereby, our new strategy and transgenic Plasmodium lines will ultimately benefit the discovery of urgently needed antimalarial drug and vaccine candidates. Generally, the ability to render RNAi-negative organisms RNAi-competent by mere introduction of two components, Ago2 and AgoshRNA, is a unique paradigm that should find broad applicability in other species.

Epigenetic reader complexes of the human malaria parasite, Plasmodium falciparum.

Epigenetic regulatory mechanisms are central to the development and survival of all eukaryotic organisms. These mechanisms critically depend on the marking of chromatin domains with distinctive histone tail modifications (PTMs) and their recognition by effector protein complexes. Here we used quantitative proteomic approaches to unveil interactions between PTMs and associated reader protein complexes of Plasmodium falciparum, a unicellular parasite causing malaria. Histone peptide pull-downs with the most prominent and/or parasite-specific PTMs revealed the binding preference for 14 putative and novel reader proteins. Amongst others, they highlighted the acetylation-level-dependent recruitment of the BDP1/BDP2 complex and identified an PhD-finger protein (PHD 1, PF3D7_1008100) that could mediate a cross-talk between H3K4me2/3 and H3K9ac marks. Tagging and interaction proteomics of 12 identified proteins unveiled the composition of 5 major epigenetic complexes, including the elusive TBP-associated-factor complex as well as two distinct GCN5/ADA2 complexes. Furthermore, it has highlighted a remarkable degree of interaction between these five (sub)complexes. Collectively, this study provides an extensive inventory of PTM-reader interactions and composition of epigenetic complexes. It will not only fuel further explorations of gene regulation amongst ancient eukaryotes, but also provides a stepping stone for exploration of PTM-reader interactions for antimalarial drug development.

What functional genomics has taught us about transcriptional regulation in malaria parasites.

Malaria parasites are characterized by a complex life cycle that is accompanied by dynamic gene expression patterns. The factors and mechanisms that regulate gene expression in these parasites have been searched for even before the advent of next generation sequencing technologies. Functional genomics approaches have substantially boosted this area of research and have yielded significant insights into the interplay between epigenetic, transcriptional and post-transcriptional mechanisms. Recently, considerable progress has been made in identifying sequence-specific transcription factors and DNA-encoded regulatory elements. Here, we review the insights obtained from these efforts including the characterization of core promoters, the involvement of sequence-specific transcription factors in life cycle progression and the mapping of gene regulatory elements. Furthermore, we discuss recent developments in the field of functional genomics and how they might contribute to further characterization of this complex gene regulatory network.

Chromatin Accessibility-Based Characterization of the Gene Regulatory Network Underlying Plasmodium falciparum Blood-Stage Development.

Underlying the development of malaria parasites within erythrocytes and the resulting pathogenicity is a hardwired program that secures proper timing of gene transcription and production of functionally relevant proteins. How stage-specific gene expression is orchestrated in vivo remains unclear. Here, using the assay for transposase accessible chromatin sequencing (ATAC-seq), we identified ∼4,000 regulatory regions in P. falciparum intraerythrocytic stages. The vast majority of these sites are located within 2 kb upstream of transcribed genes and their chromatin accessibility pattern correlates positively with abundance of the respective mRNA transcript. Importantly, these regions are sufficient to drive stage-specific reporter gene expression and DNA motifs enriched in stage-specific sets of regulatory regions interact with members of the P. falciparum AP2 transcription factor family. Collectively, this study provides initial insights into the in vivo gene regulatory network of P. falciparum intraerythrocytic stages and should serve as a valuable resource for future studies.

Comparative Heterochromatin Profiling Reveals Conserved and Unique Epigenome Signatures Linked to Adaptation and Development of Malaria Parasites.

Heterochromatin-dependent gene silencing is central to the adaptation and survival of Plasmodium falciparum malaria parasites, allowing clonally variant gene expression during blood infection in humans. By assessing genome-wide heterochromatin protein 1 (HP1) occupancy, we present a comprehensive analysis of heterochromatin landscapes across different Plasmodium species, strains, and life cycle stages. Common targets of epigenetic silencing include fast-evolving multi-gene families encoding surface antigens and a small set of conserved HP1-associated genes with regulatory potential. Many P. falciparum heterochromatic genes are marked in a strain-specific manner, increasing the parasite's adaptive capacity. Whereas heterochromatin is strictly maintained during mitotic proliferation of asexual blood stage parasites, substantial heterochromatin reorganization occurs in differentiating gametocytes and appears crucial for the activation of key gametocyte-specific genes and adaptation of erythrocyte remodeling machinery. Collectively, these findings provide a catalog of heterochromatic genes and reveal conserved and specialized features of epigenetic control across the genus Plasmodium.

GDV1 induces sexual commitment of malaria parasites by antagonizing HP1-dependent gene silencing.

Malaria is caused by Plasmodium parasites that proliferate in the bloodstream. During each replication cycle, some parasites differentiate into gametocytes, the only forms able to infect the mosquito vector and transmit malaria. Sexual commitment is triggered by activation of AP2-G, the master transcriptional regulator of gametocytogenesis. Heterochromatin protein 1 (HP1)-dependent silencing of ap2-g prevents sexual conversion in proliferating parasites. In this study, we identified Plasmodium falciparum gametocyte development 1 (GDV1) as an upstream activator of sexual commitment. We found that GDV1 targeted heterochromatin and triggered HP1 eviction, thus derepressing ap2-g Expression of GDV1 was responsive to environmental triggers of sexual conversion and controlled via a gdv1 antisense RNA. Hence, GDV1 appears to act as an effector protein that induces sexual differentiation by antagonizing HP1-dependent gene silencing.

Characterization of the Nucleosome Landscape by Micrococcal Nuclease-Sequencing (MNase-seq).

MNase-seq allows the genome-wide examination of the nucleosome landscape by determination of nucleosome positioning and occupancy. Typically, native or formaldehyde fixed chromatin is subjected to digestion by micrococcal nuclease (MNase), which degrades linker DNA and yields mainly mono-nucleosomes. The resulting material can be processed directly or can be subjected to an optional chromatin immunoprecipitation step (MNase-ChIP-seq). De-crosslinked and purified DNA is then subjected to next-generation sequencing. The protocol presented here has been tailored for the analysis of nucleosome landscape in the malaria parasite, Plasmodium falciparum, but most steps are directly applicable to other cell types. We also discuss general considerations for experimental design and computational analysis, which are crucial for accurate investigation of the nucleosome landscape.

Malaria parasites possess a telomere repeat-binding protein that shares ancestry with transcription factor IIIA.

Telomere repeat-binding factors (TRFs) are essential components of the molecular machinery that regulates telomere function. TRFs are widely conserved across eukaryotes and bind duplex telomere repeats via a characteristic MYB-type domain. Here, we identified the telomere repeat-binding protein PfTRZ in the malaria parasite Plasmodium falciparum, a member of the Alveolate phylum for which TRFs have not been described so far. PfTRZ lacks an MYB domain and binds telomere repeats via a C2H2-type zinc finger domain instead. In vivo, PfTRZ binds with high specificity to the telomeric tract and to interstitial telomere repeats upstream of subtelomeric virulence genes. Conditional depletion experiments revealed that PfTRZ regulates telomere length homeostasis and is required for efficient cell cycle progression. Intriguingly, we found that PfTRZ also binds to and regulates the expression of 5S rDNA genes. Combined with detailed phylogenetic analyses, our findings identified PfTRZ as a remote functional homologue of the basic transcription factor TFIIIA, which acquired a new function in telomere maintenance early in the apicomplexan lineage. Our work sheds unexpected new light on the evolution of telomere repeat-binding proteins and paves the way for dissecting the presumably divergent mechanisms regulating telomere functionality in one of the most deadly human pathogens.

H3.3 demarcates GC-rich coding and subtelomeric regions and serves as potential memory mark for virulence gene expression in Plasmodium falciparum.

Histones, by packaging and organizing the DNA into chromatin, serve as essential building blocks for eukaryotic life. The basic structure of the chromatin is established by four canonical histones (H2A, H2B, H3 and H4), while histone variants are more commonly utilized to alter the properties of specific chromatin domains. H3.3, a variant of histone H3, was found to have diverse localization patterns and functions across species but has been rather poorly studied in protists. Here we present the first genome-wide analysis of H3.3 in the malaria-causing, apicomplexan parasite, P. falciparum, which revealed a complex occupancy profile consisting of conserved and parasite-specific features. In contrast to other histone variants, PfH3.3 primarily demarcates euchromatic coding and subtelomeric repetitive sequences. Stable occupancy of PfH3.3 in these regions is largely uncoupled from the transcriptional activity and appears to be primarily dependent on the GC-content of the underlying DNA. Importantly, PfH3.3 specifically marks the promoter region of an active and poised, but not inactive antigenic variation (var) gene, thereby potentially contributing to immune evasion. Collectively, our data suggest that PfH3.3, together with other histone variants, indexes the P. falciparum genome to functionally distinct domains and contribute to a key survival strategy of this deadly pathogen.

Integrated transcriptomic and proteomic analyses of P. falciparum gametocytes: molecular insight into sex-specific processes and translational repression.

Sexual differentiation of malaria parasites into gametocytes in the vertebrate host and subsequent gamete fertilization in mosquitoes is essential for the spreading of the disease. The molecular processes orchestrating these transitions are far from fully understood. Here, we report the first transcriptome analysis of male and female Plasmodium falciparum gametocytes coupled with a comprehensive proteome analysis. In male gametocytes there is an enrichment of proteins involved in the formation of flagellated gametes; proteins involved in DNA replication, chromatin organization and axoneme formation. On the other hand, female gametocytes are enriched in proteins required for zygote formation and functions after fertilization; protein-, lipid- and energy-metabolism. Integration of transcriptome and proteome data revealed 512 highly expressed maternal transcripts without corresponding protein expression indicating large scale translational repression in P. falciparum female gametocytes for the first time. Despite a high degree of conservation between Plasmodium species, 260 of these 'repressed transcripts' have not been previously described. Moreover, for some of these genes, protein expression is only reported in oocysts and sporozoites indicating that repressed transcripts can be partitioned into short- and long-term storage. Finally, these data sets provide an essential resource for identification of vaccine/drug targets and for further mechanistic studies.

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.