Highly cooperative chimeric super-SOX induces naive pluripotency across species.

Our understanding of pluripotency remains limited: iPSC generation has only been established for a few model species, pluripotent stem cell lines exhibit inconsistent developmental potential, and germline transmission has only been demonstrated for mice and rats. By swapping structural elements between Sox2 and Sox17, we built a chimeric super-SOX factor, Sox2-17, that enhanced iPSC generation in five tested species: mouse, human, cynomolgus monkey, cow, and pig. A swap of alanine to valine at the interface between Sox2 and Oct4 delivered a gain of function by stabilizing Sox2/Oct4 dimerization on DNA, enabling generation of high-quality OSKM iPSCs capable of supporting the development of healthy all-iPSC mice. Sox2/Oct4 dimerization emerged as the core driver of naive pluripotency with its levels diminished upon priming. Transient overexpression of the SK cocktail (Sox+Klf4) restored the dimerization and boosted the developmental potential of pluripotent stem cells across species, providing a universal method for naive reset in mammals.

An engineered Sox17 induces somatic to neural stem cell fate transitions independently from pluripotency reprogramming.

Advanced strategies to interconvert cell types provide promising avenues to model cellular pathologies and to develop therapies for neurological disorders. Yet, methods to directly transdifferentiate somatic cells into multipotent induced neural stem cells (iNSCs) are slow and inefficient, and it is unclear whether cells pass through a pluripotent state with full epigenetic reset. We report iNSC reprogramming from embryonic and aged mouse fibroblasts as well as from human blood using an engineered Sox17 (eSox17FNV). eSox17FNV efficiently drives iNSC reprogramming while Sox2 or Sox17 fail. eSox17FNV acquires the capacity to bind different protein partners on regulatory DNA to scan the genome more efficiently and has a more potent transactivation domain than Sox2. Lineage tracing and time-resolved transcriptomics show that emerging iNSCs do not transit through a pluripotent state. Our work distinguishes lineage from pluripotency reprogramming with the potential to generate more authentic cell models for aging-associated neurodegenerative diseases.

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.

Restricting epigenetic activity promotes the reprogramming of transformed cells to pluripotency in a line-specific manner.

Somatic cell reprogramming and oncogenic transformation share surprisingly similar features, yet transformed cells are resistant to reprogramming. Epigenetic barriers must block transformed cells from reprogramming, but the nature of those barriers is unclear. In this study, we generated a systematic panel of transformed mouse embryonic fibroblasts (MEFs) using oncogenic transgenes and discovered transformed cell lines compatible with reprogramming when transfected with Oct4/Sox2/Klf4/Myc. By comparing the reprogramming-capable and incapable transformed lines we identified multiple stages of failure in the reprogramming process. Some transformed lines failed at an early stage, whilst other lines seemed to progress through a conventional reprogramming process. Finally, we show that MEK inhibition overcomes one critical reprogramming barrier by indirectly suppressing a hyperacetylated active epigenetic state. This study reveals that diverse epigenetic barriers underly resistance to reprogramming of transformed cells.

The chemical reprogramming of unipotent adult germ cells towards authentic pluripotency and de novo establishment of imprinting.

Although somatic cells can be reprogrammed to pluripotent stem cells (PSCs) with pure chemicals, authentic pluripotency of chemically induced pluripotent stem cells (CiPSCs) has never been achieved through tetraploid complementation assay. Spontaneous reprogramming of spermatogonial stem cells (SSCs) was another non-transgenic way to obtain PSCs, but this process lacks mechanistic explanation. Here, we reconstructed the trajectory of mouse SSC reprogramming and developed a five-chemical combination, boosting the reprogramming efficiency by nearly 80- to 100-folds. More importantly, chemical induced germline-derived PSCs (5C-gPSCs), but not gPSCs and chemical induced pluripotent stem cells, had authentic pluripotency, as determined by tetraploid complementation. Mechanistically, SSCs traversed through an inverted pathway of in vivo germ cell development, exhibiting the expression signatures and DNA methylation dynamics from spermatogonia to primordial germ cells and further to epiblasts. Besides, SSC-specific imprinting control regions switched from biallelic methylated states to monoallelic methylated states by imprinting demethylation and then re-methylation on one of the two alleles in 5C-gPSCs, which was apparently distinct with the imprinting reprogramming in vivo as DNA methylation simultaneously occurred on both alleles. Our work sheds light on the unique regulatory network underpinning SSC reprogramming, providing insights to understand generic mechanisms for cell-fate decision and epigenetic-related disorders in regenerative medicine.

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.

Evolutionary Landscape of SOX Genes to Inform Genotype-to-Phenotype Relationships.

The SOX transcription factor family is pivotal in controlling aspects of development. To identify genotype-phenotype relationships of SOX proteins, we performed a non-biased study of SOX using 1890 open-reading frame and 6667 amino acid sequences in combination with structural dynamics to interpret 3999 gnomAD, 485 ClinVar, 1174 Geno2MP, and 4313 COSMIC human variants. We identified, within the HMG (High Mobility Group)- box, twenty-seven amino acids with changes in multiple SOX proteins annotated to clinical pathologies. These sites were screened through Geno2MP medical phenotypes, revealing novel SOX15 R104G associated with musculature abnormality and SOX8 R159G with intellectual disability. Within gnomAD, SOX18 E137K (rs201931544), found within the HMG box of ~0.8% of Latinx individuals, is associated with seizures and neurological complications, potentially through blood-brain barrier alterations. A total of 56 highly conserved variants were found at sites outside the HMG-box, including several within the SOX2 HMG-box-flanking region with neurological associations, several in the SOX9 dimerization region associated with Campomelic Dysplasia, SOX14 K88R (rs199932938) flanking the HMG box associated with cardiovascular complications within European populations, and SOX7 A379V (rs143587868) within an SOXF conserved far C-terminal domain heterozygous in 0.716% of African individuals with associated eye phenotypes. This SOX data compilation builds a robust genotype-to-phenotype association for a gene family through more robust ortholog data integration.

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.

Expanding DSD Phenotypes Associated with Variants in the DEAH-Box RNA Helicase DHX37.

Missense variants in the RNA-helicase DHX37 are associated with either 46,XY gonadal dysgenesis or 46,XY testicular regression syndrome (TRS). DHX37 is required for ribosome biogenesis, and this subgroup of XY DSD is a new human ribosomopathy. In a cohort of 140 individuals with 46,XY DSD, we identified 7 children with either 46,XY complete gonadal dysgenesis or 46,XY TRS carrying rare or novel DHX37 variants. A novel p.R390H variant within the RecA1 domain was identified in a girl with complete gonadal dysgenesis. A paternally inherited p.R487H variant, previously associated with a recessive congenital developmental syndrome, was carried by a boy with a syndromic form of 46,XY DSD. His phenotype may be explained in part by a novel homozygous loss-of-function variant in the NGLY1 gene, which causes a congenital disorder of deglycosylation. Remarkably, a homozygous p.T477H variant was identified in a boy with TRS. His fertile father had unilateral testicular regression with typical male genital development. This expands the DSD phenotypes associated with DHX37. Structural analysis of all variants predicted deleterious effects on helicase function. Similar to all other known ribosomopathies, the mechanism of pathogenesis is unknown.

The chromatin accessibility landscape reveals distinct transcriptional regulation in the induction of human primordial germ cell-like cells from pluripotent stem cells.

In vitro induction of human primordial germ cell-like cells (hPGCLCs) provides an ideal platform to recapitulate hPGC development. However, the detailed molecular mechanisms regulating the induction of hPGCLCs remain largely uncharacterized. Here, we profiled the chromatin accessibility and transcriptome dynamics throughout the process of hPGCLC induction. Genetic ablation of SOX15 indicated the crucial roles of SOX15 in the maintenance of hPGCLCs. Mechanistically, SOX15 exerted its roles via suppressing somatic gene expression and sustaining latent pluripotency. Notably, ETV5, a downstream regulator of SOX15, was also uncovered to be essential for hPGCLC maintenance. Finally, a stepwise switch of OCT4/SOX2, OCT4/SOX17, and OCT4/SOX15 binding motifs were found to be enriched in closed-to-open regions of human embryonic stem cells, and early- and late-stage hPGCLCs, respectively. Collectively, our data characterized the chromatin accessibility and transcriptome landscapes throughout hPGCLC induction and defined the SOX15-mediated regulatory networks underlying this process.

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.

JMJD3 acts in tandem with KLF4 to facilitate reprogramming to pluripotency.

The interplay between the Yamanaka factors (OCT4, SOX2, KLF4 and c-MYC) and transcriptional/epigenetic co-regulators in somatic cell reprogramming is incompletely understood. Here, we demonstrate that the histone H3 lysine 27 trimethylation (H3K27me3) demethylase JMJD3 plays conflicting roles in mouse reprogramming. On one side, JMJD3 induces the pro-senescence factor Ink4a and degrades the pluripotency regulator PHF20 in a reprogramming factor-independent manner. On the other side, JMJD3 is specifically recruited by KLF4 to reduce H3K27me3 at both enhancers and promoters of epithelial and pluripotency genes. JMJD3 also promotes enhancer-promoter looping through the cohesin loading factor NIPBL and ultimately transcriptional elongation. This competition of forces can be shifted towards improved reprogramming by using early passage fibroblasts or boosting JMJD3's catalytic activity with vitamin C. Our work, thus, establishes a multifaceted role for JMJD3, placing it as a key partner of KLF4 and a scaffold that assists chromatin interactions and activates gene transcription.

FOXA1 Mutations Reveal Distinct Chromatin Profiles and Influence Therapeutic Response in Breast Cancer.

Mutations in the pioneer transcription factor FOXA1 are a hallmark of estrogen receptor-positive (ER+) breast cancers. Examining FOXA1 in ∼5,000 breast cancer patients identifies several hotspot mutations in the Wing2 region and a breast cancer-specific mutation SY242CS, located in the third ß strand. Using a clinico-genomically curated cohort, together with breast cancer models, we find that FOXA1 mutations associate with a lower response to aromatase inhibitors. Mechanistically, Wing2 mutations display increased chromatin binding at ER loci upon estrogen stimulation, and an enhanced ER-mediated transcription without changes in chromatin accessibility. In contrast, SY242CS shows neomorphic properties that include the ability to open distinct chromatin regions and activate an alternative cistrome and transcriptome. Structural modeling predicts that SY242CS confers a conformational change that mediates stable binding to a non-canonical DNA motif. Taken together, our results provide insights into how FOXA1 mutations perturb its function to dictate cancer progression and therapeutic response.

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.

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.