snakePipes: facilitating flexible, scalable and integrative epigenomic analysisSUMMARY: Due to the rapidly increasing scale and diversity of epigenomic data, modular and scalable analysis workflows are of wide interest. Here we present snakePipes, a workflow package for processing and downstream analysis of data from common epigenomic assays: ChIP-seq, RNA-seq, Bisulfite-seq, ATAC-seq, Hi-C and single-cell RNA-seq. snakePipes enables users to assemble variants of each workflow and to easily install and upgrade the underlying tools, via its simple command-line wrappers and yaml files. AVAILABILITY AND IMPLEMENTATION: snakePipes can be installed via conda: `conda install -c mpi-ie -c bioconda -c conda-forge snakePipes'. Source code (https://github.com/maxplanck-ie/snakepipes) and documentation (https://snakepipes.readthedocs.io/en/latest/) are available online. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
Histone variant H2A.Z regulates zygotic genome activationDuring embryogenesis, the genome shifts from transcriptionally quiescent to extensively active in a process known as Zygotic Genome Activation (ZGA). In Drosophila, the pioneer factor Zelda is known to be essential for the progression of development; still, it regulates the activation of only a small subset of genes at ZGA. However, thousands of genes do not require Zelda, suggesting that other mechanisms exist. By conducting GRO-seq, HiC and ChIP-seq in Drosophila embryos, we demonstrate that up to 65% of zygotically activated genes are enriched for the histone variant H2A.Z. H2A.Z enrichment precedes ZGA and RNA Polymerase II loading onto chromatin. In vivo knockdown of maternally contributed Domino, a histone chaperone and ATPase, reduces H2A.Z deposition at transcription start sites, causes global downregulation of housekeeping genes at ZGA, and compromises the establishment of the 3D chromatin structure. We infer that H2A.Z is essential for the de novo establishment of transcriptional programs during ZGA via chromatin reorganization.
ELAV and FNE Determine Neuronal Transcript Signatures through EXon-Activated RescueInheritance of H3K9 methylation regulates genome architecture in Drosophila early embryosConstitutive heterochromatin is essential for transcriptional silencing and genome integrity. The establishment of constitutive heterochromatin in early embryos and its role in early fruitfly development are unknown. Lysine 9 trimethylation of histone H3 (H3K9me3) and recruitment of its epigenetic reader, heterochromatin protein 1a (HP1a), are hallmarks of constitutive heterochromatin. Here, we show that H3K9me3 is transmitted from the maternal germline to the next generation. Maternally inherited H3K9me3, and the histone methyltransferases (HMT) depositing it, are required for the organization of constitutive heterochromatin: early embryos lacking H3K9 methylation display de-condensation of pericentromeric regions, centromere-centromere de-clustering, mitotic defects, and nuclear shape irregularities, resulting in embryo lethality. Unexpectedly, quantitative CUT&Tag and 4D microscopy measurements of HP1a coupled with biophysical modeling revealed that H3K9me2/3 is largely dispensable for HP1a recruitment. Instead, the main function of H3K9me2/3 at this developmental stage is to drive HP1a clustering and subsequent heterochromatin compaction. Our results show that HP1a binding to constitutive heterochromatin in the absence of H3K9me2/3 is not sufficient to promote proper embryo development and heterochromatin formation. The loss of H3K9 HMTs and H3K9 methylation alters genome organization and hinders embryonic development.
ELAV mediates circular RNA biogenesis in neuronsCircular RNAs (circRNAs) arise from back-splicing of precursor RNAs and accumulate in the nervous systems of animals, where they are thought to regulate gene expression and synaptic function. Here, we show that neuronal circRNA biosynthesis is mediated by the pan-neuronal RNA-binding protein ELAV. In Drosophila embryos, we characterized the circRNA landscape in normal and elav mutant neurons. We found that neuronal circRNAs are globally (>75%) depleted upon ELAV knockout, and induction of ELAV expression drives ectopic RNA circularization. In brain tissue, ELAV binds to pre-mRNA introns flanking putative circRNAs and decreases efficiency of linear splicing in favor of intron pairing at reverse complementary matches, inducing circularization. Together, our data demonstrate that ELAV directly modulates splicing decisions to generate the neuronal circRNA landscape.