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.

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.

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.

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.

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.

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.

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.

Reactivation of the pluripotency program precedes formation of the cranial neural crest.

During development, cells progress from a pluripotent state to a more restricted fate within a particular germ layer. However, cranial neural crest cells (CNCCs), a transient cell population that generates most of the craniofacial skeleton, have much broader differentiation potential than their ectodermal lineage of origin. Here, we identify a neuroepithelial precursor population characterized by expression of canonical pluripotency transcription factors that gives rise to CNCCs and is essential for craniofacial development. Pluripotency factor Oct4 is transiently reactivated in CNCCs and is required for the subsequent formation of ectomesenchyme. Furthermore, open chromatin landscapes of Oct4+ CNCC precursors resemble those of epiblast stem cells, with additional features suggestive of priming for mesenchymal programs. We propose that CNCCs expand their developmental potential through a transient reacquisition of molecular signatures of pluripotency.

Insights into the genetic architecture of the human face.

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.

Loss of Extreme Long-Range Enhancers in Human Neural Crest Drives a Craniofacial Disorder.

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.

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.

Zscan4 binds nucleosomal microsatellite DNA and protects mouse two-cell embryos from DNA damage.

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.

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.

Systematic perturbation of retroviral LTRs reveals widespread long-range effects on human gene regulation.

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.

Genome-wide mapping of global-to-local genetic effects on human facial shape.

Genome-wide association scans of complex multipartite traits like the human face typically use preselected phenotypic measures. Here we report a data-driven approach to phenotyping facial shape at multiple levels of organization, allowing for an open-ended description of facial variation while preserving statistical power. In a sample of 2,329 persons of European ancestry, we identified 38 loci, 15 of which replicated in an independent European sample (n = 1,719). Four loci were completely new. For the others, additional support (n = 9) or pleiotropic effects (n = 2) were found in the literature, but the results reported here were further refined. All 15 replicated loci highlighted distinctive patterns of global-to-local genetic effects on facial shape and showed enrichment for active chromatin elements in human cranial neural crest cells, suggesting an early developmental origin of the facial variation captured. These results have implications for studies of facial genetics and other complex morphological traits.

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.

Transcription-coupled changes in nuclear mobility of mammalian cis-regulatory elements.

To achieve guide RNA (gRNA) multiplexing and an efficient delivery of tens of distinct gRNAs into single cells, we developed a molecular assembly strategy termed chimeric array of gRNA oligonucleotides (CARGO). We coupled CARGO with dCas9 (catalytically dead Cas9) imaging to quantitatively measure the movement of enhancers and promoters that undergo differentiation-associated activity changes in live embryonic stem cells. Whereas all examined functional elements exhibited subdiffusive behavior, their relative mobility increased concurrently with transcriptional activation. Furthermore, acute perturbation of RNA polymerase II activity can reverse these activity-linked increases in loci mobility. Through quantitative CARGO-dCas9 imaging, we provide direct measurements of cis-regulatory element dynamics in living cells and distinct cellular and activity states and uncover an intrinsic connection between cis-regulatory element mobility and transcription.

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.

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.