Oct4 redox sensitivity potentiates reprogramming and differentiation.

The transcription factor Oct4/Pou5f1 is a component of the regulatory circuitry governing pluripotency and is widely used to induce pluripotency from somatic cells. Here we used domain swapping and mutagenesis to study Oct4's reprogramming ability, identifying a redox-sensitive DNA binding domain, cysteine residue (Cys48), as a key determinant of reprogramming and differentiation. Oct4 Cys48 sensitizes the protein to oxidative inhibition of DNA binding activity and promotes oxidation-mediated protein ubiquitylation. Pou5f1 C48S point mutation has little effect on undifferentiated embryonic stem cells (ESCs) but upon retinoic acid (RA) treatment causes retention of Oct4 expression, deregulated gene expression, and aberrant differentiation. Pou5f1 C48S ESCs also form less differentiated teratomas and contribute poorly to adult somatic tissues. Finally, we describe Pou5f1 C48S (Janky) mice, which in the homozygous condition are severely developmentally restricted after E4.5. Rare animals bypassing this restriction appear normal at birth but are sterile. Collectively, these findings uncover a novel Oct4 redox mechanism involved in both entry into and exit from pluripotency.

Single-cell analysis of the developing human ovary defines distinct insights into ovarian somatic and germline progenitors.

Formation of either an ovary or a testis during human embryonic life is one of the most important sex-specific events leading to the emergence of secondary sexual characteristics and sex assignment of babies at birth. Our study focused on the sex-specific and sex-indifferent characteristics of the prenatal ovarian stromal cells, cortical cords, and germline, with the discovery that the ovarian mesenchymal cells of the stroma are transcriptionally indistinguishable from the mesenchymal cells of the testicular interstitium. We found that first-wave pre-granulosa cells emerge at week 7 from early supporting gonadal cells with stromal identity and are spatially defined by KRT19 levels. We also identified rare transient state f0 spermatogonia cells within the ovarian cords between weeks 10 and 16. Taken together, our work illustrates a unique plasticity of the embryonic ovary during human development.

Chromatin Reprogramming of In Vitro Fertilized and Somatic Cell Nuclear Transfer Bovine Embryos During Embryonic Genome Activation.

Reprogramming of the gamete into a developmentally competent embryo identity is a fundamental aspect of preimplantation development. One of the most important processes of this reprogramming is the transcriptional awakening during embryonic genome activation (EGA), which robustly occurs in fertilized embryos but is defective in most somatic cell nuclear transfer (SCNT) embryos. However, little is known about the genome-wide underlying chromatin landscape during EGA in SCNT embryos and how it differs from a fertilized embryo. By profiling open chromatin genome-wide in both types of bovine embryos, we find that SCNT embryos fail to reprogram a subset of the EGA gene targets that are normally activated in fertilized embryos. Importantly, a small number of transcription factor (TF) motifs explain most chromatin regions that fail to open in SCNT embryos suggesting that over-expression of a limited number of TFs may provide more robust reprogramming. One such TF, the zygotically-expressed bovine gene DUXC which is a homologue of EGA factors DUX/DUX4 in mouse/human, is alone capable of activating ∻84% of all EGA transcripts that fail to activate normally in SCNT embryos. Additionally, single-cell chromatin profiling revealed low intra-embryo heterogeneity but high inter-embryo heterogeneity in SCNT embryos and an uncoupling of cell division and open chromatin reprogramming during EGA. Surprisingly, our data also indicate that transcriptional defects may arise downstream of promoter chromatin opening in SCNT embryos, suggesting additional mechanistic insights into how and why transcription at EGA is dysregulated. We anticipate that our work will lead to altered SCNT protocols to increase the developmental competency of bovine SCNT embryos.

Multiple repeat regions within mouse DUX recruit chromatin regulators to facilitate an embryonic gene expression program.

The embryonic transcription factor DUX regulates chromatin opening and gene expression in totipotent cleavage-stage mouse embryos, and its expression in embryonic stem cells promotes their conversion to 2-cell embryo-like cells (2CLCs) with extraembryonic potential. However, little is known regarding which domains within mouse DUX interact with particular chromatin and transcription regulators. Here, we reveal that the C-terminus of mouse DUX contains five uncharacterized ~100 amino acid (aa) repeats followed by an acidic 14 amino acid tail. Unexpectedly, structure-function approaches classify two repeats as 'active' and three as 'inactive' in cleavage/2CLC transcription program enhancement, with differences narrowed to a key 6 amino acid section. Our proximity dependent biotin ligation (BioID) approach identified factors selectively associated with active DUX repeat derivatives (including the 14aa 'tail'), including transcription and chromatin factors such as SWI/SNF (BAF) complex, as well as nucleolar factors that have been previously implicated in regulating the Dux locus. Finally, our mechanistic studies reveal cooperativity between DUX active repeats and the acidic tail in cofactor recruitment, DUX target opening, and transcription. Taken together, we provide several new insights into DUX structure-function, and mechanisms of chromatin and gene regulation.

Single-cell analysis of human testis aging and correlation with elevated body mass index.

Aging men display reduced reproductive health; however, testis aging is poorly understood at the molecular and genomic levels. Here, we utilized single-cell RNA-seq to profile over 44,000 cells from both young and older men and examined age-related changes in germline development and in the testicular somatic cells. Age-related changes in spermatogonial stem cells appeared modest, whereas age-related dysregulation of spermatogenesis and somatic cells ranged from moderate to severe. Altered pathways included signaling and inflammation in multiple cell types, metabolic signaling in Sertoli cells, hedgehog signaling and testosterone production in Leydig cells, cell death and growth in testicular peritubular cells, and possible developmental regression in both Leydig and peritubular cells. Remarkably, the extent of dysregulation correlated with body mass index in older but not in younger men. Collectively, we reveal candidate molecular mechanisms underlying the complex testicular changes conferred by aging and their possible exacerbation by concurrent chronic conditions such as obesity.

GFI1 Cooperates with IKZF1/IKAROS to Activate Gene Expression in T-cell Acute Lymphoblastic Leukemia.

Growth factor independence-1 (GFI1) is a transcriptional repressor and master regulator of normal and malignant hematopoiesis. Repression by GFI1 is attributable to recruitment of LSD1-containing protein complexes via its SNAG domain. However, the full complement of GFI1 partners in transcriptional control is not known. We show that in T-acute lymphoblastic leukemia (ALL) cells, GFI1 and IKAROS are transcriptional partners that co-occupy regulatory regions of hallmark T-cell development genes. Transcriptional profiling reveals a subset of genes directly transactivated through the GFI1-IKAROS partnership. Among these is NOTCH3, a key factor in T-ALL pathogenesis. Surprisingly, NOTCH3 expression by GFI1 and IKAROS requires the GFI1 SNAG domain but occurs independent of SNAG-LSD1 binding. GFI1 variants deficient in LSD1 binding fail to activate NOTCH3, but conversely, small molecules that disrupt the SNAG-LSD1 interaction while leaving the SNAG primary structure intact stimulate NOTCH3 expression. These results identify a noncanonical transcriptional control mechanism in T-ALL which supports GFI1-mediated transactivation in partnership with IKAROS and suggest competition between LSD1-containing repressive complexes and others favoring transactivation. IMPLICATIONS: Combinatorial diversity and cooperation between DNA binding proteins and complexes assembled by them can direct context-dependent transcriptional outputs to control cell fate and may offer new insights for therapeutic targeting in cancer.

Establishment of developmental gene silencing by ordered polycomb complex recruitment in early zebrafish embryos.

Vertebrate embryos achieve developmental competency during zygotic genome activation (ZGA) by establishing chromatin states that silence yet poise developmental genes for subsequent lineage-specific activation. Here, we reveal the order of chromatin states in establishing developmental gene poising in preZGA zebrafish embryos. Poising is established at promoters and enhancers that initially contain open/permissive chromatin with 'Placeholder' nucleosomes (bearing H2A.Z, H3K4me1, and H3K27ac), and DNA hypomethylation. Silencing is initiated by the recruitment of polycomb repressive complex 1 (PRC1), and H2Aub1 deposition by catalytic Rnf2 during preZGA and ZGA stages. During postZGA, H2Aub1 enables Aebp2-containing PRC2 recruitment and H3K27me3 deposition. Notably, preventing H2Aub1 (via Rnf2 inhibition) eliminates recruitment of Aebp2-PRC2 and H3K27me3, and elicits transcriptional upregulation of certain developmental genes during ZGA. However, upregulation is independent of H3K27me3 - establishing H2Aub1 as the critical silencing modification at ZGA. Taken together, we reveal the logic and mechanism for establishing poised/silent developmental genes in early vertebrate embryos.

Dissecting mammalian spermatogenesis using spatial transcriptomics.

Single-cell RNA sequencing has revealed extensive molecular diversity in gene programs governing mammalian spermatogenesis but fails to delineate their dynamics in the native context of seminiferous tubules, the spatially confined functional units of spermatogenesis. Here, we use Slide-seq, a spatial transcriptomics technology, to generate an atlas that captures the spatial gene expression patterns at near-single-cell resolution in the mouse and human testis. Using Slide-seq data, we devise a computational framework that accurately localizes testicular cell types in individual seminiferous tubules. Unbiased analysis systematically identifies spatially patterned genes and gene programs. Combining Slide-seq with targeted in situ RNA sequencing, we demonstrate significant differences in the cellular compositions of spermatogonial microenvironment between mouse and human testes. Finally, a comparison of the spatial atlas generated from the wild-type and diabetic mouse testis reveals a disruption in the spatial cellular organization of seminiferous tubules as a potential mechanism of diabetes-induced male infertility.

p53 convergently activates Dux/DUX4 in embryonic stem cells and in facioscapulohumeral muscular dystrophy cell models.

In mammalian embryos, proper zygotic genome activation (ZGA) underlies totipotent development. Double homeobox (DUX)-family factors participate in ZGA, and mouse Dux is required for forming cultured two-cell (2C)-like cells. Remarkably, in mouse embryonic stem cells, Dux is activated by the tumor suppressor p53, and Dux expression promotes differentiation into expanded-fate cell types. Long-read sequencing and assembly of the mouse Dux locus reveals its complex chromatin regulation including putative positive and negative feedback loops. We show that the p53-DUX/DUX4 regulatory axis is conserved in humans. Furthermore, we demonstrate that cells derived from patients with facioscapulohumeral muscular dystrophy (FSHD) activate human DUX4 during p53 signaling via a p53-binding site in a primate-specific subtelomeric long terminal repeat (LTR)10C element. In summary, our work shows that p53 activation convergently evolved to couple p53 to Dux/DUX4 activation in embryonic stem cells, embryos and cells from patients with FSHD, potentially uniting the developmental and disease regulation of DUX-family factors and identifying evidence-based therapeutic opportunities for FSHD.

A Role for SMARCB1 in Synovial Sarcomagenesis Reveals That SS18-SSX Induces Canonical BAF Destruction.

Reduced protein levels of SMARCB1 (also known as BAF47, INI1, SNF5) have long been observed in synovial sarcoma. Here, we show that combined Smarcb1 genetic loss with SS18-SSX expression in mice synergized to produce aggressive tumors with histomorphology, transcriptomes, and genome-wide BAF-family complex distributions distinct from SS18-SSX alone, indicating a defining role for SMARCB1 in synovial sarcoma. Smarcb1 silencing alone in mesenchyme modeled epithelioid sarcomagenesis. In mouse and human synovial sarcoma cells, SMARCB1 was identified within PBAF and canonical BAF (CBAF) complexes, coincorporated with SS18-SSX in the latter. Recombinant expression of CBAF components in human cells reconstituted CBAF subcomplexes that contained equal levels of SMARCB1 regardless of SS18 or SS18-SSX inclusion. In vivo, SS18-SSX expression led to whole-complex CBAF degradation, rendering increases in the relative prevalence of other BAF-family subtypes, PBAF and GBAF complexes, over time. Thus, SS18-SSX alters BAF subtypes levels/balance and genome distribution, driving synovial sarcomagenesis. SIGNIFICANCE: The protein level of BAF component SMARCB1 is reduced in synovial sarcoma but plays a defining role, incorporating into PBAF and SS18-SSX-containing canonical BAF complexes. Reduced levels of SMARCB1 derive from whole-complex degradation of canonical BAF driven by SS18-SSX, with relative increases in the abundance of other BAF-family subtypes.See related commentary by Maxwell and Hargreaves, p. 2375.This article is highlighted in the In This Issue feature, p. 2355.

Chromatin architecture transitions from zebrafish sperm through early embryogenesis.

Chromatin architecture mapping in 3D formats has increased our understanding of how regulatory sequences and gene expression are connected and regulated in a genome. The 3D chromatin genome shows extensive remodeling during embryonic development, and although the cleavage-stage embryos of most species lack structure before zygotic genome activation (pre-ZGA), zebrafish has been reported to have structure. Here, we aimed to determine the chromosomal architecture in paternal/sperm zebrafish gamete cells to discern whether it either resembles or informs early pre-ZGA zebrafish embryo chromatin architecture. First, we assessed the higher-order architecture through advanced low-cell in situ Hi-C. The structure of zebrafish sperm, packaged by histones, lacks topological associated domains and instead displays "hinge-like" domains of ∼150 kb that repeat every 1-2 Mbs, suggesting a condensed repeating structure resembling mitotic chromosomes. The pre-ZGA embryos lacked chromosomal structure, in contrast to prior work, and only developed structure post-ZGA. During post-ZGA, we find chromatin architecture beginning to form at small contact domains of a median length of ∼90 kb. These small contact domains are established at enhancers, including super-enhancers, and chemical inhibition of Ep300a (p300) and Crebbpa (CBP) activity, lowering histone H3K27ac, but not transcription inhibition, diminishes these contacts. Together, this study reveals hinge-like domains in histone-packaged zebrafish sperm chromatin and determines that the initial formation of high-order chromatin architecture in zebrafish embryos occurs after ZGA primarily at enhancers bearing high H3K27ac.

CTCF looping is established during gastrulation in medaka embryos.

Chromatin looping plays an important role in genome regulation. However, because ChIP-seq and loop-resolution Hi-C (DNA-DNA proximity ligation) are extremely challenging in mammalian early embryos, the developmental stage at which cohesin-mediated loops form remains unknown. Here, we study early development in medaka (the Japanese killifish, Oryzias latipes) at 12 time points before, during, and after gastrulation (the onset of cell differentiation) and characterize transcription, protein binding, and genome architecture. We find that gastrulation is associated with drastic changes in genome architecture, including the formation of the first loops between sites bound by the insulator protein CTCF and a large increase in the size of contact domains. In contrast, the binding of the CTCF is fixed throughout embryogenesis. Loops form long after genome-wide transcriptional activation, and long after domain formation seen in mouse embryos. These results suggest that, although loops may play a role in differentiation, they are not required for zygotic transcription. When we repeated our experiments in zebrafish, loops did not emerge until gastrulation, that is, well after zygotic genome activation. We observe that loop positions are highly conserved in synteny blocks of medaka and zebrafish, indicating that the 3D genome architecture has been maintained for >110-200 million years of evolution.

PANDORA-seq expands the repertoire of regulatory small RNAs by overcoming RNA modifications.

Although high-throughput RNA sequencing (RNA-seq) has greatly advanced small non-coding RNA (sncRNA) discovery, the currently widely used complementary DNA library construction protocol generates biased sequencing results. This is partially due to RNA modifications that interfere with adapter ligation and reverse transcription processes, which prevent the detection of sncRNAs bearing these modifications. Here, we present PANDORA-seq (panoramic RNA display by overcoming RNA modification aborted sequencing), employing a combinatorial enzymatic treatment to remove key RNA modifications that block adapter ligation and reverse transcription. PANDORA-seq identified abundant modified sncRNAs-mostly transfer RNA-derived small RNAs (tsRNAs) and ribosomal RNA-derived small RNAs (rsRNAs)-that were previously undetected, exhibiting tissue-specific expression across mouse brain, liver, spleen and sperm, as well as cell-specific expression across embryonic stem cells (ESCs) and HeLa cells. Using PANDORA-seq, we revealed unprecedented landscapes of microRNA, tsRNA and rsRNA dynamics during the generation of induced pluripotent stem cells. Importantly, tsRNAs and rsRNAs that are downregulated during somatic cell reprogramming impact cellular translation in ESCs, suggesting a role in lineage differentiation.

Author Correction: PANDORA-seq expands the repertoire of regulatory small RNAs by overcoming RNA modifications.

A Correction to this paper has been published: https://doi.org/10.1038/s41556-021-00687-w.

Single-cell analysis of the developing human testis reveals somatic niche cell specification and fetal germline stem cell establishment.

Human testis development in prenatal life involves complex changes in germline and somatic cell identity. To better understand, we profiled and analyzed ∼32,500 single-cell transcriptomes of testicular cells from embryonic, fetal, and infant stages. Our data show that at 6-7 weeks postfertilization, as the testicular cords are established, the Sertoli and interstitial cells originate from a common heterogeneous progenitor pool, which then resolves into fetal Sertoli cells (expressing tube-forming genes) or interstitial cells (including Leydig-lineage cells expressing steroidogenesis genes). Almost 10 weeks later, beginning at 14-16 weeks postfertilization, the male primordial germ cells exit mitosis, downregulate pluripotent transcription factors, and transition into cells that strongly resemble the state 0 spermatogonia originally defined in the infant and adult testes. Therefore, we called these fetal spermatogonia "state f0." Overall, we reveal multiple insights into the coordinated and temporal development of the embryonic, fetal, and postnatal male germline together with the somatic niche.

Cancer-Associated Gain-of-Function Mutations Activate a SWI/SNF-Family Regulatory Hub.

SWI/SNF-family remodelers (BAF/PBAF in mammals) are essential chromatin regulators, and mutations in human BAF/PBAF components are associated with ∼20% of cancers. Cancer-associated missense mutations in human BRG1 (encoding the catalytic ATPase) have been characterized previously as conferring loss-of-function. Here, we show that cancer-associated missense mutations in BRG1, when placed into the orthologous Sth1 ATPase of the yeast RSC remodeler, separate into two categories: loss-of-function enzymes, or instead, gain-of-function enzymes that greatly improve DNA translocation efficiency and nucleosome remodeling in vitro. Our work identifies a structural "hub," formed by the association of several Sth1 domains, that regulates ATPase activity and DNA translocation efficiency. Remarkably, all gain-of-function cancer-associated mutations and all loss-of-function mutations physically localize to distinct adjacent regions in the hub, which specifically regulate and implement DNA translocation, respectively. In vivo, only gain-of-function cancer-associated mutations conferred precocious chromatin accessibility. Taken together, we provide a structure-function mechanistic basis for cancer-associated hyperactivity.

NRF2 loss recapitulates heritable impacts of paternal cigarette smoke exposure.

Paternal cigarette smoke (CS) exposure is associated with increased risk of behavioral disorders and cancer in offspring, but the mechanism has not been identified. Here we use mouse models to investigate mechanisms and impacts of paternal CS exposure. We demonstrate that CS exposure induces sperm DNAme changes that are partially corrected within 28 days of removal from CS exposure. Additionally, paternal smoking is associated with changes in prefrontal cortex DNAme and gene expression patterns in offspring. Remarkably, the epigenetic and transcriptional effects of CS exposure that we observed in wild type mice are partially recapitulated in Nrf2-/- mice and their offspring, independent of smoking status. Nrf2 is a central regulator of antioxidant gene transcription, and mice lacking Nrf2 consequently display elevated oxidative stress, suggesting that oxidative stress may underlie CS-induced heritable epigenetic changes. Importantly, paternal sperm DNAme changes do not overlap with DNAme changes measured in offspring prefrontal cortex, indicating that the observed DNAme changes in sperm are not directly inherited. Additionally, the changes in sperm DNAme associated with CS exposure were not observed in sperm of unexposed offspring, suggesting the effects are likely not maintained across multiple generations.

Specialization of the chromatin remodeler RSC to mobilize partially-unwrapped nucleosomes.

SWI/SNF-family chromatin remodeling complexes, such as S. cerevisiae RSC, slide and eject nucleosomes to regulate transcription. Within nucleosomes, stiff DNA sequences confer spontaneous partial unwrapping, prompting whether and how SWI/SNF-family remodelers are specialized to remodel partially-unwrapped nucleosomes. RSC1 and RSC2 are orthologs of mammalian PBRM1 (polybromo) which define two separate RSC sub-complexes. Remarkably, in vitro the Rsc1-containing complex remodels partially-unwrapped nucleosomes much better than does the Rsc2-containing complex. Moreover, a rsc1Δ mutation, but not rsc2Δ, is lethal with histone mutations that confer partial unwrapping. Rsc1/2 isoforms both cooperate with the DNA-binding proteins Rsc3/30 and the HMG protein, Hmo1, to remodel partially-unwrapped nucleosomes, but show differential reliance on these factors. Notably, genetic impairment of these factors strongly reduces the expression of genes with wide nucleosome-deficient regions (e.g., ribosomal protein genes), known to harbor partially-unwrapped nucleosomes. Taken together, Rsc1/2 isoforms are specialized through composition and interactions to manage and remodel partially-unwrapped nucleosomes.

Developmentally Programmed Tankyrase Activity Upregulates ß-Catenin and Licenses Progression of Embryonic Genome Activation.

Embryonic genome activation (EGA) is orchestrated by an intrinsic developmental program initiated during oocyte maturation with translation of stored maternal mRNAs. Here, we show that tankyrase, a poly(ADP-ribosyl) polymerase that regulates ß-catenin levels, undergoes programmed translation during oocyte maturation and serves an essential role in mouse EGA. Newly translated TNKS triggers proteasomal degradation of axin, reducing targeted destruction of ß-catenin and promoting ß-catenin-mediated transcription of target genes, including Myc. MYC mediates ribosomal RNA transcription in 2-cell embryos, supporting global protein synthesis. Suppression of tankyrase activity using knockdown or chemical inhibition causes loss of nuclear ß-catenin and global reductions in transcription and histone H3 acetylation. Chromatin and transcriptional profiling indicate that development arrests prior to the mid-2-cell stage, mediated in part by reductions in ß-catenin and MYC. These findings indicate that post-transcriptional regulation of tankyrase serves as a ligand-independent developmental mechanism for post-translational ß-catenin activation and is required to complete EGA.