Prdm14 is a sequence-specific transcriptional regulator of embryonic stem cell (ESC) pluripotency and primordial germ cell (PGC) formation. It exerts its function, at least in part, through repressing genes associated with epigenetic modification and cell differentiation. Here, we show that this repressive function is mediated through an ETO-family co-repressor Mtgr1, which tightly binds to the pre-SET/SET domains of Prdm14 and co-occupies its genomic targets in mouse ESCs. We generated two monobodies, synthetic binding proteins, targeting the Prdm14 SET domain and demonstrate their utility, respectively, in facilitating crystallization and structure determination of the Prdm14-Mtgr1 complex, or as genetically encoded inhibitor of the Prdm14-Mtgr1 interaction. Structure-guided point mutants and the monobody abrogated the Prdm14-Mtgr1 association and disrupted Prdm14's function in mESC gene expression and PGC formation in vitro. Altogether, our work uncovers the molecular mechanism underlying Prdm14-mediated repression and provides renewable reagents for studying and controlling Prdm14 functions.

Eukaryotic cells possess a variety of signaling pathways that prevent accumulation of unfolded and misfolded proteins. Chief among these is the heat shock response (HSR), which is assumed to respond to unfolded proteins in the cytosol and nucleus alike. In this study, we probe this axiom further using engineered proteins called 'destabilizing domains', whose folding state we control with a small molecule. The sudden appearance of unfolded protein in mammalian cells elicits a robust transcriptional response, which is distinct from the HSR and other known pathways that respond to unfolded proteins. The cellular response to unfolded protein is strikingly different in the nucleus and the cytosol, although unfolded protein in either compartment engages the p53 network. This response provides cross-protection during subsequent proteotoxic stress, suggesting that it is a central component of protein quality control networks, and like the HSR, is likely to influence the initiation and progression of human pathologies.

Conflicts between transcription and replication are a potent source of DNA damage. Co-transcriptional R-loops could aggravate such conflicts by creating an additional barrier to replication fork progression. Here, we use a defined episomal system to investigate how conflict orientation and R-loop formation influence genome stability in human cells. R-loops, but not normal transcription complexes, induce DNA breaks and orientation-specific DNA damage responses during conflicts with replication forks. Unexpectedly, the replisome acts as an orientation-dependent regulator of R-loop levels, reducing R-loops in the co-directional (CD) orientation but promoting their formation in the head-on (HO) orientation. Replication stress and deregulated origin firing increase the number of HO collisions leading to genome-destabilizing R-loops. Our findings connect DNA replication to R-loop homeostasis and suggest a mechanistic basis for genome instability resulting from deregulated DNA replication, observed in cancer and other disease states.

Histone H3 lysine 4 (K4) methylation has been linked to the transcriptional activation in a variety of eukaryotic species. Here we show that a common component of MLL1, MLL2, and hSet1 H3 K4 methyltransferase complexes, the WD40-repeat protein WDR5, directly associates with histone H3 di- and trimethylated at K4 and with H3-K4-dimethylated nucleosomes. WDR5 is required for binding of the methyltransferase complex to the K4-dimethylated H3 tail as well as for global H3 K4 trimethylation and HOX gene activation in human cells. WDR5 is essential for vertebrate development, in that WDR5-depleted X. laevis tadpoles exhibit a variety of developmental defects and abnormal spatial Hox gene expression. Our results are the first demonstration that a WD40-repeat protein acts as a module for recognition of a specific histone modification and suggest a mechanism for reading and writing an epigenetic mark for gene activation.

Transposable elements, also known as transposons, are now recognized not only as parasitic DNA, the spread of which in the genome must be controlled by the host, but also as major players in genome evolution and regulation. Long interspersed element-1 (LINE-1, also known as L1), the only currently autonomous mobile transposon in humans, occupies 17% of the genome and generates inter- and intra-individual genetic variation, in some cases resulting in disease. However, how L1 activity is controlled and the function of L1s in host gene regulation are not completely understood. Here we use CRISPR-Cas9 screening strategies in two distinct human cell lines to provide a genome-wide survey of genes involved in the control of L1 retrotransposition. We identify functionally diverse genes that either promote or restrict L1 retrotransposition. These genes, which are often associated with human diseases, control the L1 life cycle at the transcriptional or the post-transcriptional level in a manner that can depend on the endogenous L1 nucleotide sequence, underscoring the complexity of L1 regulation. We further investigate the restriction of L1 by the protein MORC2 and by the human silencing hub (HUSH) complex subunits MPP8 and TASOR. HUSH and MORC2 can selectively bind evolutionarily young, full-length L1s located within transcriptionally permissive euchromatic environments, and promote deposition of histone H3 Lys9 trimethylation (H3K9me3) for transcriptional silencing. Notably, these silencing events often occur within introns of transcriptionally active genes, and lead to the downregulation of host gene expression in a HUSH-, MORC2-, and L1-dependent manner. Together, these results provide a rich resource for studies of L1 retrotransposition, elucidate a novel L1 restriction pathway and illustrate how epigenetic silencing of transposable elements rewires host gene expression programs.

The infectious cycle of primate lentiviruses is intimately linked to interactions between cells of the immune system. Nef, a potent virulence factor, alters cellular environments to increase lentiviral replication in the host, yet the mechanisms underlying these effects have remained elusive. Since Nef likely functions as an adaptor protein, we exploited a proteomic approach to directly identify molecules that Nef targets to subvert the signaling machinery in T cells. We purified to near homogeneity a major Nef-associated protein complex from T cells and identified by mass spectroscopy its subunits as DOCK2-ELMO1, a key activator of Rac in antigen- and chemokine-initiated signaling pathways, and Rac. We show that Nef activates Rac in T cell lines and in primary T cells following infection with HIV-1 in the absence of antigenic stimuli. Nef activates Rac by binding the DOCK2-ELMO1 complex, and this interaction is linked to the abilities of Nef to inhibit chemotaxis and promote T cell activation. Our data indicate that Nef targets a critical switch that regulates Rac GTPases downstream of chemokine- and antigen-initiated signaling pathways. This interaction enables Nef to influence multiple aspects of T cell function and thus provides an important mechanism by which Nef impacts pathogenesis by primate lentiviruses.

The human face is complex and multipartite, and characterization of its genetic architecture remains challenging. Using a multivariate genome-wide association study meta-analysis of 8,246 European individuals, we identified 203 genome-wide-significant signals (120 also study-wide significant) associated with normal-range facial variation. Follow-up analyses indicate that the regions surrounding these signals are enriched for enhancer activity in cranial neural crest cells and craniofacial tissues, several regions harbor multiple signals with associations to different facial phenotypes, and there is evidence for potential coordinated actions of variants. In summary, our analyses provide insights into the understanding of how complex morphological traits are shaped by both individual and coordinated genetic actions.

Facial morphology is highly variable, both within and among human populations, and a sizable portion of this variation is attributable to genetics. Previous genome scans have revealed more than 100 genetic loci associated with different aspects of normal-range facial variation. Most of these loci have been detected in Europeans, with few studies focusing on other ancestral groups. Consequently, the degree to which facial traits share a common genetic basis across diverse sets of humans remains largely unknown. We therefore investigated the genetic basis of facial morphology in an East African cohort. We applied an open-ended data-driven phenotyping approach to a sample of 2,595 3D facial images collected on Tanzanian children. This approach segments the face into hierarchically arranged, multivariate features that capture the shape variation after adjusting for age, sex, height, weight, facial size and population stratification. Genome scans of these multivariate shape phenotypes revealed significant (p < 2.5 × 10-8) signals at 20 loci, which were enriched for active chromatin elements in human cranial neural crest cells and embryonic craniofacial tissue, consistent with an early developmental origin of the facial variation. Two of these associations were in highly conserved regions showing craniofacial-specific enhancer activity during embryological development (5q31.1 and 12q21.31). Six of the 20 loci surpassed a stricter threshold accounting for multiple phenotypes with study-wide significance (p < 6.25 × 10-10). Cross-population comparisons indicated 10 association signals were shared with Europeans (seven sharing the same associated SNP), and facilitated fine-mapping of causal variants at previously reported loci. Taken together, these results may point to both shared and population-specific components to the genetic architecture of facial variation.

Transcriptional regulation exhibits extensive robustness, but human genetics indicates sensitivity to transcription factor (TF) dosage. Reconciling such observations requires quantitative studies of TF dosage effects at trait-relevant ranges, largely lacking so far. TFs play central roles in both normal-range and disease-associated variation in craniofacial morphology; we therefore developed an approach to precisely modulate TF levels in human facial progenitor cells and applied it to SOX9, a TF associated with craniofacial variation and disease (Pierre Robin sequence (PRS)). Most SOX9-dependent regulatory elements (REs) are buffered against small decreases in SOX9 dosage, but REs directly and primarily regulated by SOX9 show heightened sensitivity to SOX9 dosage; these RE responses partially predict gene expression responses. Sensitive REs and genes preferentially affect functional chondrogenesis and PRS-like craniofacial shape variation. We propose that such REs and genes underlie the sensitivity of specific phenotypes to TF dosage, while buffering of other genes leads to robust, nonlinear dosage-to-phenotype relationships.

Recent work suggests extensive adaptation of transposable elements (TEs) for host gene regulation. However, high numbers of integrations typical of TEs, coupled with sequence divergence within families, have made systematic interrogation of the regulatory contributions of TEs challenging. Here, we employ CARGO, our recent method for CRISPR gRNA multiplexing, to facilitate targeting of LTR5HS, an ape-specific class of HERVK (HML-2) LTRs that is active during early development and present in ~700 copies throughout the human genome. We combine CARGO with CRISPR activation or interference to, respectively, induce or silence LTR5HS en masse, and demonstrate that this system robustly targets the vast majority of LTR5HS insertions. Remarkably, activation/silencing of LTR5HS is associated with reciprocal up- and down-regulation of hundreds of human genes. These effects require the presence of retroviral sequences, but occur over long genomic distances, consistent with a pervasive function of LTR5HS elements as early embryonic enhancers in apes.

Monomethylation of histone H3 at lysine 4 (H3K4me1) and acetylation of histone H3 at lysine 27 (H3K27ac) are correlated with transcriptionally engaged enhancer elements, but the functional impact of these modifications on enhancer activity is not well understood. Here we used CRISPR/Cas9 genome editing to separate catalytic activity-dependent and independent functions of Mll3 (Kmt2c) and Mll4 (Kmt2d, Mll2), the major enhancer H3K4 monomethyltransferases. Loss of H3K4me1 from enhancers in Mll3/4 catalytically deficient cells causes partial reduction of H3K27ac, but has surprisingly minor effects on transcription from either enhancers or promoters. In contrast, loss of Mll3/4 proteins leads to strong depletion of enhancer Pol II occupancy and eRNA synthesis, concomitant with downregulation of target genes. Interestingly, downregulated genes exhibit reduced polymerase levels in gene bodies, but not at promoters, suggestive of pause-release defects. Altogether, our results suggest that enhancer H3K4me1 provides only a minor contribution to the long-range coactivator function of Mll3/4.

The developmental disorder Floating-Harbor syndrome (FHS) is caused by heterozygous truncating mutations in SRCAP, a gene encoding a chromatin remodeler mediating incorporation of histone variant H2A.Z. Here, we demonstrate that FHS-associated mutations result in loss of SRCAP nuclear localization, alter neural crest gene programs in human in vitro models and Xenopus embryos, and cause craniofacial defects. These defects are mediated by one of two H2A.Z subtypes, H2A.Z.2, whose knockdown mimics and whose overexpression rescues the FHS phenotype. Selective rescue by H2A.Z.2 is conferred by one of the three amino acid differences between the H2A.Z subtypes, S38/T38. We further show that H2A.Z.1 and H2A.Z.2 genomic occupancy patterns are qualitatively similar, but quantitatively distinct, and H2A.Z.2 incorporation at AT-rich enhancers and expression of their associated genes are both sensitized to SRCAP truncations. Altogether, our results illuminate the mechanism underlying a human syndrome and uncover selective functions of H2A.Z subtypes during development.

Zinc finger protein Zscan4 is selectively expressed in mouse two-cell (2C) embryos undergoing zygotic genome activation (ZGA) and in a rare subpopulation of embryonic stem cells with 2C-like features. Here, we show that Zscan4 specifically recognizes a subset of (CA)n microsatellites, repeat sequences prone to genomic instability. Zscan4-associated microsatellite regions are characterized by low nuclease sensitivity and high histone occupancy. In vitro, Zscan4 binds nucleosomes and protects them from disassembly upon torsional strain. Furthermore, Zscan4 depletion leads to elevated DNA damage in 2C mouse embryos in a transcription-dependent manner. Together, our results identify Zscan4 as a DNA sequence-dependent microsatellite binding factor and suggest a developmentally regulated mechanism, which protects fragile genomic regions from DNA damage at a time of embryogenesis associated with high transcriptional burden and genomic stress.

Non-coding mutations at the far end of a large gene desert surrounding the SOX9 gene result in a human craniofacial disorder called Pierre Robin sequence (PRS). Leveraging a human stem cell differentiation model, we identify two clusters of enhancers within the PRS-associated region that regulate SOX9 expression during a restricted window of facial progenitor development at distances up to 1.45 Mb. Enhancers within the 1.45 Mb cluster exhibit highly synergistic activity that is dependent on the Coordinator motif. Using mouse models, we demonstrate that PRS phenotypic specificity arises from the convergence of two mechanisms: confinement of Sox9 dosage perturbation to developing facial structures through context-specific enhancer activity and heightened sensitivity of the lower jaw to Sox9 expression reduction. Overall, we characterize the longest-range human enhancers involved in congenital malformations, directly demonstrate that PRS is an enhanceropathy, and illustrate how small changes in gene expression can lead to morphological variation.

Diffuse intrinsic pontine glioma (DIPG) is a fatal pediatric cancer with limited therapeutic options. The majority of cases of DIPG exhibit a mutation in histone-3 (H3K27M) that results in oncogenic transcriptional aberrancies. We show here that DIPG is vulnerable to transcriptional disruption using bromodomain inhibition or CDK7 blockade. Targeting oncogenic transcription through either of these methods synergizes with HDAC inhibition, and DIPG cells resistant to HDAC inhibitor therapy retain sensitivity to CDK7 blockade. Identification of super-enhancers in DIPG provides insights toward the cell of origin, highlighting oligodendroglial lineage genes, and reveals unexpected mechanisms mediating tumor viability and invasion, including potassium channel function and EPH receptor signaling. The findings presented demonstrate transcriptional vulnerabilities and elucidate previously unknown mechanisms of DIPG pathobiology.

Enhancer clusters overlapping disease-associated mutations in Pierre Robin sequence (PRS) patients regulate SOX9 expression at genomic distances over 1.25 Mb. We applied optical reconstruction of chromatin architecture (ORCA) imaging to trace 3D locus topology during PRS-enhancer activation. We observed pronounced changes in locus topology between cell types. Subsequent analysis of single-chromatin fiber traces revealed that these ensemble-average differences arise through changes in the frequency of commonly sampled topologies. We further identified two CTCF-bound elements, internal to the SOX9 topologically associating domain, which promote stripe formation, are positioned near the domain's 3D geometric center, and bridge enhancer-promoter contacts in a series of chromatin loops. Ablation of these elements results in diminished SOX9 expression and altered domain-wide contacts. Polymer models with uniform loading across the domain and frequent cohesin collisions recapitulate this multi-loop, centrally clustered geometry. Together, we provide mechanistic insights into architectural stripe formation and gene regulation over ultra-long genomic ranges.

R-loops are dynamic, co-transcriptional nucleic acid structures that facilitate physiological processes but can also cause DNA damage in certain contexts. Perturbations of transcription or R-loop resolution are expected to change their genomic distribution. Next-generation sequencing approaches to map RNA-DNA hybrids, a component of R-loops, have so far not allowed quantitative comparisons between such conditions. Here, we describe quantitative differential DNA-RNA immunoprecipitation (qDRIP), a method combining synthetic RNA-DNA-hybrid internal standards with high-resolution, strand-specific sequencing. We show that qDRIP avoids biases inherent to read-count normalization by accurately profiling signal in regions unaffected by transcription inhibition in human cells, and by facilitating accurate differential peak calling between conditions. We also use these quantitative comparisons to make the first estimates of the absolute count of RNA-DNA hybrids per cell and their half-lives genome-wide. Finally, we identify a subset of RNA-DNA hybrids with high GC skew which are partially resistant to RNase H. Overall, qDRIP allows for accurate normalization in conditions where R-loops are perturbed and for quantitative measurements that provide previously unattainable biological insights.

Transcription factors (TFs) can define distinct cellular identities despite nearly identical DNA-binding specificities. One mechanism for achieving regulatory specificity is DNA-guided TF cooperativity. Although in vitro studies suggest that it may be common, examples of such cooperativity remain scarce in cellular contexts. Here, we demonstrate how "Coordinator," a long DNA motif composed of common motifs bound by many basic helix-loop-helix (bHLH) and homeodomain (HD) TFs, uniquely defines the regulatory regions of embryonic face and limb mesenchyme. Coordinator guides cooperative and selective binding between the bHLH family mesenchymal regulator TWIST1 and a collective of HD factors associated with regional identities in the face and limb. TWIST1 is required for HD binding and open chromatin at Coordinator sites, whereas HD factors stabilize TWIST1 occupancy at Coordinator and titrate it away from HD-independent sites. This cooperativity results in the shared regulation of genes involved in cell-type and positional identities and ultimately shapes facial morphology and evolution.

The hormone estrogen (E2) binds the estrogen receptor to promote transcription of E2-responsive genes in the breast and other tissues. E2 also has links to genomic instability, and elevated E2 levels are tied to breast cancer. Here, we show that E2 stimulation causes a rapid, global increase in the formation of R-loops, co-transcriptional RNA-DNA products, which in some instances have been linked to DNA damage. We show that E2-dependent R-loop formation and breast cancer rearrangements are highly enriched at E2-responsive genomic loci and that E2 induces DNA replication-dependent double-strand breaks (DSBs). Strikingly, many DSBs that accumulate in response to E2 are R-loop dependent. Thus, R-loops resulting from the E2 transcriptional response are a significant source of DNA damage. This work reveals a novel mechanism by which E2 stimulation leads to genomic instability and highlights how transcriptional programs play an important role in shaping the genomic landscape of DNA damage susceptibility.

CD47 is a cell surface molecule that inhibits phagocytosis of cells that express it by binding to its receptor, SIRPα, on macrophages and other immune cells. CD47 is expressed at different levels by neoplastic and normal cells. Here, to reveal mechanisms by which different neoplastic cells generate this dominant 'don't eat me' signal, we analyse the CD47 regulatory genomic landscape. We identify two distinct super-enhancers (SEs) associated with CD47 in certain cancer cell types. We show that a set of active constituent enhancers, located within the two CD47 SEs, regulate CD47 expression in different cancer cell types and that disruption of CD47 SEs reduces CD47 gene expression. Finally we report that the TNF-NFKB1 signalling pathway directly regulates CD47 by interacting with a constituent enhancer located within a CD47-associated SE specific to breast cancer. These results suggest that cancers can evolve SE to drive CD47 overexpression to escape immune surveillance.