Chromosome segregation during spermatogenesis occurs through a unique center-kinetic mechanism in holocentric moth species.

Precise regulation of chromosome dynamics in the germline is essential for reproductive success across species. Yet, the mechanisms underlying meiotic chromosomal events such as homolog pairing and chromosome segregation are not fully understood in many species. Here, we employ Oligopaint DNA FISH to investigate mechanisms of meiotic homolog pairing and chromosome segregation in the holocentric pantry moth, Plodia interpunctella, and compare our findings to new and previous studies in the silkworm moth, Bombyx mori, which diverged from P. interpunctella over 100 million years ago. We find that pairing in both Bombyx and Plodia spermatogenesis is initiated at gene-rich chromosome ends. Additionally, both species form rod shaped cruciform-like bivalents at metaphase I. However, unlike the telomere-oriented chromosome segregation mechanism observed in Bombyx, Plodia can orient bivalents in multiple different ways at metaphase I. Surprisingly, in both species we find that kinetochores consistently assemble at non-telomeric loci toward the center of chromosomes regardless of where chromosome centers are located in the bivalent. Additionally, sister kinetochores do not seem to be paired in these species. Instead, four distinct kinetochores are easily observed at metaphase I. Despite this, we find clear end-on microtubule attachments and not lateral microtubule attachments co-orienting these separated kinetochores. These findings challenge the classical view of segregation where paired, poleward-facing kinetochores are required for accurate homolog separation in meiosis I. Our studies here highlight the importance of exploring fundamental processes in non-model systems, as employing novel organisms can lead to the discovery of novel biology.

Dosage compensation in Bombyx mori is achieved by partial repression of both Z chromosomes in males.

SignificanceGenes on sex chromosomes (i.e. human chX) are regulated differently in males and females to balance gene expression levels between sexes (XY vs. XX). This sex-specific regulation is called dosage compensation (DC). DC is achieved by altering the shape and compaction of sex chromosomes specifically in one sex. In this study, we use Oligopaints to examine DC in silkworms. This study visualizes this phenomenon in a species with ZW sex chromosomes, which evolved independently of XY. Our data support a long-standing model for how DC mechanisms evolved across species, and we show potential similarity between DC in silkworms and nematodes, suggesting that this type of DC may have emerged multiple independent times throughout evolution.

Transposable element landscapes in aging Drosophila.

Genetic mechanisms that repress transposable elements (TEs) in young animals decline during aging, as reflected by increased TE expression in aged animals. Does increased TE expression during aging lead to more genomic TE copies in older animals? To address this question, we quantified TE Landscapes (TLs) via whole genome sequencing of young and aged Drosophila strains of wild-type and mutant backgrounds. We quantified TLs in whole flies and dissected brains and validated the feasibility of our approach in detecting new TE insertions in aging Drosophila genomes when small RNA and RNA interference (RNAi) pathways are compromised. We also describe improved sequencing methods to quantify extra-chromosomal DNA circles (eccDNAs) in Drosophila as an additional source of TE copies that accumulate during aging. Lastly, to combat the natural progression of aging-associated TE expression, we show that knocking down PAF1, a conserved transcription elongation factor that antagonizes RNAi pathways, may bolster suppression of TEs during aging and extend lifespan. Our study suggests that in addition to a possible influence by different genetic backgrounds, small RNA and RNAi mechanisms may mitigate genomic TL expansion despite the increase in TE transcripts during aging.

NURF301 contributes to gypsy chromatin insulator-mediated nuclear organization.

Chromatin insulators are DNA-protein complexes that can prevent the spread of repressive chromatin and block communication between enhancers and promoters to regulate gene expression. In Drosophila, the gypsy chromatin insulator complex consists of three core proteins: CP190, Su(Hw), and Mod(mdg4)67.2. These factors concentrate at nuclear foci termed insulator bodies, and changes in insulator body localization have been observed in mutants defective for insulator function. Here, we identified NURF301/E(bx), a nucleosome remodeling factor, as a novel regulator of gypsy insulator body localization through a high-throughput RNAi imaging screen. NURF301 promotes gypsy-dependent insulator barrier activity and physically interacts with gypsy insulator proteins. Using ChIP-seq, we found that NURF301 co-localizes with insulator proteins genome-wide, and NURF301 promotes chromatin association of Su(Hw) and CP190 at gypsy insulator binding sites. These effects correlate with NURF301-dependent nucleosome repositioning. At the same time, CP190 and Su(Hw) both facilitate recruitment of NURF301 to chromatin. Finally, Oligopaint FISH combined with immunofluorescence revealed that NURF301 promotes 3D contact between insulator bodies and gypsy insulator DNA binding sites, and NURF301 is required for proper nuclear positioning of gypsy binding sites. Our data provide new insights into how a nucleosome remodeling factor and insulator proteins cooperatively contribute to nuclear organization.

Oligopaint DNA FISH reveals telomere-based meiotic pairing dynamics in the silkworm, Bombyx mori.

Accurate chromosome segregation during meiosis is essential for reproductive success. Yet, many fundamental aspects of meiosis remain unclear, including the mechanisms regulating homolog pairing across species. This gap is partially due to our inability to visualize individual chromosomes during meiosis. Here, we employ Oligopaint FISH to investigate homolog pairing and compaction of meiotic chromosomes and resurrect a classical model system, the silkworm Bombyx mori. Our Oligopaint design combines multiplexed barcoding with secondary oligo labeling for high flexibility and low cost. These studies illustrate that Oligopaints are highly specific in whole-mount gonads and on meiotic squashes. We show that meiotic pairing is robust in both males and females and that pairing can occur through numerous partially paired intermediate structures. We also show that pairing in male meiosis occurs asynchronously and seemingly in a transcription-biased manner. Further, we reveal that meiotic bivalent formation in B. mori males is highly similar to bivalent formation in C. elegans, with both of these pathways ultimately resulting in the pairing of chromosome ends with non-paired ends facing the spindle pole. Additionally, microtubule recruitment in both C. elegans and B. mori is likely dependent on kinetochore proteins but independent of the centromere-specifying histone CENP-A. Finally, using super-resolution microscopy in the female germline, we show that homologous chromosomes remain associated at telomere domains in the absence of chiasma and after breakdown and modification to the synaptonemal complex in pachytene. These studies reveal novel insights into mechanisms of meiotic homolog pairing both with or without recombination.

DiffChIPL: a differential peak analysis method for high-throughput sequencing data with biological replicates based on limma.

MOTIVATION: ChIP-seq detects protein-DNA interactions within chromatin, such as that of chromatin structural components and transcription machinery. ChIP-seq profiles are often noisy and variable across replicates, posing a challenge to the development of effective algorithms to accurately detect differential peaks. Methods have recently been designed for this purpose but sometimes yield conflicting results that are inconsistent with the underlying biology. Most existing algorithms perform well on limited datasets. To improve differential analysis of ChIP-seq, we present a novel Differential analysis method for ChIP-seq based on Limma (DiffChIPL). RESULTS: DiffChIPL is adaptive to asymmetrical or symmetrical data and can accurately report global differences. We used simulated and real datasets for transcription factors (TFs) and histone modification marks to validate and benchmark our algorithm. DiffChIPL shows superior performance in sensitivity and false positive rate in different simulations and control datasets. DiffChIPL also performs well on real ChIP-seq, CUT&RUN, CUT&Tag and ATAC-seq datasets. DiffChIPL is an accurate and robust method, exhibiting better performance in differential analysis for a variety of applications including TF binding, histone modifications and chromatin accessibility. AVAILABILITY AND IMPLEMENTATION: https://github.com/yancychy/DiffChIPL. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Widespread rearrangement of 3D chromatin organization underlies polycomb-mediated stress-induced silencing.

Chromosomes of metazoan organisms are partitioned in the interphase nucleus into discrete topologically associating domains (TADs). Borders between TADs are formed in regions containing active genes and clusters of architectural protein binding sites. The transcription of most genes is repressed after temperature stress in Drosophila. Here we show that temperature stress induces relocalization of architectural proteins from TAD borders to inside TADs, and this is accompanied by a dramatic rearrangement in the 3D organization of the nucleus. TAD border strength declines, allowing for an increase in long-distance inter-TAD interactions. Similar but quantitatively weaker effects are observed upon inhibition of transcription or depletion of individual architectural proteins. Heat shock-induced inter-TAD interactions result in increased contacts among enhancers and promoters of silenced genes, which recruit Pc and form Pc bodies in the nucleolus. These results suggest that the TAD organization of metazoan genomes is plastic and can be reconfigured quickly.

The zinc-finger protein CLAMP promotes gypsy chromatin insulator function in Drosophila.

Chromatin insulators are DNA-protein complexes that establish independent higher-order DNA domains to influence transcription. Insulators are functionally defined by two properties: they can block communication between an enhancer and a promoter, and also act as a barrier between heterochromatin and euchromatin. In Drosophila, the gypsy insulator complex contains three core components; Su(Hw), CP190 and Mod(mdg4)67.2. Here, we identify a novel role for Chromatin-linked adaptor for MSL proteins (CLAMP) in promoting gypsy chromatin insulator function. When clamp is knocked down, gypsy-dependent enhancer-blocking and barrier activities are strongly reduced. CLAMP associates physically with the core gypsy insulator complex, and ChIP-seq analysis reveals extensive overlap, particularly with promoter-bound CP190 on chromatin. Depletion of CLAMP disrupts CP190 binding at a minority of shared sites, whereas depletion of CP190 results in extensive loss of CLAMP chromatin association. Finally, reduction of CLAMP disrupts CP190 localization within the nucleus. Our results support a positive functional relationship between CLAMP and CP190 to promote gypsy chromatin insulator activity.

Shep regulates Drosophila neuronal remodeling by controlling transcription of its chromatin targets.

Neuronal remodeling is crucial for formation of the mature nervous system and disruption of this process can lead to neuropsychiatric diseases. Global gene expression changes in neurons during remodeling as well as the factors that regulate these changes remain poorly defined. To elucidate this process, we performed RNA-seq on isolated Drosophila larval and pupal neurons and found upregulated synaptic signaling and downregulated gene expression regulators as a result of normal neuronal metamorphosis. We further tested the role of alan shepard (shep), which encodes an evolutionarily conserved RNA-binding protein required for proper neuronal remodeling. Depletion of shep in neurons prevents the execution of metamorphic gene expression patterns, and shep-regulated genes correspond to Shep chromatin and/or RNA-binding targets. Reduced expression of a Shep-inhibited target gene that we identified, brat, is sufficient to rescue neuronal remodeling defects of shep knockdown flies. Our results reveal direct regulation of transcriptional programs by Shep to regulate neuronal remodeling during metamorphosis.

Argonaute2 and LaminB modulate gene expression by controlling chromatin topology.

Drosophila Argonaute2 (AGO2) has been shown to regulate expression of certain loci in an RNA interference (RNAi)-independent manner, but its genome-wide function on chromatin remains unknown. Here, we identified the nuclear scaffolding protein LaminB as a novel interactor of AGO2. When either AGO2 or LaminB are depleted in Kc cells, similar transcription changes are observed genome-wide. In particular, changes in expression occur mainly in active or potentially active chromatin, both inside and outside LaminB-associated domains (LADs). Furthermore, we identified a somatic target of AGO2 transcriptional repression, no hitter (nht), which is immersed in a LAD located within a repressive topologically-associated domain (TAD). Null mutation but not catalytic inactivation of AGO2 leads to ectopic expression of nht and downstream spermatogenesis genes. Depletion of either AGO2 or LaminB results in reduced looping interactions within the nht TAD as well as ectopic inter-TAD interactions, as detected by 4C-seq analysis. Overall, our findings reveal coordination of AGO2 and LaminB function to dictate genome architecture and thereby regulate gene expression.

A compendium of RNA-binding motifs for decoding gene regulation.

RNA-binding proteins are key regulators of gene expression, yet only a small fraction have been functionally characterized. Here we report a systematic analysis of the RNA motifs recognized by RNA-binding proteins, encompassing 205 distinct genes from 24 diverse eukaryotes. The sequence specificities of RNA-binding proteins display deep evolutionary conservation, and the recognition preferences for a large fraction of metazoan RNA-binding proteins can thus be inferred from their RNA-binding domain sequence. The motifs that we identify in vitro correlate well with in vivo RNA-binding data. Moreover, we can associate them with distinct functional roles in diverse types of post-transcriptional regulation, enabling new insights into the functions of RNA-binding proteins both in normal physiology and in human disease. These data provide an unprecedented overview of RNA-binding proteins and their targets, and constitute an invaluable resource for determining post-transcriptional regulatory mechanisms in eukaryotes.

Different enhancer classes in Drosophila bind distinct architectural proteins and mediate unique chromatin interactions and 3D architecture.

Eukaryotic gene expression is regulated by enhancer-promoter interactions but the molecular mechanisms that govern specificity have remained elusive. Genome-wide studies utilizing STARR-seq identified two enhancer classes in Drosophila that interact with different core promoters: housekeeping enhancers (hkCP) and developmental enhancers (dCP). We hypothesized that the two enhancer classes are occupied by distinct architectural proteins, affecting their enhancer-promoter contacts. By evaluating ChIP-seq occupancy of architectural proteins, typical enhancer-associated proteins, and histone modifications, we determine that both enhancer classes are enriched for RNA Polymerase II, CBP, and architectural proteins but there are also distinctions. hkCP enhancers contain H3K4me3 and exclusively bind Cap-H2, Chromator, DREF and Z4, whereas dCP enhancers contain H3K4me1 and are more enriched for Rad21 and Fs(1)h-L. Additionally, we map the interactions of each enhancer class utilizing a Hi-C dataset with <1 kb resolution. Results suggest that hkCP enhancers are more likely to form multi-TSS interaction networks and be associated with topologically associating domain (TAD) borders, while dCP enhancers are more often bound to one or two TSSs and are enriched at chromatin loop anchors. The data support a model suggesting that the unique architectural protein occupancy within enhancers is one contributor to enhancer-promoter interaction specificity.

RNAi-independent role for Argonaute2 in CTCF/CP190 chromatin insulator function.

A major role of the RNAi pathway in Schizosaccharomyces pombe is to nucleate heterochromatin, but it remains unclear whether this mechanism is conserved. To address this question in Drosophila, we performed genome-wide localization of Argonaute2 (AGO2) by chromatin immunoprecipitation (ChIP)-seq in two different embryonic cell lines and found that AGO2 localizes to euchromatin but not heterochromatin. This localization pattern is further supported by immunofluorescence staining of polytene chromosomes and cell lines, and these studies also indicate that a substantial fraction of AGO2 resides in the nucleus. Intriguingly, AGO2 colocalizes extensively with CTCF/CP190 chromatin insulators but not with genomic regions corresponding to endogenous siRNA production. Moreover, AGO2, but not its catalytic activity or Dicer-2, is required for CTCF/CP190-dependent Fab-8 insulator function. AGO2 interacts physically with CTCF and CP190, and depletion of either CTCF or CP190 results in genome-wide loss of AGO2 chromatin association. Finally, mutation of CTCF, CP190, or AGO2 leads to reduction of chromosomal looping interactions, thereby altering gene expression. We propose that RNAi-independent recruitment of AGO2 to chromatin by insulator proteins promotes the definition of transcriptional domains throughout the genome.

Drosophila Argonaute2 turnover is regulated by the ubiquitin proteasome pathway.

Argonaute (AGO) proteins play a central role in the RNA interference (RNAi) pathway, which is a cytoplasmic mechanism important for post-transcriptional regulation of gene expression. In Drosophila, AGO2 also functions in the nucleus to regulate chromatin insulator activity and transcription. Although there are a number of studies focused on AGO2 function, the regulation of AGO2 turnover is not well understood. We found that mutation of T1149 or R1158 in the conserved PIWI domain causes AGO2 protein instability, but only T1149 affects RNAi activity. Mass spec analysis shows that several proteasome components co-purify with both wildtype and mutant AGO2, and knockdown of two proteasome pathway components results in AGO2 protein accumulation. Finally, AGO2 protein levels increase after treatment with the proteasome inhibitor MG132. Our results indicate that the ubiquitin-proteasome pathway is involved in AGO2 protein turnover.

The RNA-binding protein Rumpelstiltskin antagonizes gypsy chromatin insulator function in a tissue-specific manner.

Chromatin insulators are DNA-protein complexes that are situated throughout the genome that are proposed to contribute to higher-order organization and demarcation into distinct transcriptional domains. Mounting evidence in different species implicates RNA and RNA-binding proteins as regulators of chromatin insulator activities. Here, we identify the Drosophila hnRNP M homolog Rumpelstiltskin (Rump) as an antagonist of gypsy chromatin insulator enhancer-blocking and barrier activities. Despite ubiquitous expression of Rump, decreasing Rump levels leads to improvement of barrier activity only in tissues outside of the central nervous system (CNS). Furthermore, rump mutants restore insulator body localization in an insulator mutant background only in non-CNS tissues. Rump associates physically with core gypsy insulator proteins, and chromatin immunoprecipitation and sequencing analysis of Rump demonstrates extensive colocalization with a subset of insulator sites across the genome. The genome-wide binding profile and tissue specificity of Rump contrast with that of Shep, a recently identified RNA-binding protein that antagonizes gypsy insulator activity primarily in the CNS. Our findings indicate parallel roles for RNA-binding proteins in mediating tissue-specific regulation of chromatin insulator activity.

metaseq: a Python package for integrative genome-wide analysis reveals relationships between chromatin insulators and associated nuclear mRNA.

Here we introduce metaseq, a software library written in Python, which enables loading multiple genomic data formats into standard Python data structures and allows flexible, customized manipulation and visualization of data from high-throughput sequencing studies. We demonstrate its practical use by analyzing multiple datasets related to chromatin insulators, which are DNA-protein complexes proposed to organize the genome into distinct transcriptional domains. Recent studies in Drosophila and mammals have implicated RNA in the regulation of chromatin insulator activities. Moreover, the Drosophila RNA-binding protein Shep has been shown to antagonize gypsy insulator activity in a tissue-specific manner, but the precise role of RNA in this process remains unclear. Better understanding of chromatin insulator regulation requires integration of multiple datasets, including those from chromatin-binding, RNA-binding, and gene expression experiments. We use metaseq to integrate RIP- and ChIP-seq data for Shep and the core gypsy insulator protein Su(Hw) in two different cell types, along with publicly available ChIP-chip and RNA-seq data. Based on the metaseq-enabled analysis presented here, we propose a model where Shep associates with chromatin cotranscriptionally, then is recruited to insulator complexes in trans where it plays a negative role in insulator activity.

Surviving an identity crisis: a revised view of chromatin insulators in the genomics era.

The control of complex, developmentally regulated loci and partitioning of the genome into active and silent domains is in part accomplished through the activity of DNA-protein complexes termed chromatin insulators. Together, the multiple, well-studied classes of insulators in Drosophila melanogaster appear to be generally functionally conserved. In this review, we discuss recent genomic-scale experiments and attempt to reconcile these newer findings in the context of previously defined insulator characteristics based on classical genetic analyses and transgenic approaches. Finally, we discuss the emerging understanding of mechanisms of chromatin insulator regulation. This article is part of a Special Issue entitled: Chromatin and epigenetic regulation of animal development.

Messenger RNA is a functional component of a chromatin insulator complex.

Chromatin insulators are DNA protein complexes situated throughout the genome capable of demarcating independent transcriptional domains. Previous studies point to an important role for RNA in gypsy chromatin insulator function in Drosophila; however, the identity of these putative insulator-associated RNAs is not currently known. Here we utilize RNA-immunoprecipitation and high throughput sequencing (RIP-seq) to isolate RNAs stably associated with gypsy insulator complexes. Strikingly, these RNAs correspond to specific sense-strand, spliced and polyadenylated mRNAs, including two insulator protein transcripts. In order to assess the functional significance of these associated mRNAs independent of their coding function, we expressed untranslatable versions of these transcripts in developing flies and observed both alteration of insulator complex nuclear localization as well as improvement of enhancer-blocking activity. Together, these data suggest a novel, noncoding mechanism by which certain mRNAs contribute to chromatin insulator function.

Genome-wide localization of exosome components to active promoters and chromatin insulators in Drosophila.

Chromatin insulators are functionally conserved DNA-protein complexes situated throughout the genome that organize independent transcriptional domains. Previous work implicated RNA as an important cofactor in chromatin insulator activity, although the precise mechanisms are not yet understood. Here we identify the exosome, the highly conserved major cellular 3' to 5' RNA degradation machinery, as a physical interactor of CP190-dependent chromatin insulator complexes in Drosophila. Genome-wide profiling of exosome by ChIP-seq in two different embryonic cell lines reveals extensive and specific overlap with the CP190, BEAF-32 and CTCF insulator proteins. Colocalization occurs mainly at promoters but also boundary elements such as Mcp, Fab-8, scs and scs', which overlaps with a promoter. Surprisingly, exosome associates primarily with promoters but not gene bodies of active genes, arguing against simple cotranscriptional recruitment to RNA substrates. Similar to insulator proteins, exosome is also significantly enriched at divergently transcribed promoters. Directed ChIP of exosome in cell lines depleted of insulator proteins shows that CTCF is required specifically for exosome association at Mcp and Fab-8 but not other sites, suggesting that alternate mechanisms must also contribute to exosome chromatin recruitment. Taken together, our results reveal a novel positive relationship between exosome and chromatin insulators throughout the genome.