The malaria parasite Plasmodium falciparum has a nonphotosynthetic plastid called the apicoplast, which contains its own genome. Regulatory mechanisms for apicoplast gene expression remain poorly understood, despite this organelle being crucial for the parasite life cycle. Here, we identify a nuclear-encoded apicoplast RNA polymerase σ subunit (sigma factor) which, along with the α subunit, appears to mediate apicoplast transcript accumulation. This has a periodicity reminiscent of parasite circadian or developmental control. Expression of the apicoplast subunit gene, apSig, together with apicoplast transcripts, increased in the presence of the blood circadian signaling hormone melatonin. Our data suggest that the host circadian rhythm is integrated with intrinsic parasite cues to coordinate apicoplast genome transcription. This evolutionarily conserved regulatory system might be a future target for malaria treatment.
De novo variants (DNVs) cause many genetic diseases. When DNVs are examined in the whole coding regions of genes in next-generation sequencing analyses, pathogenic DNVs often cluster in a specific region. One such region is the last exon and the last 50 bp of the penultimate exon, where truncating DNVs cause escape from nonsense-mediated mRNA decay [NMD(-) region]. Such variants can have dominant-negative or gain-of-function effects. Here, we first developed a resource of rates of truncating DNVs in NMD(-) regions under the null model of DNVs. Utilizing this resource, we performed enrichment analysis of truncating DNVs in NMD(-) regions in 346 developmental and epileptic encephalopathy (DEE) trios. We observed statistically significant enrichment of truncating DNVs in semaphorin 6B (SEMA6B) (p value: 2.8 × 10-8; exome-wide threshold: 2.5 × 10-6). The initial analysis of the 346 individuals and additional screening of 1,406 and 4,293 independent individuals affected by DEE and developmental disorders collectively identified four truncating DNVs in the SEMA6B NMD(-) region in five individuals who came from unrelated families (p value: 1.9 × 10-13) and consistently showed progressive myoclonic epilepsy. RNA analysis of lymphoblastoid cells established from an affected individual showed that the mutant allele escaped NMD, indicating stable production of the truncated protein. Importantly, heterozygous truncating variants in the NMD(+) region of SEMA6B are observed in general populations, and SEMA6B is most likely loss-of-function tolerant. Zebrafish expressing truncating variants in the NMD(-) region of SEMA6B orthologs displayed defective development of brain neurons and enhanced pentylenetetrazole-induced seizure behavior. In summary, we show that truncating DNVs in the final exon of SEMA6B cause progressive myoclonic epilepsy.