39 40 Malaria pathogenesis results from the asexual replication of Plasmodium falciparum 41 within human red blood cells, which relies on a precisely timed cascade of gene 42 expression over a 48-hour life cycle. Although substantial post-transcriptional 43 regulation of this hardwired program has been observed, it remains unclear how 44 these processes are mediated on a transcriptome-wide level. To this end, we 45 identified mRNA modifications in the P. falciparum transcriptome and performed a 46 comprehensive characterization of N 6 -methyladenosine (m 6 A) over the course of 47 blood stage development. Using mass spectrometry and m 6 A RNA sequencing, we 48 demonstrate that m 6 A is highly developmentally regulated, exceeding m 6 A levels 49 known in any other eukaryote. We identify an evolutionarily conserved m 6 A writer 50 complex and show that knockdown of the putative m 6 A methyltransferase by 51 CRISPR interference leads to increased levels of transcripts that normally contain 52 m 6 A. In accordance, we find an inverse correlation between m 6 A status and mRNA 53 stability or translational efficiency. Our data reveal the crucial role of extensive m 6 A 54 mRNA methylation in dynamically fine-tuning the transcriptional program of a 55 unicellular eukaryote as well as a new 'epitranscriptomic' layer of gene regulation in 56 malaria parasites. 57 58 59 Keywords 60 61 Plasmodium falciparum development; mRNA methylation; N 6 -methyladenosine; 62 mRNA degradation; translational efficiency; CRISPR interference 63 64 65 66 Malaria, a mosquito-borne human disease caused by the unicellular apicomplexan 67 parasite Plasmodium falciparum, remains a major global health threat (World Health 68 Organization, 2018). Human pathogenesis results from the 48-hour intra-erythrocytic 69 developmental cycle (IDC), during which each parasite undergoes schizogony within 70 red blood cells (RBC) to create up to 32 new daughter cells (Fig. 1a). Each step of 71 the IDC -RBC invasion (0-hour post infection [hpi]), host cell remodeling, genome 72replication, and eventually RBC egress -is effected by a precisely timed cascade of 73 gene expression. Throughout this process, mRNAs from the majority of genes 74 expressed reach peak abundance only once, which is thought to correspond to the 75 time point when its encoded product is most required 1 .
76At the transcriptional level of gene regulation, silencing is achieved via 77 heterochromatinization of genes encoding variant surface antigens 2 and the inducer 78 of sexual stage development 3 . Heterochromatin, however, is restricted to 79 subtelomeric and several central chromosomal regions 4 , and the large majority of 80 the P. falciparum genome remains in a euchromatic and transcriptionally permissive 81 state throughout the IDC 5 . Recent studies of nucleosome occupancy have shown 82 that dynamic chromatin accessibility in promoter regions strongly correlates with 83 mRNA abundance, and that these sequences can be recognized by transcription 84 factors in vitro 6-8 . Intriguingly, ...