Plasticity in patterns of histone modifications and chromosomal proteins in Drosophila heterochromatin

Three-Dimensional Genome Organization and Function in Drosophila

Understanding how the metazoan genome is used during development and cell differentiation is one of the major challenges in the postgenomic era. Early studies in Drosophila suggested that three-dimensional (3D) chromosome organization plays important regulatory roles in this process and recent technological advances started to reveal connections at the molecular level. Here we will consider general features of the architectural organization of the Drosophila genome, providing historical perspective and insights from recent work. We will compare the linear and spatial segmentation of the fly genome and focus on the two key regulators of genome architecture: insulator components and Polycomb group proteins. With its unique set of genetic tools and a compact, well annotated genome, Drosophila is poised to remain a model system of choice for rapid progress in understanding principles of genome organization and to serve as a proving ground for development of 3D genome-engineering techniques.

A single double-strand break system reveals repair dynamics and mechanisms in heterochromatin and euchromatin

Janssen et al. developed an in vivo single double-strand break (DSB) system for both heterochromatic and euchromatic loci in Drosophila melanogaster. Live imaging and sequence analysis of repair products reveal that DSBs in euchromatin and heterochromatin are repaired with similar kinetics, employ both NHEJ and HR, and can use homologous chromosomes as an HR template.

Repair of DNA double-strand breaks (DSBs) must be properly orchestrated in diverse chromatin regions to maintain genome stability. The choice between two main DSB repair pathways, nonhomologous end-joining (NHEJ) and homologous recombination (HR), is regulated by the cell cycle as well as chromatin context.

Pericentromeric heterochromatin forms a distinct nuclear domain that is enriched for repetitive DNA sequences that pose significant challenges for genome stability. Heterochromatic DSBs display specialized temporal and spatial dynamics that differ from euchromatic DSBs. Although HR is thought to be the main pathway used to repair heterochromatic DSBs, direct tests of this hypothesis are lacking. Here, we developed an in vivo single DSB system for both heterochromatic and euchromatic loci in Drosophila melanogaster. Live imaging of single DSBs in larval imaginal discs recapitulates the spatio–temporal dynamics observed for irradiation (IR)-induced breaks in cell culture. Importantly, live imaging and sequence analysis of repair products reveal that DSBs in euchromatin and heterochromatin are repaired with similar kinetics, employ both NHEJ and HR, and can use homologous chromosomes as an HR template. This direct analysis reveals important insights into heterochromatin DSB repair in animal tissues and provides a foundation for further explorations of repair mechanisms in different chromatin domains.

Regulatory functions and chromatin loading dynamics of linker histone H1 during endoreplication in Drosophila

Eukaryotic DNA replicates asynchronously, with discrete genomic loci replicating during different stages of S phase. Drosophila larval tissues undergo endoreplication without cell division, and the latest replicating regions occasionally fail to complete endoreplication, resulting in underreplicated domains of polytene chromosomes. Here we show that linker histone H1 is required for the underreplication (UR) phenomenon in Drosophila salivary glands. H1 directly interacts with the Suppressor of UR (SUUR) protein and is required for SUUR binding to chromatin in vivo. These observations implicate H1 as a critical factor in the formation of underreplicated regions and an upstream effector of SUUR. We also demonstrate that the localization of H1 in chromatin changes profoundly during the endocycle. At the onset of endocycle S (endo-S) phase, H1 is heavily and specifically loaded into late replicating genomic regions and is then redistributed during the course of endoreplication. Our data suggest that cell cycle-dependent chromosome occupancy of H1 is governed by several independent processes. In addition to the ubiquitous replication-related disassembly and reassembly of chromatin, H1 is deposited into chromatin through a novel pathway that is replication-independent, rapid, and locus-specific. This cell cycle-directed dynamic localization of H1 in chromatin may play an important role in the regulation of DNA replication timing.

The Drosophila speciation factor HMR localizes to genomic insulator sites

Hybrid incompatibility between Drosophila melanogaster and D. simulans is caused by a lethal interaction of the proteins encoded by the Hmr and Lhr genes. In D. melanogaster the loss of HMR results in mitotic defects, an increase in transcription of transposable elements and a deregulation of heterochromatic genes. To better understand the molecular mechanisms that mediate HMR’s function, we measured genome-wide localization of HMR in D. melanogaster tissue culture cells by chromatin immunoprecipitation. Interestingly, we find HMR localizing to genomic insulator sites that can be classified into two groups. One group belongs to gypsy insulators and another one borders HP1a bound regions at active genes. The transcription of the latter group genes is strongly affected in larvae and ovaries of Hmr mutant flies. Our data suggest a novel link between HMR and insulator proteins, a finding that implicates a potential role for genome organization in the formation of species.

upSET, the Drosophila homologue of SET3, Is Required for Viability and the Proper Balance of Active and Repressive Chromatin Marks

Chromatin plays a critical role in faithful implementation of gene expression programs. Different post-translational modifications (PTMs) of histone proteins reflect the underlying state of gene activity, and many chromatin proteins write, erase, bind, or are repelled by, these histone marks. One such protein is UpSET, the Drosophila homolog of yeast Set3 and mammalian KMT2E (MLL5). Here, we show that UpSET is necessary for the proper balance between active and repressed states. Using CRISPR/Cas-9 editing, we generated S2 cells that are mutant for upSET We found that loss of UpSET is tolerated in S2 cells, but that heterochromatin is misregulated, as evidenced by a strong decrease in H3K9me2 levels assessed by bulk histone PTM quantification. To test whether this finding was consistent in the whole organism, we deleted the upSET coding sequence using CRISPR/Cas-9, which we found to be lethal in both sexes in flies. We were able to rescue this lethality using a tagged upSET transgene, and found that UpSET protein localizes to transcriptional start sites (TSS) of active genes throughout the genome. Misregulated heterochromatin is apparent by suppressed position effect variegation of the wm4 allele in heterozygous upSET-deleted flies. Using nascent-RNA sequencing in the upSET-mutant S2 lines, we show that this result applies to heterochromatin genes generally. Our findings support a critical role for UpSET in maintaining heterochromatin, perhaps by delimiting the active chromatin environment.

Nuclear Dynamics of Heterochromatin Repair

Repairing double-strand breaks (DSBs) is particularly challenging in pericentromeric heterochromatin, where the abundance of repeated sequences exacerbates the risk of ectopic recombination and chromosome rearrangements. Recent studies in Drosophila cells revealed that faithful homologous recombination (HR) repair of heterochromatic DSBs relies on the relocalization of DSBs to the nuclear periphery before Rad51 recruitment. We summarize here the exciting progress in understanding this pathway, including conserved responses in mammalian cells and surprising similarities with mechanisms in yeast that deal with DSBs in distinct sites that are difficult to repair, including other repeated sequences. We will also point out some of the most important open questions in the field and emerging evidence suggesting that deregulating these pathways might have dramatic consequences for human health.

Role of Chromatin Modifications in Drosophila Germline Stem Cell Differentiation

During Drosophila oogenesis, germline stem cells (GSCs) self-renew and differentiate to give rise to a mature egg. Self-renewal and differentiation of GSCs are regulated by both intrinsic mechanisms such as regulation of gene expression in the germ line and extrinsic signaling pathways from the surrounding somatic niche. Epigenetic mechanisms, including histone-modifying proteins, nucleosome remodeling complexes, and histone variants, play a critical role in regulating intrinsic gene expression and extrinsic signaling cues from the somatic niche. In the GSCs, intrinsic epigenetic modifiers are required to maintain a stem cell fate by promoting expression of self-renewal factors and repressing the differentiation program. Subsequently, in the GSC daughters, epigenetic regulators activate the differentiation program to promote GSC differentiation. During differentiation, the GSC daughter undergoes meiosis to give rise to the developing egg, containing a compacted chromatin architecture called the karyosome. Epigenetic modifiers control the attachment of chromosomes to the nuclear lamina to aid in meiotic recombination and the release from the lamina for karyosome formation. The germ line is in close contact with the soma for the entirety of this developmental process. This proximity facilitates signaling from the somatic niche to the developing germ line. Epigenetic modifiers play a critical role in the somatic niche, modulating signaling pathways in order to coordinate the transition of GSC to an egg. Together, intrinsic and extrinsic epigenetic mechanisms modulate this exquisitely balanced program.

α-Synuclein enhances histone H3 lysine-9 dimethylation and H3K9me2-dependent transcriptional responses

α-Synuclein (αS) is a protein linked to Parkinson’s disease (PD) and related neurodegenerative disorders. It is mostly localized within synapses, but αS has also been suggested to play a role in the nucleus. We used transgenic Drosophila and inducible SH-SY5Y neuroblastoma cells to investigate the effects of αS on chromatin with a particular focus on histone modifications. Overexpression of αS in male flies as well as in retinoic acid pre-treated neuroblastoma cells led to an elevation of histone H3K9 methylations, mostly mono- (H3K9me1) and di- (H3K9me2). The transient increase of H3K9 methylation in αS-induced SH-SY5Y cells was preceded by mRNA induction of the euchromatic histone lysine N-methyltransferase 2 (EHMT2). EHMT2 and H3K9me2 can function within the REST complex. Chromatin immunoprecipitation (ChIP) analyses of selected candidate, REST regulated genes showed significantly increased H3K9me2 promoter occupancy of genes encoding the L1CAM cell adhesion molecule and the synaptosomal-associated protein SNAP25, whose reduced expression levels were confirmed by RT-qPCR in αS induced cells. Treatment with EHMT inhibitor UNC0638 restored the mRNA levels of L1CAM and SNAP25. Thus, αS overexpression enhances H3K9 methylations via ΕΗΜΤ2 resulting in elevated H3K9me2 at the SNAP25 promoter, possibly affecting SNARE complex assembly and hence synaptic vesicle fusion events regulated by αS.

A decade of understanding spatio-temporal regulation of DNA repair by the nuclear architecture

The nucleus is a hub for gene expression and is a highly organized entity. The nucleoplasm is heterogeneous, owing to the preferential localization of specific metabolic factors, which lead to the definition of nuclear compartments or bodies. The genome is organized into chromosome territories, as well as heterochromatin and euchromatin domains. Recent observations have indicated that nuclear organization is important for maintaining genomic stability. For example, nuclear organization has been implicated in stabilizing damaged DNA, repair-pathway choice, and in preventing chromosomal rearrangements. Over the past decade, several studies have revealed that dynamic changes in the nuclear architecture are important during double-strand break repair. Stemming from work in yeast, relocation of a damaged site prior to repair appears to be at least partially conserved in multicellular eukaryotes. In this review, we will discuss genome and nucleoplasm architecture, particularly the importance of the nuclear periphery in genome stability. We will also discuss how the site of relocation regulates repair-pathway choice.

Antagonistic roles of Drosophila Tctp and Brahma in chromatin remodelling and stabilizing repeated sequences

Genome stability is essential for all organisms. Translationally controlled tumour protein (TCTP) is a conserved protein associated with cancers. TCTP is involved in multiple intracellular functions, but its role in transcription and genome stability is poorly understood. Here, we demonstrate new functions of Drosophila TCTP (Tctp) in transcription and the stability of repeated sequences (rDNA and pericentromeric heterochromatin). Tctp binds Brahma (Brm) chromatin remodeler to negatively modulate its activity. Tctp mutants show abnormally high levels of transcription in a large set of genes and transposons. These defects are ameliorated by brm mutations. Furthermore, Tctp promotes the stability of repeated sequences by opposing the Brm function. Additional regulation of pericentromeric heterochromatin by Tctp is mediated by su(var)3-9 transcriptional regulation. Altogether, Tctp regulates transcription and the stability of repeated sequences by antagonizing excess Brm activity. This study provides insights into broader nuclear TCTP functions for the maintenance of genome stability.

Genome stability is important for normal cellular function. Here, Hong and Choi show that translationally controlled tumour protein (TCTP) in Drosophila regulates pericentromeric chromatin remodelling and transcription via negatively regulating a chromatin remodeler Brahma.

Cervantes and Quijote protect heterochromatin from aberrant recombination and lead the way to the nuclear periphery

Repairing double-strand breaks (DSBs) is particularly challenging in heterochromatin, where the abundance of repeated sequences exacerbates the risk of ectopic recombination and chromosome rearrangements. In Drosophila cells, faithful homologous recombination (HR) repair of heterochromatic DSBs relies on a specialized pathway that relocalizes repair sites to the nuclear periphery before Rad51 recruitment. Here we show that HR progression is initially blocked inside the heterochromatin domain by SUMOylation and the coordinated activity of two distinct Nse2 SUMO E3 ligases: Quijote (Qjt) and Cervantes (Cerv). In addition, the SUMO-targeted ubiquitin ligase (STUbL) Dgrn, but not its partner dRad60, is recruited to heterochromatic DSBs at early stages of repair and mediates relocalization. However, Dgrn is not required to prevent Rad51 recruitment inside the heterochromatin domain, suggesting that the block to HR progression inside the domain and relocalization to the nuclear periphery are genetically separable pathways. Further, SUMOylation defects affect relocalization without blocking heterochromatin expansion, revealing that expansion is not required for relocalization. Finally, nuclear pores and inner nuclear membrane proteins (INMPs) anchor STUbL/RENi components and repair sites to the nuclear periphery, where repair continues. Together, these studies reveal a critical role of SUMOylation and nuclear architecture in the spatial and temporal regulation of heterochromatin repair and the protection of genome integrity.

Chromatin-modifying genetic interventions suppress age-associated transposable element activation and extend life span in Drosophila

Transposable elements (TEs) are mobile genetic elements, highly enriched in heterochromatin, that constitute a large percentage of the DNA content of eukaryotic genomes. Aging in Drosophila melanogaster is characterized by loss of repressive heterochromatin structure and loss of silencing of reporter genes in constitutive heterochromatin regions. Using next-generation sequencing, we found that transcripts of many genes native to heterochromatic regions and TEs increased with age in fly heads and fat bodies. A dietary restriction regimen, known to extend life span, repressed the age-related increased expression of genes located in heterochromatin, as well as TEs. We also observed a corresponding age-associated increase in TE transposition in fly fat body cells that was delayed by dietary restriction. Furthermore, we found that manipulating genes known to affect heterochromatin structure, including overexpression of Sir2, Su(var)3-9, and Dicer-2, as well as decreased expression of Adar, mitigated age-related increases in expression of TEs. Increasing expression of either Su(var)3-9 or Dicer-2 also led to an increase in life span. Mutation of Dicer-2 led to an increase in DNA double-strand breaks. Treatment with the reverse transcriptase inhibitor 3TC resulted in decreased TE transposition as well as increased life span in TE-sensitized Dicer-2 mutants. Together, these data support the retrotransposon theory of aging, which hypothesizes that epigenetically silenced TEs become deleteriously activated as cellular defense and surveillance mechanisms break down with age. Furthermore, interventions that maintain repressive heterochromatin and preserve TE silencing may prove key to preventing damage caused by TE activation and extending healthy life span.

Direct interrogation of the role of H3K9 in metazoan heterochromatin function

A defining feature of heterochromatin is methylation of Lys9 of histone H3 (H3K9me), a binding site for heterochromatin protein 1 (HP1). Penke et al. generated and analyzed H3K9R mutant flies, separating the functions of H3K9 and nonhistone substrates of H3K9 methyltransferases.

A defining feature of heterochromatin is methylation of Lys9 of histone H3 (H3K9me), a binding site for heterochromatin protein 1 (HP1). Although H3K9 methyltransferases and HP1 are necessary for proper heterochromatin structure, the specific contribution of H3K9 to heterochromatin function and animal development is unknown. Using our recently developed platform to engineer histone genes in Drosophila, we generated H3K9R mutant flies, separating the functions of H3K9 and nonhistone substrates of H3K9 methyltransferases. Nucleosome occupancy and HP1a binding at pericentromeric heterochromatin are markedly decreased in H3K9R mutants. Despite these changes in chromosome architecture, a small percentage of H3K9R mutants complete development. Consistent with this result, expression of most protein-coding genes, including those within heterochromatin, is similar between H3K9R and controls. In contrast, H3K9R mutants exhibit increased open chromatin and transcription from piRNA clusters and transposons, resulting in transposon mobilization. Hence, transposon silencing is a major developmental function of H3K9.

Comparative Genomic Analyses Provide New Insights into the Evolutionary Dynamics of Heterochromatin in Drosophila

The term heterochromatin has been long considered synonymous with gene silencing, but it is now clear that the presence of transcribed genes embedded in pericentromeric heterochromatin is a conserved feature in the evolution of eukaryotic genomes. Several studies have addressed the epigenetic changes that enable the expression of genes in pericentric heterochromatin, yet little is known about the evolutionary processes through which this has occurred. By combining genome annotation analysis and high-resolution cytology, we have identified and mapped 53 orthologs of D. melanogaster heterochromatic genes in the genomes of two evolutionarily distant species, D. pseudoobscura and D. virilis. Our results show that the orthologs of the D. melanogaster heterochromatic genes are clustered at three main genomic regions in D. virilis and D. pseudoobscura. In D. virilis, the clusters lie in the middle of euchromatin, while those in D. pseudoobscura are located in the proximal portion of the chromosome arms. Some orthologs map to the corresponding Muller C element in D. pseudoobscura and D. virilis, while others localize on the Muller B element, suggesting that chromosomal rearrangements that have been instrumental in the fusion of two separate elements involved the progenitors of genes currently located in D. melanogaster heterochromatin. These results demonstrate an evolutionary repositioning of gene clusters from ancestral locations in euchromatin to the pericentromeric heterochromatin of descendent D. melanogaster chromosomes. Remarkably, in both D. virilis and D. pseudoobscura the gene clusters show a conserved association with the HP1a protein, one of the most highly evolutionarily conserved epigenetic marks. In light of these results, we suggest a new scenario whereby ancestral HP1-like proteins (and possibly other epigenetic marks) may have contributed to the evolutionary repositioning of gene clusters into heterochromatin.

CHANGE POINT ANALYSIS OF HISTONE MODIFICATIONS REVEALS EPIGENETIC BLOCKS LINKING TO PHYSICAL DOMAINS

Histone modification is a vital epigenetic mechanism for transcriptional control in eukaryotes. High-throughput techniques have enabled whole-genome analysis of histone modifications in recent years. However, most studies assume one combination of histone modification invariantly translates to one transcriptional output regardless of local chromatin environment. In this study we hypothesize that, the genome is organized into local domains that manifest similar enrichment pattern of histone modification, which leads to orchestrated regulation of expression of genes with relevant biological functions. We propose a multivariate Bayesian Change Point (BCP) model to segment the Drosophila melanogaster genome into consecutive blocks on the basis of combinatorial patterns of histone marks. By modeling the sparse distribution of histone marks with a zero-inflated Gaussian mixture, our partitions capture local BLOCKs that manifest relatively homogeneous enrichment pattern of histone marks. We further characterized BLOCKs by their transcription levels, distribution of genes, degree of co-regulation and GO enrichment. Our results demonstrate that these BLOCKs, although inferred merely from histone modifications, reveal strong relevance with physical domains, which suggests their important roles in chromatin organization and coordinated gene regulation.

Conservation and divergence of the histone code in nucleomorphs

BACKGROUND

Nucleomorphs, the remnant nuclei of photosynthetic algae that have become endosymbionts to other eukaryotes, represent a unique example of convergent reductive genome evolution in eukaryotes, having evolved independently on two separate occasions in chlorarachniophytes and cryptophytes. The nucleomorphs of the two groups have evolved in a remarkably convergent manner, with numerous very similar features. Chief among them is the extreme reduction and compaction of nucleomorph genomes, with very small chromosomes and extremely short or even completely absent intergenic spaces. These characteristics pose a number of intriguing questions regarding the mechanisms of transcription and gene regulation in such a crowded genomic context, in particular in terms of the functioning of the histone code, which is common to almost all eukaryotes and plays a central role in chromatin biology.

RESULTS

This study examines the sequences of nucleomorph histone proteins in order to address these issues. Remarkably, all classical transcription- and repression-related components of the histone code seem to be missing from chlorarachniophyte nucleomorphs. Cryptophyte nucleomorph histones are generally more similar to the conventional eukaryotic state; however, they also display significant deviations from the typical histone code. Based on the analysis of specific components of the code, we discuss the state of chromatin and the transcriptional machinery in these nuclei.

CONCLUSIONS

The results presented here shed new light on the mechanisms of nucleomorph transcription and gene regulation and provide a foundation for future studies of nucleomorph chromatin and transcriptional biology.

Meiotic recombination is suppressed near the histone-defined border of euchromatin and heterochromatin on chromosome 2L of Drosophila melanogaster

In Drosophila melanogaster, the borders between pericentric heterochromatin and euchromatin on the major chromosome arms have been defined in various ways, including chromatin-specific histone modifications, the binding patterns of heterochromatin-enriched chromosomal proteins, and various cytogenetic techniques. Elucidation of the genetic properties that independently define the different chromatin states associated with heterochromatin and euchromatin should help refine the boundary. Since meiotic recombination is present in euchromatin, but absent in heterochromatin, it constitutes a key genetic property that can be observed transitioning between chromatin states. Using P element insertion lines marked with a su(Hw) insulated mini-white gene, meiotic recombination was found to transition in a region consistent with the H3K9me2 transition observed in ovaries.

Rapid evolution of a Y-chromosome heterochromatin protein underlies sex chromosome meiotic drive

Sex chromosome meiotic drive, the non-Mendelian transmission of sex chromosomes, is the expression of an intragenomic conflict that can have extreme evolutionary consequences. However, the molecular bases of such conflicts remain poorly understood. Here, we show that a young and rapidly evolving X-linked heterochromatin protein 1 (HP1) gene, HP1D2, plays a key role in the classical Paris sex-ratio (SR) meiotic drive occurring in Drosophila simulans Driver HP1D2 alleles prevent the segregation of the Y chromatids during meiosis II, causing female-biased sex ratio in progeny. HP1D2 accumulates on the heterochromatic Y chromosome in male germ cells, strongly suggesting that it controls the segregation of sister chromatids through heterochromatin modification. We show that Paris SR drive is a consequence of dysfunctional HP1D2 alleles that fail to prepare the Y chromosome for meiosis, thus providing evidence that the rapid evolution of genes controlling the heterochromatin structure can be a significant source of intragenomic conflicts.

Spatially Resolved Quantification of Chromatin Condensation through Differential Local Rheology in Cell Nuclei Fluorescence Lifetime Imaging

The linear sequence of DNA encodes access to the complete set of proteins that carry out cellular functions. Yet, much of the functionality appropriate for each cell is nested within layers of dynamic regulation and organization, including a hierarchy of chromatin structural states and spatial arrangement within the nucleus. There remain limitations in our understanding of gene expression within the context of nuclear organization from an inability to characterize hierarchical chromatin organization in situ. Here we demonstrate the use of fluorescence lifetime imaging microscopy (FLIM) to quantify and spatially resolve chromatin condensation state using cell-permeable, DNA-binding dyes (Hoechst 33342 and PicoGreen). Through in vitro and in situ experiments we demonstrate the sensitivity of fluorescence lifetime to condensation state through the mechanical effects that accompany the structural changes and are reflected through altered viscosity. The establishment of FLIM for resolving and quantifying chromatin condensation state opens the door for single-measurement mechanical studies of the nucleus and for characterizing the role of genome structure and organization in nuclear processes that accompany physiological and pathological changes.

Targeting of P-Element Reporters to Heterochromatic Domains by Transposable Element 1360 in Drosophila melanogaster

Heterochromatin is a common DNA packaging form employed by eukaryotes to constitutively silence transposable elements. Determining which sequences to package as heterochromatin is vital for an organism. Here, we use Drosophila melanogaster to study heterochromatin formation, exploiting position-effect variegation, a process whereby a transgene is silenced stochastically if inserted in proximity to heterochromatin, leading to a variegating phenotype. Previous studies identified the transposable element 1360 as a target for heterochromatin formation. We use transgene reporters with either one or four copies of 1360 to determine if increasing local repeat density can alter the fraction of the genome supporting heterochromatin formation. We find that including 1360 in the reporter increases the frequency with which variegating phenotypes are observed. This increase is due to a greater recovery of insertions at the telomere-associated sequences (∼50% of variegating inserts). In contrast to variegating insertions elsewhere, the phenotype of telomere-associated sequence insertions is largely independent of the presence of 1360 in the reporter. We find that variegating and fully expressed transgenes are located in different types of chromatin and that variegating reporters in the telomere-associated sequences differ from those in pericentric heterochromatin. Indeed, chromatin marks at the transgene insertion site can be used to predict the eye phenotype. Our analysis reveals that increasing the local repeat density (via the transgene reporter) does not enlarge the fraction of the genome supporting heterochromatin formation. Rather, additional copies of 1360 appear to target the reporter to the telomere-associated sequences with greater efficiency, thus leading to an increased recovery of variegating insertions.

Heterochromatin remodeling by CDK12 contributes to learning in Drosophila

Dynamic regulation of chromatin structure is required to modulate the transcription of genes in eukaryotes. However, the factors that contribute to the plasticity of heterochromatin structure are elusive. Here, we report that cyclin-dependent kinase 12 (CDK12), a transcription elongation-associated RNA polymerase II (RNAPII) kinase, antagonizes heterochromatin enrichment in Drosophila chromosomes. Notably, loss of CDK12 induces the ectopic accumulation of heterochromatin protein 1 (HP1) on euchromatic arms, with a prominent enrichment on the X chromosome. Furthermore, ChIP and sequencing analysis reveals that the heterochromatin enrichment on the X chromosome mainly occurs within long genes involved in neuronal functions. Consequently, heterochromatin enrichment reduces the transcription of neuronal genes in the adult brain and results in a defect in Drosophila courtship learning. Taken together, these results define a previously unidentified role of CDK12 in controlling the epigenetic transition between euchromatin and heterochromatin and suggest a chromatin regulatory mechanism in neuronal behaviors.

Heterochromatic breaks move to the nuclear periphery to continue recombinational repair

Heterochromatin mostly comprises repeated sequences prone to harmful ectopic recombination during double-strand break (DSB) repair. In Drosophila cells, ‘safe’ homologous recombination (HR) repair of heterochromatic breaks relies on a specialized pathway that relocalizes damaged sequences away from the heterochromatin domain before strand invasion. Here we show that heterochromatic DSBs move to the nuclear periphery to continue HR repair. Relocalization depends on nuclear pore and inner nuclear membrane proteins (INMPs) that anchor repair sites to the nuclear periphery via the Smc5/6-interacting proteins STUbL/RENi. Both the initial block to HR progression inside the heterochromatin domain, and the targeting of repair sites to the nuclear periphery, rely on SUMO and SUMO E3 ligases. This study reveals a critical role for SUMOylation in the spatial and temporal regulation of HR repair in heterochromatin, and identifies the nuclear periphery as a specialized site for heterochromatin repair in a multicellular eukaryote.

Modulation of Heterochromatin by Male Specific Lethal Proteins and roX RNA in Drosophila melanogaster Males

The ribonucleoprotein Male Specific Lethal (MSL) complex is required for X chromosome dosage compensation in Drosophila melanogaster males. Beginning at 3 h of development the MSL complex binds transcribed X-linked genes and modifies chromatin. A subset of MSL complex proteins, including MSL1 and MSL3, is also necessary for full expression of autosomal heterochromatic genes in males, but not females. Loss of the non-coding roX RNAs, essential components of the MSL complex, lowers the expression of heterochromatic genes and suppresses position effect variegation (PEV) only in males, revealing a sex-limited disruption of heterochromatin. To explore the molecular basis of this observation we examined additional proteins that participate in compensation and found that MLE, but not Jil-1 kinase, contributes to heterochromatic gene expression. To determine if identical regions of roX RNA are required for dosage compensation and heterochromatic silencing, we tested a panel of roX1 transgenes and deletions and find that the X chromosome and heterochromatin functions are separable by some mutations. Chromatin immunoprecipitation of staged embryos revealed widespread autosomal binding of MSL3 before and after localization of the MSL complex to the X chromosome at 3 h AEL. Autosomal MSL3 binding was dependent on MSL1, supporting the idea that a subset of MSL proteins associates with chromatin throughout the genome during early development. The broad localization of these proteins early in embryogenesis supports the idea of direct action at autosomal sites. We postulate that this may contribute to the sex-specific differences in heterochromatin that we, and others, have noted.

Identification of chromatin marks at TERRA promoter and encoding region

TERRA is a long non-coding RNA that is essential for telomere integrity. Although it is transcribed from subtelomeres and telomeres, how it is expressed in heterochromatic region is currently unknown. In this study, we focused our analysis on TERRA-encoding region TelBam3.4 and TelBam3.4-like sequences, and determined their transcription start sites, as well as enrichment of RNA polymerase II and histone modifications. We found that H3K4me3 and H3K9me3 are present at TERRA promoters, whereas H3K27ac and H3K9me3 are present at telomeric repeats. Consistently, we show that presence of active histone modifications H3K4me3 and H3K27ac are correlated to TERRA expression. These results mark an important step towards understanding telomere maintenance and transcription.

Tools for Targeted Genome Engineering of Established Drosophila Cell Lines

We describe an adaptation of φC31 integrase-mediated targeted cassette exchange for use in Drosophila cell lines. Single copies of an attP-bounded docking platform carrying a GFP-expression marker, with or without insulator elements flanking the attP sites, were inserted by P-element transformation into the Kc167 and Sg4 cell lines; each of the resulting docking-site lines carries a single mapped copy of one of the docking platforms. Vectors for targeted substitution contain a cloning cassette flanked by attB sites. Targeted substitution occurs by integrase-mediated substitution between the attP sites (integrated) and the attB sites (vector). We describe procedures for isolating cells carrying the substitutions and for eliminating the products of secondary off-target events. We demonstrate the technology by integrating a cassette containing a Cu(2+)-inducible mCherry marker, and we report the expression properties of those lines. When compared with clonal lines made by traditional transformation methods, which lead to the illegitimate insertion of tandem arrays, targeted insertion lines give more uniform expression, lower basal expression, and higher induction ratios. Targeted substitution, though intricate, affords results that should greatly improve comparative expression assays-a major emphasis of cell-based studies.

Genome architecture: from linear organisation of chromatin to the 3D assembly in the nucleus

The genetic information is stored in the eukaryotic nucleus in the form of chromatin. This is a macromolecular entity that includes genomic DNA and histone proteins that form nucleosomes, plus a large variety of chromatin-associated non-histone proteins. Chromatin is structurally and functionally organised at various levels. One reveals the linear topography of DNA, histones and their post-translational modifications and non-histone proteins along each chromosome. This level provides regulatory information about the association of genomic elements with particular signatures that have been used to define chromatin states. Importantly, these chromatin states correlate with structural and functional genomic features. Another regulatory layer is established at the level of the 3D organisation of chromatin within the nucleus, which has been revealed clearly as non-random. Instead, a variety of intra- and inter-chromosomal genomic domains with specific epigenetic and functional properties has been identified. In this review, we discuss how the recent advances in genomic approaches have contributed to our understanding of these two levels of genome architecture. We have emphasised our analysis with the aim of integrating information available for yeast, Arabidopsis, Drosophila, and mammalian cells. We consider that this comparative study helps define common and unique features in each system, providing a basis to better understand the complexity of genome organisation.

Lessons from modENCODE

The modENCODE (Model Organism Encyclopedia of DNA Elements) Consortium aimed to map functional elements-including transcripts, chromatin marks, regulatory factor binding sites, and origins of DNA replication-in the model organisms Drosophila melanogaster and Caenorhabditis elegans. During its five-year span, the consortium conducted more than 2,000 genome-wide assays in developmentally staged animals, dissected tissues, and homogeneous cell lines. Analysis of these data sets provided foundational insights into genome, epigenome, and transcriptome structure and the evolutionary turnover of regulatory pathways. These studies facilitated a comparative analysis with similar data types produced by the ENCODE Consortium for human cells. Genome organization differs drastically in these distant species, and yet quantitative relationships among chromatin state, transcription, and cotranscriptional RNA processing are deeply conserved. Of the many biological discoveries of the modENCODE Consortium, we highlight insights that emerged from integrative studies. We focus on operational and scientific lessons that may aid future projects of similar scale or aims in other, emerging model systems.

Ancestral Chromatin Configuration Constrains Chromatin Evolution on Differentiating Sex Chromosomes in Drosophila

Sex chromosomes evolve distinctive types of chromatin from a pair of ancestral autosomes that are usually euchromatic. In Drosophila, the dosage-compensated X becomes enriched for hyperactive chromatin in males (mediated by H4K16ac), while the Y chromosome acquires silencing heterochromatin (enriched for H3K9me2/3). Drosophila autosomes are typically mostly euchromatic but the small dot chromosome has evolved a heterochromatin-like milieu (enriched for H3K9me2/3) that permits the normal expression of dot-linked genes, but which is different from typical pericentric heterochromatin. In Drosophila busckii, the dot chromosomes have fused to the ancestral sex chromosomes, creating a pair of 'neo-sex' chromosomes. Here we collect genomic, transcriptomic and epigenomic data from D. busckii, to investigate the evolutionary trajectory of sex chromosomes from a largely heterochromatic ancestor. We show that the neo-sex chromosomes formed <1 million years ago, but nearly 60% of neo-Y linked genes have already become non-functional. Expression levels are generally lower for the neo-Y alleles relative to their neo-X homologs, and the silencing heterochromatin mark H3K9me2, but not H3K9me3, is significantly enriched on silenced neo-Y genes. Despite rampant neo-Y degeneration, we find that the neo-X is deficient for the canonical histone modification mark of dosage compensation (H4K16ac), relative to autosomes or the compensated ancestral X chromosome, possibly reflecting constraints imposed on evolving hyperactive chromatin in an originally heterochromatic environment. Yet, neo-X genes are transcriptionally more active in males, relative to females, suggesting the evolution of incipient dosage compensation on the neo-X. Our data show that Y degeneration proceeds quickly after sex chromosomes become established through genomic and epigenetic changes, and are consistent with the idea that the evolution of sex-linked chromatin is influenced by its ancestral configuration.

Chromatin structure analysis enables detection of DNA insertions into the mammalian nuclear genome

Background

Genetically modified organisms (GMOs) have numerous biomedical, agricultural and environmental applications. Development of accurate methods for the detection of GMOs is a prerequisite for the identification and control of authorized and unauthorized release of these engineered organisms into the environment and into the food chain. Current detection methods are unable to detect uncharacterized GMOs, since either the DNA sequence of the transgene or the amino acid sequence of the protein must be known for DNA-based or immunological-based detection, respectively.

Methods

Here we describe the application of an epigenetics-based approach for the detection of mammalian GMOs via analysis of chromatin structural changes occurring in the host nucleus upon the insertion of foreign or endogenous DNA.

Results

Immunological methods combined with DNA next generation sequencing enabled direct interrogation of chromatin structure and identification of insertions of various size foreign (human or viral) DNA sequences, DNA sequences often used as genome modification tools (e.g. viral sequences, transposon elements), or endogenous DNA sequences into the nuclear genome of a model animal organism.

Conclusions

The results provide a proof-of-concept that epigenetic approaches can be used to detect the insertion of endogenous and exogenous sequences into the genome of higher organisms where the method of genetic modification, the sequence of inserted DNA, and the exact genomic insertion site(s) are unknown.

General significance

Measurement of chromatin dynamics as a sensor for detection of genomic manipulation and, more broadly, organism exposure to environmental or other factors affecting the epigenomic landscape are discussed.

•

Insertion of DNA sequences into a host genome causes chromatin structure remodeling.

•

ChIP-seq identifies molecular signatures of DNA insertion into the mammalian genome.

•

Focus on epigenetic marks limits sequencing data amount required for GMO detection.

•

Proof-of-concept for use of chromatin dynamics as a sensor of genomic manipulation.

The Role of piRNA-Mediated Epigenetic Silencing in the Population Dynamics of Transposable Elements in Drosophila melanogaster

The piwi-interacting RNAs (piRNA) are small RNAs that target selfish transposable elements (TEs) in many animal genomes. Until now, piRNAs’ role in TE population dynamics has only been discussed in the context of their suppression of TE transposition, which alone is not sufficient to account for the skewed frequency spectrum and stable containment of TEs. On the other hand, euchromatic TEs can be epigenetically silenced via piRNA-dependent heterochromatin formation and, similar to the widely known “Position-effect variegation”, heterochromatin induced by TEs can “spread” into nearby genes. We hypothesized that the piRNA-mediated spread of heterochromatin from TEs into adjacent genes has deleterious functional effects and leads to selection against individual TEs. Unlike previously identified deleterious effects of TEs due to the physical disruption of DNA, the functional effect we investigated here is mediated through the epigenetic influences of TEs. We found that the repressive chromatin mark, H3K9me, is elevated in sequences adjacent to euchromatic TEs at multiple developmental stages in Drosophila melanogaster. Furthermore, the heterochromatic states of genes depend not only on the number of and distance from adjacent TEs, but also on the likelihood that their nearest TEs are targeted by piRNAs. These variations in chromatin status probably have functional consequences, causing genes near TEs to have lower expression. Importantly, we found stronger selection against TEs that lead to higher H3K9me enrichment of adjacent genes, demonstrating the pervasive evolutionary consequences of TE-induced epigenetic silencing. Because of the intrinsic biological mechanism of piRNA amplification, spread of TE heterochromatin could result in the theoretically required synergistic deleterious effects of TE insertions for stable containment of TE copy number. The indirect deleterious impact of piRNA-mediated epigenetic silencing of TEs is a previously unexplored, yet important, element for the evolutionary dynamics of TEs.

The piwi-interacting RNAs (piRNAs) are small RNAs that can suppress the expression of selfish transposable elements (TEs) in many animal genomes. One mechanism by which piRNAs silence TEs is through the formation of heterochromatin, which is condensed chromatin and generally associated with repressed gene expression. Several functional studies have demonstrated that piRNA-mediated heterochromatin of TEs can spread to adjacent genes. We hypothesized that this spread of TE-induced heterochromatin influences the function of adjacent genes, ultimately resulting in selection against individual TEs. Consistent with our hypothesis, we found that sequences and genes adjacent to TEs are enriched in heterochromatic marks. We determine that this TE-induced variation in epigenetic status is probably piRNA-dependent and that this change in chromatin state influences the expression levels of adjacent genes. Importantly, TEs that lead to higher heterochromatin enrichment of adjacent genes are more strongly selected against, demonstrating the evolutionary consequences of TE-induced epigenetic silencing. In contrast to previously studied deleterious impacts of TEs, which depend on TEs’ physical disruptions of DNAs, our proposed functional effect of TEs is mediated through their epigenetic influence. Our study suggests that the piRNA-dependent epigenetic impact of TEs may play an important role in the evolutionary dynamics of TEs.

A balancing act: heterochromatin protein 1a and the Polycomb group coordinate their levels to silence chromatin in Drosophila

BACKGROUND

The small non-histone protein Heterochromatin protein 1a (HP1a) plays a vital role in packaging chromatin, most notably in forming constitutive heterochromatin at the centromeres and telomeres. A second major chromatin regulating system is that of the Polycomb/trithorax groups of genes which, respectively, maintain the repressed/activated state of euchromatin. Recent analyses suggest they affect the expression of a multitude of genes, beyond the homeotics whose alteration in expression lead to their initial discovery.

RESULTS

Our data suggest that early in Drosophila development, HP1a collaborates with the Polycomb/trithorax groups of proteins to regulate gene expression and that the two chromatin systems do not act separately as convention describes. HP1a affects the levels of both the Polycomb complexes and RNA polymerase II at promoters, as assayed by chromatin immunoprecipitation analysis. Deposition of both the repressive (H3K27me3) and activating (H3K4me3) marks promoted by the Polycomb/trithorax group genes at gene promoters is affected. Additionally, depending on which parent contributes the null mutation of the HP1a gene, the levels of the H3K27me3 and H3K9me3 silencing marks at both promoters and heterochromatin are different. Changes in levels of the H3K27me3 and H3K9me3 repressive marks show a mostly reciprocal nature. The time around the mid-blastula transition, when the zygotic genome begins to be actively transcribed, appears to be a transition/decision point for setting the levels.

CONCLUSIONS

We find that HP1a, which is normally critical for the formation of constitutive heterochromatin, also affects the generation of the epigenetic marks of the Polycomb/trithorax groups of proteins, chromatin modifiers which are key to maintaining gene expression in euchromatin. At gene promoters, deposition of both the repressive H3K27me3 and activating H3K4me3 marks of histone modifications shows a dependence on HP1a. Around the mid-blastula transition, when the zygotic genome begins to be actively transcribed, a pivotal decision for the level of silencing appears to take place. This is also when the embryo organizes its genome into heterochromatin and euchromatin. A balance between the HP1a and Polycomb group silencing systems appears to be set for the chromatin types that each system will primarily regulate.

Drosophila muller f elements maintain a distinct set of genomic properties over 40 million years of evolution

The Muller F element (4.2 Mb, ~80 protein-coding genes) is an unusual autosome of Drosophila melanogaster; it is mostly heterochromatic with a low recombination rate. To investigate how these properties impact the evolution of repeats and genes, we manually improved the sequence and annotated the genes on the D. erecta, D. mojavensis, and D. grimshawi F elements and euchromatic domains from the Muller D element. We find that F elements have greater transposon density (25–50%) than euchromatic reference regions (3–11%). Among the F elements, D. grimshawi has the lowest transposon density (particularly DINE-1: 2% vs. 11–27%). F element genes have larger coding spans, more coding exons, larger introns, and lower codon bias. Comparison of the Effective Number of Codons with the Codon Adaptation Index shows that, in contrast to the other species, codon bias in D. grimshawi F element genes can be attributed primarily to selection instead of mutational biases, suggesting that density and types of transposons affect the degree of local heterochromatin formation. F element genes have lower estimated DNA melting temperatures than D element genes, potentially facilitating transcription through heterochromatin. Most F element genes (~90%) have remained on that element, but the F element has smaller syntenic blocks than genome averages (3.4–3.6 vs. 8.4–8.8 genes per block), indicating greater rates of inversion despite lower rates of recombination. Overall, the F element has maintained characteristics that are distinct from other autosomes in the Drosophila lineage, illuminating the constraints imposed by a heterochromatic milieu.

DNA replication and transcription programs respond to the same chromatin cues

DNA replication is a dynamic process that occurs in a temporal order along each of the chromosomes. A consequence of the temporally coordinated activation of replication origins is the establishment of broad domains (>100 kb) that replicate either early or late in S phase. This partitioning of the genome into early and late replication domains is important for maintaining genome stability, gene dosage, and epigenetic inheritance; however, the molecular mechanisms that define and establish these domains are poorly understood. The modENCODE Project provided an opportunity to investigate the chromatin features that define the Drosophila replication timing program in multiple cell lines. The majority of early and late replicating domains in the Drosophila genome were static across all cell lines; however, a small subset of domains was dynamic and exhibited differences in replication timing between the cell lines. Both origin selection and activation contribute to defining the DNA replication program. Our results suggest that static early and late replicating domains were defined at the level of origin selection (ORC binding) and likely mediated by chromatin accessibility. In contrast, dynamic domains exhibited low ORC densities in both cell types, suggesting that origin activation and not origin selection governs the plasticity of the DNA replication program. Finally, we show that the male-specific early replication of the X chromosome is dependent on the dosage compensation complex (DCC), suggesting that the transcription and replication programs respond to the same chromatin cues. Specifically, MOF-mediated hyperacetylation of H4K16 on the X chromosome promotes both the up-regulation of male-specific transcription and origin activation.

The Release 6 reference sequence of the Drosophila melanogaster genome

Drosophila melanogaster plays an important role in molecular, genetic, and genomic studies of heredity, development, metabolism, behavior, and human disease. The initial reference genome sequence reported more than a decade ago had a profound impact on progress in Drosophila research, and improving the accuracy and completeness of this sequence continues to be important to further progress. We previously described improvement of the 117-Mb sequence in the euchromatic portion of the genome and 21 Mb in the heterochromatic portion, using a whole-genome shotgun assembly, BAC physical mapping, and clone-based finishing. Here, we report an improved reference sequence of the single-copy and middle-repetitive regions of the genome, produced using cytogenetic mapping to mitotic and polytene chromosomes, clone-based finishing and BAC fingerprint verification, ordering of scaffolds by alignment to cDNA sequences, incorporation of other map and sequence data, and validation by whole-genome optical restriction mapping. These data substantially improve the accuracy and completeness of the reference sequence and the order and orientation of sequence scaffolds into chromosome arm assemblies. Representation of the Y chromosome and other heterochromatic regions is particularly improved. The new 143.9-Mb reference sequence, designated Release 6, effectively exhausts clone-based technologies for mapping and sequencing. Highly repeat-rich regions, including large satellite blocks and functional elements such as the ribosomal RNA genes and the centromeres, are largely inaccessible to current sequencing and assembly methods and remain poorly represented. Further significant improvements will require sequencing technologies that do not depend on molecular cloning and that produce very long reads.

Non-coding roX RNAs prevent the binding of the MSL-complex to heterochromatic regions

Long non-coding RNAs contribute to dosage compensation in both mammals and Drosophila by inducing changes in the chromatin structure of the X-chromosome. In Drosophila melanogaster, roX1 and roX2 are long non-coding RNAs that together with proteins form the male-specific lethal (MSL) complex, which coats the entire male X-chromosome and mediates dosage compensation by increasing its transcriptional output. Studies on polytene chromosomes have demonstrated that when both roX1 and roX2 are absent, the MSL-complex becomes less abundant on the male X-chromosome and is relocated to the chromocenter and the 4^th chromosome. Here we address the role of roX RNAs in MSL-complex targeting and the evolution of dosage compensation in Drosophila. We performed ChIP-seq experiments which showed that MSL-complex recruitment to high affinity sites (HAS) on the X-chromosome is independent of roX and that the HAS sequence motif is conserved in D. simulans. Additionally, a complete and enzymatically active MSL-complex is recruited to six specific genes on the 4^th chromosome. Interestingly, our sequence analysis showed that in the absence of roX RNAs, the MSL-complex has an affinity for regions enriched in Hoppel transposable elements and repeats in general. We hypothesize that roX mutants reveal the ancient targeting of the MSL-complex and propose that the role of roX RNAs is to prevent the binding of the MSL-complex to heterochromatin.

In both fruit flies and humans, males and females have different sets of sex chromosomes. This generates differences in gene dosage that must be compensated for by adjusting the transcriptional output of most genes located on the X-chromosome. The specific recognition and targeting of the X-chromosome is essential for such dosage compensation. In fruit flies, dosage compensation is mediated by the male-specific lethal (MSL) complex, which upregulates gene transcription on the male X-chromosome. The MSL-complex consists of five proteins and two non-coding RNAs, roX1 and roX2. While non-coding RNAs are known to be critical for dosage compensation in both flies and mammals, their precise functions remain elusive. Here we present a study on the targeting and function of the MSL-complex in the absence of roX RNAs. The results obtained suggest that the dosage compensating MSL-complex has an intrinsic tendency to target repeat-rich regions and that the function of roX RNAs is to prevent its binding to such targets. Our findings reveal an ancient targeting and regulatory function of the MSL-complex that has been adapted for use in dosage compensation and modified by the rapidly evolving noncoding roX RNAs.

cisMEP: an integrated repository of genomic epigenetic profiles and cis-regulatory modules in Drosophila

BACKGROUND

Cis-regulatory modules (CRMs), or the DNA sequences required for regulating gene expression, play the central role in biological researches on transcriptional regulation in metazoan species. Nowadays, the systematic understanding of CRMs still mainly resorts to computational methods due to the time-consuming and small-scale nature of experimental methods. But the accuracy and reliability of different CRM prediction tools are still unclear. Without comparative cross-analysis of the results and combinatorial consideration with extra experimental information, there is no easy way to assess the confidence of the predicted CRMs. This limits the genome-wide understanding of CRMs.

DESCRIPTION

It is known that transcription factor binding and epigenetic profiles tend to determine functions of CRMs in gene transcriptional regulation. Thus integration of the genome-wide epigenetic profiles with systematically predicted CRMs can greatly help researchers evaluate and decipher the prediction confidence and possible transcriptional regulatory functions of these potential CRMs. However, these data are still fragmentary in the literatures. Here we performed the computational genome-wide screening for potential CRMs using different prediction tools and constructed the pioneer database, cisMEP (cis-regulatory module epigenetic profile database), to integrate these computationally identified CRMs with genomic epigenetic profile data. cisMEP collects the literature-curated TFBS location data and nine genres of epigenetic data for assessing the confidence of these potential CRMs and deciphering the possible CRM functionality.

CONCLUSIONS

cisMEP aims to provide a user-friendly interface for researchers to assess the confidence of different potential CRMs and to understand the functions of CRMs through experimentally-identified epigenetic profiles. The deposited potential CRMs and experimental epigenetic profiles for confidence assessment provide experimentally testable hypotheses for the molecular mechanisms of metazoan gene regulation. We believe that the information deposited in cisMEP will greatly facilitate the comparative usage of different CRM prediction tools and will help biologists to study the modular regulatory mechanisms between different TFs and their target genes.

Comparative analysis of metazoan chromatin organization

Genome function is dynamically regulated in part by chromatin, which consists of the histones, non-histone proteins and RNA molecules that package DNA. Studies in Caenorhabditis elegans and Drosophila melanogaster have contributed substantially to our understanding of molecular mechanisms of genome function in humans, and have revealed conservation of chromatin components and mechanisms. Nevertheless, the three organisms have markedly different genome sizes, chromosome architecture and gene organization. On human and fly chromosomes, for example, pericentric heterochromatin flanks single centromeres, whereas worm chromosomes have dispersed heterochromatin-like regions enriched in the distal chromosomal 'arms', and centromeres distributed along their lengths. To systematically investigate chromatin organization and associated gene regulation across species, we generated and analysed a large collection of genome-wide chromatin data sets from cell lines and developmental stages in worm, fly and human. Here we present over 800 new data sets from our ENCODE and modENCODE consortia, bringing the total to over 1,400. Comparison of combinatorial patterns of histone modifications, nuclear lamina-associated domains, organization of large-scale topological domains, chromatin environment at promoters and enhancers, nucleosome positioning, and DNA replication patterns reveals many conserved features of chromatin organization among the three organisms. We also find notable differences in the composition and locations of repressive chromatin. These data sets and analyses provide a rich resource for comparative and species-specific investigations of chromatin composition, organization and function.

pENCODE: a plant encyclopedia of DNA elements

ENCODE projects exist for many eukaryotes, including humans, but as of yet no defined project exists for plants. A plant ENCODE would be invaluable to the research community and could be more readily produced than its metazoan equivalents by capitalizing on the preexisting infrastructure provided from similar projects. Collecting and normalizing plant epigenomic data for a range of species will facilitate hypothesis generation, cross-species comparisons, annotation of genomes, and an understanding of epigenomic functions throughout plant evolution. Here, we discuss the need for such a project, outline the challenges it faces, and suggest ways forward to build a plant ENCODE.

Insights into the epigenomic landscape of the human malaria vector Anopheles gambiae

The epigenome of the human malaria vector Anopheles gambiae was characterized in midgut cells by mapping the distribution and levels of two post-translational histone modifications, H3K27ac and H3K27me3. These histone profiles were then correlated with levels of gene expression obtained by RNA-seq. Analysis of the transcriptome of A. gambiae midguts and salivary glands led to the discovery of 13,898 new transcripts not present in the most recent genome assembly. A subset of these transcripts is differentially expressed between midgut and salivary glands. The enrichment profiles of H3K27ac and H3K27me3 are mutually exclusive and associate with high and low levels of transcription, respectively. This distribution agrees with previous findings in Drosophila showing association of these two histone modifications with either active or inactive transcriptional states, including Polycomb-associated domains in silenced genes. This study provides a mosquito epigenomics platform for future comparative studies in other mosquito species, opening future investigations into the role of epigenetic processes in vector-borne systems of medical and economic importance.

Heterochromatin and sibling species of Simulium praelargum s.l. (Diptera: Simuliidae)

Cytological descriptions are given for two sympatric sibling species, Simulium praelargum IIIL-st and Simulium praelargum IIIL-1.2, from Darjeeling, West Bengal, India. Sibling IIIL-1.2 differs from IIIL-st by a two step fixed included inversion (IIIL-1.2) found on the long arm of chromosome III. Both siblings possess heterochromatinized IIIS polytene chromosome ends, while IIIL-1.2 is unique within the Simuliidae in that it exhibits a very large enhanced chromocentre that persists in normal somatic tissue, contrary to other chromocentre-containing taxa. In IIIL-1.2, the chromocentre occurs as a positively allocyclic heteropycnotic body in normal mitotic interphase cells of neuroblast ganglia, oogonia, and spermatogonia. In mitotic metaphase chromosomes, the chromocentre forms large pronounced primary (centromere) constrictions and appears to be associated with nonhomologous pairing. Initial cytological studies on taxa within the feuerborni group have revealed heterochromatic chromosomal polymorphisms. Progressive fixation of these polymorphisms in different taxa within the group opens up the possibility of using comparative molecular or genomic approaches to begin to define the functional and structural aspects of the epigenome and to further characterize mitosis and meiosis in S. praelargum IIIL-1.2.

Nucleosomes shape DNA polymorphism and divergence

An estimated 80% of genomic DNA in eukaryotes is packaged as nucleosomes, which, together with the remaining interstitial linker regions, generate higher order chromatin structures [1]. Nucleosome sequences isolated from diverse organisms exhibit ∼10 bp periodic variations in AA, TT and GC dinucleotide frequencies. These sequence elements generate intrinsically curved DNA and help establish the histone-DNA interface. We investigated an important unanswered question concerning the interplay between chromatin organization and genome evolution: do the DNA sequence preferences inherent to the highly conserved histone core exert detectable natural selection on genomic divergence and polymorphism? To address this hypothesis, we isolated nucleosomal DNA sequences from Drosophila melanogaster embryos and examined the underlying genomic variation within and between species. We found that divergence along the D. melanogaster lineage is periodic across nucleosome regions with base changes following preferred nucleotides, providing new evidence for systematic evolutionary forces in the generation and maintenance of nucleosome-associated dinucleotide periodicities. Further, Single Nucleotide Polymorphism (SNP) frequency spectra show striking periodicities across nucleosomal regions, paralleling divergence patterns. Preferred alleles occur at higher frequencies in natural populations, consistent with a central role for natural selection. These patterns are stronger for nucleosomes in introns than in intergenic regions, suggesting selection is stronger in transcribed regions where nucleosomes undergo more displacement, remodeling and functional modification. In addition, we observe a large-scale (∼180 bp) periodic enrichment of AA/TT dinucleotides associated with nucleosome occupancy, while GC dinucleotide frequency peaks in linker regions. Divergence and polymorphism data also support a role for natural selection in the generation and maintenance of these super-nucleosomal patterns. Our results demonstrate that nucleosome-associated sequence periodicities are under selective pressure, implying that structural interactions between nucleosomes and DNA sequence shape sequence evolution, particularly in introns.

In eukaryotic cells, the majority of DNA is packaged in nucleosomes comprised of ∼147 bp of DNA wound tightly around the highly conserved histone octamer. Nucleosomal DNA from diverse organisms shows an anti-correlated ∼10 bp periodicity of AT-rich and GC-rich dinucleotides. These sequence features influence DNA bending and shape, facilitating structural interactions. We asked whether natural selection mediated through the periodic sequence preferences of nucleosomes shapes the evolution of non-protein-coding regions of D. melanogaster by examining the inter- and intra-species genomic variation relative to these fundamental chromatin building blocks. The sequence changes across nucleosome-bound regions on the melanogaster lineage mirror the observed nucleosome dinucleotide periodicities. Importantly, we show that the frequencies of polymorphisms in natural populations vary across these regions, paralleling divergence, with higher frequencies of preferred alleles. These patterns are most evident for intronic regions and indicate that non-protein coding regions are evolving toward sequences that facilitate the canonical association with the histone core. This result is consistent with the hypothesis that interactions between DNA and the core have systematic impacts on function that are subject to natural selection and are not solely due to mutational bias. These ubiquitous interactions with the histone core partially account for the evolutionary constraint observed in unannotated genomic regions, and may drive broad changes in base composition.

The chromatin landscape of the moss Physcomitrella patens and its dynamics during development and drought stress

The moss Physcomitrella patens is an important model organism for evo-devo studies. Here, we determined the genome-wide chromatin landscape of five important histone three (H3) modifications (H3K4me3, H3K27me3, H3K27Ac, H3K9Ac and H3K9me2) and describe the changes to these histone marks in two contrasted situations, developmental transition and abiotic (drought) stress. Integrative analysis of these histone H3 modifications revealed their preferential association into 15 chromatin states (CS) in genic regions of the P. patens genome. Synergistic relationships that influence expression levels were revealed for the three activating marks H3K4me3, H3K27Ac and H3K9Ac, while an antagonistic relationship was found between CS containing the H3K27me3 and H3K27Ac marks, suggesting that H3K27 is a key indexing residue regarding transcriptional output. Concerning the alteration of histone marks in response to developmental transition (juvenile to adult) and drought stress, the three activating marks H3K4me3, H3K27Ac and H3K9Ac show significant changes in both situations. However, changes to H3K27me3 are central only for genes differentially expressed during development. Interestingly, genes induced during drought stress show significant histone mark toggling during developmental transition. This situation suggests that drought induced adult (gametophore expressed) genes are primed to respond to this stress during the juvenile to adult transition.

Cell lines

We review the properties and uses of cell lines in Drosophila research, emphasizing the variety of lines, the large body of genomic and transcriptional data available for many of the lines, and the variety of ways the lines have been used to provide tools for and insights into the developmental, molecular, and cell biology of Drosophila and mammals.

The chromatin landscape of Drosophila: comparisons between species, sexes, and chromosomes

The chromatin landscape is key for gene regulation, but little is known about how it differs between sexes or between species. Here, we study the sex-specific chromatin landscape of Drosophila miranda, a species with young sex chromosomes, and compare it with Drosophila melanogaster. We analyze six histone modifications in male and female larvae of D. miranda (H3K4me1, H3K4me3, H3K36me3, H4K16ac, H3K27me3, and H3K9me2), and define seven biologically meaningful chromatin states that show different enrichments for transcribed and silent genes, repetitive elements, housekeeping, and tissue-specific genes. The genome-wide distribution of both active and repressive chromatin states differs between males and females. In males, active chromatin is enriched on the X, relative to females, due to dosage compensation of the hemizygous X. Furthermore, a smaller fraction of the euchromatic portion of the genome is in a repressive chromatin state in males relative to females. However, sex-specific chromatin states appear not to explain sex-biased expression of genes. Overall, conservation of chromatin states between male and female D. miranda is comparable to conservation between D. miranda and D. melanogaster, which diverged >30 MY ago. Active chromatin states are more highly conserved across species, while heterochromatin shows very low levels of conservation. Divergence in chromatin profiles contributes to expression divergence between species, with ∼26% of genes in different chromatin states in the two species showing species-specific or species-biased expression, an enrichment of approximately threefold over null expectation. Our data suggest that heteromorphic sex chromosomes in males (that is, a hypertranscribed X and an inactivated Y) may contribute to global redistribution of active and repressive chromatin marks between chromosomes and sexes.

The Little Fly that Could: Wizardry and Artistry of Drosophila Genomics

For more than 100 years now, the fruit fly Drosophila melanogaster has been at the forefront of our endeavors to unlock the secrets of the genome. From the pioneering studies of chromosomes and heredity by Morgan and his colleagues, to the generation of fly models for human disease, Drosophila research has been at the forefront of genetics and genomics. We present a broad overview of some of the most powerful genomics tools that keep Drosophila research at the cutting edge of modern biomedical research.

Impact of nuclear Piwi elimination on chromatin state in Drosophila melanogaster ovaries

The Piwi-interacting RNA (piRNA)-interacting Piwi protein is involved in transcriptional silencing of transposable elements in ovaries of Drosophila melanogaster. Here we characterized the genome-wide effect of nuclear Piwi elimination on the presence of the heterochromatic H3K9me3 mark and HP1a, as well as on the transcription-associated mark H3K4me2. Our results demonstrate that a significant increase in the H3K4me2 level upon nuclear Piwi loss is not accompanied by the alterations in H3K9me3 and HP1a levels for several germline-expressed transposons, suggesting that in this case Piwi prevents transcription by a mechanism distinct from H3K9 methylation. We found that the targets of Piwi-dependent chromatin repression are mainly related to the elements that display a higher level of H3K4me2 modification in the absence of silencing, i.e. most actively transcribed elements. We also show that Piwi-guided silencing does not significantly influence the chromatin state of dual-strand piRNA-producing clusters. In addition, host protein-coding gene expression is essentially not affected due to the nuclear Piwi elimination, but we noted an increase in small nuclear spliceosomal RNAs abundance and propose Piwi involvement in their post-transcriptional regulation. Our work reveals new aspects of transposon silencing in Drosophila, indicating that transcription of transposons can underpin their Piwi dependent silencing, while canonical heterochromatin marks are not obligatory for their repression.