Psoriasis genetic background depends on polygenic and multifactorial mode of inheritance. As in other complex disorders, the estimation of the disease risk based on individual genetic variants is impossible. For this reason, recent investigations have been focused on combinations of known psoriasis susceptibility markers in order to improve the disease risk evaluation. Our aim was to compare psoriasis genetic risk score (GRS) for five susceptibility loci involved in the immunological response (HLA-C, ERAP1, ZAP70) and in the skin barrier function (LCE3, CSTA) between patients with chronic plaque psoriasis (n = 148) and the control group (n = 146). A significantly higher number of predisposing alleles was observed in patients with psoriasis in comparison to healthy individuals (6.1 vs. 5.2, respectively; P = 8.8×10-7). The statistical significance was even more profound when GRS weighted by logarithm odds ratios was evaluated (P = 9.9×10-14). Our results demonstrate the developed panel of five susceptibility loci to be more efficient in predicting psoriasis risk in the Polish population and to possess higher sensitivity and specificity for the disease than any of the markers analyzed separately, including the most informative HLA-C*06 allele.

Naive and primed pluripotency is characterized by distinct signaling requirements, transcriptomes, and developmental properties, but both cellular states share key transcriptional regulators: Oct4, Sox2, and Nanog. Here, we demonstrate that transition between these two pluripotent states is associated with widespread Oct4 relocalization, mirrored by global rearrangement of enhancer chromatin landscapes. Our genomic and biochemical analyses identified candidate mediators of primed state-specific Oct4 binding, including Otx2 and Zic2/3. Even when differentiation cues are blocked, premature Otx2 overexpression is sufficient to exit the naive state, induce transcription of a substantial subset of primed pluripotency-associated genes, and redirect Oct4 to previously inaccessible enhancer sites. However, the ability of Otx2 to engage new enhancer regions is determined by its levels, cis-encoded properties of the sites, and the signaling environment. Our results illuminate regulatory mechanisms underlying pluripotency and suggest that the capacity of transcription factors such as Otx2 and Oct4 to pioneer new enhancer sites is highly context dependent.

Following implantation, mouse epiblast cells transit from a naive to a primed state in which they are competent for both somatic and primordial germ cell (PGC) specification. Using mouse embryonic stem cells as an in vitro model to study the transcriptional regulatory principles orchestrating peri-implantation development, here we show that the transcription factor Foxd3 is necessary for exit from naive pluripotency and progression to a primed pluripotent state. During this transition, Foxd3 acts as a repressor that dismantles a significant fraction of the naive pluripotency expression program through decommissioning of active enhancers associated with key naive pluripotency and early germline genes. Subsequently, Foxd3 needs to be silenced in primed pluripotent cells to allow re-activation of relevant genes required for proper PGC specification. Our findings therefore uncover a cycle of activation and deactivation of Foxd3 required for exit from naive pluripotency and subsequent PGC specification.

The association of guttate psoriasis with a streptococcal throat infection and HLA-Cw*06 allele is well established in different populations. Nevertheless, only few studies on this form of disease have been performed in the Polish population.

Emergence of form and function during embryogenesis arises in large part through cell-type- and cell-state-specific variation in gene expression patterns, mediated by specialized cis-regulatory elements called enhancers. Recent large-scale epigenomic mapping revealed unexpected complexity and dynamics of enhancer utilization patterns, with 400,000 putative human enhancers annotated by the ENCODE project alone. These large-scale efforts were largely enabled through the understanding that enhancers share certain stereotypical chromatin features. However, an important question still lingers: what is the functional significance of enhancer chromatin modification? Here we give an overview of enhancer-associated modifications of histones and DNA and discuss enzymatic activities involved in their dynamic deposition and removal. We describe potential downstream effectors of these marks and propose models for exploring functions of chromatin modification in regulating enhancer activity during development.

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.

The human herpes simplex virus (HSV) host cell factor HCF-1 is a transcriptional coregulator that associates with both histone methyl- and acetyltransferases, and a histone deacetylase and regulates cell proliferation and division. In HSV-infected cells, HCF-1 associates with the viral protein VP16 to promote formation of a multiprotein-DNA transcriptional activator complex. The ability of HCF proteins to stabilize this VP16-induced complex has been conserved in diverse animal species including Drosophila melanogaster and Caenorhabditis elegans suggesting that VP16 targets a conserved cellular function of HCF-1.

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.

Inappropriate activation of the p53 transcription factor contributes to numerous developmental syndromes characterized by distinct constellations of phenotypes. How p53 drives exquisitely specific sets of symptoms in diverse syndromes, however, remains enigmatic. Here, we deconvolute the basis of p53-driven developmental syndromes by leveraging an array of mouse strains to modulate the spatial expression pattern, temporal profile, and magnitude of p53 activation during embryogenesis. We demonstrate that inappropriate p53 activation in the neural crest, facial ectoderm, anterior heart field, and endothelium induces distinct spectra of phenotypes. Moreover, altering the timing and degree of p53 hyperactivation substantially affects the phenotypic outcomes. Phenotypes are associated with p53-driven cell-cycle arrest or apoptosis, depending on the cell type, with gene expression programs, rather than extent of mitochondrial priming, largely governing the specific response. Together, our findings provide a critical framework for decoding the role of p53 as a mediator of diverse developmental syndromes.

Gene regulatory divergence is thought to play a central role in determining human-specific traits. However, our ability to link divergent regulation to divergent phenotypes is limited. Here, we utilized human-chimpanzee hybrid induced pluripotent stem cells to study gene expression separating these species. The tetraploid hybrid cells allowed us to separate cis- from trans-regulatory effects, and to control for nongenetic confounding factors. We differentiated these cells into cranial neural crest cells, the primary cell type giving rise to the face. We discovered evidence of lineage-specific selection on the hedgehog signaling pathway, including a human-specific sixfold down-regulation of EVC2 (LIMBIN), a key hedgehog gene. Inducing a similar down-regulation of EVC2 substantially reduced hedgehog signaling output. Mice and humans lacking functional EVC2 show striking phenotypic parallels to human-chimpanzee craniofacial differences, suggesting that the regulatory divergence of hedgehog signaling may have contributed to the unique craniofacial morphology of humans.

The analysis of contemporary genomic data typically operates on one-dimensional phenotypic measurements (e.g. standing height). Here we report on a data-driven, family-informed strategy to facial phenotyping that searches for biologically relevant traits and reduces multivariate 3D facial shape variability into amendable univariate measurements, while preserving its structurally complex nature. We performed a biometric identification of siblings in a sample of 424 children, defining 1,048 sib-shared facial traits. Subsequent quantification and analyses in an independent European cohort (n = 8,246) demonstrated significant heritability for a subset of traits (0.17-0.53) and highlighted 218 genome-wide significant loci (38 also study-wide) associated with facial variation shared by siblings. These loci showed preferential enrichment for active chromatin marks in cranial neural crest cells and embryonic craniofacial tissues and several regions harbor putative craniofacial genes, thereby enhancing our knowledge on the genetic architecture of normal-range facial variation.

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.

Many genes are regulated by multiple enhancers that often simultaneously activate their target gene. However, how individual enhancers collaborate to activate transcription is not well understood. Here, we dissect the functions and interdependencies of five enhancer elements that together activate Fgf5 expression during exit from naive murine pluripotency. Four intergenic elements form a super-enhancer, and most of the elements contribute to Fgf5 induction at distinct time points. A fifth, poised enhancer located in the first intron contributes to Fgf5 expression at every time point by amplifying overall Fgf5 expression levels. Despite low individual enhancer activity, together these elements strongly induce Fgf5 expression in a super-additive fashion that involves strong accumulation of RNA polymerase II at the intronic enhancer. Finally, we observe a strong anti-correlation between RNA polymerase II levels at enhancers and their distance to the closest promoter, and we identify candidate elements with properties similar to the intronic enhancer.

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.

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.

Chromosome instability and aneuploidies occur very frequently in human embryos, impairing proper embryogenesis and leading to cell cycle arrest, loss of cell viability, and developmental failures in 50-80% of cleavage-stage embryos. This high frequency of cellular extinction events represents a significant experimental obstacle challenging analyses of individual cells isolated from human preimplantation embryos. We carried out single cell expression profiling of 241 individual cells recovered from 32 human embryos during the early and late stages of viable human blastocyst (VHB) differentiation. Classification of embryonic cells was performed solely based on expression patterns of human pluripotency-associated transcripts (HPAT), which represent a family of primate-specific transposable element-derived lincRNAs highly expressed in human embryonic stem cells and regulating nuclear reprogramming and pluripotency induction. We then validated our findings by analyzing transcriptomes of 1,708 individual cells recovered from more than 100 human embryos and 259 mouse cells from more than 40 mouse embryos at different stages of preimplantation embryogenesis. HPAT's expression-guided spatiotemporal reconstruction of human embryonic development inferred from single-cell expression analysis of VHB differentiation enabled identification of telomerase-positive embryonic cells co-expressing key pluripotency regulatory genes and genetic markers of three major lineages. Follow-up validation analyses confirmed the emergence in human embryos prior to lineage segregation of telomerase-positive cells co-expressing genetic markers of multiple lineages. Observations reported in this contribution support the hypothesis of a developmental pathway of creation embryonic lineages and extraembryonic tissues from telomerase-positive pre-lineage cells manifesting multi-lineage precursor phenotype.

We report the identification of 67 previously undescribed histone modifications, increasing the current number of known histone marks by about 70%. We further investigated one of the marks, lysine crotonylation (Kcr), confirming that it represents an evolutionarily-conserved histone posttranslational modification. The unique structure and genomic localization of histone Kcr suggest that it is mechanistically and functionally different from histone lysine acetylation (Kac). Specifically, in both human somatic and mouse male germ cell genomes, histone Kcr marks either active promoters or potential enhancers. In male germinal cells immediately following meiosis, Kcr is enriched on sex chromosomes and specifically marks testis-specific genes, including a significant proportion of X-linked genes that escape sex chromosome inactivation in haploid cells. These results therefore dramatically extend the repertoire of histone PTM sites and designate Kcr as a specific mark of active sex chromosome-linked genes in postmeiotic male germ cells.

RNA functions at enhancers remain mysterious. Here we show that the 7SK small nuclear RNA (snRNA) inhibits enhancer transcription by modulating nucleosome position. 7SK occupies enhancers and super enhancers genome wide in mouse and human cells, and it is required to limit enhancer-RNA initiation and synthesis in a manner distinct from promoter pausing. Clustered elements at super enhancers uniquely require 7SK to prevent convergent transcription and DNA-damage signaling. 7SK physically interacts with the BAF chromatin-remodeling complex, recruits BAF to enhancers and inhibits enhancer transcription by modulating chromatin structure. In turn, 7SK occupancy at enhancers coincides with that of Brd4 and is exquisitely sensitive to the bromodomain inhibitor JQ1. Thus, 7SK uses distinct mechanisms to counteract the diverse consequences of pervasive transcription that distinguish super enhancers, enhancers and promoters.

Eukaryotic genomes are transcribed into numerous regulatory long non-coding RNAs (lncRNAs). Compared to mRNAs, lncRNAs display higher developmental stage-, tissue-, and cell-subtype-specificity of expression, and are generally less abundant in a population of cells. Despite the progress in single-cell-focused research, the origins of low population-level expression of lncRNAs in homogeneous populations of cells are poorly understood. Here, we identify HIPSTR (Heterogeneously expressed from the Intronic Plus Strand of the TFAP2A-locus RNA), a novel lncRNA gene in the developmentally regulated TFAP2A locus. HIPSTR has evolutionarily conserved expression patterns, its promoter is most active in undifferentiated cells, and depletion of HIPSTR in HEK293 and in pluripotent H1BP cells predominantly affects the genes involved in early organismal development and cell differentiation. Most importantly, we find that HIPSTR is specifically induced and heterogeneously expressed in the 8-cell-stage human embryos during the major wave of embryonic genome activation. We systematically explore the phenomenon of cell-to-cell variation of gene expression and link it to low population-level expression of lncRNAs, showing that, similar to HIPSTR, the expression of thousands of lncRNAs is more highly heterogeneous than the expression of mRNAs in the individual, otherwise indistinguishable cells of totipotent human embryos, primordial germ cells, and stable cell lines.

During preimplantation development, the embryo must establish totipotency and enact the earliest differentiation choices, processes that involve extensive chromatin modification. To identify novel developmental regulators, we screened for genes that are preferentially transcribed in the pluripotent inner cell mass (ICM) of the mouse blastocyst. Genes that encode chromatin remodeling factors were prominently represented in the ICM, including Chd1l, a member of the Snf2 gene family. Chd1l is developmentally regulated and expressed in embryonic stem (ES) cells, but its role in development has not been investigated. Here we show that inhibiting Chd1l protein production by microinjection of antisense morpholinos causes arrest prior to the blastocyst stage. Despite this important function in vivo, Chd1l is non-essential for cultured ES cell survival, pluripotency, or differentiation, suggesting that Chd1l is vital for events in embryos that are distinct from events in ES cells. Our data reveal a novel role for the chromatin remodeling factor Chd1l in the earliest cell divisions of mammalian development.