Genetics is a major determinant of susceptibility to autoimmune disorders. Here, we examined whether genome organization provides resilience or susceptibility to sequence variations, and how this would contribute to the molecular etiology of an autoimmune disease. We generated high-resolution maps of linear and 3D genome organization in thymocytes of NOD mice, a model of type 1 diabetes (T1D), and the diabetes-resistant C57BL/6 mice. Multi-enhancer interactions formed at genomic regions harboring genes with prominent roles in T cell development in both strains. However, diabetes risk-conferring loci coalesced enhancers and promoters in NOD, but not C57BL/6 thymocytes. 3D genome mapping of NODxC57BL/6 F1 thymocytes revealed that genomic misfolding in NOD mice is mediated in cis. Moreover, immune cells infiltrating the pancreas of humans with T1D exhibited increased expression of genes located on misfolded loci in mice. Thus, genetic variation leads to altered 3D chromatin architecture and associated changes in gene expression that may underlie autoimmune pathology.

A defining feature of early embryogenesis is the transition from maternal to zygotic control. This transition requires embryo-wide zygotic genome activation (ZGA), but the extent of spatiotemporal coordination of ZGA between individual cells is unknown. Multiple interrelated parameters, including elapsed time, completed cycles of cell division, and cell size may impact ZGA onset; however, the principal determinant of ZGA during vertebrate embryogenesis is debated. Here, we perform single-cell imaging of large-scale ZGA in whole-mount Xenopus embryos. We find a striking new spatiotemporal pattern of ZGA whose onset tightly correlates with cell size but not with elapsed time or number of cell divisions. Further, reducing cell size induces premature ZGA, dose dependently. We conclude that large-scale ZGA is not spatially uniform and that its onset is determined at the single-cell level, primarily by cell size. Our study suggests that spatial patterns of ZGA onset may be a common feature of embryonic systems.

Cellular ability to mount an enhanced transcriptional response upon repeated exposure to external cues is termed transcriptional memory, which can be maintained epigenetically through cell divisions and can depend on a nuclear pore component Nup98. The majority of mechanistic knowledge on transcriptional memory has been derived from bulk molecular assays. To gain additional perspective on the mechanism and contribution of Nup98 to memory, we used single-molecule RNA FISH (smFISH) to examine the dynamics of transcription in Drosophila cells upon repeated exposure to the steroid hormone ecdysone. We combined smFISH with mathematical modeling and found that upon hormone exposure, cells rapidly activate a low-level transcriptional response, but simultaneously begin a slow transition into a specialized memory state characterized by a high rate of expression. Strikingly, our modeling predicted that this transition between non-memory and memory states is independent of the transcription stemming from initial activation. We confirmed this prediction experimentally by showing that inhibiting transcription during initial ecdysone exposure did not interfere with memory establishment. Together, our findings reveal that Nup98's role in transcriptional memory is to stabilize the forward rate of conversion from low to high expressing state, and that induced genes engage in two separate behaviors - transcription itself and the establishment of epigenetically propagated transcriptional memory.

Heterodimeric TGF-β ligands outperform homodimers in a variety of developmental, cell culture, and therapeutic contexts; however, the mechanisms underlying this increased potency remain uncharacterized. Here, we use dorsal-ventral axial patterning of the zebrafish embryo to interrogate the BMP2/7 heterodimer signaling mechanism. We demonstrate that differential interactions with BMP antagonists do not account for the reduced signaling ability of homodimers. Instead, we find that while overexpressed BMP2 homodimers can signal, they require two nonredundant type I receptors, one from the Acvr1 subfamily and one from the Bmpr1 subfamily. This implies that all BMP signaling within the zebrafish gastrula, even BMP2 homodimer signaling, requires Acvr1. This is particularly surprising as BMP2 homodimers do not bind Acvr1 in vitro. Furthermore, we find that the roles of the two type I receptors are subfunctionalized within the heterodimer signaling complex, with the kinase activity of Acvr1 being essential, while that of Bmpr1 is not. These results suggest that the potency of the Bmp2/7 heterodimer arises from the ability to recruit both Acvr1 and Bmpr1 into the same signaling complex.

Nuclear speckles are prominent nuclear bodies that contain proteins and RNA involved in gene expression. Although links between nuclear speckles and gene activation are emerging, the mechanisms regulating association of genes with speckles are unclear. We find that speckle association of p53 target genes is driven by the p53 transcription factor. Focusing on p21, a key p53 target, we demonstrate that speckle association boosts expression by elevating nascent RNA amounts. p53-regulated speckle association did not depend on p53 transactivation functions but required an intact proline-rich domain and direct DNA binding, providing mechanisms within p53 for regulating gene-speckle association. Beyond p21, a substantial subset of p53 targets have p53-regulated speckle association. Strikingly, speckle-associating p53 targets are more robustly activated and occupy a distinct niche of p53 biology compared with non-speckle-associating p53 targets. Together, our findings illuminate regulated speckle association as a mechanism used by a transcription factor to boost gene expression.

The Bicoid morphogen gradient directs the patterning of cell fates along the anterior-posterior axis of the syncytial Drosophila embryo and serves as a paradigm of morphogen-mediated patterning. The simplest models of gradient formation rely on constant protein synthesis and diffusion from anteriorly localized source mRNA, coupled with uniform protein degradation. However, currently such models cannot account for all known gradient characteristics. Recent work has proposed that bicoid mRNA spatial distribution is sufficient to produce the observed protein gradient, minimizing the role of protein transport. Here, we adapt a novel method of fluorescent in situ hybridization to quantify the global spatio-temporal dynamics of bicoid mRNA particles. We determine that >90% of all bicoid mRNA is continuously present within the anterior 20% of the embryo. bicoid mRNA distribution along the body axis remains nearly unchanged despite dynamic mRNA translocation from the embryo core to the cortex. To evaluate the impact of mRNA distribution on protein gradient dynamics, we provide detailed quantitative measurements of nuclear Bicoid levels during the formation of the protein gradient. We find that gradient establishment begins 45 minutes after fertilization and that the gradient requires about 50 minutes to reach peak levels. In numerical simulations of gradient formation, we find that incorporating the actual bicoid mRNA distribution yields a closer prediction of the observed protein dynamics compared to modeling protein production from a point source at the anterior pole. We conclude that the spatial distribution of bicoid mRNA contributes to, but cannot account for, protein gradient formation, and therefore that protein movement, either active or passive, is required for gradient formation.

Patterning the embryonic dorsoventral axis of both vertebrates and invertebrates requires signalling through bone morphogenetic proteins (BMPs). Although a well-studied process, the identity of the physiologically relevant BMP signalling complex in the Drosophila melanogaster embryo is controversial, is generally inferred from cell culture studies and has not been investigated in vertebrates. Here, we demonstrate that dorsoventral patterning in zebrafish, Danio rerio, requires two classes of non-redundant type I BMP receptors, Alk3/6 and Alk8 (activin-like kinases 3/6 and 8). We show, under physiological conditions in the embryo, that these two type I receptor classes form a complex in a manner that depends on Bmp2 and Bmp7. We found that both Bmp2-7 heterodimers, as well as Bmp2 and Bmp7 homodimers, form in the embryo. However, only recombinant ligand heterodimers can activate BMP signalling in the early embryo, whereas a combination of Bmp2 and Bmp7 homodimers cannot. We propose that only heterodimers, signalling through two distinct classes of type I receptor, possess sufficient receptor affinity in an environment of extracellular antagonists to elicit the signalling response required for dorsoventral patterning.

Dosage compensation in Drosophila melanogaster involves a 2-fold transcriptional upregulation of the male X chromosome, which relies on the X-chromosome-binding males-specific lethal (MSL) complex. However, how such 2-fold precision is accomplished remains unclear. Here, we show that a nuclear pore component, Mtor, is involved in setting the correct levels of transcription from the male X chromosome. Using larval tissues, we demonstrate that the depletion of Mtor results in selective upregulation at MSL targets of the male X, beyond the required 2-fold. Mtor and MSL components interact genetically, and depletion of Mtor can rescue the male lethality phenotype of MSL components. Using RNA fluorescence in situ hybridization (FISH) analysis and nascent transcript sequencing, we find that the effect of Mtor is not due to defects in mRNA export but occurs at the level of nascent transcription. These findings demonstrate a physiological role for Mtor in the process of dosage compensation, as a transcriptional attenuator of X chromosome gene expression.

Nuclear pore complexes have emerged in recent years as chromatin-binding nuclear scaffolds, able to influence target gene expression. However, how nucleoporins (Nups) exert this control remains poorly understood. Here we show that ectopically tethering Drosophila Nups, especially Sec13, to chromatin is sufficient to induce chromatin decondensation. This decondensation is mediated through chromatin-remodeling complex PBAP, as PBAP is both robustly recruited by Sec13 and required for Sec13-induced decondensation. This phenomenon is not correlated with localization of the target locus to the nuclear periphery, but is correlated with robust recruitment of Nup Elys. Furthermore, we identified a biochemical interaction between endogenous Sec13 and Elys with PBAP, and a role for endogenous Elys in global as well as gene-specific chromatin decompaction. Together, these findings reveal a functional role and mechanism for specific nuclear pore components in promoting an open chromatin state.

Chromatin loops enable transcription-factor-bound distal enhancers to interact with their target promoters to regulate transcriptional programs. Although developmental transcription factors such as active forms of Notch can directly stimulate transcription by activating enhancers, the effect of their oncogenic subversion on the 3D organization of cancer genomes is largely undetermined. By mapping chromatin looping genome-wide in Notch-dependent triple-negative breast cancer and B cell lymphoma, we show that beyond the well-characterized role of Notch as an activator of distal enhancers, Notch regulates its direct target genes by instructing enhancer repositioning. Moreover, a large fraction of Notch-instructed regulatory loops form highly interacting enhancer and promoter spatial clusters termed "3D cliques." Loss- and gain-of-function experiments show that Notch preferentially targets hyperconnected 3D cliques that regulate the expression of crucial proto-oncogenes. Our observations suggest that oncogenic hijacking of developmental transcription factors can dysregulate transcription through widespread effects on the spatial organization of cancer genomes.