DNA-guided transcription factor cooperativity shapes face and limb mesenchyme.

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.

Single Amino Acid Change Underlies Distinct Roles of H2A.Z Subtypes in Human Syndrome.

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.

Transcription-Replication Conflict Orientation Modulates R-Loop Levels and Activates Distinct DNA Damage Responses.

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.

Structural elements promote architectural stripe formation and facilitate ultra-long-range gene regulation at a human disease locus.

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.

Co-transcriptional genome surveillance by HUSH is coupled to termination machinery.

The HUSH complex recognizes and silences foreign DNA such as viruses, transposons, and transgenes without prior exposure to its targets. Here, we show that endogenous targets of the HUSH complex fall into two distinct classes based on the presence or absence of H3K9me3. These classes are further distinguished by their transposon content and differential response to the loss of HUSH. A de novo genomic rearrangement at the Sox2 locus induces a switch from H3K9me3-independent to H3K9me3-associated HUSH targeting, resulting in silencing. We further demonstrate that HUSH interacts with the termination factor WDR82 and-via its component MPP8-with nascent RNA. HUSH accumulates at sites of high RNAPII occupancy including long exons and transcription termination sites in a manner dependent on WDR82 and CPSF. Together, our results uncover the functional diversity of HUSH targets and show that this vertebrate-specific complex exploits evolutionarily ancient transcription termination machinery for co-transcriptional chromatin targeting and genome surveillance.

Enhancer divergence and cis-regulatory evolution in the human and chimp neural crest.

cis-regulatory changes play a central role in morphological divergence, yet the regulatory principles underlying emergence of human traits remain poorly understood. Here, we use epigenomic profiling from human and chimpanzee cranial neural crest cells to systematically and quantitatively annotate divergence of craniofacial cis-regulatory landscapes. Epigenomic divergence is often attributable to genetic variation within TF motifs at orthologous enhancers, with a novel motif being most predictive of activity biases. We explore properties of this cis-regulatory change, revealing the role of particular retroelements, uncovering broad clusters of species-biased enhancers near genes associated with human facial variation, and demonstrating that cis-regulatory divergence is linked to quantitative expression differences of crucial neural crest regulators. Our work provides a wealth of candidates for future evolutionary studies and demonstrates the value of "cellular anthropology," a strategy of using in-vitro-derived embryonic cell types to elucidate both fundamental and evolving mechanisms underlying morphological variation in higher primates.

Mll3 and Mll4 Facilitate Enhancer RNA Synthesis and Transcription from Promoters Independently of H3K4 Monomethylation.

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.

Fallen immortals.

The memory of somatic cell gene expression is reset in the germline in a process that is accompanied by dramatic changes in chromatin modifications. In this issue, Katz et al. (2009) show that the histone demethylase Lsd1/Spr-5 may participate in this resetting process in the worm, thereby preventing a decline in germ cell epigenetic stability and viability over ensuing generations.

Jarid2/Jumonji coordinates control of PRC2 enzymatic activity and target gene occupancy in pluripotent cells.

Polycomb Repressive Complex 2 (PRC2) regulates key developmental genes in embryonic stem (ES) cells and during development. Here we show that Jarid2/Jumonji, a protein enriched in pluripotent cells and a founding member of the Jumonji C (JmjC) domain protein family, is a PRC2 subunit in ES cells. Genome-wide ChIP-seq analyses of Jarid2, Ezh2, and Suz12 binding reveal that Jarid2 and PRC2 occupy the same genomic regions. We further show that Jarid2 promotes PRC2 recruitment to the target genes while inhibiting PRC2 histone methyltransferase activity, suggesting that it acts as a "molecular rheostat" that finely calibrates PRC2 functions at developmental genes. Using Xenopus laevis as a model we demonstrate that Jarid2 knockdown impairs the induction of gastrulation genes in blastula embryos and results in failure of differentiation. Our findings illuminate a mechanism of histone methylation regulation in pluripotent cells and during early cell-fate transitions.

Selective silencing of euchromatic L1s revealed by genome-wide screens for L1 regulators.

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.

Tissue-selective effects of nucleolar stress and rDNA damage in developmental disorders.

Many craniofacial disorders are caused by heterozygous mutations in general regulators of housekeeping cellular functions such as transcription or ribosome biogenesis. Although it is understood that many of these malformations are a consequence of defects in cranial neural crest cells, a cell type that gives rise to most of the facial structures during embryogenesis, the mechanism underlying cell-type selectivity of these defects remains largely unknown. By exploring molecular functions of DDX21, a DEAD-box RNA helicase involved in control of both RNA polymerase (Pol) I- and II-dependent transcriptional arms of ribosome biogenesis, we uncovered a previously unappreciated mechanism linking nucleolar dysfunction, ribosomal DNA (rDNA) damage, and craniofacial malformations. Here we demonstrate that genetic perturbations associated with Treacher Collins syndrome, a craniofacial disorder caused by heterozygous mutations in components of the Pol I transcriptional machinery or its cofactor TCOF1 (ref. 1), lead to relocalization of DDX21 from the nucleolus to the nucleoplasm, its loss from the chromatin targets, as well as inhibition of rRNA processing and downregulation of ribosomal protein gene transcription. These effects are cell-type-selective, cell-autonomous, and involve activation of p53 tumour-suppressor protein. We further show that cranial neural crest cells are sensitized to p53-mediated apoptosis, but blocking DDX21 loss from the nucleolus and chromatin rescues both the susceptibility to apoptosis and the craniofacial phenotypes associated with Treacher Collins syndrome. This mechanism is not restricted to cranial neural crest cells, as blood formation is also hypersensitive to loss of DDX21 functions. Accordingly, ribosomal gene perturbations associated with Diamond-Blackfan anaemia disrupt DDX21 localization. At the molecular level, we demonstrate that impaired rRNA synthesis elicits a DNA damage response, and that rDNA damage results in tissue-selective and dosage-dependent effects on craniofacial development. Taken together, our findings illustrate how disruption in general regulators that compromise nucleolar homeostasis can result in tissue-selective malformations.

Genome scans of facial features in East Africans and cross-population comparisons reveal novel associations.

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.

FaceBase 3: analytical tools and FAIR resources for craniofacial and dental research.

The FaceBase Consortium was established by the National Institute of Dental and Craniofacial Research in 2009 as a 'big data' resource for the craniofacial research community. Over the past decade, researchers have deposited hundreds of annotated and curated datasets on both normal and disordered craniofacial development in FaceBase, all freely available to the research community on the FaceBase Hub website. The Hub has developed numerous visualization and analysis tools designed to promote integration of multidisciplinary data while remaining dedicated to the FAIR principles of data management (findability, accessibility, interoperability and reusability) and providing a faceted search infrastructure for locating desired data efficiently. Summaries of the datasets generated by the FaceBase projects from 2014 to 2019 are provided here. FaceBase 3 now welcomes contributions of data on craniofacial and dental development in humans, model organisms and cell lines. Collectively, the FaceBase Consortium, along with other NIH-supported data resources, provide a continuously growing, dynamic and current resource for the scientific community while improving data reproducibility and fulfilling data sharing requirements.

qDRIP: a method to quantitatively assess RNA-DNA hybrid formation genome-wide.

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.

A unique chromatin signature uncovers early developmental enhancers in humans.

Cell-fate transitions involve the integration of genomic information encoded by regulatory elements, such as enhancers, with the cellular environment. However, identification of genomic sequences that control human embryonic development represents a formidable challenge. Here we show that in human embryonic stem cells (hESCs), unique chromatin signatures identify two distinct classes of genomic elements, both of which are marked by the presence of chromatin regulators p300 and BRG1, monomethylation of histone H3 at lysine 4 (H3K4me1), and low nucleosomal density. In addition, elements of the first class are distinguished by the acetylation of histone H3 at lysine 27 (H3K27ac), overlap with previously characterized hESC enhancers, and are located proximally to genes expressed in hESCs and the epiblast. In contrast, elements of the second class, which we term 'poised enhancers', are distinguished by the absence of H3K27ac, enrichment of histone H3 lysine 27 trimethylation (H3K27me3), and are linked to genes inactive in hESCs and instead are involved in orchestrating early steps in embryogenesis, such as gastrulation, mesoderm formation and neurulation. Consistent with the poised identity, during differentiation of hESCs to neuroepithelium, a neuroectoderm-specific subset of poised enhancers acquires a chromatin signature associated with active enhancers. When assayed in zebrafish embryos, poised enhancers are able to direct cell-type and stage-specific expression characteristic of their proximal developmental gene, even in the absence of sequence conservation in the fish genome. Our data demonstrate that early developmental enhancers are epigenetically pre-marked in hESCs and indicate an unappreciated role of H3K27me3 at distal regulatory elements. Moreover, the wealth of new regulatory sequences identified here provides an invaluable resource for studies and isolation of transient, rare cell populations representing early stages of human embryogenesis.

CHD7 cooperates with PBAF to control multipotent neural crest formation.

Heterozygous mutations in the gene encoding the CHD (chromodomain helicase DNA-binding domain) member CHD7, an ATP-dependent chromatin remodeller homologous to the Drosophila trithorax-group protein Kismet, result in a complex constellation of congenital anomalies called CHARGE syndrome, which is a sporadic, autosomal dominant disorder characterized by malformations of the craniofacial structures, peripheral nervous system, ears, eyes and heart. Although it was postulated 25 years ago that CHARGE syndrome results from the abnormal development of the neural crest, this hypothesis remained untested. Here we show that, in both humans and Xenopus, CHD7 is essential for the formation of multipotent migratory neural crest (NC), a transient cell population that is ectodermal in origin but undergoes a major transcriptional reprogramming event to acquire a remarkably broad differentiation potential and ability to migrate throughout the body, giving rise to craniofacial bones and cartilages, the peripheral nervous system, pigmentation and cardiac structures. We demonstrate that CHD7 is essential for activation of the NC transcriptional circuitry, including Sox9, Twist and Slug. In Xenopus embryos, knockdown of Chd7 or overexpression of its catalytically inactive form recapitulates all major features of CHARGE syndrome. In human NC cells CHD7 associates with PBAF (polybromo- and BRG1-associated factor-containing complex) and both remodellers occupy a NC-specific distal SOX9 enhancer and a conserved genomic element located upstream of the TWIST1 gene. Consistently, during embryogenesis CHD7 and PBAF cooperate to promote NC gene expression and cell migration. Our work identifies an evolutionarily conserved role for CHD7 in orchestrating NC gene expression programs, provides insights into the synergistic control of distal elements by chromatin remodellers, illuminates the patho-embryology of CHARGE syndrome, and suggests a broader function for CHD7 in the regulation of cell motility.

Chromatin remodelling factor Mll1 is essential for neurogenesis from postnatal neural stem cells.

Epigenetic mechanisms that maintain neurogenesis throughout adult life remain poorly understood. Trithorax group (trxG) and Polycomb group (PcG) gene products are part of an evolutionarily conserved chromatin remodelling system that activate or silence gene expression, respectively. Although PcG member Bmi1 has been shown to be required for postnatal neural stem cell self-renewal, the role of trxG genes remains unknown. Here we show that the trxG member Mll1 (mixed-lineage leukaemia 1) is required for neurogenesis in the mouse postnatal brain. Mll1-deficient subventricular zone neural stem cells survive, proliferate and efficiently differentiate into glial lineages; however, neuronal differentiation is severely impaired. In Mll1-deficient cells, early proneural Mash1 (also known as Ascl1) and gliogenic Olig2 expression are preserved, but Dlx2, a key downstream regulator of subventricular zone neurogenesis, is not expressed. Overexpression of Dlx2 can rescue neurogenesis in Mll1-deficient cells. Chromatin immunoprecipitation demonstrates that Dlx2 is a direct target of MLL in subventricular zone cells. In differentiating wild-type subventricular zone cells, Mash1, Olig2 and Dlx2 loci have high levels of histone 3 trimethylated at lysine 4 (H3K4me3), consistent with their transcription. In contrast, in Mll1-deficient subventricular zone cells, chromatin at Dlx2 is bivalently marked by both H3K4me3 and histone 3 trimethylated at lysine 27 (H3K27me3), and the Dlx2 gene fails to properly activate. These data support a model in which Mll1 is required to resolve key silenced bivalent loci in postnatal neural precursors to the actively transcribed state for the induction of neurogenesis, but not for gliogenesis.

A common cis-regulatory variant impacts normal-range and disease-associated human facial shape through regulation of PKDCC during chondrogenesis.

Genome-wide association studies (GWAS) identified thousands of genetic variants linked to phenotypic traits and disease risk. However, mechanistic understanding of how GWAS variants influence complex morphological traits and can, in certain cases, simultaneously confer normal-range phenotypic variation and disease predisposition, is still largely lacking. Here, we focus on rs6740960, a single nucleotide polymorphism (SNP) at the 2p21 locus, which in GWAS studies has been associated both with normal-range variation in jaw shape and with an increased risk of non-syndromic orofacial clefting. Using in vitro derived embryonic cell types relevant for human facial morphogenesis, we show that this SNP resides in an enhancer that regulates chondrocytic expression of PKDCC - a gene encoding a tyrosine kinase involved in chondrogenesis and skeletal development. In agreement, we demonstrate that the rs6740960 SNP is sufficient to confer chondrocyte-specific differences in PKDCC expression. By deploying dense landmark morphometric analysis of skull elements in mice, we show that changes in Pkdcc dosage are associated with quantitative changes in the maxilla, mandible, and palatine bone shape that are concordant with the facial phenotypes and disease predisposition seen in humans. We further demonstrate that the frequency of the rs6740960 variant strongly deviated among different human populations, and that the activity of its cognate enhancer diverged in hominids. Our study provides a mechanistic explanation of how a common SNP can mediate normal-range and disease-associated morphological variation, with implications for the evolution of human facial features.

A genome-wide genetic screen uncovers determinants of human pigmentation.

Skin color, one of the most diverse human traits, is determined by the quantity, type, and distribution of melanin. In this study, we leveraged the light-scattering properties of melanin to conduct a genome-wide screen for regulators of melanogenesis. We identified 169 functionally diverse genes that converge on melanosome biogenesis, endosomal transport, and gene regulation, of which 135 represented previously unknown associations with pigmentation. In agreement with their melanin-promoting function, the majority of screen hits were up-regulated in melanocytes from darkly pigmented individuals. We further unraveled functions of KLF6 as a transcription factor that regulates melanosome maturation and pigmentation in vivo, and of the endosomal trafficking protein COMMD3 in modulating melanosomal pH. Our study reveals a plethora of melanin-promoting genes, with broad implications for human variation, cell biology, and medicine.

Precise modulation of transcription factor levels identifies features underlying dosage sensitivity.

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.