Abstract
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.
eLife assessment
This valuable study underscores the significance of PfMORC in shaping chromatin and guiding transitions in the malaria parasite, Plasmodium falciparum, that are essential for its survival. Solid evidence reveals PfMORC's influence on genes related to antigenic variation and the parasite's lifecycle, marking PfMORC as a key regulator of parasite heterochromatin.
Introduction
Malaria is a mosquito-borne infectious disease that is caused by protozoan parasites of the genus Plasmodium. Among the five human-infecting species, Plasmodium falciparum is the deadliest, with over 600,000 deaths in 20211. To adapt to extreme environmental challenges, P. falciparum possesses unique strategies that direct the tightly coordinated changes in gene expression and control transcriptional switching in genes encoded by multigene families to ensure antigenic variation and immune evasion. Changes in gene expression throughout the parasite life cycle are controlled by a surprisingly low repertoire of transcription factors (TFs) encoded in the Plasmodium genome2–9. The 27 apicomplexan APETALA2 (ApiAP2) DNA-binding proteins are the only well documented parasite TFs known to contribute to the modulation of gene expression throughout various stages of the parasite’s development 3,5–8,10–13. Since the discovery of the AP2 gene family 3, evidence has alluded to their role as master regulators of gene expression throughout transitory phases of the parasite life cycle. The most well documented are the gametocyte specific TFs (AP2-G, AP2-G2) 7,14 with subsequent discovery of those associated with transition to sexual differentiation 13 as well as sporozoite (AP2-SP) and liver stages (AP2-L) 6,10. The Plasmodium genome, however, encompasses well over 5,000 protein-coding genes suggesting there are most likely other molecular components responsible for gene expression. It is believed that, to offset the relatively low TF range, Plasmodium has evolved additional mechanisms regulating gene expression including mechanisms that use epigenetics factors, RNA binding proteins or regulate chromatin structure15–20. All together these components work in combination to regulate the dynamic organization of DNA and the cascade of gene expression required for the parasite life cycle progression 21–28. Despite their importance, the identification and functional characterization of regulatory players controlling chromatin remains a challenge in this intractable organism. However, with the advent of sensitive technologies capable of capturing important chromatin associated regulatory complexes such as chromatin immunoprecipitation (ChIP), chromatin conformation capture Hi-C technologies, and chromatin enrichment for proteomics (ChEP), we can now solve important facets of molecular components controlling epigenetics and chromatin structure.
Microrchidia (MORC) belong to a highly conserved nuclear protein superfamily with widespread domain architectures that link MORCs with signaling-dependent chromatin remodeling and epigenetic regulation across plant and animal kingdoms including the apicomplexan parasites 29–32. In all organisms, MORC contains several domains that form a catalytically active ATPase. Apicomplexan MORC ATPases are encircled by Kelch-repeat β- propellers as well as a CW-type zinc finger domain functioning as a histone reader 33. In higher eukaryotes, MORCs were first identified as epigenetic regulators and chromatin remodelers in germ cell development. Currently, these proteins are shown to be involved in various human diseases including cancers and are expected to serve as important biomarkers for diagnosis and treatment 34. In the apicomplexan parasite, Toxoplasma gondii, TgMORC was shown to recruit the histone deacetylase HDAC3 to particular genome loci to regulate chromatin accessibility, restricting sexual commitment 33,35. TgMORC-depleted cells also resulted in change in gene expression with up regulation of secondary AP2 factors and a significant shift from asexual to sexual differentiation 32,33,36–38. Recently, ChEP experiments done in P. falciparum identified PfMORC as one of the highest enriched chromatin-bound proteins at different stages of the parasite intraerythrocytic development cycle (IDC) 28. PfMORC was also detected at a relatively high level throughout the parasite life cycle including sporozoites and liver stages 33,39 and has been identified in several protein pull down experiments as targeting both AP2 TFs and epigenetic factors 32,33. Immunoprecipitation experiments demonstrated that PfMORC seems to interact with AP2-G2 32, a TF that plays a critical role in the maturation of gametocytes. However, the genome wide distribution and exact function of this protein throughout the parasite life cycle remained elusive.
In this study, we apply CRISPR/Cas9 genomic editing technologies to determine the function, genome distribution, and indispensability of PfMORC throughout the parasite’s IDC. Immunoprecipitation of an HA-tagged parasite line validates the role of PfMORC in heterochromatin structure maintenance. Using downregulation of PfMORC induced through TetR-DOZI system, we demonstrate the functional significance of PfMORC in heterochromatin stability and gene repression. Immunofluorescence based assays and ChIP-seq experiments show that PfMORC localizes to heterochromatin clusters at or near var genes with significant overlap with a well-known signature of the parasite heterochromatin post-translational modification (PTM) H3K9-trimethylation (H3K9me3) mark, confirming a possible role in gene repression. Our findings also shed light on the function of PfMORC as a potential chromatin insulator near transcriptionally active genes required for cell cycle progression. Finally, Hi-C analyses demonstrate that downregulation of PfMORC results in significant dysregulation of chromatin architectural stability resulting in parasite death. All together our work provides significant insight into the role of PfMORC in the maintenance of the pathogens’ chromosomal architectural integrity and validates this protein as a promising target for novel therapeutic interventions.
Results
Generation of PfMORC-HA transgenic line shows nuclear localization in heterochromatin clusters
To characterize the role of PfMORC in P. falciparum, we applied CRISPR/Cas9 genome editing tool to add a 3X-HA tag at the C-terminal coding region of the Pfmorc locus (PF3D7_1468100) in an NF54 line 40 (Fig. 1a, Sup. data 1). Recovered transgenic parasites were cloned using limiting dilution and the correct incorporation of the tag into the parasite genome was validated via PCR and whole genome sequencing (WGS) (Fig. 1b, Sup. data 2). WGS analysis confirmed the presence of the HA tag at the expected locus with no obvious off-target effect but uncovered an additional nonsense mutation in the gametocyte development protein 1 (gdv1) gene (PF3D7_0935400) at amino acid position 561 out of 599 total. Mutations in gdv1 have been detected in the past and seem to emerge relatively often indicating a fitness benefit in the in vitro culture system used in the laboratory 8. The involvement of GDV1 as essential in sexual commitment is well substantiated and further study involving the role of PfMORC in sexual commitment could not be fully addressed 7,8,14,41–43 (Sup. data 2). Despite this finding, we were able to design experiments to determine the role of PfMORC throughout the IDC as no other major detrimental mutations were detected in the genome. We first validated PfMORC protein expression through Western blot analysis (Fig. S1). Regardless of potential protein degradation, results showed an expected band size within the ∼290 kDa range that was absent in NF54 control compared to our tagged line (Fig. S1). Immunofluorescence (IFA) analysis of intracellular parasites revealed that PfMORC is localized in the nucleus throughout the parasite IDC in punctate patterns (Fig. 1c). A single foci per parasite was detected in the nuclear periphery at the ring stage. An increased number of foci with a more diffuse signal could be detected as the number of DNA copies increased at the schizont stage. Moreover, punctuates were shown to have strong colocalization signals with histone H3K9me3 marks present at the different parasite developmental stages analyzed. This was most evident during the early ring stage of the parasite IDC where the chromatin organization is more compact with a single DNA copy. Altogether, our findings indicated that PfMORC is a true nuclear protein that is most likely associated with P. falciparum heterochromatin cluster(s).
![](https://prod--epp.elifesciences.org/iiif/2/92499%2Fv1%2Fcontent%2F557253v1_fig1.tif/full/max/0/default.jpg)
PfMORC-HA is associated with heterochromatin.
(a) Illustration of the PfMORC containing domains including Kelch type b propeller and ATPase domains using InterProScan. (b) Design strategy applied for PfMORC C-terminal HA tagging. PCR amplification of the genomic C-terminus end of Pfmorc region extending towards the 3’UTR (F1) as well as extension from the C-terminus towards the HA flanking sequence (F2) verifies correct insertion site. NF54 genomic DNA was used as negative control. (c) IFA experiment: PfMORC foci 22 expressing co-localization with H3K9me3 marks (red). Cell nuclei are stained with DAPI (blue). BF: brightfield (d) Protein immunoprecipitation: Significance plot representing PfMORC interactome recovered through immunoprecipitation followed by mass spectrometry (IP-MS). Graph list MORC (orange) bindings partners of highest affinity associated with TF regulation 22 22 and chromatin remodelers, erasers, and writers (purple). Proteins enriched in the PfMORC-HA samples compared with controls were filtered with log2 FC ≥ 2 and Z statistic > 5. (e) Gene ontology enrichment analysis of the significantly enriched proteins. The top 2 terms of Cellular Component (top) and top 3 terms of Molecular Function (bottom) are represented as -log10 (adjusted P-value) (Fisher’s exact test with Bonferroni adjustment).
PfMORC interacting partners include proteins involved in heterochromatin maintenance, transcription regulation and chromatin remodeling
To confirm the association of PfMORC with proteins involved in heterochromatin cluster(s), we performed immunoprecipitation (IP) followed by mass-spectrometry (MS) analysis using mature stages of PfMORC-HA and parental line as control. Our IP-MS analysis matched 1109 spectra to 105 unique peptide sequences corresponding to the PfMORC protein (on average) consistent with HA-affinity enrichment of PfMORC. To define a set of PfMORC associated proteins, we used the QPROT statistical framework 44 to compare proteins detected in PfMORC samples and controls. We identified 73 PfMORC associated proteins (log2 FC ≥ 2 and Z statistic > 5) (Fig. 1d, Sup. data 3). Gene Ontology (GO) enrichment analysis indicated that these proteins are nuclear and associated with heterochromatin, DNA binding, and transcription factor activity (Fig. 1e). Detailed analysis showed the detection of 8 AP2 transcription factors, including AP2-G2, AP2-O2, AP2-G5, and SIP2. The latter is involved in heterochromatin formation and chromosome end regulation 45. Among the other AP2s detected, 5 were previously identified by ChIP-seq as Plasmodium falciparum AP2 Heterochromatin-associated Factors (PfAP2-HFs) 46. Additional chromatin remodeling proteins were detected such as ISWI chromatin-remodeling complex ATPase (PF3D7_0624600) and SWIB/MDM2 domain-containing protein (PF3D7_0611400), as well as the epigenetic readers BDP3 (PF3D7_0110500) and PHD2 (PF3D7_1433400) 47. We also detected a significant enrichment of histone erasers and writers such as a N-acetyltransferase (PF3D7_1020700), HDAC1 (PF3D7_0925700), and SET3 (PF3D7_0827800). These data indicated that PfMORC interacts with multiple ApiAP2 TFs, chromatin remodelers, and epigenetic players associated with heterochromatin regions.
Genome-wide distribution of PfMORC indicates dual functionality
Given the multifunctional role of MORC proteins throughout apicomplexan parasites, we further explored the association of PfMORC with chromatin accessibility and stage-specific gene regulation by assessing the genome-wide distribution of PfMORC using chromatin immunoprecipitation followed by deep sequencing (ChIP–seq) in duplicate at the ring, trophozoite and schizont stages. Inspection of PfMORC binding sites across each chromosome revealed a strong signal for telomeric and subtelomeric regions of the parasite genome, as well as the internal clusters of genes responsible for antigenic variation such as var and rifin (Fig. 2a-b, Sup. data 4). We found that while only 47% of reads mapped to antigenic genes in the ring stage, 95% of reads mapped to the same genes during the trophozoite stage and 82% during the schizont stage. Although a vast majority of the genome is actively transcribed during the trophozoite stage, these antigenic gene families are under tight regulatory control to ensure mutually exclusive expression and participate in immune evasion 48,49. This could explain the disparity in the level of PfMORC binding of the coding regions between the ring and trophozoite stages (Fig. 2b, Fig. S2c).
![](https://prod--epp.elifesciences.org/iiif/2/92499%2Fv1%2Fcontent%2F557253v1_fig2.tif/full/max/0/default.jpg)
Genome wide distribution of PfMORC proteins.
(a) Chromosome distribution plots of PfMORC binding showing a predisposition for subtelomeric and internal var gene regions (red). Each track is input subtracted, and per-million read count normalized before normalizing the track height to allow for direct comparison between stages. Grey boxes indicate position of centromeres. (b) Profile plots showing PfMORC coverage from 0.5 kb 5’ of the transcription start site (TSS) to 1.5 kb 3’ of the TSS in the ring, trophozoite, and schizont stages. Each plot includes per-million read count normalized coverage at 1 bp resolution for all genes within the var and rifin gene families as well as gametocyte-specific genes. (c) PfMORC coverage of the gametocyte-specific transcription factor ap2-g.
We also evaluated PfMORC binding of stage-specific gene families including gametocyte-related genes, and merozoite surface proteins (msp) 50. Although most gametocyte genes contain low levels of PfMORC binding compared to antigenic gene families (Fig. 2b), the transcription factor AP2-G, known to be repressed during the IDC and required for sexual commitment, contains high levels of PfMORC throughout the gene body (Fig. 2d). For some of the msp genes that are usually expressed later during the IDC to prepare the parasite for egress and invasion of new erythrocytes, PfMORC binding was detected in the gene bodies at the trophozoite stage. We also detected a small switch in PfMORC binding sites, moving from their gene bodies to their intergenic and promoter regions at the trophozoite stage (Fig. S2d). This data implicates PfMORC as a potential insulator for some of these genes in addition to its role as a transcriptional repressor. PfMORC most likely moved away from the gene body to the regulatory regions surrounding the transcription start site (TSS) to guide RNA Polymerase and transcription factors and aid in activating expression of genes at the schizont stage that are crucial for egress and invasion. These results provide strong evidence for the direct effects that PfMORC binding has on tightly controlled antigenic and stage-specific genes including crucial gametocyte genes.
PfMORC is essential for P. falciparum survival
We next sought to confirm the functional relevance of PfMORC protein in parasite survival. While partial PfMORC knockdown using the glmS-ribozyme system had previously shown to have no significant effect on the parasite survival 39, another study using transposon mutagenesis (piggyBac) identified PfMORC as likely essential 51. To resolve these conflicting results, we applied a complementary approach using the CRISPR/Cas9 gene editing strategy to incorporate an inducible TetR-DOZI system to knockdown (KD) PfMORC 52 through administration of anhydrotetracycline (aTC) (Fig. 3a) 40,53. The protein was also modified to include a C-terminal 3x-HA tag. Parental and transgenic clones were validated via PCR (Fig. 3b) and WGS to confirm the correct insertion of the inducible system (Fig. S3) as well as the absence of major mutation that could explain some of the phenotypes observed. While results from our WGS validated our editing strategy without the detection of any obvious deleterious off-target effect, we identified a nonsense mutation in AP2-G gene (PF3D7_1222600) explaining our inability to obtain mature gametocytes. Protein expression and successful KD was confirmed via western blot with ∼ 20% and ∼ 60% decrease in PfMORC expression at 24h and 36h after aTC removal, respectively (Fig. S3). This observed downregulation of PfMORC was slightly above the level observed by Singh and colleagues 39.
![](https://prod--epp.elifesciences.org/iiif/2/92499%2Fv1%2Fcontent%2F557253v1_fig3.tif/full/max/0/default.jpg)
PfMORC is essential for cell survival.
(a) Diagram representation of PfMORC-HA-TetR-DOZI plasmid. (b). PCR amplification is used to verify genomic insertion using primers sets targeting 1.5 kbp of WT Pfmorc genome locus absent in transgenic line (MORC construct) (F3) as well as verification of HA insertion (F4) and TetR-DOZI system extending along 3’ UTR of the construct (F5). (c) Phenotypic and quantitative (d) analysis of parasite cell progression after aTC withdrawal at the ring stage (0-6 hpi) (2-way ANOVA, n=3, p≤0.0001). (e) Phenotypic and quantitative (f) analysis of parasite cell progression after aTC removal at the trophozoite stage of cell cycle progression (24 hpi) (2-way ANOVA, n=3, p≤0.0001).
We then performed a phenotypic assay on (+/-) aTC PfMORC-HA-TetR-DOZI parasite cultures. The assay was conducted in replicates of synchronized cultures that were then split into (+) aTC and (-) aTC conditions either at the ring or trophozoite stages on PfMORC-HA-TetR-DOZI clones and WT control lines (Fig. 3c-f, Sup. data 5). Sequential morphological screens were performed through Giemsa-stained smears and monitored by microscope imaging. As opposed to what was observed previously with the glmS-ribozyme system 39, aTC removal at the ring stage induced clear signs of stress and cell cycle arrest in early trophozoite and schizont stages of the first intraerythrocytic cycle (Fig. 3c-d) compared to (+) aTC PfMORC and WT parasites. Our data showed an ∼53%, 77%, and 84% (n=3, p≤0.0001) drop in parasitemia at 48-, 72-and 120-hours post invasion (hpi), respectively, when compared to control conditions. aTC supplemented and WT cultures showed unperturbed cell cycle progression, reinvasion, and morphological development. Interestingly, when aTC was withheld at the late trophozoite stages (24 to 30 hpi), parasites could complete the first cycle and successfully reinvading into new red blood cells (RBCs) (Fig. 3e-f). The evident hindered phenotypic response may either be associated with the delay time for the protein to be completely depleted or may suggest a more dispensable role of PfMORC in the schizont stage. However, clear indications of stress and cell cycle arrest were ultimately detected at trophozoite and early schizont stages of the second cell cycle. Quantitative analysis of (-) aTC PfMORC cultures revealed significantly decreased parasitemia of ∼37% and 66% (n=3, p≤0.0001) at 96 hpi and 120 hpi, respectively, compared to (+) aTC PfMORC and WT control lines (Fig. 3e-f, Sup. data 5), confirming the importance of the protein at the trophozoite and early schizont stages.
Effect of PfMORC Knockdown on parasite transcriptome
To define the effects of PfMORC KD on transcription, we conducted detailed time-course measurements of mRNA levels using RNA-sequencing (RNA-seq) throughout the parasite asexual cycle on (+/-) aTC PfMORC lines. Parasites were first synchronized and aTC was removed from one set of samples. Total RNA was then extracted at the trophozoite (24 hpi) and the schizont (36 hpi) stages to allow for detection of early and late changes in gene expression after aTC removal. Pairwise correlations analysis (Fig. S4a) between (+/-) aTC PfMORC treated lines at the two different time points confirmed appropriate reproducibility of our RNA-seq experiments.
Our results identified a relatively low number of differentially expressed genes at 24 hpi with 96 and 93 upregulated and downregulated genes, respectively (FDR < 0.05 and log2 FC > 0.5 or < -0.5) (Fig 4a, Sup. data 6). This data is in agreement with the absence of major phenotypic changes observed in (-) aTC PfMORC parasite cultures after 24 hours. GO enrichment analysis indicated upregulation of genes in response to xenobiotic stimulus including several phosphatases, hydrolases and heat shock proteins suggesting stress sensing of the parasites at this stage (Fig. 4c). At 36 hpi the number of genes that exhibit changes in gene expression were significantly higher with 1319 upregulated and 1150 downregulated (FDR < 0.05 and log2 FC > 0.5 or < -0.5) in (-) aTC PfMORC conditions, here also confirming the significant changes observed at the phenotypic level (Fig. 4b). These results emphasize the cascading effects of PfMORC KD as the parasite progressed throughout its cell cycle. GO analysis indicated an enrichment of upregulated genes involved in multiple pathways including ATP metabolic process, mitochondria, translation, food vacuole and protein folding (Fig. 4d); characteristic of not only stress, but also a clear signal of cell cycle arrest in (-) aTC PfMORC samples. Among upregulated genes, we also observed several var genes and genes exported at the surface of the red blood cell that could be linked to significant decreased in PfMORC binding of these genes families as well as a disorganization of the heterochromatin cluster(s) at the trophozoite and schizont stages. GO enrichment analysis for genes that were downregulated included many genes involved in DNA replication, chromosome organization and mitotic spindle further emphasizing a strong cell cycle arrest and absence of potential major compensatory mechanisms for cell division and chromatin organization (Fig. 4d). Additional downregulated genes were required for invasion such as several merozoite surface proteins. Our results clearly indicate stress sensing and an arrest in cell cycle progression at the trophozoite and schizont stages (Sup. data 5) 54.
![](https://prod--epp.elifesciences.org/iiif/2/92499%2Fv1%2Fcontent%2F557253v1_fig4.tif/full/max/0/default.jpg)
PfMORC KD on parasite transcriptome.
Volcano plots denoting upregulated (red), and downregulated (blue) genes discovered through differential expression analysis following PfMORC knockdown at (a) 24 hpi (b) 36 hpi. Gene ontology enrichment analysis for upregulated (red) and downregulated (blue) genes at (c) 24 hpi and (d) 36 hpi. (e) Overlap of differentially expressed genes and genes containing significant peaks called by PfMORC ChIP-seq analysis at the trophozoite and schizont stages.
Although these results are compelling, the small overlap observed between the ChIP-seq signals and RNA-seq results indicate that major changes observed in gene expression at the schizont stages are not simply the result of reduced PfMORC binding in targeted gene bodies but a combination of direct and indirect effects of the degradation of PfMORC that leads to cell cycle arrest and potential collapse of the heterochromatin (Fig. 4e).
Knockdown of PfMORC expression results in the loss of tightly regulated heterochromatin structures
To better understand the effect that downregulation of PfMORC has on the chromatin structure and investigate whether changes in chromatin accessibility may explain large changes in gene expression observed using RNA-seq, we performed Hi-C on (+/-) aTC PfMORC cultures at 24 hpi (trophozoite) and 36 hpi (schizont). Biological replicates for each sample were collected and used to generate Hi-C libraries with >37 million reads per replicate/sample, ensuring comprehensive coverage of both intrachromosomal and interchromosomal interactions. After processing (pairing, mapping, and quality filtering) the raw sequences via HiC-Pro 55, there were approximately 15 million (σ = ∼10 million) high quality interaction pairs per sample. Due to the high number of reads and relatively small size of the P. falciparum genome (23.3 Mb) compared to higher eukaryotes, we elected to bin our Hi-C data at 10-kb resolution, which allowed the identification of genome-wide patterns while not introducing much noise by binning at too high of a resolution. A high stratum-adjusted correlation (Fig. S4), especially at the 24 hpi time point, suggested that the chromatin structure was consistent between biological replicates, therefore we chose to combine replicates for downstream analyses and visualization.
Because of variation in the number of valid interaction pairs between (+/-) aTC PfMORC samples, 100 iterations of random sampling were performed on samples with higher read count to obtain ∼35 million and ∼9 million interaction pairs at 24 hpi and 36 hpi, respectively. ICED normalized intrachromosomal heatmaps displayed a high proportion of interactions at distances less than 10% the length of each chromosome as well as strong subtelomeric interactions and internal regions containing genes involved in antigenic variation (Fig. 5a-b, Fig. S6-S9) This pattern was similar to those observed previously at various stages of the life cycle 27,55,56 and confirm that genes involved in antigenic variation are usually confined within dense heterochromatin rich regions to aid the tight control of var genes necessary for their mutually exclusive expression and immune evasion 48,49. Data generated at 36 hpi in (-) aTC PfMORC parasite cultures were however tumultuous, with a decrease in defined heterochromatin borders and fewer long-range intrachromosomal interactions across the genome. This may indicate a significant loss of chromatin maintenance in (-) aTC PfMORC parasites (Fig. 5a-b).
![](https://prod--epp.elifesciences.org/iiif/2/92499%2Fv1%2Fcontent%2F557253v1_fig5.tif/full/max/0/default.jpg)
Loss of PfMORC expression correlates with heterochromatin expansion. Intrachromosomal interaction heatmaps of (+/-) aTC PfMORC for chromosome 4 at
(a) 24 hpi and (b) 36 hpi displaying heterochromatin clustering within antigenic (var, rifin and stevor) gene-dense regions (red). Differential interaction heatmaps highlight changes in chromatin structure following removal of aTC and subsequent PfMORC knockdown at (c) 24 hpi and (d) 36 hpi. (e) Whole-genome 3D models of the chromatin structure at both time points (24 hpi and 36 hpi) and (+/-) aTC.
To pinpoint which regions were most strongly affected by knockdown of PfMORC, we used Selfish 57 to identify differential intrachromosomal and interchromosomal interactions in the (+/-) aTC PfMORC samples. Although the 24 hpi and 36 hpi time points shared similarities in loci most highly affected by PfMORC KD on many chromosomes, the changes at 36 hpi were more significant with a higher log2 FC (FDR < 0.05) at most loci (Fig. 5c-d, Fig. S10-S11). We also observed consistent loss of interactions between most regions containing var genes. These results indicate a failed attempt of the (-) aTC PfMORC parasites to maintain their overall chromatin structure with a significant weakening of the tightly controlled heterochromatin cluster.
We further validated our results using PASTIS to generate coordinate matrices and subsequently visualize a consensus three-dimensional model of the chromosome folding and their overall organization within the nucleus 58. This can indicate changes in spatial distance that best describe the interaction data. The overall change in the chromatin 3D structure was clearly shown by our models of (+/-) aTC PfMORC samples (Fig. 5e). While the co-localization of centromeres and telomeres in different regions of the nucleus were conserved, the 3D chromatin structure of (-) aTC PfMORC at 36 hpi displayed clear opening of the chromatin and loss of interactions in most regions. This could explain the large changes in gene expression, including increased expression of all var and rifin genes in the (-) aTC PfMORC line (Sup. data 6). While heterochromatin maintenance may be essential for tight control of var gene expression, preservation of the overall structure of the chromatin may be necessary to regulate the accurate expression of most P. falciparum transcripts.
Discussion
Despite intensive investigations on the dynamic nature of the Plasmodium chromatin, progress has been slow in capturing the regulatory factors that maintain and control chromatin structure and gene expression throughout parasite development. Preliminary studies done by our lab and others have identified PfMORC as one of the most abundant chromatin-bound proteins 35,59. MORC proteins from plant to animal cells have a broad range of binding sites within the genome and participate in either DNA methylation establishment or heterochromatin formation 30,60–65. They have also been shown to co-localize with TFs in unmethylated promoter regions to alter chromatin accessibility and regulate TF binding and gene expression 62–65. Protein pull-down experiments in Plasmodium spp. and T. gondii studies have detected MORC in complexes with several AP2 TFs. In T. gondii, TgMORC was characterized as an upstream transcriptional repressor of sexual commitment 35. In Plasmodium, the role of PfMORC was still poorly understood. Conflicting results on the essentiality of PfMORC using either KD or KO strategies remained unresolved 39,51. Using the CRISPR/Cas9 genome editing tool, we have determined the nuclear localization, genome-wide distribution, and regulatory impacts of PfMORC on the parasite chromatin, transcriptome, and cell cycle development. We now have powerful evidence demonstrating that PfMORC is not only critical for parasite cell cycle progression and its survival but has a direct role in heterochromatin formation and gene silencing, including regulating immune evasion through antigenic variation. Disparities observed in the essentiality of PfMORC for the parasite survival in previous studies 39 are most likely the results of weak protein disruption highlighting the need for a significant downregulation of PfMORC for true functional analysis. Using IFA and protein pull downs, we have confirmed that PfMORC localizes to the nucleus and interacts with multiple ApiAP2 TFs, chromatin remodelers, and epigenetic players associated with heterochromatin regions including H3K9me3, SIP2 and the ISWI chromatin-remodeling complex (SWIB/MDM2)46 (Fig. 1). Using a combination of ChIP-seq, protein knock down, RNA-seq and Hi-C experiments, we have demonstrated that the MORC protein is essential for the tight regulation of gene expression through chromatin compaction, preventing access to gene promoters from TFs and the general transcriptional machinery in a stage specific manner. It is, however, important to recognize that the many pathways affected at the transcriptional level throughout the asexual stages in (-) aTC PfMORC lines are not only the direct results of PfMORC downregulation and decrease of its targeted DNA binding site. They are most likely a combination of direct and indirect effects of PfMORC KD warranted by the cell cycle arrest observed at the phenotypic level and the collapse of the chromatin organization confirmed using chromatin conformation capture experiment. These direct and indirect effects should be carefully considered, and a combination of functional genomic studies should be completed when interpreting changes in gene expression in mutant generated lines.
We also demonstrated that in addition to its direct role in heterochromatin formation and antigenic variation, PfMORC may act as a potential repressor and insulator to control a set of parasite specific genes including genes involved in parasite egress and invasion between the trophozoite to schizont stage transition 46. While MORC amino acid sequences reveal some conservation across species in their N-terminal portions including the tripartite GHKL ATPase domains, there is less conservation in their carboxy-terminal domains. As our data confirm the importance of PfMORC and the parasite specificity of several interacting partners, it is tempting to speculate that drugs targeting these protein complexes could lead to novel antiparasitic strategies.
Methods
Asexual parasite Culture and maintenance
Asexual P. falciparum strain NF54 parasites (MRA-1000) were propagated in 5% of human O+ erythrocytes and 10 mL of RPMI-1640 medium containing 0.5% Albumax II (Invitrogen), 2 mM L-glutamine, 50 mg/L hypoxanthine, 25 mM HEPES, 0.225% NaHCO3 and 10 mg/mL gentamicin. They were maintained at 37°C and gassed with a sterile mixture of 5% O2, 5% CO2 and 90% N2 66,67. Parasite synchronization was achieved using two 5% D-sorbitol treatments 8 hours apart 66,67.
Plasmids Construction
Flagging of the P. falciparum MORC (Pf3D7_1468100) gene spanning position (Ch14: 2785386 - 2797934 (-)) was performed using a two-plasmid design to insert a 3x HA tag. The pCasG-Cas9-sgRNA plasmid vector (gifted by Dr. Sean Prigge) contains the site to express the sgRNA, along with the yDHODH gene as the positive selection marker. The gRNA oligos (Sup. data 1) were ligated after digestion of the plasmid with BsaI.The pDC2-cam-Cas9-U6 plasmid (gifted by Dr. Marcus Lee) was digested with BamHI and Apal to remove the eGFP tag from the backbone. 453 bp of C-terminal region of PfMORC and 469 bp of 3’UTR were amplified from P. falciparum genomic DNA with their respective primers (Sup. data 1). A 3x HA-tag was fused to the C-terminal region of the amplified product along with the 3’UTR formed through Gibson assembly master mix (NEB, E2611S). To generate the PfMORC-HA knockdown constructs, a pKDPfAUBL plasmid (gifted by Dr. Sean Prigge) 68 was digested with AscI and AatII. Homology arms HA1 and HA2 of PfMORC were amplified with respective primers (Sup. data 1) having 20 bp overhang and inserted into the digested plasmid using Gibson assembly mix. The resulting pKD-MORC-HA-TetR-DOZI construct was linearized with EcoRV enzyme prior to transfection. All constructs were confirmed through restriction enzyme digestion and Sanger sequencing.
Plasmids were isolated from 250 mL cultures of Escherichia coli (XL10-Gold Ultracompetent Cells, Agilent Cat. 200314) and 60 µg of pDC2-MORC-HA or linearized pKD-MORC-HA-TetR-DOZI, were used with 60 µg of pCasG-plasmid containing gRNA to transfect 200 µl of fresh red blood cells (RBCs) infected with 3-5% early ring stage parasites. After one erythrocytic cycle, transfected cultures were supplemented with 1.5 µM WR99210 (2.6nM) (provided by the (Jacobus Pharmaceuticals, Princeton, NJ) and 2.5 μg/mL Blasticidin (RPI Corp B12150-0.1). PfMORC-HA-TetR-DOZI transfected parasites were maintained with 500 nM anhydrotetracycline (aTC) 40,53. Media and drug selection was replenished every 24 hours for 7 consecutive days after which DSM-1 drug selection was halted Once parasites were detected by microscopy, integration of the insert was confirmed by PCR amplification. To generate genetically homogenous parasite lines, the transfected parasites were serially diluted to approximately 0.5% parasite/well, into 96 well plates.
Molecular analysis of the transgenic lines
Genomic DNA (gDNA) was extracted and purified using the DNeasy Blood & Tissue kit (Qiagen) following instructions from the manufacturer. The diagnostic PCR analysis was used to genotype the transfected lines using the primers listed in Supplementary data 1. The PCR amplification was conducted using KAPA HiFi HotStart ReadyMix (Roche) and amplicons were analyzed by gel electrophoresis followed by sequencing. For whole genome sequencing, genomic DNA was fragmented using a Covaris S220 ultrasonicator and libraries were generated using KAPA LTP Library Preparation Kit (Roche, KK8230). To verify that the insertion was present in the genome at the correct location in both transfected lines, reads were mapped using Bowtie2 (v2.4.4) to the P. falciparum 3D7 reference genome (PlasmoDB, v48), edited to include the insertion sequence in the intended location. Integrative Genomic Viewer (IGV, Broad Institute) was used to verify that reads aligned to the modified sequence.
Variant analysis by genome-wide sequencing
To call variants (SNPs/indels) in the transfected lines compared to a previously sequenced control 3D7 line, genomic DNA reads were first trimmed to adapters and aligned to the Homo sapiens genome (assembly GRCh38) to remove human-mapped reads. Remaining reads were aligned to the P. falciparum 3D7 genome using bwa (version 0.7.17) and PCR duplicates were removed using PicardTools (Broad Institute). GATK HaplotypeCaller (https://gatk.broadinstitute.org/hc/en-us) was used to call variants between the sample and the 3D7 reference genome for both the transfected lines and the NF54 control. Only variants that were present in both transfected lines but not the NF54 control line were kept. We examined only coding-region variants and removed those that were synonymous variants or were located in var, rifin, or stevor genes. Quality control of variants was done by hard filtering using GATK guidelines.
Immunofluorescence assays
A double-staining immunofluorescence assay was used on NF54 control and transgenic parasite lines of mixed parasite population growth stages. Parasites were washed in incomplete medium prior to fixing onto coverslips with 4% paraformaldehyde for 20 min at RT under darkness. After fixation, samples were washed 3-5 times with 1x PBS followed by permeabilization with 0.5%Triton X-100 in PBS for 25 min at RT. Subsequently, samples were subjected to PBS washes and then incubated overnight at 4°C in blocking buffer (2mg/ml BSA solution in PBS containing 0.05% Tween-20 solution). Following overnight blocking, samples were washed and incubated for 1 hr at RT with Anti-HA Rb Ab (1:500, Abcam, ab9110) in blocking buffer. After primary Ab incubation, samples were subject to 3x washes with wash buffer (1x PBS containing 0.05% Tween-20) followed by incubation with anti-rabbit DyLight 550 (Abcam ab98489; 1: 500) secondary antibody for 1 h at room temperature. After incubation and series of washes with wash buffer, slides were incubated with anti-H3K9me3 antibody, Alexa Fluor 488 conjugate (Millipore 07–442-AF488; 1:100) for 1 hr at RT. Slides are then washed and mounted in Vectashield Antifade Mounting Medium with DAPI (Vector Laboratories, H-1200). Images were acquired using a Keyence BZ-X810 Fluorescence Microscope and were processed through ImageJ.
Western blotting
Parasites were synchronized twice at 8 hr. intervals between synchronizations. After one cycle, culture was washed 3x with complete media to remove residual aTC followed by dividing the cultures into two equal conditions. Parasites were grown with or without 500 nM aTC for 24 and 36 hrs. after which RBCs were lysed using 0.15 % saponin and parasites were collected after being washed with ice cold 1X PBS. Proteins were recovered from the lysed parasites after 30 min of incubation in lysis buffer (150 mM NaCl, 0.5 % NP40, 50 mM Tris pH 8, 1 mM EDTA and protease inhibitors) and 10 sec of sonication. Proteins were quantified with Pierce™ BCA Protein Assay Kit (Thermo Fisher 23227). 40 µg of proteins were loaded onto 3–8% precast polyacrylamide gel (Bio-Rad, 3450129). After migration, proteins were transferred onto a PVDF membrane and the membrane was blocked and then probed overnight with anti-HA tag antibody (1:2500, Abcam, ab9110). After primary Ab incubation, blots were subsequently washed 3x with a washing buffer followed by HRP-labeled Goat anti-Rabbit IgG (H + L) (1:10,000, NovexTM, A16104). ClarityTM Western ECL Substrate (Bio-Rad, 1705060) was applied to reveal the blots. Relative abundance of (+/-) aTC PfMORC was calculated by pixel density using ImageJ.
Immunoprecipitation followed by MudPIT mass spectrometry
Mid-to late-stage asexual parasites were collected following saponin treatment and purified samples were then resuspended into fresh IP buffer (50 mM Tris-HCl pH 7.5, 300 mM NaCl, 0.5 mM EDTA, 0.5 mM EGTA, 2 mM AEBSF 0.5 % Triton X-100, and EDTA-free protease inhibitor cocktail (Roche)). Post cell lysis solution was homogenized via sonication for 6-9 rounds. The soluble extracts were centrifuged at 13,000 x g for 15 min at 4°C. The lysates were precleared with Dynabeads™ Protein A (Invitrogen) for 1h at 4°C. Anti-HA tag antibody (1:2500, Abcam, ab9110) are added to control and HA tagged PfMORC precleared protein extract samples for 1hr at 4°C under constant rotation followed by the addition of fresh Dynabeads™ Protein A beads to each sample and incubated overnight at 4°C. Dynabeads were washed 3 times with 500 µL of buffer (PBS, 0.05% Tween-20). Proteins were eluted into an elution buffer (50 mM Tris-HCl pH 6.7, 100 mM DTT and 2% SDS). The eluent was subsequently precipitated overnight in 20% TCA followed by cold acetone washes. The urea-denatured, reduced, alkylated, and digested proteins were analyzed by Multidimensional Protein Identification Technology (MudPIT) on a Orbitrap Elite mass spectrometer coupled to an Agilent 1260 series HPLC, as described previously 69.
Proteomics data processing and analysis
Tandem mass (MS/MS) spectra were interpreted using ProluCID v.1.3.3 70 against a database consisting of 5527 non-redundant (NR) Plasmodium falciparum 3D7 proteins (PlasmoDB, v42), 36661 NR human proteins (NCBI, 2018-03-30 release), 419 common contaminants (human keratins, IgGs, and proteolytic enzymes), together with shuffled versions of all of these sequences. DTASelect v.1.9 71 and swallow v.0.0.1, an in-house developed software (https://github.com/tzw-wen/kite) were used to control FDRs resulting in protein FDRs less than 1.86%. All datasets were contrasted against their merged data set, respectively, using Contrast v1.9 71 and in-house developed sandmartin v.0.0.1 (https://github.com/tzw-wen/kite/tree/master/ kitelinux). Our in-house developed software, NSAF7 v.0.0.1 (https://github.com/tzw-wen/kite/tree/master/windowsapp/NSAF7x64), was used to generate spectral count-based label free quantitation results 72. QPROT 73 was used to calculate values of log2 fold change and Z-statistic to compare two replicate PfMORC affinity purifications to two negative controls. Proteins enriched in the PfMORC-HA samples with values of log2 fold change ≥ 2 and Z-statistic ≥ 5 were considered significantly enriched.
ChIP assay
PfMORC-HA asexual stage parasites were supplemented with and without aTC along with NF54 parasites (as a control) and were harvested at 24 and 36 HPS and cross linked with 1% formaldehyde for 10 min at 37°C followed by quenching with 150 mM glycine and 3 x washing with 1 x PBS. The pellets were resuspended in 1 mL of nuclear extraction buffer (10 mM HEPES, 10 mM KCl, 0.1 mM EDTA, 0.1 mM EGTA, 1 mM DTT, 0.5 mM AEBSF and 1x Protease inhibitor cocktail), and incubated on ice for 30 min before addition of NP-40/Igepal to a final of 0.25%. After lysis, the parasites are subject to homogenization by passing through a 26.5G needle/syringe to burst the nucleus. After centrifugation for 20 min at 2,500 x g at 4°C, samples were resuspended in 130 µL of shearing buffer (0.1% SDS, 1 mM EDTA, 10 mM Tris-HCl pH 7.5 and 1X Protease inhibitor cocktail) and transferred to a 130 µL Covaris sonication microtube. The samples were then sonicated using a Covaris S220 Ultrasonicator for 8 min (Duty cycle: 5%, intensity peak power:140, cycles per burst: 200, bath temperature: 6°C). The samples were transferred to ChIP dilution buffer (30 mM Tris-HCl pH 8.0, 3 mM EDTA, 0.1% SDS, 30 mM NaCl, 1.8% Triton X-100, 1X protease inhibitor tablet, 1X phosphatase inhibitor tablet) and centrifuged for 10 min at 13,000 rpm at 4°C, retaining the supernatant. For each sample, 13 μL of protein A agarose/salmon sperm DNA beads were washed three times with 500 µL ChIP dilution buffer by centrifuging for 1 min at 1,000 rpm at room temperature. For pre-clearing, the diluted chromatin samples were added to the beads and incubated for 1 hour at 4°C with rotation, then pelleted by centrifugation for 1 min at 1,000 rpm. Before adding antibody, ∼10% of each sample was taken as input. 2 µg of anti-HA tag antibodies (1:2500, Abcam, ab9110) were added to the sample and incubated overnight at 4°C with rotation. For each sample, rabbit IgG antibody (Cat.689 No. 10500c, Invitrogen) was used as a negative-control library. Per sample, 25 µL of protein A agarose/salmon sperm DNA beads were washed with ChIP dilution buffer (no inhibitors), blocked with 1 mg/mL BSA for 1 hour at 4°C, and then washed three more times with buffer. 25 µL of washed and blocked beads were added to the sample and incubated for 1 hour at 4°C with continuous mixing to collect the antibody/protein complex. Beads were pelleted by centrifugation for 1 min at 1,000 rpm at 4°C. The bead/antibody/protein complex was then washed with rotation using 1 mL of each buffer twice; low salt immune complex wash buffer (1% SDS, 1% Triton X-100, 2 mM EDTA, 20 mM Tris-HCl pH 8.0, 150 mM NaCl), high salt immune complex wash buffer (1% SDS, 1% Triton X-100, 2 mM EDTA, 20 mM Tris-HCl pH 8.0, 500 mM NaCl), and TE wash buffer (10 mM Tris-HCl pH 8.0, 1 mM EDTA). Complexes were eluted from antibodies by adding 250 μL of freshly prepared elution buffer (1% SDS, 0.1 M sodium bicarbonate). 5 M NaCl were added to the eluate and cross-linking was reversed by heating at 45°C overnight followed by addition of 15 μL of 20 mg/mL RNase A with 30 min incubation at 37°C. After this, 10 μL of 0.5 M EDTA, 20 μL of 1 M Tris-HCl pH 7.5, and 2 μL of 20 mg/mL proteinase K were added to the eluate and incubated for 2 hours at 45°C. DNA was recovered by phenol/chloroform extraction and ethanol precipitation. DNA was purified using Agencourt AMPure XP beads. Libraries were then prepared from this DNA using KAPA LTP Library Preparation Kit (Roche, KK8230) and sequenced on a NovaSeq 6000 machine.
ChIP-seq analysis
FastQC (version 0.11.8, https://www.bioinformatics.babraham.ac.uk/projects/fastqc/) was used to analyze raw read quality. Any adapter sequences were removed using Trimmomatic (version 0.39, http://www.usadellab.org/cms/?page=trimmomatic). Bases with Phred quality scores below 20 were trimmed using Sickle (version 1.33, https://asexual stage parasitesb.com/najoshi/sickle). The resulting reads were mapped against the P. falciparum genome (v48) using Bowtie2 (version 2.4.4). Using Samtools (version 1.11), only properly paired reads with mapping quality 40 or higher were retained and reads marked as PCR duplicates were removed by PicardTools MarkDuplicates (version 2.18.0, Broad Institute). Genome-wide read counts per nucleotide were normalized by dividing millions of mapped reads for each sample (for all samples including input) and subtracting IgG read counts from the anti-HA IP counts. Peak calling was performed using MACS3 (version 3.0.0a7) 74 with the upper and lower limits for fold enrichment set to 2 and 50, respectively.
Phenotypic analyses
P. falciparum NF54 line along with transgenic PfMORC-HA-TetR-DOZI line were synchronized, grown for one cycle post-synchronization and subject to 3 consecutive washes with 1X PBS before being supplemented either with or without (aTC) for 48 hrs. or 72 hrs. (depending on which stage the parasites were washed). Culture media were replaced daily with or without aTC supplementation.
Quantitative Growth assay
Synchronous cultures either ring or trophozoite stage at 0.5 % parasitemia were grown with or without aTC (500 nM) into a 96 well plate in triplicates for each time point. Parasites were collected at every 24 hrs. time for 3 cycles and subjected to SYBR green assay (Thermo Fisher, S7523).
RNA-seq library preparation
Parasites at ring, trophozoite or schizont stage were extracted following saponin treatment before flash freezing. Two independent biological replicates were generated for each time point, culture condition and line. Total RNA was extracted with TRIzol® LS Reagent (Invitrogen) followed by incubation for 1 hr with 4 units of DNase I (NEB) at 37°C. RNA samples were visualized by RNA electrophoresis and quantified on Synergy™ HT (BioTek). mRNA was then purified using NEBNext® Poly(A) mRNA Magnetic Isolation Module (NEB, E7490) according to the manufacturer’s instructions. Libraries were prepared using NEBNext® Ultra™ Directional RNA Library Prep Kit (NEB, E7420) and amplified by PCR with KAPA HiFi HotStart Ready Mix (Roche). PCR conditions consisted of 15 min at 37°C followed by 12 cycles of 98°C for 30 sec, 55°C for 10 s and 62°C for 1 min 15 sec, then finally one cycle for 5 min at 62°C. The quantity and quality of the final libraries were assessed using a Bioanalyzer (Agilent Technology Inc). Libraries were sequenced using a NovaSeq 6000 DNA sequencer (Illumina), producing paired-end 100-bp reads.
RNA-seq data processing and differential expression analysis
FastQC [https://www.bioinformatics.babraham.ac.uk/projects/fastqc/] was used to analyze raw read quality and thus 11 bp of each read and any adapter sequences were removed using Trimmomatic (v0.39) [http://www.usadellab.org/cms/?page=trimmomatic]. Bases were trimmed from reads using Sickle with a Phred quality threshold of 25 (v1.33) [https://github.com/najoshi/sickle] and reads shorter than 18 bp were removed. The resulting reads were mapped against the P. falciparum 3D7 genome (PlasmoDB, v53) using HISAT2 (v2-2.2.1) with default parameters. Uniquely mapped, properly paired reads with mapping quality of 40 or higher were retained using SAMtools (v1.11) [http://samtools.sourceforge.net/]. Raw read counts were determined for each gene in the P. falciparum genome using BedTools [https://bedtools.readthedocs.io/en/latest/#] to intersect the aligned reads with the genome annotation. Differential expression analysis was performed using DESeq2 to call up-and down-regulated genes (FDR < 0.05 and log2 FC > 0.5). Volcano plots were made using the R package EnhancedVolcano.
Hi-C library preparation
(+/-) aTC parasites at 24 hpi and 36 hpi were cross linked using 1.25% formaldehyde for 15 minutes at 37°C in 10 mL total volume followed by quenching with 150 mM glycine. Parasites were then washed three times with chilled 1x PBS on a rocking platform. Parasite nuclei were released using lysis buffer (10 mM Tris-HCl, pH 8.0, 10 mM NaCl, 2 mM AEBSF, 0.25% Igepal CA-630, and 1X EDTA-free protease inhibitor cocktail (Roche)) and 15 syringe passages through a 26.5-gauge needle. After releasing the crosslinked chromatin with 0.1% SDS, chromatin was digested using MboI restriction enzyme overnight at 37°C. Thereafter Hi-C libraries were prepared as previously described 27,75.
Hi-C data processing, differential interaction analysis and generation of 3D models
Paired-end Hi-C libraries were processed (pairing, mapping, quality filtering, binning, and normalization) using HiC-Pro 55. A mapping quality cutoff of 30 was set while aligning to the P. falciparum genome (PlasmoDB, v58) and the resulting reads were binned at 10 kb resolution and ICED-normalized. Stratum-adjusted correlation coefficients were calculated using HiCRep 76, then replicates were merged to generate a single consensus sample for each time point and condition. The normalized matrices were per-million read count normalized and maximum values for the heatmap color scale on each chromosome was set to the maximum value for all samples to allow for direct comparisons between each condition and time point. Furthermore, due to the overrepresentation of short-range interactions, the maximum values for each heatmap were also set to the 90th percentile of each chromosome to slightly compress highly interacting regions and enhance visualization of interactions. Significant intrachromosomal interactions were identified with FitHiC 77 to identify the log-linear relationship between contact probability and genomic distance. Differential interchromosomal and intrachromosomal interactions between (+) aTC and (-) aTC were identified using Selfish 57 with default parameters (FDR < 0.05). The maximum and minimum values for the color scales in the differential heatmaps were set to the absolute value of the largest log2 fold change within each chromosome. PASTIS 58 was used to generate coordinate matrices from the raw read count matrices output by HiC-Pro, and then visualized as 3D chromatin models in ChimeraX 78, highlighting regions containing var genes, telomeres, and centromeres.
Statistical analyses
Parasitemia and proportion of asexual stages were analyzed using a two-way ANOVA with Tukey’s test for multiple comparisons. Significant differences were indicated as following: * for p < 0.05; ** for p < 0.01, *** for p < 0.001 and **** for p < 0.0001. Statistical tests were performed with GraphPad Prism. Figures were generated with GraphPad Prism and BioRender.
Data Availability
Sequence reads for all sequencing experiments have been deposited in the NCBI Sequence Read Archive with accession PRJNA994684. Original data underlying this manuscript generated at the Stowers Institute can be accessed from the Stowers Original Data Repository at http://www.stowers.org/research/publications/LIBPB-2404. The mass spectrometry dataset generated for this study is available from the MassIVE data repository (ftp://MSV000092353@massive.ucsd.edu) using the identifier MSV000092353. All other data are available in the main text, methods, and supplementary data.
Acknowledgements
This work was supported by NIH grants to KGLR (nos. 1R01 AI136511 and R21 AI142506-01) and by the University of California, Riverside to KGLR (no. NIFA-Hatch-225935). This publication includes data generated at the UC San Diego IGM Genomics Center utilizing an Illumina NovaSeq 6000 that was purchased with funding from a National Institutes of Health SIG grant (#S10 OD026929).
Declarations of Interest
The authors declare no competing interests.
Supplementary Data
Supplementary data 1. PfMORC plasmid constructs to generate transgenic lines. All constructs were performed through listed primers, restriction sites and gRNAs.
Supplementary data 2. Mapped Whole Genome sequencing results of PfMORC transfectants.
Supplementary data 3. (a) Proteins identified by MudPIT analysis after MORC-HA immunoprecipitation. PfMORC (PF3D7_1468100) is highlighted in yellow and significantly purified proteins are in red. The filters used are QPROT log2 fold change > 2 and Z statistic > 5. Detailed protein list is provided in an additional sheet. (b) Raw data showing proteins identified by MudPIT analysis after PfMORC-HA immunoprecipitation.
Supplementary data 4. Peak calling result of PfMORC ChIP-seq. Results obtained at (a) ring stages (b) trophozoite stages and (c) schizont stages of cell progression.
Supplementary data 5. Parasite Survival Assay illustrating effect of PfMORC down-regulation on parasites at ring stage (top) or trophozoite stage (bottom) cell cycle. Experiments were conducted in triplicates with parasitemia quantified via Relative Fluorescence Units (RFUs) obtained by SYBR green assay for each time point. Significance derived through 2-way ANOVA.
Supplementary data 6. (a) RNA-seq read counts of PfMORC at 24 and 36 hours. DEseq2 analysis of differentially expressed genes of PfMORC at (b) 24 hours and (c) 36 hrs.
Figure S1. Expression of PfMORC-HA in WT control and PfMORC-HA lines. Western blot image revealed expression of PfMORC-HA protein in cultures using anti-HA antibody.
Figure S2. Correlation, peak calling, and gene family specific coverage of PfMORC ChIP-seq. (a) Heatmap indicating correlation between time points and replicates. (b) Overlap of called peaks between time points. (c) PfMORC coverage of select var genes from 500 bp 5’ of the TSS to 500 bp 3’ of the coding region. (d) PfMORC coverage of three highly expressed msp genes from 500 bp 5’ of the TSS to 500 bp 3’ of the coding region.
Figure S3. Expression of PfMORC-HA in PfMORC-HA-TetR-DOZI lines. Western blot image revealed expression of PfMORC-HA protein in samples extracted at 24 HPS or 36 HPS in +/-aTC conditions using anti-HA antibody.
Figure S4. RNA-seq correlation heatmap. (a) Heatmap indicating correlation between (+/-) aTC condition, time points and replicates.
Figure S5. Hi-C correlation analyses. (a) Stratum-adjusted correlation between (+/-) aTC condition, time points and replicates. (b) Negative log-linear relationship between contact probability and genomic distance.
Figure S6-9. Hi-C interaction heatmaps binned at 10kb resolution. Each sample includes intrachromosomal contact count heatmaps for all 14 chromosomes within the P. falciparum genome and a genome-wide interchromosomal interaction heatmap. Data for each heatmap is ICED and per-million read count normalized, and the maximum y-value for each heatmap is set to the highest value within the given chromosome for all samples. Antigenic gene-containing regions (red) and centromeres (gray) are highlighted.
Figure S10-S11. Differential Hi-C interaction heatmaps binned at 10kb resolution. Heatmaps generated from differential interaction matrices, identifying regions with a positive (red) or negative (blue) log2 fold change between (-) aTC and (+) aTC. Data for each heatmap is ICED and per-million read count normalized, and the maximum y-value for each heatmap is set to the highest value within the given chromosome for all samples. Antigenic gene-containing regions (red) and centromeres (gray) are highlighted.
References
- 1WHO. World malaria report 2021. (2021).
- 2Comparative genomics of transcriptional control in the human malaria parasite Plasmodium falciparumGenome Res 14:1548–1554https://doi.org/10.1101/gr.2218604
- 3Discovery of the principal specific transcription factors of Apicomplexa and their implication for the evolution of the AP2-integrase DNA binding domainsNucleic Acids Res 33:3994–4006https://doi.org/10.1093/nar/gki709
- 4Specific DNA-binding by apicomplexan AP2 transcription factorsProc Natl Acad Sci U S A 105:8393–8398https://doi.org/10.1073/pnas.0801993105
- 5Identification of a transcription factor in the mosquito-invasive stage of malaria parasitesMol Microbiol 71:1402–1414https://doi.org/10.1111/j.1365-2958.2009.06609.x
- 6Identification of an AP2-family protein that is critical for malaria liver stage developmentPLoS One 7https://doi.org/10.1371/journal.pone.0047557
- 7A cascade of DNA-binding proteins for sexual commitment and development in PlasmodiumNature 507:253–257https://doi.org/10.1038/nature12970
- 8A transcriptional switch underlies commitment to sexual development in malaria parasitesNature 507:248–252https://doi.org/10.1038/nature12920
- 9Cooperative binding of ApiAP2 transcription factors is crucial for the expression of virulence genes in Toxoplasma gondiiNucleic Acids Res 46:6057–6068https://doi.org/10.1093/nar/gky373
- 10Transcription factor AP2-Sp and its target genes in malarial sporozoitesMol Microbiol 75:854–863https://doi.org/10.1111/j.1365-2958.2009.07005.x
- 11A Knockout Screen of ApiAP2 Genes Reveals Networks of Interacting Transcriptional Regulators Controlling the Plasmodium Life CycleCell Host Microbe 21:11–22https://doi.org/10.1016/j.chom.2016.12.003
- 12Multiple regulatory roles of AP2/ERF transcription factor in angiospermBot Stud 58https://doi.org/10.1186/s40529-016-0159-1
- 13Female-specific gene regulation in malaria parasites by an AP2-family transcription factorMol Microbiol 113:40–51https://doi.org/10.1111/mmi.14334
- 14Global transcriptional repression: An initial and essential step for Plasmodium sexual developmentProc Natl Acad Sci U S A 112:12824–12829https://doi.org/10.1073/pnas.1504389112
- 15Expression profiling of the schizont and trophozoite stages of Plasmodium falciparum with a long-oligonucleotide microarrayGenome Biol 4https://doi.org/10.1186/gb-2003-4-2-r9
- 16Discovery of gene function by expression profiling of the malaria parasite life cycleScience 301:1503–1508https://doi.org/10.1126/science.1087025
- 17The transcriptome of the intraerythrocytic developmental cycle of Plasmodium falciparumPLoS Biol 1https://doi.org/10.1371/journal.pbio.0000005
- 18Analysis of the Plasmodium and Anopheles transcriptomes during oocyst differentiationJ Biol Chem 279:5581–5587https://doi.org/10.1074/jbc.M307587200
- 19Genome-wide identification of genes upregulated at the onset of gametocytogenesis in Plasmodium falciparumMol Biochem Parasitol 143:100–110https://doi.org/10.1016/j.molbiopara.2005.04.015
- 20Full-Length Transcriptome Analysis of Plasmodium falciparum by Single-Molecule Long-Read SequencingFront Cell Infect Microbiol 11https://doi.org/10.3389/fcimb.2021.631545
- 21Capturing chromosome conformationScience 295:1306–1311https://doi.org/10.1126/science.1067799
- 22Chromosome Conformation Capture Carbon Copy (5C): a massively parallel solution for mapping interactions between genomic elementsGenome Res 16:1299–1309https://doi.org/10.1101/gr.5571506
- 23Histone lysine methyltransferases and demethylases in Plasmodium falciparumInt J Parasitol 38:1083–1097https://doi.org/10.1016/j.ijpara.2008.01.002
- 24Topological domains in mammalian genomes identified by analysis of chromatin interactionsNature 485:376–380https://doi.org/10.1038/nature11082
- 25Persistence and activation of malaria hypnozoites in long-term primary hepatocyte culturesNat Med 20:307–312https://doi.org/10.1038/nm.3461
- 26Bipartite structure of the inactive mouse X chromosomeGenome Biol 16https://doi.org/10.1186/s13059-015-0728-8
- 27Comparative 3D genome organization in apicomplexan parasitesProc Natl Acad Sci U S A 116:3183–3192https://doi.org/10.1073/pnas.1810815116
- 28The chromatin bound proteome of the human malaria parasiteMicrob Genom 6https://doi.org/10.1099/mgen.0.000327
- 29MORC proteins and epigenetic regulationPlant Signal Behav 7:1561–1565https://doi.org/10.4161/psb.22460
- 30The MORC family: new epigenetic regulators of transcription and DNA damage responseEpigenetics 8:685–693https://doi.org/10.4161/epi.24976
- 31MORC-1 Integrates Nuclear RNAi and Transgenerational Chromatin Architecture to Promote Germline ImmortalityDev Cell 41:408–423https://doi.org/10.1016/j.devcel.2017.04.023
- 32The PfAP2-G2 transcription factor is a critical regulator of gametocyte maturationMol Microbiol 115:1005–1024https://doi.org/10.1111/mmi.14676
- 33A MORC-driven transcriptional switch controls Toxoplasma developmental trajectories and sexual commitmentNat Microbiol 5:570–583https://doi.org/10.1038/s41564-020-0674-4
- 34MORC protein family-related signature within human disease and cancerCell Death Dis 12https://doi.org/10.1038/s41419-021-04393-1
- 35MORC proteins regulate transcription factor binding by mediating chromatin compaction in active chromatin regionsGenome Biol 24https://doi.org/10.1186/s13059-023-02939-4
- 36Histone-modifying complexes regulate gene expression pertinent to the differentiation of the protozoan parasite Toxoplasma gondiiMol Cell Biol 25:10301–10314https://doi.org/10.1128/MCB.25.23.10301-10314.2005
- 37Drug inhibition of HDAC3 and epigenetic control of differentiation in Apicomplexa parasitesJ Exp Med 206:953–966https://doi.org/10.1084/jem.20082826
- 38Landscape of the Plasmodium Interactome Reveals Both Conserved and Species-Specific FunctionalityCell Rep 28:1635–1647https://doi.org/10.1016/j.celrep.2019.07.019
- 39PfMORC confers melatonin dependent synchrony of the human malaria parasite P. falciparum in the asexual stageSci Rep 11https://doi.org/10.1038/s41598-021-81235-2
- 40Synthetic RNA-protein modules integrated with native translation mechanisms to control gene expression in malaria parasitesNat Commun 7https://doi.org/10.1038/ncomms10727
- 41GDV1 induces sexual commitment of malaria parasites by antagonizing HP1-dependent gene silencingScience 359:1259–1263https://doi.org/10.1126/science.aan6042
- 42Inducible developmental reprogramming redefines commitment to sexual development in the malaria parasite Plasmodium bergheiNat Microbiol 3:1206–1213https://doi.org/10.1038/s41564-018-0223-6
- 43Mechanisms of triggering malaria gametocytogenesis by AP2-GParasitol Int 84https://doi.org/10.1016/j.parint.2021.102403
- 44QPROT: Statistical method for testing differential expression using protein-level intensity data in label-free quantitative proteomicsJ Proteomics 129:121–126https://doi.org/10.1016/j.jprot.2015.07.036
- 45A major role for the Plasmodium falciparum ApiAP2 protein PfSIP2 in chromosome end biologyPLoS Pathog 6https://doi.org/10.1371/journal.ppat.1000784
- 46Genome-wide landscape of ApiAP2 transcription factors reveals a heterochromatin-associated regulatory network during Plasmodium falciparum blood-stage developmentNucleic Acids Res 50:3413–3431https://doi.org/10.1093/nar/gkac176
- 47Epigenetic reader complexes of the human malaria parasite, Plasmodium falciparumNucleic Acids Res 47:11574–11588https://doi.org/10.1093/nar/gkz1044
- 48Developmental selection of var gene expression in Plasmodium falciparumNature 394:392–395https://doi.org/10.1038/28660
- 49Antigenic variation in malaria: in situ switching, relaxed and mutually exclusive transcription of var genes during intra-erythrocytic development in Plasmodium falciparumEMBO J 17:5418–5426https://doi.org/10.1093/emboj/17.18.5418
- 50Probabilistic data integration identifies reliable gametocyte-specific proteins and transcripts in malaria parasitesSci Rep 8https://doi.org/10.1038/s41598-017-18840-7
- 51Uncovering the essential genes of the human malaria parasite Plasmodium falciparum by saturation mutagenesisScience 360https://doi.org/10.1126/science.aap7847
- 52Versatile control of Plasmodium falciparum gene expression with an inducible protein-RNA interactionNat Commun 5https://doi.org/10.1038/ncomms6329
- 53An integrated platform for genome engineering and gene expression perturbation in Plasmodium falciparumSci Rep 11https://doi.org/10.1038/s41598-020-77644-4
- 54Small nucleolar RNAs: an abundant group of noncoding RNAs with diverse cellular functionsCell 109:145–148https://doi.org/10.1016/s0092-8674(02)00718-3
- 55HiC-Pro: an optimized and flexible pipeline for Hi-C data processingGenome Biol 16https://doi.org/10.1186/s13059-015-0831-x
- 56Changes in genome organization of parasite-specific gene families during the Plasmodium transmission stagesNat Commun 9https://doi.org/10.1038/s41467-018-04295-5
- 57Selfish: discovery of differential chromatin interactions via a self-similarity measureBioinformatics 35:i145–i153https://doi.org/10.1093/bioinformatics/btz362
- 58A statistical approach for inferring the 3D structure of the genomeBioinformatics 30:i26–33https://doi.org/10.1093/bioinformatics/btu268
- 59An apicomplexan bromodomain proteinTgBD 1https://doi.org/10.1128/mbio.03573-22
- 60Molecular mechanism of the MORC4 ATPase activationNat Commun 11https://doi.org/10.1038/s41467-020-19278-8
- 61Involvement of histone deacetylation in MORC2-mediated down-regulation of carbonic anhydrase IXNucleic Acids Res 38:2813–2824https://doi.org/10.1093/nar/gkq006
- 62Structural insight into the zinc finger CW domain as a histone modification readerStructure 18:1127–1139https://doi.org/10.1016/j.str.2010.06.012
- 63MORC3 Forms Nuclear Condensates through Phase SeparationiScience 17:182–189https://doi.org/10.1016/j.isci.2019.06.030
- 64MORC4 Promotes Chemoresistance of Luminal A/B Breast Cancer via STAT3-Mediated MID2 UpregulationOnco Targets Ther 13:6795–6803https://doi.org/10.2147/OTT.S260509
- 65Two-step colocalization of MORC3 with PML nuclear bodiesJ Cell Sci 123:2014–2024https://doi.org/10.1242/jcs.063586
- 66BHuman malaria parasites in continuous culture. 1976. J Parasitol 91, 484-486 https://doi.org/10.1645/0022-3395(2005)091[0484:HMPICC]2.0.CO;2
- 67Cycloguanil and its parent compound proguanil demonstrate distinct activities against Plasmodium falciparum malaria parasites transformed with human dihydrofolate reductaseMol Pharmacol 54:1140–1147https://doi.org/10.1124/mol.54.6.1140
- 68Redesigned TetR-Aptamer System To Control Gene Expression in Plasmodium falciparummSphere 5https://doi.org/10.1128/mSphere.00457-20
- 69Proteomic analysis by multidimensional protein identification technologyMethods Mol Biol 328:159–175https://doi.org/10.1385/1-59745-026-X:159
- 70ProLuCID: An improved SEQUEST-like algorithm with enhanced sensitivity and specificityJ Proteomics 129:16–24https://doi.org/10.1016/j.jprot.2015.07.001
- 713rd. DTASelect and Contrast: tools for assembling and comparing protein identifications from shotgun proteomicsJ Proteome Res 1:21–26https://doi.org/10.1021/pr015504q
- 72Refinements to label free proteome quantitation: how to deal with peptides shared by multiple proteinsAnal Chem 82:2272–2281https://doi.org/10.1021/ac9023999
- 73Significance analysis of spectral count data in label-free shotgun proteomicsMol Cell Proteomics 7:2373–2385https://doi.org/10.1074/mcp.M800203-MCP200
- 74Model-based analysis of ChIP-Seq (MACS)Genome Biol 9https://doi.org/10.1186/gb-2008-9-9-r137
- 75Chromosomes Conformation Capture Coupled with Next-Generation Sequencing (Hi-C) in Plasmodium falciparumMethods Mol Biol 2369 :15–25https://doi.org/10.1007/978-1-0716-1681-9_2
- 76HiCRep: assessing the reproducibility of Hi-C data using a stratum-adjusted correlation coefficientGenome Res 27:1939–1949https://doi.org/10.1101/gr.220640.117
- 77Statistical confidence estimation for Hi-C data reveals regulatory chromatin contactsGenome Res 24:999–1011https://doi.org/10.1101/gr.160374.113
- 78UCSF ChimeraX: Meeting modern challenges in visualization and analysisProtein Sci 27:14–25https://doi.org/10.1002/pro.3235
Article and author information
Version history
- Preprint posted:
- Sent for peer review:
- Reviewed Preprint version 1:
Copyright
© 2023, Chahine et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
- views
- 300
- downloads
- 13
- citation
- 1
Views, downloads and citations are aggregated across all versions of this paper published by eLife.