The homeodomain of Oct4 is a dimeric binder of methylated CpG elements.

Oct4 is essential to maintain pluripotency and has a pivotal role in establishing the germline. Its DNA-binding POU domain was recently found to bind motifs with methylated CpG elements normally associated with epigenetic silencing. However, the mode of binding and the consequences of this capability has remained unclear. Here, we show that Oct4 binds to a compact palindromic DNA element with a methylated CpG core (CpGpal) in alternative states of pluripotency and during cellular reprogramming towards induced pluripotent stem cells (iPSCs). During cellular reprogramming, typical Oct4 bound enhancers are uniformly demethylated, with the prominent exception of the CpGpal sites where DNA methylation is often maintained. We demonstrate that Oct4 cooperatively binds the CpGpal element as a homodimer, which contrasts with the ectoderm-expressed POU factor Brn2. Indeed, binding to CpGpal is Oct4-specific as other POU factors expressed in somatic cells avoid this element. Binding assays combined with structural analyses and molecular dynamic simulations show that dimeric Oct4-binding to CpGpal is driven by the POU-homeodomain whilst the POU-specific domain is detached from DNA. Collectively, we report that Oct4 exerts parts of its regulatory function in the context of methylated DNA through a DNA recognition mechanism that solely relies on its homeodomain.

Evaluation of the determinants for improved pluripotency induction and maintenance by engineered SOX17.

An engineered SOX17 variant with point mutations within its DNA binding domain termed SOX17FNV is a more potent pluripotency inducer than SOX2, yet the underlying mechanism remains unclear. Although wild-type SOX17 was incapable of inducing pluripotency, SOX17FNV outperformed SOX2 in mouse and human pluripotency reprogramming. In embryonic stem cells, SOX17FNV could replace SOX2 to maintain pluripotency despite considerable sequence differences and upregulated genes expressed in cleavage-stage embryos. Mechanistically, SOX17FNV co-bound OCT4 more cooperatively than SOX2 in the context of the canonical SoxOct DNA element. SOX2, SOX17, and SOX17FNV were all able to bind nucleosome core particles in vitro, which is a prerequisite for pioneer transcription factors. Experiments using purified proteins and in cellular contexts showed that SOX17 variants phase-separated more efficiently than SOX2, suggesting an enhanced ability to self-organise. Systematic deletion analyses showed that the N-terminus of SOX17FNV was dispensable for its reprogramming activity. However, the C-terminus encodes essential domains indicating multivalent interactions that drive transactivation and reprogramming. We defined a minimal SOX17FNV (miniSOX) that can support reprogramming with high activity, reducing the payload of reprogramming cassettes. This study uncovers the mechanisms behind SOX17FNV-induced pluripotency and establishes engineered SOX factors as powerful cell engineering tools.

SOX9 and SOX10 control fluid homeostasis in the inner ear for hearing through independent and cooperative mechanisms.

The in vivo mechanisms underlying dominant syndromes caused by mutations in SRY-Box Transcription Factor 9 (SOX9) and SOX10 (SOXE) transcription factors, when they either are expressed alone or are coexpressed, are ill-defined. We created a mouse model for the campomelic dysplasia SOX9Y440X mutation, which truncates the transactivation domain but leaves DNA binding and dimerization intact. Here, we find that SOX9Y440X causes deafness via distinct mechanisms in the endolymphatic sac (ES)/duct and cochlea. By contrast, conditional heterozygous Sox9-null mice are normal. During the ES development of Sox9Y440X/+ heterozygotes, Sox10 and genes important for ionic homeostasis are down-regulated, and there is developmental persistence of progenitors, resulting in fewer mature cells. Sox10 heterozygous null mutants also display persistence of ES/duct progenitors. By contrast, SOX10 retains its expression in the early Sox9Y440X/+ mutant cochlea. Later, in the postnatal stria vascularis, dominant interference by SOX9Y440X is implicated in impairing the normal cooperation of SOX9 and SOX10 in repressing the expression of the water channel Aquaporin 3, thereby contributing to endolymphatic hydrops. Our study shows that for a functioning endolymphatic system in the inner ear, SOX9 regulates Sox10, and depending on the cell type and target gene, it works either independently of or cooperatively with SOX10. SOX9Y440X can interfere with the activity of both SOXE factors, exerting effects that can be classified as haploinsufficient/hypomorphic or dominant negative depending on the cell/gene context. This model of disruption of transcription factor partnerships may be applicable to congenital deafness, which affects ∼0.3% of newborns, and other syndromic disorders.

OCT4 interprets and enhances nucleosome flexibility.

Pioneer transcription factors are proteins that induce cellular identity transitions by binding to inaccessible regions of DNA in nuclear chromatin. They contribute to chromatin opening and recruit other factors to regulatory DNA elements. The structural features and dynamics modulating their interaction with nucleosomes are still unresolved. From a combination of experiments and molecular simulations, we reveal here how the pioneer factor and master regulator of pluripotency, Oct4, interprets and enhances nucleosome structural flexibility. The magnitude of Oct4's impact on nucleosome dynamics depends on the binding site position and the mobility of the unstructured tails of nucleosomal histone proteins. Oct4 uses both its DNA binding domains to propagate and stabilize open nucleosome conformations, one for specific sequence recognition and the other for nonspecific interactions with nearby regions of DNA. Our findings provide a structural basis for the versatility of transcription factors in engaging with nucleosomes and have implications for understanding how pioneer factors induce chromatin dynamics.

Transposable element sequence fragments incorporated into coding and noncoding transcripts modulate the transcriptome of human pluripotent stem cells.

Transposable elements (TEs) occupy nearly 40% of mammalian genomes and, whilst most are fragmentary and no longer capable of transposition, they can nevertheless contribute to cell function. TEs within genes transcribed by RNA polymerase II can be copied as parts of primary transcripts; however, their full contribution to mature transcript sequences remains unresolved. Here, using long and short read (LR and SR) RNA sequencing data, we show that 26% of coding and 65% of noncoding transcripts in human pluripotent stem cells (hPSCs) contain TE-derived sequences. Different TE families are incorporated into RNAs in unique patterns, with consequences to transcript structure and function. The presence of TE sequences within a transcript is correlated with TE-type specific changes in its subcellular distribution, alterations in steady-state levels and half-life, and differential association with RNA Binding Proteins (RBPs). We identify hPSC-specific incorporation of endogenous retroviruses (ERVs) and LINE:L1 into protein-coding mRNAs, which generate TE sequence-derived peptides. Finally, single cell RNA-seq reveals that hPSCs express ERV-containing transcripts, whilst differentiating subpopulations lack ERVs and express SINE and LINE-containing transcripts. Overall, our comprehensive analysis demonstrates that the incorporation of TE sequences into the RNAs of hPSCs is more widespread and has a greater impact than previously appreciated.

Directed Evolution of an Enhanced POU Reprogramming Factor for Cell Fate Engineering.

Transcription factor-driven cell fate engineering in pluripotency induction, transdifferentiation, and forward reprogramming requires efficiency, speed, and maturity for widespread adoption and clinical translation. Here, we used Oct4, Sox2, Klf4, and c-Myc driven pluripotency reprogramming to evaluate methods for enhancing and tailoring cell fate transitions, through directed evolution with iterative screening of pooled mutant libraries and phenotypic selection. We identified an artificially evolved and enhanced POU factor (ePOU) that substantially outperforms wild-type Oct4 in terms of reprogramming speed and efficiency. In contrast to Oct4, not only can ePOU induce pluripotency with Sox2 alone, but it can also do so in the absence of Sox2 in a three-factor ePOU/Klf4/c-Myc cocktail. Biochemical assays combined with genome-wide analyses showed that ePOU possesses a new preference to dimerize on palindromic DNA elements. Yet, the moderate capacity of Oct4 to function as a pioneer factor, its preference to bind octamer DNA and its capability to dimerize with Sox2 and Sox17 proteins remain unchanged in ePOU. Compared with Oct4, ePOU is thermodynamically stabilized and persists longer in reprogramming cells. In consequence, ePOU: 1) differentially activates several genes hitherto not implicated in reprogramming, 2) reveals an unappreciated role of thyrotropin-releasing hormone signaling, and 3) binds a distinct class of retrotransposons. Collectively, these features enable ePOU to accelerate the establishment of the pluripotency network. This demonstrates that the phenotypic selection of novel factor variants from mammalian cells with desired properties is key to advancing cell fate conversions with artificially evolved biomolecules.

IL10RA modulates crizotinib sensitivity in NPM1-ALK+ anaplastic large cell lymphoma.

Anaplastic large cell lymphoma (ALCL) is a T-cell malignancy predominantly driven by a hyperactive anaplastic lymphoma kinase (ALK) fusion protein. ALK inhibitors, such as crizotinib, provide alternatives to standard chemotherapy with reduced toxicity and side effects. Children with lymphomas driven by nucleophosmin 1 (NPM1)-ALK fusion proteins achieved an objective response rate to ALK inhibition therapy of 54% to 90% in clinical trials; however, a subset of patients progressed within the first 3 months of treatment. The mechanism for the development of ALK inhibitor resistance is unknown. Through genome-wide clustered regularly interspaced short palindromic repeats (CRISPR) activation and knockout screens in ALCL cell lines, combined with RNA sequencing data derived from ALK inhibitor-relapsed patient tumors, we show that resistance to ALK inhibition by crizotinib in ALCL can be driven by aberrant upregulation of interleukin 10 receptor subunit alpha (IL10RA). Elevated IL10RA expression rewires the STAT3 signaling pathway, bypassing otherwise critical phosphorylation by NPM1-ALK. IL-10RA expression does not correlate with response to standard chemotherapy in pediatric patients, suggesting that a combination of crizotinib and chemotherapy could prevent ALK inhibitor resistance-specific relapse.

Concurrent binding to DNA and RNA facilitates the pluripotency reprogramming activity of Sox2.

Some transcription factors that specifically bind double-stranded DNA appear to also function as RNA-binding proteins. Here, we demonstrate that the transcription factor Sox2 is able to directly bind RNA in vitro as well as in mouse and human cells. Sox2 targets RNA via a 60-amino-acid RNA binding motif (RBM) positioned C-terminally of the DNA binding high mobility group (HMG) box. Sox2 can associate with RNA and DNA simultaneously to form ternary RNA/Sox2/DNA complexes. Deletion of the RBM does not affect selection of target genes but mitigates binding to pluripotency related transcripts, switches exon usage and impairs the reprogramming of somatic cells to a pluripotent state. Our findings designate Sox2 as a multi-functional factor that associates with RNA whilst binding to cognate DNA sequences, suggesting that it may co-transcriptionally regulate RNA metabolism during somatic cell reprogramming.

SOX17 in cellular reprogramming and cancer.

SOX17 is a transcription factor directing the specification and development of the primitive endoderm, primitive germ cells, definitive endoderm and, subsequently, is involved in the cardiovascular system and several endoderm-derived organs. The analysis of cancer genome sequencing data classified SOX17 as mutated cancer driver gene in endometrial cancer. These studies identified hotspot missense mutations within its DNA binding and transactivation domains. In somatic cell reprogramming, structure-based protein re-engineering showed a single missense mutation in SOX17 can change the DNA dependent heterodimer formation with OCT4 and enables the replacement of SOX2 with SOX17 mutants to induce pluripotency. This reveals the profound impact of specific missense mutations on gene function and regulatory activity. Here, we review the roles of SOX17 in cancer and discuss its cross-talk with the WNT/ß-catenin pathway, potentially reconciling its activity as re-engineered reprogramming factor and mutated cancer driver gene.

Capturing the interactome of newly transcribed RNA.

We combine the labeling of newly transcribed RNAs with 5-ethynyluridine with the characterization of bound proteins. This approach, named capture of the newly transcribed RNA interactome using click chemistry (RICK), systematically captures proteins bound to a wide range of RNAs, including nascent RNAs and traditionally neglected nonpolyadenylated RNAs. RICK has identified mitotic regulators amongst other novel RNA-binding proteins with preferential affinity for nonpolyadenylated RNAs, revealed a link between metabolic enzymes/factors and nascent RNAs, and expanded the known RNA-bound proteome of mouse embryonic stem cells. RICK will facilitate an in-depth interrogation of the total RNA-bound proteome in different cells and systems.

Synergistic co-regulation and competition by a SOX9-GLI-FOXA phasic transcriptional network coordinate chondrocyte differentiation transitions.

The growth plate mediates bone growth where SOX9 and GLI factors control chondrocyte proliferation, differentiation and entry into hypertrophy. FOXA factors regulate hypertrophic chondrocyte maturation. How these factors integrate into a Gene Regulatory Network (GRN) controlling these differentiation transitions is incompletely understood. We adopted a genome-wide whole tissue approach to establish a Growth Plate Differential Gene Expression Library (GP-DGEL) for fractionated proliferating, pre-hypertrophic, early and late hypertrophic chondrocytes, as an overarching resource for discovery of pathways and disease candidates. De novo motif discovery revealed the enrichment of SOX9 and GLI binding sites in the genes preferentially expressed in proliferating and prehypertrophic chondrocytes, suggesting the potential cooperation between SOX9 and GLI proteins. We integrated the analyses of the transcriptome, SOX9, GLI1 and GLI3 ChIP-seq datasets, with functional validation by transactivation assays and mouse mutants. We identified new SOX9 targets and showed SOX9-GLI directly and cooperatively regulate many genes such as Trps1, Sox9, Sox5, Sox6, Col2a1, Ptch1, Gli1 and Gli2. Further, FOXA2 competes with SOX9 for the transactivation of target genes. The data support a model of SOX9-GLI-FOXA phasic GRN in chondrocyte development. Together, SOX9-GLI auto-regulate and cooperate to activate and repress genes in proliferating chondrocytes. Upon hypertrophy, FOXA competes with SOX9, and control toward terminal differentiation passes to FOXA, RUNX, AP1 and MEF2 factors.

Mutations involving the SRY-related gene SOX8 are associated with a spectrum of human reproductive anomalies.

SOX8 is an HMG-box transcription factor closely related to SRY and SOX9. Deletion of the gene encoding Sox8 in mice causes reproductive dysfunction but the role of SOX8 in humans is unknown. Here, we show that SOX8 is expressed in the somatic cells of the early developing gonad in the human and influences human sex determination. We identified two individuals with 46, XY disorders/differences in sex development (DSD) and chromosomal rearrangements encompassing the SOX8 locus and a third individual with 46, XY DSD and a missense mutation in the HMG-box of SOX8. In vitro functional assays indicate that this mutation alters the biological activity of the protein. As an emerging body of evidence suggests that DSDs and infertility can have common etiologies, we also analysed SOX8 in a cohort of infertile men (n = 274) and two independent cohorts of women with primary ovarian insufficiency (POI; n = 153 and n = 104). SOX8 mutations were found at increased frequency in oligozoospermic men (3.5%; P < 0.05) and POI (5.06%; P = 4.5 × 10-5) as compared with fertile/normospermic control populations (0.74%). The mutant proteins identified altered SOX8 biological activity as compared with the wild-type protein. These data demonstrate that SOX8 plays an important role in human reproduction and SOX8 mutations contribute to a spectrum of phenotypes including 46, XY DSD, male infertility and 46, XX POI.

Pathogenic variants in the DEAH-box RNA helicase DHX37 are a frequent cause of 46,XY gonadal dysgenesis and 46,XY testicular regression syndrome.

PURPOSE: XY individuals with disorders/differences of sex development (DSD) are characterized by reduced androgenization caused, in some children, by gonadal dysgenesis or testis regression during fetal development. The genetic etiology for most patients with 46,XY gonadal dysgenesis and for all patients with testicular regression syndrome (TRS) is unknown. METHODS: We performed exome and/or Sanger sequencing in 145 individuals with 46,XY DSD of unknown etiology including gonadal dysgenesis and TRS. RESULTS: Thirteen children carried heterozygous missense pathogenic variants involving the RNA helicase DHX37, which is essential for ribosome biogenesis. Enrichment of rare/novel DHX37 missense variants in 46,XY DSD is highly significant compared with controls (P value = 5.8 × 10-10). Five variants are de novo (P value = 1.5 × 10-5). Twelve variants are clustered in two highly conserved functional domains and were specifically associated with gonadal dysgenesis and TRS. Consistent with a role in early testis development, DHX37 is expressed specifically in somatic cells of the developing human and mouse testis. CONCLUSION: DHX37 pathogenic variants are a new cause of an autosomal dominant form of 46,XY DSD, including gonadal dysgenesis and TRS, showing that these conditions are part of a clinical spectrum. This raises the possibility that some forms of DSD may be a ribosomopathy.

DNA-mediated dimerization on a compact sequence signature controls enhancer engagement and regulation by FOXA1.

FOXA1 is a transcription factor capable to bind silenced chromatin to direct context-dependent cell fate conversion. Here, we demonstrate that a compact palindromic DNA element (termed 'DIV' for its diverging half-sites) induces the homodimerization of FOXA1 with strongly positive cooperativity. Alternative structural models are consistent with either an indirect DNA-mediated cooperativity or a direct protein-protein interaction. The cooperative homodimer formation is strictly constrained by precise half-site spacing. Re-analysis of chromatin immunoprecipitation sequencing data indicates that the DIV is effectively targeted by FOXA1 in the context of chromatin. Reporter assays show that FOXA1-dependent transcriptional activity declines when homodimeric binding is disrupted. In response to phosphatidylinositol-3 kinase inhibition DIV sites pre-bound by FOXA1 such as at the PVT1/MYC locus exhibit a strong increase in accessibility suggesting a role of the DIV configuration in the chromatin closed-open dynamics. Moreover, several disease-associated single nucleotide polymorphisms map to DIV elements and show allelic differences in FOXA1 homodimerization, reporter gene expression and are annotated as quantitative trait loci. This includes the rs541455835 variant at the MAPT locus encoding the Tau protein associated with Parkinson's disease. Collectively, the DIV guides chromatin engagement and regulation by FOXA1 and its perturbation could be linked to disease etiologies.

Molecular basis for the genome engagement by Sox proteins.

The Sox transcription factor family consists of 20 members in the human genome. Many of them are key determinants of cellular identities and possess the capacity to reprogram cell fates by pioneering the epigenetic remodeling of the genome. This activity is intimately tied to their ability to specifically bind and bend DNA alone or with other proteins. Here we discuss our current knowledge on how Sox transcription factors such as Sox2, Sox17, Sox18 and Sox9 'read' the genome to find and regulate their target genes and highlight the roles of partner factors including Pax6, Nanog, Oct4 and Brn2. We integrate insights from structural and biochemical studies as well as high-throughput assays to probe DNA specificity in vitro as well as in cells and tissues.

Changing POU dimerization preferences converts Oct6 into a pluripotency inducer.

The transcription factor Oct4 is a core component of molecular cocktails inducing pluripotent stem cells (iPSCs), while other members of the POU family cannot replace Oct4 with comparable efficiency. Rather, group III POU factors such as Oct6 induce neural lineages. Here, we sought to identify molecular features determining the differential DNA-binding and reprogramming activity of Oct4 and Oct6. In enhancers of pluripotency genes, Oct4 cooperates with Sox2 on heterodimeric SoxOct elements. By re-analyzing ChIP-Seq data and performing dimerization assays, we found that Oct6 homodimerizes on palindromic OctOct more cooperatively and more stably than Oct4. Using structural and biochemical analyses, we identified a single amino acid directing binding to the respective DNA elements. A change in this amino acid decreases the ability of Oct4 to generate iPSCs, while the reverse mutation in Oct6 does not augment its reprogramming activity. Yet, with two additional amino acid exchanges, Oct6 acquires the ability to generate iPSCs and maintain pluripotency. Together, we demonstrate that cell type-specific POU factor function is determined by select residues that affect DNA-dependent dimerization.

Diversity among POU transcription factors in chromatin recognition and cell fate reprogramming.

The POU (Pit-Oct-Unc) protein family is an evolutionary ancient group of transcription factors (TFs) that bind specific DNA sequences to direct gene expression programs. The fundamental importance of POU TFs to orchestrate embryonic development and to direct cellular fate decisions is well established, but the molecular basis for this activity is insufficiently understood. POU TFs possess a bipartite 'two-in-one' DNA binding domain consisting of two independently folding structural units connected by a poorly conserved and flexible linker. Therefore, they represent a paradigmatic example to study the molecular basis for the functional versatility of TFs. Their modular architecture endows POU TFs with the capacity to accommodate alternative composite DNA sequences by adopting different quaternary structures. Moreover, associations with partner proteins crucially influence the selection of their DNA binding sites. The plentitude of DNA binding modes confers the ability to POU TFs to regulate distinct genes in the context of different cellular environments. Likewise, different binding modes of POU proteins to DNA could trigger alternative regulatory responses in the context of different genomic locations of the same cell. Prominent POU TFs such as Oct4, Brn2, Oct6 and Brn4 are not only essential regulators of development but have also been successfully employed to reprogram somatic cells to pluripotency and neural lineages. Here we review biochemical, structural, genomic and cellular reprogramming studies to examine how the ability of POU TFs to select regulatory DNA, alone or with partner factors, is tied to their capacity to epigenetically remodel chromatin and drive specific regulatory programs that give cells their identities.

Identification of a Wells-Dawson polyoxometalate-based AP-2γ inhibitor with pro-apoptotic activity.

AP-2 gamma (AP-2γ) is a transcription factor that plays pivotal roles in breast cancer biology. To search for small molecule inhibitors of AP-2γ, we performed a high-throughput fluorescence anisotropy screen and identified a polyoxometalate compound with Wells-Dawson structure K6[P2Mo18O62] (Dawson-POM) that blocks the DNA-binding activity of AP-2γ. We showed that this blocking activity is due to the direct binding of Dawson-POM to AP-2γ. We also provided evidence to show that Dawson-POM decreases AP-2γ-dependent transcription similar to silencing the gene. Finally, we demonstrated that Dawson-POM contains anti-proliferative and pro-apoptotic effects in breast cancer cells. In summary, we identified the first small molecule inhibitor of AP-2γ and showed Dawson-POM-mediated inhibition of AP-2γ as a potential avenue for cancer therapy.

Quantitative profiling of selective Sox/POU pairing on hundreds of sequences in parallel by Coop-seq.

Cooperative binding of transcription factors is known to be important in the regulation of gene expression programs conferring cellular identities. However, current methods to measure cooperativity parameters have been laborious and therefore limited to studying only a few sequence variants at a time. We developed Coop-seq (cooperativity by sequencing) that is capable of efficiently and accurately determining the cooperativity parameters for hundreds of different DNA sequences in a single experiment. We apply Coop-seq to 12 dimer pairs from the Sox and POU families of transcription factors using 324 unique sequences with changed half-site orientation, altered spacing and discrete randomization within the binding elements. The study reveals specific dimerization profiles of different Sox factors with Oct4. By contrast, Oct4 and the three neural class III POU factors Brn2, Brn4 and Oct6 assemble with Sox2 in a surprisingly indistinguishable manner. Two novel half-site configurations can support functional Sox/Oct dimerization in addition to known composite motifs. Moreover, Coop-seq uncovers a nucleotide switch within the POU half-site when spacing is altered, which is mirrored in genomic loci bound by Sox2/Oct4 complexes.

Structure and decoy-mediated inhibition of the SOX18/Prox1-DNA interaction.

The transcription factor (TF) SOX18 drives lymphatic vessel development in both embryogenesis and tumour-induced neo-lymphangiogenesis. Genetic disruption of Sox18 in a mouse model protects from tumour metastasis and established the SOX18 protein as a molecular target. Here, we report the crystal structure of the SOX18 DNA binding high-mobility group (HMG) box bound to a DNA element regulating Prox1 transcription. The crystals diffracted to 1.75Å presenting the highest resolution structure of a SOX/DNA complex presently available revealing water structure, structural adjustments at the DNA contact interface and non-canonical conformations of the DNA backbone. To explore alternatives to challenging small molecule approaches for targeting the DNA-binding activity of SOX18, we designed a set of five decoys based on modified Prox1-DNA. Four decoys potently inhibited DNA binding of SOX18 in vitro and did not interact with non-SOX TFs. Serum stability, nuclease resistance and thermal denaturation assays demonstrated that a decoy circularized with a hexaethylene glycol linker and terminal phosphorothioate modifications is most stable. This SOX decoy also interfered with the expression of a luciferase reporter under control of a SOX18-dependent VCAM1 promoter in COS7 cells. Collectively, we propose SOX decoys as potential strategy for inhibiting SOX18 activity to disrupt tumour-induced neo-lymphangiogenesis.