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

Chromatin Landscape Is Associated With Sex-Biased Expression and Drosophila-Like Dosage Compensation of the Z Chromosome in Artemia franciscana

The males and females of the brine shrimp Artemia franciscana are highly dimorphic, and this dimorphism is associated with substantial sex-biased gene expression in heads and gonads. How these sex-specific patterns of expression are regulated at the molecular level is unknown. A. franciscana also has differentiated ZW sex chromosomes, with complete dosage compensation, but the molecular mechanism through which compensation is achieved is unknown. Here, we conducted CUT&TAG assays targeting 7 post-translational histone modifications (H3K27me3, H3K9me2, H3K9me3, H3K36me3, H3K27ac, H3K4me3, and H4K16ac) in heads and gonads of A. franciscana, allowing us to divide the genome into 12 chromatin states. We further defined functional chromatin signatures for all genes, which were correlated with transcript level abundances. Differences in the occupancy of the profiled epigenetic marks between sexes were associated with differential gene expression between males and females. Finally, we found a significant enrichment of the permissive H4K16ac histone mark in the Z-specific region in both tissues of females but not males, supporting the role of this histone mark in mediating dosage compensation of the Z chromosome.

Polymorphic transposable elements contribute to variation in recombination landscapes

Meiotic recombination is a prominent force shaping genome evolution, and understanding why recombination rates vary within and between species has remained a central, though challenging, question. Variation in recombination is widely thought to influence the efficacy of selection in purging transposable elements (TEs), prevalent selfish genetic elements, leading to widely observed negative correlations between TE abundance and recombination rates across taxa. However, accumulating evidence suggests that TEs could instead be the cause rather than the consequence of this relationship. To test this prediction, we formally investigated the influence of polymorphic, putatively active TEs on recombination rates. We developed and benchmarked an approach that uses PacBio long-read sequencing to efficiently, accurately, and cost-effectively identify crossovers (COs), a key recombination product, among large numbers of pooled recombinant individuals. By applying this approach to Drosophila strains with distinct TE insertion profiles, we found that polymorphic TEs, especially RNA-based TEs and TEs with local enrichment of repressive marks, reduce the occurrence of COs. Such an effect leads to different CO frequencies between homologous sequences with and without TEs, contributing to varying CO maps between individuals. The suppressive effect of TEs on CO is further supported by two orthogonal approaches-analyzing the distributions of COs in panels of recombinant inbred lines in relation to TE polymorphism and applying marker-assisted estimations of CO frequencies to isogenic strains with and without transgenically inserted TEs. Our investigations reveal how the constantly changing TE landscape can actively modify recombination, shaping genome evolution within and between species.

HP1a promotes chromatin liquidity and drives spontaneous heterochromatin compartmentalization

Compartmentalization of the nucleus into heterochromatin and euchromatin is highly conserved across eukaryotes. Constitutive heterochromatin (C-Het) constitutes a liquid-like condensate that packages the repetitive regions of the genome through the enrichment of histone modification H3K9me3 and recruitment of its cognate reader protein Heterochromatin Protein-1 (HP1a). The ability for well-ordered nucleosome arrays and HP1a to independently form biomolecular condensates suggests that the emergent material properties of C-Het compartments may contribute to its functions such as force-buffering, dosage-dependent gene silencing, and selective permeability. Using an in vitro reconstitution system we directly assess the contributions of H3K9me3 and HP1a on the biophysical properties of C-Het. In the presence of HP1a, H3K9me3 (Me-) and unmodified (U-) chromatin form co-condensates composed of distinct, immiscible domains. These chromatin domains form spontaneously and are reversible. Independently of HP1a, H3K9me3 modifications are sufficient to increase linker-DNA length within chromatin arrays and slow chromatin condensate growth. HP1a increases the liquidity of chromatin condensates while dramatically differentiating the viscoelastic properties of Me-chromatin versus U-chromatin. Mutating key residues in HP1a show that HP1a interactions with itself and chromatin determine the relative interfacial tension between chromatin compartments, however the formation of condensates is driven by the underlying chromatin. These direct measurements map the energetic landscape that determines C-Het compartmentalization, demonstrating that nuclear compartmentalization is a spontaneous and energetically favorable process in which HP1a plays a critical role in establishing a hierarchy of affinities between H3K9me3-chromatin and unmodified-chromatin.

STAR-RL: Spatial-Temporal Hierarchical Reinforcement Learning for Interpretable Pathology Image Super-Resolution

Pathology image are essential for accurately interpreting lesion cells in cytopathology screening, but acquiring high-resolution digital slides requires specialized equipment and long scanning times. Though super-resolution (SR) techniques can alleviate this problem, existing deep learning models recover pathology image in a black-box manner, which can lead to untruthful biological details and misdiagnosis. Additionally, current methods allocate the same computational resources to recover each pixel of pathology image, leading to the sub-optimal recovery issue due to the large variation of pathology image. In this paper, we propose the first hierarchical reinforcement learning framework named Spatial-Temporal hierARchical Reinforcement Learning (STAR-RL), mainly for addressing the aforementioned issues in pathology image super-resolution problem. We reformulate the SR problem as a Markov decision process of interpretable operations and adopt the hierarchical recovery mechanism in patch level, to avoid sub-optimal recovery. Specifically, the higher-level spatial manager is proposed to pick out the most corrupted patch for the lower-level patch worker. Moreover, the higher-level temporal manager is advanced to evaluate the selected patch and determine whether the optimization should be stopped earlier, thereby avoiding the over-processed problem. Under the guidance of spatial-temporal managers, the lower-level patch worker processes the selected patch with pixel-wise interpretable actions at each time step. Experimental results on medical images degraded by different kernels show the effectiveness of STAR-RL. Furthermore, STAR-RL validates the promotion in tumor diagnosis with a large margin and shows generalizability under various degradations. The source code is available at https://github.com/CUHK-AIM-Group/STAR-RL.

DNA double-strand break movement in heterochromatin depends on the histone acetyltransferase dGcn5

Cells employ diverse strategies to repair double-strand breaks (DSBs), a dangerous form of DNA damage that threatens genome integrity. Eukaryotic nuclei consist of different chromatin environments, each displaying distinct molecular and biophysical properties that can significantly influence the DSB-repair process. DSBs arising in the compact and silenced heterochromatin domains have been found to move to the heterochromatin periphery in mouse and Drosophila to prevent aberrant recombination events. However, it is poorly understood how chromatin components, such as histone post-translational modifications, contribute to these DSB movements within heterochromatin. Using irradiation as well as locus-specific DSB induction in Drosophila tissues and cultured cells, we find enrichment of histone H3 lysine 9 acetylation (H3K9ac) at DSBs in heterochromatin but not euchromatin. We find this increase is mediated by the histone acetyltransferase dGcn5, which rapidly localizes to heterochromatic DSBs. Moreover, we demonstrate that in the absence of dGcn5, heterochromatic DSBs display impaired recruitment of the SUMO E3 ligase Nse2/Qjt and fail to relocate to the heterochromatin periphery to complete repair. In summary, our results reveal a previously unidentified role for dGcn5 and H3K9ac in heterochromatic DSB repair and underscore the importance of differential chromatin responses at heterochromatic and euchromatic DSBs to promote safe repair.

Epigenetic inheritance and gene expression regulation in early Drosophila embryos

Precise spatiotemporal regulation of gene expression is of paramount importance for eukaryotic development. The maternal-to-zygotic transition (MZT) during early embryogenesis in Drosophila involves the gradual replacement of maternally contributed mRNAs and proteins by zygotic gene products. The zygotic genome is transcriptionally activated during the first 3 hours of development, in a process known as "zygotic genome activation" (ZGA), by the orchestrated activities of a few pioneer factors. Their decisive role during ZGA has been characterized in detail, whereas the contribution of chromatin factors to this process has been historically overlooked. In this review, we aim to summarize the current knowledge of how chromatin regulation impacts the first stages of Drosophila embryonic development. In particular, we will address the following questions: how chromatin factors affect ZGA and transcriptional silencing, and how genome architecture promotes the integration of these processes early during development. Remarkably, certain chromatin marks can be intergenerationally inherited, and their presence in the early embryo becomes critical for the regulation of gene expression at later stages. Finally, we speculate on the possible roles of these chromatin marks as carriers of epialleles during transgenerational epigenetic inheritance (TEI).

Varying recombination landscapes between individuals are driven by polymorphic transposable elements

Meiotic recombination is a prominent force shaping genome evolution, and understanding the causes for varying recombination landscapes within and between species has remained a central, though challenging, question. Recombination rates are widely observed to negatively associate with the abundance of transposable elements (TEs), selfish genetic elements that move between genomic locations. While such associations are usually interpreted as recombination influencing the efficacy of selection at removing TEs, accumulating findings suggest that TEs could instead be the cause rather than the consequence. To test this prediction, we formally investigated the influence of polymorphic, putatively active TEs on recombination rates. We developed and benchmarked a novel approach that uses PacBio long-read sequencing to efficiently, accurately, and cost-effectively identify crossovers (COs), a key recombination product, among large numbers of pooled recombinant individuals. By applying this approach to Drosophila strains with distinct TE insertion profiles, we found that polymorphic TEs, especially RNA-based TEs and TEs with local enrichment of repressive marks, reduce the occurrence of COs. Such an effect leads to different CO frequencies between homologous sequences with and without TEs, contributing to varying CO maps between individuals. The suppressive effect of TEs on CO is further supported by two orthogonal approaches-analyzing the distributions of COs in panels of recombinant inbred lines in relation to TE polymorphism and applying marker-assisted estimations of CO frequencies to isogenic strains with and without transgenically inserted TEs. Our investigations reveal how the constantly changing mobilome can actively modify recombination landscapes, shaping genome evolution within and between species.

3D chromatin structures associated with ncRNA roX2 for hyperactivation and coactivation across the entire X chromosome

The three-dimensional (3D) organization of chromatin within the nucleus is crucial for gene regulation. However, the 3D architectural features that coordinate the activation of an entire chromosome remain largely unknown. We introduce an omics method, RNA-associated chromatin DNA-DNA interactions, that integrates RNA polymerase II (RNAPII)-mediated regulome with stochastic optical reconstruction microscopy to investigate the landscape of noncoding RNA roX2-associated chromatin topology for gene equalization to achieve dosage compensation. Our findings reveal that roX2 anchors to the target gene transcription end sites (TESs) and spreads in a distinctive boot-shaped configuration, promoting a more open chromatin state for hyperactivation. Furthermore, roX2 arches TES to transcription start sites to enhance transcriptional loops, potentially facilitating RNAPII convoying and connecting proximal promoter-promoter transcriptional hubs for synergistic gene regulation. These TESs cluster as roX2 compartments, surrounded by inactive domains for coactivation of multiple genes within the roX2 territory. In addition, roX2 structures gradually form and scaffold for stepwise coactivation in dosage compensation.

"Off-pore" nucleoporins relocalize heterochromatic breaks through phase separation

Phase separation forms membraneless compartments in the nuclei, including by establishing heterochromatin "domains" and repair foci. Pericentromeric heterochromatin mostly comprises repeated sequences prone to aberrant recombination, and "safe" homologous recombination (HR) repair of these sequences requires the movement of repair sites to the nuclear periphery before Rad51 recruitment and strand invasion. How this mobilization initiates is unknown, and the contribution of phase separation to these dynamics is unclear. Here, we show that Nup98 nucleoporin is recruited to heterochromatic repair sites before relocalization through Sec13 or Nup88 nucleoporins, and downstream from the Smc5/6 complex and SUMOylation. Remarkably, the phase separation properties of Nup98 are required and sufficient to mobilize repair sites and exclude Rad51, thus preventing aberrant recombination while promoting HR repair. Disrupting this pathway results in heterochromatin repair defects and widespread chromosome rearrangements, revealing a novel "off-pore" role for nucleoporins and phase separation in nuclear dynamics and genome integrity in a multicellular eukaryote.

Rapid evolution of piRNA clusters in the Drosophila melanogaster ovary

The piRNA pathway is a highly conserved mechanism to repress transposable element (TE) activity in the animal germline via a specialized class of small RNAs called piwi-interacting RNAs (piRNAs). piRNAs are produced from discrete genomic regions called piRNA clusters (piCs). Although the molecular processes by which piCs function are relatively well understood in Drosophila melanogaster, much less is known about the origin and evolution of piCs in this or any other species. To investigate piC origin and evolution, we use a population genomic approach to compare piC activity and sequence composition across eight geographically distant strains of D. melanogaster with high-quality long-read genome assemblies. We perform annotations of ovary piCs and genome-wide TE content in each strain. Our analysis uncovers extensive variation in piC activity across strains and signatures of rapid birth and death of piCs. Most TEs inferred to be recently active show an enrichment of insertions into old and large piCs, consistent with the previously proposed "trap" model of piC evolution. In contrast, a small subset of active LTR families is enriched for the formation of new piCs, suggesting that these TEs have higher proclivity to form piCs. Thus, our findings uncover processes leading to the origin of piCs. We propose that piC evolution begins with the emergence of piRNAs from individual insertions of a few select TE families prone to seed new piCs that subsequently expand by accretion of insertions from most other TE families during evolution to form larger "trap" clusters. Our study shows that TEs themselves are the major force driving the rapid evolution of piCs.

Heterochromatic 3D genome organization is directed by HP1a- and H3K9-dependent and independent mechanisms

Whether and how histone post-translational modifications and the proteins that bind them drive 3D genome organization remains unanswered. Here, we evaluate the contribution of H3K9-methylated constitutive heterochromatin to 3D genome organization in Drosophila tissues. We find that the predominant organizational feature of wild-type tissues is the segregation of euchromatic chromosome arms from heterochromatic pericentromeres. Reciprocal perturbation of HP1a⋅H3K9me binding, using a point mutation in the HP1a chromodomain or replacement of the replication-dependent histone H3 with H3K9R mutant histones, revealed that HP1a binding to methylated H3K9 in constitutive heterochromatin is required to limit contact frequency between pericentromeres and chromosome arms and regulate the distance between arm and pericentromeric regions. Surprisingly, the self-association of pericentromeric regions is largely preserved despite the loss of H3K9 methylation and HP1a occupancy. Thus, the HP1a⋅H3K9 interaction contributes to but does not solely drive the segregation of euchromatin and heterochromatin inside the nucleus.

Nup153 is not required for anchoring heterochromatic DSBs to the nuclear periphery

Pericentromeric heterochromatin mostly comprises repeated DNA sequences prone to ectopic recombination. In Drosophila cells, 'safe' homologous recombination repair requires relocalization of heterochromatic repair sites to the nuclear periphery before Rad51 recruitment and strand invasion. DSBs are anchored to the nuclear periphery through the Nup107/160 nucleoporin complex. Previous studies suggested that the nuclear pore 'basket' protein Nup153 could also mediate anchoring, but Nup153 RNAi depletion also affects Nup107 association with the pores, preventing a direct assessment of Nup153 role. Using a separation of function mutant, here we show that Nup153 is not required for anchoring heterochromatic DSBs to the nuclear periphery.

Alternative end-joining results in smaller deletions in heterochromatin relative to euchromatin

Pericentromeric heterochromatin is highly enriched for repetitive sequences prone to aberrant recombination. Previous studies showed that homologous recombination (HR) repair is uniquely regulated in this domain to enable 'safe' repair while preventing aberrant recombination. In Drosophila cells, DNA double-strand breaks (DSBs) relocalize to the nuclear periphery through nuclear actin-driven directed motions before recruiting the strand invasion protein Rad51 and completing HR repair. End-joining (EJ) repair also occurs with high frequency in heterochromatin of fly tissues, but how alternative EJ (alt-EJ) pathways operate in heterochromatin remains largely uncharacterized. Here, we induce DSBs in single euchromatic and heterochromatic sites using a new system that combines the DR- white reporter and I-SceI expression in spermatogonia of flies. Using this approach, we detect higher frequency of HR repair in heterochromatin, relative to euchromatin. Further, sequencing of mutagenic repair junctions reveals the preferential use of different EJ pathways across distinct euchromatic and heterochromatic sites. Interestingly, synthesis-dependent microhomology-mediated end joining (SD-MMEJ) appears differentially regulated in the two domains, with a preferential use of motifs close to the cut site in heterochromatin relative to euchromatin, resulting in smaller deletions. Together, these studies establish a new approach to study repair outcomes in fly tissues, and support the conclusion that heterochromatin uses more HR and less mutagenic EJ repair relative to euchromatin.

Rapid evolution of piRNA clusters in the Drosophila melanogaster ovary

Animal genomes are parasitized by a horde of transposable elements (TEs) whose mutagenic activity can have catastrophic consequences. The piRNA pathway is a conserved mechanism to repress TE activity in the germline via a specialized class of small RNAs associated with effector Piwi proteins called piwi-associated RNAs (piRNAs). piRNAs are produced from discrete genomic regions called piRNA clusters (piCs). While piCs are generally enriched for TE sequences and the molecular processes by which they are transcribed and regulated are relatively well understood in Drosophila melanogaster, much less is known about the origin and evolution of piCs in this or any other species. To investigate piC evolution, we use a population genomics approach to compare piC activity and sequence composition across 8 geographically distant strains of D. melanogaster with high quality long-read genome assemblies. We perform extensive annotations of ovary piCs and TE content in each strain and test predictions of two proposed models of piC evolution. The 'de novo' model posits that individual TE insertions can spontaneously attain the status of a small piC to generate piRNAs silencing the entire TE family. The 'trap' model envisions large and evolutionary stable genomic clusters where TEs tend to accumulate and serves as a long-term "memory" of ancient TE invasions and produce a great variety of piRNAs protecting against related TEs entering the genome. It remains unclear which model best describes the evolution of piCs. Our analysis uncovers extensive variation in piC activity across strains and signatures of rapid birth and death of piCs in natural populations. Most TE families inferred to be recently or currently active show an enrichment of strain-specific insertions into large piCs, consistent with the trap model. By contrast, only a small subset of active LTR retrotransposon families is enriched for the formation of strain-specific piCs, suggesting that these families have an inherent proclivity to form de novo piCs. Thus, our findings support aspects of both 'de novo' and 'trap' models of piC evolution. We propose that these two models represent two extreme stages along an evolutionary continuum, which begins with the emergence of piCs de novo from a few specific LTR retrotransposon insertions that subsequently expand by accretion of other TE insertions during evolution to form larger 'trap' clusters. Our study shows that piCs are evolutionarily labile and that TEs themselves are the major force driving the formation and evolution of piCs.

Genetic variation in chromatin state across multiple tissues in Drosophila melanogaster

We use ATAC-seq to examine chromatin accessibility for four different tissues in Drosophila melanogaster: adult female brain, ovaries, and both wing and eye-antennal imaginal discs from males. Each tissue is assayed in eight different inbred strain genetic backgrounds, seven associated with a reference quality genome assembly. We develop a method for the quantile normalization of ATAC-seq fragments and test for differences in coverage among genotypes, tissues, and their interaction at 44099 peaks throughout the euchromatic genome. For the strains with reference quality genome assemblies, we correct ATAC-seq profiles for read mis-mapping due to nearby polymorphic structural variants (SVs). Comparing coverage among genotypes without accounting for SVs results in a highly elevated rate (55%) of identifying false positive differences in chromatin state between genotypes. After SV correction, we identify 1050, 30383, and 4508 regions whose peak heights are polymorphic among genotypes, among tissues, or exhibit genotype-by-tissue interactions, respectively. Finally, we identify 3988 candidate causative variants that explain at least 80% of the variance in chromatin state at nearby ATAC-seq peaks.

Actin-related protein 6 facilitates proneural protein-induced gene activation for rapid neural differentiation

Neurogenesis is initiated by basic helix-loop-helix proneural proteins. Here, we show that Actin-related protein 6 (Arp6), a core component of the H2A.Z exchange complex SWR1, interacts with proneural proteins and is crucial for efficient onset of proneural protein target gene expression. Arp6 mutants exhibit reduced transcription in sensory organ precursors (SOPs) downstream of the proneural protein patterning event. This leads to retarded differentiation and division of SOPs and smaller sensory organs. These phenotypes are also observed in proneural gene hypomorphic mutants. Proneural protein expression is not reduced in Arp6 mutants. Enhanced proneural gene expression fails to rescue retarded differentiation in Arp6 mutants, suggesting that Arp6 acts downstream of or in parallel with proneural proteins. H2A.Z mutants display Arp6-like retardation in SOPs. Transcriptomic analyses demonstrate that loss of Arp6 and H2A.Z preferentially decreases expression of proneural protein-activated genes. H2A.Z enrichment in nucleosomes around the transcription start site before neurogenesis correlates highly with greater activation of proneural protein target genes by H2A.Z. We propose that upon proneural protein binding to E-box sites, H2A.Z incorporation around the transcription start site allows rapid and efficient activation of target genes, promoting rapid neural differentiation.

Epigenetic Silencing of P-Element Reporter Genes Induced by Transcriptionally Active Domains of Constitutive Heterochromatin in Drosophila melanogaster

Reporter genes inserted via P-element integration into different locations of the Drosophila melanogaster genome have been routinely used to monitor the functional state of chromatin domains. It is commonly thought that P-element-derived reporter genes are subjected to position effect variegation (PEV) when transposed into constitutive heterochromatin because they acquire heterochromatin-like epigenetic modifications that promote silencing. However, sequencing and annotation of the D. melanogaster genome have shown that constitutive heterochromatin is a genetically and molecularly heterogeneous compartment. In fact, in addition to repetitive DNAs, it harbors hundreds of functional genes, together accounting for a significant fraction of its entire genomic territory. Notably, most of these genes are actively transcribed in different developmental stages and tissues, irrespective of their location in heterochromatin. An open question in the genetic and molecular studies on PEV in D. melanogaster is whether functional heterochromatin domains, i.e., heterochromatin harboring active genes, are able to silence reporter genes therein transposed or, on the contrary, can drive their expression. In this work, we provide experimental evidence showing that strong silencing of the Pw⁺ reporters is induced even when they are integrated within or near actively transcribed loci in the pericentric regions of chromosome 2. Interestingly, some Pw⁺ reporters were found insensitive to the action of a known PEV suppressor. Two of them are inserted within Yeti, a gene expressed in the deep heterochromatin of chromosome 2 which carries active chromatin marks. The difference sensitivity to suppressors-exhibited Pw⁺ reporters supports the view that different epigenetic regulators or mechanisms control different regions of heterochromatin. Together, our results suggest that there may be more complexity regarding the molecular mechanisms underlying PEV.

Genetic variation in P-element dysgenic sterility is associated with double-strand break repair and alternative splicing of TE transcripts

The germline mobilization of transposable elements (TEs) by small RNA mediated silencing pathways is conserved across eukaryotes and critical for ensuring the integrity of gamete genomes. However, genomes are recurrently invaded by novel TEs through horizontal transfer. These invading TEs are not targeted by host small RNAs, and their unregulated activity can cause DNA damage in germline cells and ultimately lead to sterility. Here we use hybrid dysgenesis-a sterility syndrome of Drosophila caused by transposition of invading P-element DNA transposons-to uncover host genetic variants that modulate dysgenic sterility. Using a panel of highly recombinant inbred lines of Drosophila melanogaster, we identified two linked quantitative trait loci (QTL) that determine the severity of dysgenic sterility in young and old females, respectively. We show that ovaries of fertile genotypes exhibit increased expression of splicing factors that suppress the production of transposase encoding transcripts, which likely reduces the transposition rate and associated DNA damage. We also show that fertile alleles are associated with decreased sensitivity to double-stranded breaks and enhanced DNA repair, explaining their ability to withstand high germline transposition rates. Together, our work reveals a diversity of mechanisms whereby host genotype modulates the cost of an invading TE, and points to genetic variants that were likely beneficial during the P-element invasion.

The chromatin factor ROW cooperates with BEAF-32 in regulating long-range inducible genes

Insulator proteins located at the boundaries of topological associated domains (TAD) are involved in higher-order chromatin organization and transcription regulation. However, it is still not clear how long-range contacts contribute to transcriptional regulation. Here, we show that relative-of-WOC (ROW) is essential for the long-range transcription regulation mediated by the boundary element-associated factor of 32kD (BEAF-32). We find that ROW physically interacts with heterochromatin proteins (HP1b and HP1c) and the insulator protein (BEAF-32). These proteins interact at TAD boundaries where ROW, through its AT-hook motifs, binds AT-rich sequences flanked by BEAF-32-binding sites and motifs. Knockdown of row downregulates genes that are long-range targets of BEAF-32 and bound indirectly by ROW (without binding motif). Analyses of high-throughput chromosome conformation capture (Hi-C) data reveal long-range interactions between promoters of housekeeping genes bound directly by ROW and promoters of developmental genes bound indirectly by ROW. Thus, our results show cooperation between BEAF-32 and the ROW complex, including HP1 proteins, to regulate the transcription of developmental and inducible genes through long-range interactions.

Light, chromatin, action: nuclear events regulating light signaling in Arabidopsis

The plant nucleus provides a major hub for environmental signal integration at the chromatin level. Multiple light signaling pathways operate and exchange information by regulating a large repertoire of gene targets that shape plant responses to a changing environment. In addition to the established role of transcription factors in triggering photoregulated changes in gene expression, there are eminent reports on the significance of chromatin regulators and nuclear scaffold dynamics in promoting light-induced plant responses. Here, we report and discuss recent advances in chromatin-regulatory mechanisms modulating plant architecture and development in response to light, including the molecular and physiological roles of key modifications such as DNA, RNA and histone methylation, and/or acetylation. The significance of the formation of biomolecular condensates of key light signaling components is discussed and potential applications to agricultural practices overviewed.

Species-specific chromatin landscape determines how transposable elements shape genome evolution

Transposable elements (TEs) are selfish genetic parasites that increase their copy number at the expense of host fitness. The ‘success’, or genome-wide abundance, of TEs differs widely between species. Deciphering the causes for this large variety in TE abundance has remained a central question in evolutionary genomics. We previously proposed that species-specific TE abundance could be driven by the inadvertent consequences of host-direct epigenetic silencing of TEs—the spreading of repressive epigenetic marks from silenced TEs into adjacent sequences. Here, we compared this TE-mediated local enrichment of repressive marks, or ‘the epigenetic effect of TEs’, in six species in the Drosophila melanogaster subgroup to dissect step-by-step the role of such effect in determining genomic TE abundance. We found that TE-mediated local enrichment of repressive marks is prevalent and substantially varies across and even within species. While this TE-mediated effect alters the epigenetic states of adjacent genes, we surprisingly discovered that the transcription of neighboring genes could reciprocally impact this spreading. Importantly, our multi-species analysis provides the power and appropriate phylogenetic resolution to connect species-specific host chromatin regulation, TE-mediated epigenetic effects, the strength of natural selection against TEs, and genomic TE abundance unique to individual species. Our findings point toward the importance of host chromatin landscapes in shaping genome evolution through the epigenetic effects of a selfish genetic parasite.

All the instructions required for life are encoded in the set of DNA present in a cell. It therefore seems natural to think that every bit of this genetic information should serve the organism. And yet most species carry parasitic ‘transposable’ sequences, or transposons, whose only purpose is to multiply and insert themselves at other positions in the genome.

It is possible for cells to suppress these selfish elements. Chemical marks can be deposited onto the DNA to temporarily ‘silence’ transposons and prevent them from being able to move and replicate. However, this sometimes comes at a cost: the repressive chemical modifications can spread to nearby genes that are essential for the organism and perturb their function.

Strangely, the prevalence of transposons varies widely across the tree of life. These sequences form the majority of the genome of certain species – in fact, they represent about half of the human genetic information. But their abundance is much lower in other organisms, forming a measly 6% of the genome of puffer fish for instance. Even amongst fruit fly species, the prevalence of transposable elements can range between 2% and 25%. What explains such differences?

Huang et al. set out to examine this question through the lens of transposon silencing, systematically comparing how this process impacts nearby regions in six species of fruit flies. This revealed variations in the strength of the side effects associated with transposon silencing, resulting in different levels of perturbation on neighbouring genes. A stronger impact was associated with the species having fewer transposons in its genome, suggesting that an evolutionary pressure is at work to keep the abundance of transposons at a low level in these species. Further analyses showed that the genes which determine how silencing marks are distributed may also be responsible for the variations in the impact of transposon silencing. They could therefore be the ones driving differences in the abundance of transposons between species.

Overall, this work sheds light on the complex mechanisms shaping the evolution of genomes, and it may help to better understand how transposons are linked to processes such as aging and cancer.

The role of insulators and transcription in 3D chromatin organization of flies

The DNA in many organisms, including humans, is shown to be organized in topologically associating domains (TADs). In Drosophila, several architectural proteins are enriched at TAD borders, but it is still unclear whether these proteins play a functional role in the formation and maintenance of TADs. Here, we show that depletion of BEAF-32, Cp190, Chro, and Dref leads to changes in TAD organization and chromatin loops. Their depletion predominantly affects TAD borders located in regions moderately enriched in repressive modifications and depleted in active ones, whereas TAD borders located in euchromatin are resilient to these knockdowns. Furthermore, transcriptomic data has revealed hundreds of genes displaying differential expression in these knockdowns and showed that the majority of differentially expressed genes are located within reorganized TADs. Our work identifies a novel and functional role for architectural proteins at TAD borders in Drosophila and a link between TAD reorganization and subsequent changes in gene expression.

An Expanding Toolkit for Heterochromatin Repair Studies

Pericentromeric heterochromatin is mostly composed of repetitive DNA sequences prone to aberrant recombination. Cells have developed highly specialized mechanisms to enable ‘safe’ homologous recombination (HR) repair while preventing aberrant recombination in this domain. Understanding heterochromatin repair responses is essential to understanding the critical mechanisms responsible for genome integrity and tumor suppression. Here, we review the tools, approaches, and methods currently available to investigate double-strand break (DSB) repair in pericentromeric regions, and also suggest how technologies recently developed for euchromatin repair studies can be adapted to characterize responses in heterochromatin. With this ever-growing toolkit, we are witnessing exciting progress in our understanding of how the ‘dark matter’ of the genome is repaired, greatly improving our understanding of genome stability mechanisms.

Synergistic epistasis of the deleterious effects of transposable elements

The replicative nature and generally deleterious effects of transposable elements (TEs) raise an outstanding question about how TE copy number is stably contained in host populations. Classic theoretical analyses predict that, when the decline in fitness due to each additional TE insertion is greater than linear, or when there is synergistic epistasis, selection against TEs can result in a stable equilibrium of TE copy number. While several mechanisms are predicted to yield synergistic deleterious effects of TEs, we lack empirical investigations of the presence of such epistatic interactions. Purifying selection with synergistic epistasis generates repulsion linkage between deleterious alleles. We investigated this population genetic signal in the likely ancestral Drosophila melanogaster population and found evidence supporting the presence of synergistic epistasis among TE insertions, especially TEs expected to exert large fitness impacts. Even though synergistic epistasis of TEs has been predicted to arise through ectopic recombination and TE-mediated epigenetic silencing mechanisms, we only found mixed support for the associated predictions. We observed signals of synergistic epistasis for a large number of TE families, which is consistent with the expectation that such epistatic interaction mainly happens among copies of the same family. Curiously, significant repulsion linkage was also found among TE insertions from different families, suggesting the possibility that synergism of TEs' deleterious fitness effects could arise above the family level and through mechanisms similar to those of simple mutations. Our findings set the stage for investigating the prevalence and importance of epistatic interactions in the evolutionary dynamics of TEs.

Patterns of piRNA Regulation in Drosophila Revealed through Transposable Element Clade Inference

Transposable elements (TEs) are self-replicating "genetic parasites" ubiquitous to eukaryotic genomes. In addition to conflict between TEs and their host genomes, TEs of the same family are in competition with each other. They compete for the same genomic niches while experiencing the same regime of copy-number selection. This suggests that competition among TEs may favor the emergence of new variants that can outcompete their ancestral forms. To investigate the sequence evolution of TEs, we developed a method to infer clades: collections of TEs that share SNP variants and represent distinct TE family lineages. We applied this method to a panel of 85 Drosophila melanogaster genomes and found that the genetic variation of several TE families shows significant population structure that arises from the population-specific expansions of single clades. We used population genetic theory to classify these clades into younger versus older clades and found that younger clades are associated with a greater abundance of sense and antisense piRNAs per copy than older ones. Further, we find that the abundance of younger, but not older clades, is positively correlated with antisense piRNA production, suggesting a general pattern where hosts preferentially produce antisense piRNAs from recently active TE variants. Together these findings suggest a pattern whereby new TE variants arise by mutation and then increase in copy number, followed by the host producing antisense piRNAs that may be used to silence these emerging variants.

Nuclear Actin Dynamics in Gene Expression, DNA Repair, and Cancer

Actin is a highly conserved protein in mammals. The actin dynamics is regulated by actin-binding proteins and actin-related proteins. Nuclear actin and these regulatory proteins participate in multiple nuclear processes, including chromosome architecture organization, chromatin remodeling, transcription machinery regulation, and DNA repair. It is well known that the dysfunctions of these processes contribute to the development of cancer. Moreover, emerging evidence has shown that the deregulated actin dynamics is also related to cancer. This chapter discusses how the deregulation of nuclear actin dynamics contributes to tumorigenesis via such various nuclear events.

OUP accepted manuscript

Interspecific hybridization is often seen as a genomic stress that may lead to new gene expression patterns and deregulation of transposable elements (TEs). The understanding of expression changes in hybrids compared with parental species is essential to disentangle their putative role in speciation processes. However, to date we ignore the detailed mechanisms involved in genomic deregulation in hybrids. We studied the ovarian transcriptome and epigenome of the Drosophila buzzatii and Drosophila koepferae species together with their F₁ hybrid females. We found a trend toward underexpression of genes and TE families in hybrids. The epigenome in hybrids was highly similar to the parental epigenomes and showed intermediate histone enrichments between parental species in most cases. Differential gene expression in hybrids was often associated only with changes in H3K4me3 enrichments, whereas differential TE family expression in hybrids may be associated with changes in H3K4me3, H3K9me3, or H3K27me3 enrichments. We identified specific genes and TE families, which their differential expression in comparison with the parental species was explained by their differential chromatin mark combination enrichment. Finally, cis–trans compensatory regulation could also contribute in some way to the hybrid deregulation. This work provides the first study of histone content in Drosophila interspecific hybrids and their effect on gene and TE expression deregulation.

H3K9me2 genome-wide distribution in the holocentric insect Spodoptera frugiperda (Lepidoptera: Noctuidae)

BACKGROUND

Eukaryotic genomes are packaged by Histone proteins in a structure called chromatin. There are different chromatin types. Euchromatin is typically associated with decondensed, transcriptionally active regions and heterochromatin to more condensed regions of the chromosomes. Methylation of Lysine 9 of Histone H3 (H3K9me) is a conserved biochemical marker of heterochromatin. In many organisms, heterochromatin is usually localized at telomeric as well as pericentromeric regions but can also be found at interstitial chromosomal loci. This distribution may vary in different species depending on their general chromosomal organization. Holocentric species such as Spodoptera frugiperda (Lepidoptera: Noctuidae) possess dispersed centromeres instead of a monocentric one and thus no observable pericentromeric compartment. To identify the localization of heterochromatin in such species we performed ChIP-Seq experiments and analyzed the distribution of the heterochromatin marker H3K9me2 in the Sf9 cell line and whole 4th instar larvae (L4) in relation to RNA-Seq data.

RESULTS

In both samples we measured an enrichment of H3K9me2 at the (sub) telomeres, rDNA loci, and satellite DNA sequences, which could represent dispersed centromeric regions. We also observed that density of H3K9me2 is positively correlated with transposable elements and protein-coding genes. But contrary to most model organisms, H3K9me2 density is not correlated with transcriptional repression.

CONCLUSION

This is the first genome-wide ChIP-Seq analysis conducted in S. frugiperda for H3K9me2. Compared to model organisms, this mark is found in expected chromosomal compartments such as rDNA and telomeres. However, it is also localized at numerous dispersed regions, instead of the well described large pericentromeric domains, indicating that H3K9me2 might not represent a classical heterochromatin marker in Lepidoptera. (242 words).

Complex Genetic Interactions between Piwi and HP1a in the Repression of Transposable Elements and Tissue-Specific Genes in the Ovarian Germline

Insertions of transposable elements (TEs) in eukaryotic genomes are usually associated with repressive chromatin, which spreads to neighbouring genomic sequences. In ovaries of Drosophila melanogaster, the Piwi-piRNA pathway plays a key role in the transcriptional silencing of TEs considered to be exerted mostly through the establishment of H3K9me3 histone marks recruiting Heterochromatin Protein 1a (HP1a). Here, using RNA-seq, we investigated the expression of TEs and the adjacent genomic regions upon Piwi and HP1a germline knockdowns sharing a similar genetic background. We found that the depletion of Piwi and HP1a led to the derepression of only partially overlapping TE sets. Several TEs were silenced predominantly by HP1a, whereas the upregulation of some other TEs was more pronounced upon Piwi knockdown and, surprisingly, was diminished upon a Piwi/HP1a double-knockdown. We revealed that HP1a loss influenced the expression of thousands of protein-coding genes mostly not adjacent to TE insertions and, in particular, downregulated a putative transcriptional factor required for TE activation. Nevertheless, our results indicate that Piwi and HP1a cooperatively exert repressive effects on the transcription of euchromatic loci flanking the insertions of some Piwi-regulated TEs. We suggest that this mechanism controls the silencing of a small set of TE-adjacent tissue-specific genes, preventing their inappropriate expression in ovaries.

Is There a Histone Code for Cellular Quiescence?

Many of the cells in our bodies are quiescent, that is, temporarily not dividing. Under certain physiological conditions such as during tissue repair and maintenance, quiescent cells receive the appropriate stimulus and are induced to enter the cell cycle. The ability of cells to successfully transition into and out of a quiescent state is crucial for many biological processes including wound healing, stem cell maintenance, and immunological responses. Across species and tissues, transcriptional, epigenetic, and chromosomal changes associated with the transition between proliferation and quiescence have been analyzed, and some consistent changes associated with quiescence have been identified. Histone modifications have been shown to play a role in chromatin packing and accessibility, nucleosome mobility, gene expression, and chromosome arrangement. In this review, we critically evaluate the role of different histone marks in these processes during quiescence entry and exit. We consider different model systems for quiescence, each of the most frequently monitored candidate histone marks, and the role of their writers, erasers and readers. We highlight data that support these marks contributing to the changes observed with quiescence. We specifically ask whether there is a quiescence histone “code,” a mechanism whereby the language encoded by specific combinations of histone marks is read and relayed downstream to modulate cell state and function. We conclude by highlighting emerging technologies that can be applied to gain greater insight into the role of a histone code for quiescence.

An explainable artificial intelligence approach for decoding the enhancer histone modifications code and identification of novel enhancers in Drosophila

BACKGROUND

Enhancers are non-coding regions of the genome that control the activity of target genes. Recent efforts to identify active enhancers experimentally and in silico have proven effective. While these tools can predict the locations of enhancers with a high degree of accuracy, the mechanisms underpinning the activity of enhancers are often unclear.

RESULTS

Using machine learning (ML) and a rule-based explainable artificial intelligence (XAI) model, we demonstrate that we can predict the location of known enhancers in Drosophila with a high degree of accuracy. Most importantly, we use the rules of the XAI model to provide insight into the underlying combinatorial histone modifications code of enhancers. In addition, we identified a large set of putative enhancers that display the same epigenetic signature as enhancers identified experimentally. These putative enhancers are enriched in nascent transcription, divergent transcription and have 3D contacts with promoters of transcribed genes. However, they display only intermediary enrichment of mediator and cohesin complexes compared to previously characterised active enhancers. We also found that 10-15% of the predicted enhancers display similar characteristics to super enhancers observed in other species.

CONCLUSIONS

Here, we applied an explainable AI model to predict enhancers with high accuracy. Most importantly, we identified that different combinations of epigenetic marks characterise different groups of enhancers. Finally, we discovered a large set of putative enhancers which display similar characteristics with previously characterised active enhancers.

The Organization of Pericentromeric Heterochromatin in Polytene Chromosome 3 of the Drosophilamelanogaster Line with the Rif11; SuURES Su(var)3-906 Mutations Suppressing Underreplication

Although heterochromatin makes up 40% of the Drosophila melanogaster genome, its organization remains little explored, especially in polytene chromosomes, as it is virtually not represented in them due to underreplication. Two all-new approaches were used in this work: (i) with the use of a newly synthesized Drosophila line that carries three mutations, Rif11, SuURES and Su(var)3-906, suppressing the underreplication of heterochromatic regions, we obtained their fullest representation in polytene chromosomes and described their structure; (ii) 20 DNA fragments with known positions on the physical map as well as molecular genetic features of the genome (gene density, histone marks, heterochromatin proteins, origin recognition complex proteins, replication timing sites and satellite DNAs) were mapped in the newly polytenized heterochromatin using FISH and bioinformatics data. The borders of the heterochromatic regions and variations in their positions on arm 3L have been determined for the first time. The newly polytenized heterochromatic material exhibits two main types of morphology: a banding pattern (locations of genes and short satellites) and reticular chromatin (locations of large blocks of satellite DNA). The locations of the banding and reticular polytene heterochromatin was determined on the physical map.

Effects of demethylase KDM4B on the biological characteristics and function of yak cumulus cells in vitro

The present study aims to investigate the expression and function of lysine-specific demethylase 4B (KDM4B) in yak cumulus cells (CCs) in order to reveal the mechanisms by which KDM4B regulates biological characteristics and function of CCs. The cellular location of KDM4B and the methylation pattern of H3K9 were detected using immunofluorescence (IF) staining in CCs. The mRNA expression levels of apoptosis-related genes (BCL-2, HAX1 and BAX) and genes related to the estrogen pathway (ESR2, CYP17 and 3B-HSD) were estimated by qRT-PCR after knockdown of KDM4B expression by siRNA in yak CCs. Then, a proliferation assay, Annexin V-FITC staining, and ELISA were utilized to explore the effects of KDM4B silencing on CCs proliferation, apoptosis, and estrogen (E₂) secretion, respectively. The results showed that KDM4B is located in the nuclei of yak CCs and is distributed in a dotted pattern. Knockdown KDM4B induced a decrease in cell proliferation, an increase in apoptotic rate and a reduction in the levels of E₂ secretion of CCs. Additionally, the methylation patterns of H3K9me2 and H3K9me3 were significantly increased in CCs transfected with KDM4B siRNA-1 (P < 0.05). The mRNA expression level of apoptosis promoting BAX genes was significantly upregulated, but 3B-HSD, ESR2 and anti-apoptotic HAX1 genes were significantly downregulated in transfected CCs (P < 0.05). Furthermore, the rate of embryos developing from the 2-cell stage to blastocysts was lower in the siRNA-1 transfection group than that of the control group (28.6 ± 2.9% vs 40.4 ± 2.4%, P < 0.05). In conclusion, our study indicates that KDM4B regulates the biological characteristics and physiological function of yak CCs mainly through changing the methylation patterns of H3K9 and related gene expression levels.

The role of H3K36 methylation and associated methyltransferases in chromosome-specific gene regulation

In Drosophila, two chromosomes require special mechanisms to balance their transcriptional output to the rest of the genome. These are the male-specific lethal complex targeting the male X chromosome and Painting of fourth targeting chromosome 4. Here, we explore the role of histone H3 methylated at lysine-36 (H3K36) and the associated methyltransferases—Set2, NSD, and Ash1—in these two chromosome-specific systems. We show that the loss of Set2 impairs the MSL complex–mediated dosage compensation; however, the effect is not recapitulated by H3K36 replacement and indicates an alternative target of Set2. Unexpectedly, balanced transcriptional output from the fourth chromosome requires intact H3K36 and depends on the additive functions of NSD and Ash1. We conclude that H3K36 methylation and the associated methyltransferases are important factors to balance transcriptional output of the male X chromosome and the fourth chromosome. Furthermore, our study highlights the pleiotropic effects of these enzymes.

Chromatin Organization and Function in Drosophila

Eukaryotic genomes are packaged into high-order chromatin structures organized in discrete territories inside the cell nucleus, which is surrounded by the nuclear envelope acting as a barrier. This chromatin organization is complex and dynamic and, thus, determining the spatial and temporal distribution and folding of chromosomes within the nucleus is critical for understanding the role of chromatin topology in genome function. Primarily focusing on the regulation of gene expression, we review here how the genome of Drosophila melanogaster is organized into the cell nucleus, from small scale histone–DNA interactions to chromosome and lamina interactions in the nuclear space.

The Drosophila HP1 family is associated with active gene expression across chromatin contexts

Drosophila Heterochromatin Protein 1a (HP1a) is essential for heterochromatin formation and is involved in transcriptional silencing. However, certain loci require HP1a to be transcribed. One model posits that HP1a acts as a transcriptional silencer within euchromatin while acting as an activator within heterochromatin. However, HP1a has been observed as an activator of a set of euchromatic genes. Therefore, it is not clear whether, or how, chromatin context informs the function of HP1 proteins. To understand the role of HP1 proteins in transcription, we examined the genome-wide binding profile of HP1a as well as two other Drosophila HP1 family members, HP1B and HP1C, to determine whether coordinated binding of these proteins is associated with specific transcriptional outcomes. We found that HP1 proteins share many of their endogenous binding targets. These genes are marked by active histone modifications and are expressed at higher levels than nontarget genes in both heterochromatin and euchromatin. In addition, HP1 binding targets displayed increased RNA polymerase pausing compared with nontarget genes. Specifically, colocalization of HP1B and HP1C was associated with the highest levels of polymerase pausing and gene expression. Analysis of HP1 null mutants suggests these proteins coordinate activity at transcription start sites to regulate transcription. Depletion of HP1B or HP1C alters expression of protein-coding genes bound by HP1 family members. Our data broaden understanding of the mechanism of transcriptional activation by HP1a and highlight the need to consider particular protein-protein interactions, rather than broader chromatin context, to predict impacts of HP1 at transcription start sites.

BREC: an R package/Shiny app for automatically identifying heterochromatin boundaries and estimating local recombination rates along chromosomes

Background

Meiotic recombination is a vital biological process playing an essential role in genome's structural and functional dynamics. Genomes exhibit highly various recombination profiles along chromosomes associated with several chromatin states. However, eu-heterochromatin boundaries are not available nor easily provided for non-model organisms, especially for newly sequenced ones. Hence, we miss accurate local recombination rates necessary to address evolutionary questions.

Results

Here, we propose an automated computational tool, based on the Marey maps method, allowing to identify heterochromatin boundaries along chromosomes and estimating local recombination rates. Our method, called BREC (heterochromatin Boundaries and RECombination rate estimates) is non-genome-specific, running even on non-model genomes as long as genetic and physical maps are available. BREC is based on pure statistics and is data-driven, implying that good input data quality remains a strong requirement. Therefore, a data pre-processing module (data quality control and cleaning) is provided. Experiments show that BREC handles different markers' density and distribution issues.

Conclusions

BREC's heterochromatin boundaries have been validated with cytological equivalents experimentally generated on the fruit fly Drosophila melanogaster genome, for which BREC returns congruent corresponding values. Also, BREC's recombination rates have been compared with previously reported estimates. Based on the promising results, we believe our tool has the potential to help bring data science into the service of genome biology and evolution. We introduce BREC within an R-package and a Shiny web-based user-friendly application yielding a fast, easy-to-use, and broadly accessible resource. The BREC R-package is available at the GitHub repository https://github.com/GenomeStructureOrganization.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04233-1.

Heterochromatin-dependent transcription of satellite DNAs in the Drosophila melanogaster female germline

Large blocks of tandemly repeated DNAs-satellite DNAs (satDNAs)-play important roles in heterochromatin formation and chromosome segregation. We know little about how satDNAs are regulated; however, their misregulation is associated with genomic instability and human diseases. We use the Drosophila melanogaster germline as a model to study the regulation of satDNA transcription and chromatin. Here we show that complex satDNAs (>100-bp repeat units) are transcribed into long noncoding RNAs and processed into piRNAs (PIWI interacting RNAs). This satDNA piRNA production depends on the Rhino-Deadlock-Cutoff complex and the transcription factor Moonshiner-a previously described non-canonical pathway that licenses heterochromatin-dependent transcription of dual-strand piRNA clusters. We show that this pathway is important for establishing heterochromatin at satDNAs. Therefore, satDNAs are regulated by piRNAs originating from their own genomic loci. This novel mechanism of satDNA regulation provides insight into the role of piRNA pathways in heterochromatin formation and genome stability.

Interactions and Feedbacks in E-Cadherin Transcriptional Regulation

Epithelial tissues rely on the adhesion between participating cells to retain their integrity. The transmembrane protein E-cadherin is the major protein that mediates homophilic adhesion between neighbouring cells and is, therefore, one of the critical components for epithelial integrity. E-cadherin downregulation has been described extensively as a prerequisite for epithelial-to-mesenchymal transition and is a hallmark in many types of cancer. Due to this clinical importance, research has been mostly focused on understanding the mechanisms leading to transcriptional repression of this adhesion molecule. However, in recent years it has become apparent that re-expression of E-cadherin is a major step in the progression of many cancers during metastasis. Here, we review the currently known molecular mechanisms of E-cadherin transcriptional activation and inhibition and highlight complex interactions between individual mechanisms. We then propose an additional mechanism, whereby the competition between adhesion complexes and heterochromatin protein-1 for binding to STAT92E fine-tunes the levels of E-cadherin expression in Drosophila but also regulates other genes promoting epithelial robustness. We base our hypothesis on both existing literature and our experimental evidence and suggest that such feedback between the cell surface and the nucleus presents a powerful paradigm for epithelial resilience.

Automated generation of context-specific gene regulatory networks with a weighted approach in Drosophila melanogaster

The regulation of gene expression is a key factor in the development and maintenance of life in all organisms. Even so, little is known at whole genome scale for most genes and contexts. We propose a method, Tool for Weighted Epigenomic Networks in Drosophila melanogaster (Fly T-WEoN), to generate context-specific gene regulatory networks starting from a reference network that contains all known gene regulations in the fly. Unlikely regulations are removed by applying a series of knowledge-based filters. Each of these filters is implemented as an independent module that considers a type of experimental evidence, including DNA methylation, chromatin accessibility, histone modifications and gene expression. Fly T-WEoN is based on heuristic rules that reflect current knowledge on gene regulation in D. melanogaster obtained from the literature. Experimental data files can be generated with several standard procedures and used solely when and if available. Fly T-WEoN is available as a Cytoscape application that permits integration with other tools and facilitates downstream network analysis. In this work, we first demonstrate the reliability of our method to then provide a relevant application case of our tool: early development of D. melanogaster. Fly T-WEoN together with its step-by-step guide is available at https://weon.readthedocs.io.

Unraveling the features of somatic transposition in the Drosophila intestine

Transposable elements (TEs) play a significant role in evolution, contributing to genetic variation. However, TE mobilization in somatic cells is not well understood. Here, we address the prevalence of transposition in a somatic tissue, exploiting the Drosophila midgut as a model. Using whole-genome sequencing of in vivo clonally expanded gut tissue, we have mapped hundreds of high-confidence somatic TE integration sites genome-wide. We show that somatic retrotransposon insertions are associated with inactivation of the tumor suppressor Notch, likely contributing to neoplasia formation. Moreover, applying Oxford Nanopore long-read sequencing technology we provide evidence for tissue-specific differences in retrotransposition. Comparing somatic TE insertional activity with transcriptomic and small RNA sequencing data, we demonstrate that transposon mobility cannot be simply predicted by whole tissue TE expression levels or by small RNA pathway activity. Finally, we reveal that somatic TE insertions in the adult fly intestine are enriched in genic regions and in transcriptionally active chromatin. Together, our findings provide clear evidence of ongoing somatic transposition in Drosophila and delineate previously unknown features underlying somatic TE mobility in vivo.

The Drosophila Y Chromosome Affects Heterochromatin Integrity Genome-Wide

The Drosophila Y chromosome is gene poor and mainly consists of silenced, repetitive DNA. Nonetheless, the Y influences expression of hundreds of genes genome-wide, possibly by sequestering key components of the heterochromatin machinery away from other positions in the genome. To test the influence of the Y chromosome on the genome-wide chromatin landscape, we assayed the genomic distribution of histone modifications associated with gene activation (H3K4me3) or heterochromatin (H3K9me2 and H3K9me3) in fruit flies with varying sex chromosome complements (X0, XY, and XYY males; XX and XXY females). Consistent with the general deficiency of active chromatin modifications on the Y, we find that Y gene dose has little influence on the genomic distribution of H3K4me3. In contrast, both the presence and the number of Y chromosomes strongly influence genome-wide enrichment patterns of repressive chromatin modifications. Highly repetitive regions such as the pericentromeres, the dot, and the Y chromosome (if present) are enriched for heterochromatic modifications in wildtype males and females, and even more strongly in X0 flies. In contrast, the additional Y chromosome in XYY males and XXY females diminishes the heterochromatic signal in these normally silenced, repeat-rich regions, which is accompanied by an increase in expression of Y-linked repeats. We find hundreds of genes that are expressed differentially between individuals with aberrant sex chromosome karyotypes, many of which also show sex-biased expression in wildtype Drosophila. Thus, Y chromosomes influence heterochromatin integrity genome-wide, and differences in the chromatin landscape of males and females may also contribute to sex-biased gene expression and sexual dimorphisms.

Tirant Stealthily Invaded Natural Drosophila melanogaster Populations during the Last Century

It was long thought that solely three different transposable elements (TEs)-the I-element, the P-element, and hobo-invaded natural Drosophila melanogaster populations within the last century. By sequencing the "living fossils" of Drosophila research, that is, D. melanogaster strains sampled from natural populations at different time points, we show that a fourth TE, Tirant, invaded D. melanogaster populations during the past century. Tirant likely spread in D. melanogaster populations around 1938, followed by the I-element, hobo, and, lastly, the P-element. In addition to the recent insertions of the canonical Tirant, D. melanogaster strains harbor degraded Tirant sequences in the heterochromatin which are likely due to an ancient invasion, likely predating the split of D. melanogaster and D. simulans. These degraded insertions produce distinct piRNAs that were unable to prevent the novel Tirant invasion. In contrast to the I-element, P-element, and hobo, we did not find that Tirant induces any hybrid dysgenesis symptoms. This absence of apparent phenotypic effects may explain the late discovery of the Tirant invasion. Recent Tirant insertions were found in all investigated natural populations. Populations from Tasmania carry distinct Tirant sequences, likely due to a founder effect. By investigating the TE composition of natural populations and strains sampled at different time points, insertion site polymorphisms, piRNAs, and phenotypic effects, we provide a comprehensive study of a natural TE invasion.

Evolutionary Dynamics of the Pericentromeric Heterochromatin in Drosophila virilis and Related Species

Pericentromeric heterochromatin in Drosophila generally consists of repetitive DNA, forming the environment associated with gene silencing. Despite the expanding knowledge of the impact of transposable elements (TEs) on the host genome, little is known about the evolution of pericentromeric heterochromatin, its structural composition, and age. During the evolution of the Drosophilidae, hundreds of genes have become embedded within pericentromeric regions yet retained activity. We investigated a pericentromeric heterochromatin fragment found in D. virilis and related species, describing the evolution of genes in this region and the age of TE invasion. Regardless of the heterochromatic environment, the amino acid composition of the genes is under purifying selection. However, the selective pressure affects parts of genes in varying degrees, resulting in expansion of gene introns due to TEs invasion. According to the divergence of TEs, the pericentromeric heterochromatin of the species of virilis group began to form more than 20 million years ago by invasions of retroelements, miniature inverted repeat transposable elements (MITEs), and Helitrons. Importantly, invasions into the heterochromatin continue to occur by TEs that fall under the scope of piRNA silencing. Thus, the pericentromeric heterochromatin, in spite of its ability to induce silencing, has the means for being dynamic, incorporating the regions of active transcription.

Live Cell Imaging of Nuclear Actin Filaments and Heterochromatic Repair foci in Drosophila and Mouse Cells

Pericentromeric heterochromatin is mostly composed of repeated DNA sequences, which are prone to aberrant recombination during double-strand break (DSB) repair. Studies in Drosophila and mouse cells revealed that 'safe' homologous recombination (HR) repair of these sequences relies on the relocalization of repair sites to outside the heterochromatin domain before Rad51 recruitment. Relocalization requires a striking network of nuclear actin filaments (F-actin) and myosins that drive directed motions. Understanding this pathway requires the detection of nuclear actin filaments that are significantly less abundant than those in the cytoplasm, and the imaging and tracking of repair sites for long time periods. Here, we describe an optimized protocol for live cell imaging of nuclear F-actin in Drosophila cells, and for repair focus tracking in mouse cells, including: imaging setup, image processing approaches, and analysis methods. We emphasize approaches that can be applied to identify the most effective fluorescent markers for live cell imaging, strategies to minimize photobleaching and phototoxicity with a DeltaVision deconvolution microscope, and image processing and analysis methods using SoftWoRx and Imaris software. These approaches enable a deeper understanding of the spatial and temporal dynamics of heterochromatin repair and have broad applicability in the fields of nuclear architecture, nuclear dynamics, and DNA repair.

A unique chromatin profile defines adaptive genomic regions in a fungal plant pathogen

Genomes store information at scales beyond the linear nucleotide sequence, which impacts genome function at the level of an individual, while influences on populations and long-term genome function remains unclear. Here, we addressed how physical and chemical DNA characteristics influence genome evolution in the plant pathogenic fungus Verticillium dahliae. We identified incomplete DNA methylation of repetitive elements, associated with specific genomic compartments originally defined as Lineage-Specific (LS) regions that contain genes involved in host adaptation. Further chromatin characterization revealed associations with features such as H3 Lys-27 methylated histones (H3K27me3) and accessible DNA. Machine learning trained on chromatin data identified twice as much LS DNA as previously recognized, which was validated through orthogonal analysis, and we propose to refer to this DNA as adaptive genomic regions. Our results provide evidence that specific chromatin profiles define adaptive genomic regions, and highlight how different epigenetic factors contribute to the organization of these regions.

Origins and Evolutionary Patterns of the 1.688 Satellite DNA Family in Drosophila Phylogeny

Satellite DNAs (satDNAs) are a ubiquitous feature of eukaryotic genomes and are usually the major components of constitutive heterochromatin. The 1.688 satDNA, also known as the 359 bp satellite, is one of the most abundant repetitive sequences in Drosophila melanogaster and has been linked to several different biological functions. We investigated the presence and evolution of the 1.688 satDNA in 16 Drosophila genomes. We find that the 1.688 satDNA family is much more ancient than previously appreciated, being shared among part of the melanogaster group that diverged from a common ancestor ∼27 Mya. We found that the 1.688 satDNA family has two major subfamilies spread throughout Drosophila phylogeny (∼360 bp and ∼190 bp). Phylogenetic analysis of ∼10,000 repeats extracted from 14 of the species revealed that the 1.688 satDNA family is present within heterochromatin and euchromatin. A high number of euchromatic repeats are gene proximal, suggesting the potential for local gene regulation. Notably, heterochromatic copies display concerted evolution and a species-specific pattern, whereas euchromatic repeats display a more typical evolutionary pattern, suggesting that chromatin domains may influence the evolution of these sequences. Overall, our data indicate the 1.688 satDNA as the most perduring satDNA family described in Drosophila phylogeny to date. Our study provides a strong foundation for future work on the functional roles of 1.688 satDNA across many Drosophila species.

Interplay of pericentromeric genome organization and chromatin landscape regulates the expression of Drosophila melanogaster heterochromatic genes

Background

Transcription of genes residing within constitutive heterochromatin is paradoxical to the tenets of epigenetic code. The regulatory mechanisms of Drosophila melanogaster heterochromatic gene transcription remain largely unknown. Emerging evidence suggests that genome organization and transcriptional regulation are inter-linked. However, the pericentromeric genome organization is relatively less studied. Therefore, we sought to characterize the pericentromeric genome organization and understand how this organization along with the pericentromeric factors influences heterochromatic gene expression.

Results

Here, we characterized the pericentromeric genome organization in Drosophila melanogaster using 5C sequencing. Heterochromatic topologically associating domains (Het TADs) correlate with distinct epigenomic domains of active and repressed heterochromatic genes at the pericentromeres. These genes are known to depend on the heterochromatic landscape for their expression. However, HP1a or Su(var)3-9 RNAi has minimal effects on heterochromatic gene expression, despite causing significant changes in the global Het TAD organization. Probing further into this observation, we report the role of two other chromatin proteins enriched at the pericentromeres-dMES-4 and dADD1 in regulating the expression of a subset of heterochromatic genes.

Conclusions

Distinct pericentromeric genome organization and chromatin landscapes maintained by the interplay of heterochromatic factors (HP1a, H3K9me3, dMES-4 and dADD1) are sufficient to support heterochromatic gene expression despite the loss of global Het TAD structure. These findings open new avenues for future investigations into the mechanisms of heterochromatic gene expression.

Establishment and evolution of heterochromatin

The eukaryotic genome is packaged into transcriptionally active euchromatin and silent heterochromatin, with most studies focused on the former encompassing the majority of protein-coding genes. The recent development of various sequencing techniques has refined this classic dichromatic partition and has better illuminated the composition, establishment, and evolution of this genomic and epigenomic "dark matter" in the context of topologically associated domains and phase-separated droplets. Heterochromatin includes genomic regions that can be densely stained by chemical dyes, which have been shown to be enriched for repetitive elements and epigenetic marks, including H3K9me2/3 and H3K27me3. Heterochromatin is usually replicated late, concentrated at the nuclear periphery or around nucleoli, and usually lacks highly expressed genes; and now it is considered to be as neither genetically inert nor developmentally static. Heterochromatin guards genome integrity against transposon activities and exerts important regulatory functions by targeting beyond its contained genes. Both its nucleotide sequences and regulatory proteins exhibit rapid coevolution between species. In addition, there are dynamic transitions between euchromatin and heterochromatin during developmental and evolutionary processes. We summarize here the ever-changing characteristics of heterochromatin and propose models and principles for the evolutionary transitions of heterochromatin that have been mainly learned from studies of Drosophila and yeast. Finally, we highlight the role of sex chromosomes in studying heterochromatin evolution.