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.

A split personality for nucleosomes.

A high-resolution look at where histones touch DNA reveals a surprisingly intricate, dynamic, and modular nucleosome. Three advances in the study by Rhee et al. include unexpected interactions between the H3 tail and linker DNA, new evidence for existence of subnucleosomal particles, and asymmetric patterns of histone modification within a single nucleosome that correspond to the direction of transcription.

Opportunistic binding of EcR to open chromatin drives tissue-specific developmental responses.

Steroid hormones perform diverse biological functions in developing and adult animals. However, the mechanistic basis for their tissue specificity remains unclear. In Drosophila, the ecdysone steroid hormone is essential for coordinating developmental timing across physically separated tissues. Ecdysone directly impacts genome function through its nuclear receptor, a heterodimer of the EcR and ultraspiracle proteins. Ligand binding to EcR triggers a transcriptional cascade, including activation of a set of primary response transcription factors. The hierarchical organization of this pathway has left the direct role of EcR in mediating ecdysone responses obscured. Here, we investigate the role of EcR in controlling tissue-specific ecdysone responses, focusing on two tissues that diverge in their response to rising ecdysone titers: the larval salivary gland, which undergoes programmed destruction, and the wing imaginal disc, which initiates morphogenesis. We find that EcR functions bimodally, with both gene repressive and activating functions, even at the same developmental stage. EcR DNA binding profiles are highly tissue-specific, and transgenic reporter analyses demonstrate that EcR plays a direct role in controlling enhancer activity. Finally, despite a strong correlation between tissue-specific EcR binding and tissue-specific open chromatin, we find that EcR does not control chromatin accessibility at genomic targets. We conclude that EcR contributes extensively to tissue-specific ecdysone responses. However, control over access to its binding sites is subordinated to other transcription factors.

Hormone-dependent control of developmental timing through regulation of chromatin accessibility.

Specification of tissue identity during development requires precise coordination of gene expression in both space and time. Spatially, master regulatory transcription factors are required to control tissue-specific gene expression programs. However, the mechanisms controlling how tissue-specific gene expression changes over time are less well understood. Here, we show that hormone-induced transcription factors control temporal gene expression by regulating the accessibility of DNA regulatory elements. Using the Drosophila wing, we demonstrate that temporal changes in gene expression are accompanied by genome-wide changes in chromatin accessibility at temporal-specific enhancers. We also uncover a temporal cascade of transcription factors following a pulse of the steroid hormone ecdysone such that different times in wing development can be defined by distinct combinations of hormone-induced transcription factors. Finally, we show that the ecdysone-induced transcription factor E93 controls temporal identity by directly regulating chromatin accessibility across the genome. Notably, we found that E93 controls enhancer activity through three different modalities, including promoting accessibility of late-acting enhancers and decreasing accessibility of early-acting enhancers. Together, this work supports a model in which an extrinsic signal triggers an intrinsic transcription factor cascade that drives development forward in time through regulation of chromatin accessibility.

Direct interrogation of the role of H3K9 in metazoan heterochromatin function.

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

Expression of E93 provides an instructive cue to control dynamic enhancer activity and chromatin accessibility during development.

How temporal cues combine with spatial inputs to control gene expression during development is poorly understood. Here, we test the hypothesis that the Drosophila transcription factor E93 controls temporal gene expression by regulating chromatin accessibility. Precocious expression of E93 early in wing development reveals that it can simultaneously activate and deactivate different target enhancers. Notably, the precocious patterns of enhancer activity resemble the wild-type patterns that occur later in development, suggesting that expression of E93 alters the competence of enhancers to respond to spatial cues. Genomic profiling reveals that precocious E93 expression is sufficient to regulate chromatin accessibility at a subset of its targets. These accessibility changes mimic those that normally occur later in development, indicating that precocious E93 accelerates the wild-type developmental program. Further, we find that target enhancers that do not respond to precocious E93 in early wings become responsive after a developmental transition, suggesting that parallel temporal pathways work alongside E93. These findings support a model wherein E93 expression functions as an instructive cue that defines a broad window of developmental time through control of chromatin accessibility.

Disentangling the developmental origins of a novel phenotype: enhancement versus reversal of environmentally induced gene expression.

Increasing evidence suggests that many novel traits might have originated via plasticity-led evolution (PLE). Yet, little is known of the developmental processes that underpin PLE, especially in its early stages. One such process is 'phenotypic accommodation', which occurs when, in response to a change in the environment, an organism experiences adjustments across variable parts of its phenotype that improve its fitness. Here, we asked if environmentally induced changes in gene expression are enhanced or reversed during phenotypic accommodation of a novel, complex phenotype in spadefoot toad tadpoles (Spea multiplicata). More genes than expected were affected by both the environment and phenotypic accommodation in the liver and brain. However, although phenotypic accommodation primarily reversed environmentally induced changes in gene expression in liver tissue, it enhanced these changes in brain tissue. Thus, depending on the tissue, phenotypic accommodation may either minimize functional disruption via reversal of gene expression patterns or promote novelty via enhancement of existing expression patterns. Our study thereby provides insights into the developmental origins of a novel phenotype and the incipient stages of PLE.

Changes in chromatin accessibility ensure robust cell cycle exit in terminally differentiated cells.

During terminal differentiation, most cells exit the cell cycle and enter into a prolonged or permanent G0 in which they are refractory to mitogenic signals. Entry into G0 is usually initiated through the repression of cell cycle gene expression by formation of a transcriptional repressor complex called dimerization partner (DP), retinoblastoma (RB)-like, E2F and MuvB (DREAM). However, when DREAM repressive function is compromised during terminal differentiation, additional unknown mechanisms act to stably repress cycling and ensure robust cell cycle exit. Here, we provide evidence that developmentally programmed, temporal changes in chromatin accessibility at a small subset of critical cell cycle genes act to enforce cell cycle exit during terminal differentiation in the Drosophila melanogaster wing. We show that during terminal differentiation, chromatin closes at a set of pupal wing enhancers for the key rate-limiting cell cycle regulators Cyclin E (cycE), E2F transcription factor 1 (e2f1), and string (stg). This closing coincides with wing cells entering a robust postmitotic state that is strongly refractory to cell cycle reactivation, and the regions that close contain known binding sites for effectors of mitogenic signaling pathways such as Yorkie and Notch. When cell cycle exit is genetically disrupted, chromatin accessibility at cell cycle genes remains unaffected, and the closing of distal enhancers at cycE, e2f1, and stg proceeds independent of the cell cycling status. Instead, disruption of cell cycle exit leads to changes in accessibility and expression of a subset of hormone-induced transcription factors involved in the progression of terminal differentiation. Our results uncover a mechanism that acts as a cell cycle-independent timer to limit the response to mitogenic signaling and aberrant cycling in terminally differentiating tissues. In addition, we provide a new molecular description of the cross talk between cell cycle exit and terminal differentiation during metamorphosis.

Memes: A motif analysis environment in R using tools from the MEME Suite.

Identification of biopolymer motifs represents a key step in the analysis of biological sequences. The MEME Suite is a widely used toolkit for comprehensive analysis of biopolymer motifs; however, these tools are poorly integrated within popular analysis frameworks like the R/Bioconductor project, creating barriers to their use. Here we present memes, an R package that provides a seamless R interface to a selection of popular MEME Suite tools. memes provides a novel "data aware" interface to these tools, enabling rapid and complex discriminative motif analysis workflows. In addition to interfacing with popular MEME Suite tools, memes leverages existing R/Bioconductor data structures to store the multidimensional data returned by MEME Suite tools for rapid data access and manipulation. Finally, memes provides data visualization capabilities to facilitate communication of results. memes is available as a Bioconductor package at https://bioconductor.org/packages/memes, and the source code can be found at github.com/snystrom/memes.

Direct and widespread role for the nuclear receptor EcR in mediating the response to ecdysone in Drosophila.

The ecdysone pathway was among the first experimental systems employed to study the impact of steroid hormones on the genome. In Drosophila and other insects, ecdysone coordinates developmental transitions, including wholesale transformation of the larva into the adult during metamorphosis. Like other hormones, ecdysone controls gene expression through a nuclear receptor, which functions as a ligand-dependent transcription factor. Although it is clear that ecdysone elicits distinct transcriptional responses within its different target tissues, the role of its receptor, EcR, in regulating target gene expression is incompletely understood. In particular, EcR initiates a cascade of transcription factor expression in response to ecdysone, making it unclear which ecdysone-responsive genes are direct EcR targets. Here, we use the larval-to-prepupal transition of developing wings to examine the role of EcR in gene regulation. Genome-wide DNA binding profiles reveal that EcR exhibits widespread binding across the genome, including at many canonical ecdysone response genes. However, the majority of its binding sites reside at genes with wing-specific functions. We also find that EcR binding is temporally dynamic, with thousands of binding sites changing over time. RNA-seq reveals that EcR acts as both a temporal gate to block precocious entry to the next developmental stage as well as a temporal trigger to promote the subsequent program. Finally, transgenic reporter analysis indicates that EcR regulates not only temporal changes in target enhancer activity but also spatial patterns. Together, these studies define EcR as a multipurpose, direct regulator of gene expression, greatly expanding its role in coordinating developmental transitions.

Lysine 27 of replication-independent histone H3.3 is required for Polycomb target gene silencing but not for gene activation.

Proper determination of cell fates depends on epigenetic information that is used to preserve memory of decisions made earlier in development. Post-translational modification of histone residues is thought to be a central means by which epigenetic information is propagated. In particular, modifications of histone H3 lysine 27 (H3K27) are strongly correlated with both gene activation and gene repression. H3K27 acetylation is found at sites of active transcription, whereas H3K27 methylation is found at loci silenced by Polycomb group proteins. The histones bearing these modifications are encoded by the replication-dependent H3 genes as well as the replication-independent H3.3 genes. Owing to differential rates of nucleosome turnover, H3K27 acetylation is enriched on replication-independent H3.3 histones at active gene loci, and H3K27 methylation is enriched on replication-dependent H3 histones across silenced gene loci. Previously, we found that modification of replication-dependent H3K27 is required for Polycomb target gene silencing, but it is not required for gene activation. However, the contribution of replication-independent H3.3K27 to these functions is unknown. Here, we used CRISPR/Cas9 to mutate the endogenous replication-independent H3.3K27 to a non-modifiable residue. Surprisingly, we find that H3.3K27 is also required for Polycomb target gene silencing despite the association of H3.3 with active transcription. However, the requirement for H3.3K27 comes at a later stage of development than that found for replication-dependent H3K27, suggesting a greater reliance on replication-independent H3.3K27 in post-mitotic cells. Notably, we find no evidence of global transcriptional defects in H3.3K27 mutants, despite the strong correlation between H3.3K27 acetylation and active transcription.

Chromatin conformation and transcriptional activity are permissive regulators of DNA replication initiation in Drosophila.

Chromatin structure has emerged as a key contributor to spatial and temporal control over the initiation of DNA replication. However, despite genome-wide correlations between early replication of gene-rich, accessible euchromatin and late replication of gene-poor, inaccessible heterochromatin, a causal relationship between chromatin structure and replication initiation remains elusive. Here, we combined histone gene engineering and whole-genome sequencing in Drosophila to determine how perturbing chromatin structure affects replication initiation. We found that most pericentric heterochromatin remains late replicating in H3K9R mutants, even though H3K9R pericentric heterochromatin is depleted of HP1a, more accessible, and transcriptionally active. These data indicate that HP1a loss, increased chromatin accessibility, and elevated transcription do not result in early replication of heterochromatin. Nevertheless, a small amount of pericentric heterochromatin with increased accessibility replicates earlier in H3K9R mutants. Transcription is de-repressed in these regions of advanced replication but not in those regions of the H3K9R mutant genome that replicate later, suggesting that transcriptional repression may contribute to late replication. We also explored relationships among chromatin, transcription, and replication in euchromatin by analyzing H4K16R mutants. In Drosophila, the X Chromosome gene expression is up-regulated twofold and replicates earlier in XY males than it does in XX females. We found that H4K16R mutation prevents normal male development and abrogates hyperexpression and earlier replication of the male X, consistent with previously established genome-wide correlations between transcription and early replication. In contrast, H4K16R females are viable and fertile, indicating that H4K16 modification is dispensable for genome replication and gene expression.

Enhancer identification and activity evaluation in the red flour beetle, Tribolium castaneum.

Evolution of cis-regulatory elements (such as enhancers) plays an important role in the production of diverse morphology. However, a mechanistic understanding is often limited by the absence of methods for studying enhancers in species other than established model systems. Here, we sought to establish methods to identify and test enhancer activity in the red flour beetle, Tribolium castaneum To identify possible enhancer regions, we first obtained genome-wide chromatin profiles from various tissues and stages of Tribolium using FAIRE (formaldehyde-assisted isolation of regulatory elements)-sequencing. Comparison of these profiles revealed a distinct set of open chromatin regions in each tissue and at each stage. In addition, comparison of the FAIRE data with sets of computationally predicted (i.e. supervised cis-regulatory module-predicted) enhancers revealed a very high overlap between the two datasets. Second, using nubbin in the wing and hunchback in the embryo as case studies, we established the first universal reporter assay system that works in various contexts in Tribolium, and in a cross-species context. Together, these advances will facilitate investigation of cis-evolution and morphological diversity in Tribolium and other insects.

Transcriptomic bases of a polyphenism.

Polyphenism-in which multiple distinct phenotypes are produced from a single genotype owing to differing environmental conditions-is commonplace, but its molecular bases are poorly understood. Here, we examine the transcriptomic bases of a polyphenism in Mexican spadefoot toads (Spea multiplicata). Depending on their environment, their tadpoles develop into either a default "omnivore" morph or a novel "carnivore" morph. We compared patterns of gene expression among sibships that exhibited high versus low production of carnivores when reared in conditions that induce the carnivore morph versus those that do not. We found that production of the novel carnivore morph actually involved changes in fewer genes than did the maintenance of the default omnivore morph in the inducing environment. However, only body samples showed this pattern; head samples showed the opposite pattern. We also found that changes to lipid metabolism (especially cholesterol biosynthesis) and peroxisome contents and function might be crucial for establishing and maintaining differences between the morphs. Thus, our findings suggest that carnivore phenotype might have originally evolved following the breakdown of robustness mechanisms that maintain the default omnivore phenotype, and that the carnivore morph is developmentally regulated by lipid metabolism and peroxisomal form, function, and/or signaling. This study also serves as a springboard for further exploration into the nature and causes of plasticity in an emerging model system.

Supramolecular assembly of the beta-catenin destruction complex and the effect of Wnt signaling on its localization, molecular size, and activity in vivo.

Wnt signaling provides a paradigm for cell-cell signals that regulate embryonic development and stem cell homeostasis and are inappropriately activated in cancers. The tumor suppressors APC and Axin form the core of the multiprotein destruction complex, which targets the Wnt-effector beta-catenin for phosphorylation, ubiquitination and destruction. Based on earlier work, we hypothesize that the destruction complex is a supramolecular entity that self-assembles by Axin and APC polymerization, and that regulating assembly and stability of the destruction complex underlie its function. We tested this hypothesis in Drosophila embryos, a premier model of Wnt signaling. Combining biochemistry, genetic tools to manipulate Axin and APC2 levels, advanced imaging and molecule counting, we defined destruction complex assembly, stoichiometry, and localization in vivo, and its downregulation in response to Wnt signaling. Our findings challenge and revise current models of destruction complex function. Endogenous Axin and APC2 proteins and their antagonist Dishevelled accumulate at roughly similar levels, suggesting competition for binding may be critical. By expressing Axin:GFP at near endogenous levels we found that in the absence of Wnt signals, Axin and APC2 co-assemble into large cytoplasmic complexes containing tens to hundreds of Axin proteins. Wnt signals trigger recruitment of these to the membrane, while cytoplasmic Axin levels increase, suggesting altered assembly/disassembly. Glycogen synthase kinase3 regulates destruction complex recruitment to the membrane and release of Armadillo/beta-catenin from the destruction complex. Manipulating Axin or APC2 levels had no effect on destruction complex activity when Wnt signals were absent, but, surprisingly, had opposite effects on the destruction complex when Wnt signals were present. Elevating Axin made the complex more resistant to inactivation, while elevating APC2 levels enhanced inactivation. Our data suggest both absolute levels and the ratio of these two core components affect destruction complex function, supporting models in which competition among Axin partners determines destruction complex activity.

An online repository of solvation thermodynamic and structural maps of SARS-CoV-2 targets.

SARS-CoV-2 recently jumped species and rapidly spread via human-to-human transmission to cause a global outbreak of COVID-19. The lack of effective vaccine combined with the severity of the disease necessitates attempts to develop small molecule drugs to combat the virus. COVID19_GIST_HSA is a freely available online repository to provide solvation thermodynamic maps of COVID-19-related protein small molecule drug targets. Grid inhomogeneous solvation theory maps were generated using AmberTools cpptraj-GIST, 3D reference interaction site model maps were created with AmberTools rism3d.snglpnt and hydration site analysis maps were created using SSTMap code. The resultant data can be applied to drug design efforts: scoring solvent displacement for docking, rational lead modification, prioritization of ligand- and protein- based pharmacophore elements, and creation of water-based pharmacophores. Herein, we demonstrate the use of the solvation thermodynamic mapping data. It is hoped that this freely provided data will aid in small molecule drug discovery efforts to defeat SARS-CoV-2.

Chromatin profiling of Drosophila CNS subpopulations identifies active transcriptional enhancers.

One of the key issues in studying transcriptional regulation during development is how to employ genome-wide assays that reveals sites of open chromatin and transcription factor binding to efficiently identify biologically relevant genes and enhancers. Analysis of Drosophila CNS midline cell development provides a useful system for studying transcriptional regulation at the genomic level due to a large, well-characterized set of midline-expressed genes and in vivo validated enhancers. In this study, FAIRE-seq on FACS-purified midline cells was performed and the midline FAIRE data were compared with whole-embryo FAIRE data. We find that regions of the genome with a strong midline FAIRE peak and weak whole-embryo FAIRE peak overlap with known midline enhancers and provide a useful predictive tool for enhancer identification. In a complementary analysis, we compared a large dataset of fragments that drive midline expression in vivo with the FAIRE data. Midline enhancer fragments with a midline FAIRE peak tend to be near midline-expressed genes, whereas midline enhancers without a midline FAIRE peak were often distant from midline-expressed genes and unlikely to drive midline transcription in vivo.

Transcription start site profiling uncovers divergent transcription and enhancer-associated RNAs in Drosophila melanogaster.

BACKGROUND: High-resolution transcription start site (TSS) mapping in D. melanogaster embryos and cell lines has revealed a rich and detailed landscape of both cis- and trans-regulatory elements and factors. However, TSS profiling has not been investigated in an orthogonal in vivo setting. Here, we present a comprehensive dataset that links TSS dynamics with nucleosome occupancy and gene expression in the wandering third instar larva, a developmental stage characterized by large-scale shifts in transcriptional programs in preparation for metamorphosis. RESULTS: The data recapitulate major regulatory classes of TSSs, based on peak width, promoter-proximal polymerase pausing, and cis-regulatory element density. We confirm the paucity of divergent transcription units in D. melanogaster, but also identify notable exceptions. Furthermore, we identify thousands of novel initiation events occurring at unannotated TSSs that can be classified into functional categories by their local density of histone modifications. Interestingly, a sub-class of these unannotated TSSs overlaps with functionally validated enhancer elements, consistent with a regulatory role for "enhancer RNAs" (eRNAs) in defining developmental transcription programs. CONCLUSIONS: High-depth TSS mapping is a powerful strategy for identifying and characterizing low-abundance and/or low-stability RNAs. Global analysis of transcription initiation patterns in a developing organism reveals a vast number of novel initiation events that identify potential eRNAs as well as other non-coding transcripts critical for animal development.

Zelda is differentially required for chromatin accessibility, transcription factor binding, and gene expression in the early Drosophila embryo.

The transition from a specified germ cell to a population of pluripotent cells occurs rapidly following fertilization. During this developmental transition, the zygotic genome is largely transcriptionally quiescent and undergoes significant chromatin remodeling. In Drosophila, the DNA-binding protein Zelda (also known as Vielfaltig) is required for this transition and for transcriptional activation of the zygotic genome. Open chromatin is associated with Zelda-bound loci, as well as more generally with regions of active transcription. Nonetheless, the extent to which Zelda influences chromatin accessibility across the genome is largely unknown. Here we used formaldehyde-assisted isolation of regulatory elements to determine the role of Zelda in regulating regions of open chromatin in the early embryo. We demonstrate that Zelda is essential for hundreds of regions of open chromatin. This Zelda-mediated chromatin accessibility facilitates transcription-factor recruitment and early gene expression. Thus, Zelda possesses some key characteristics of a pioneer factor. Unexpectedly, chromatin at a large subset of Zelda-bound regions remains open even in the absence of Zelda. The GAGA factor-binding motif and embryonic GAGA factor binding are specifically enriched in these regions. We propose that both Zelda and GAGA factor function to specify sites of open chromatin and together facilitate the remodeling of the early embryonic genome.