Transcription factor binding dynamics during human ES cell differentiation

A primate-specific endogenous retroviral envelope protein sequesters SFRP2 to regulate human cardiomyocyte development

Endogenous retroviruses (ERVs) occupy a significant part of the human genome, with some encoding proteins that influence the immune system or regulate cell-cell fusion in early extra-embryonic development. However, whether ERV-derived proteins regulate somatic development is unknown. Here, we report a somatic developmental function for the primate-specific ERVH48-1 (SUPYN/Suppressyn). ERVH48-1 encodes a fragment of a viral envelope that is expressed during early embryonic development. Loss of ERVH48-1 led to impaired mesoderm and cardiomyocyte commitment and diverted cells to an ectoderm-like fate. Mechanistically, ERVH48-1 is localized to sub-cellular membrane compartments through a functional N-terminal signal peptide and binds to the WNT antagonist SFRP2 to promote its polyubiquitination and degradation, thus limiting SFRP2 secretion and blocking repression of WNT/β-catenin signaling. Knockdown of SFRP2 or expression of a chimeric SFRP2 with the ERVH48-1 signal peptide rescued cardiomyocyte differentiation. This study demonstrates how ERVH48-1 modulates WNT/β-catenin signaling and cell type commitment in somatic development.

TGF-β modulates cell fate in human ES cell-derived foregut endoderm by inhibiting Wnt and BMP signaling

Genetic differences between pluripotent stem cell lines cause variable activity of extracellular signaling pathways, limiting reproducibility of directed differentiation protocols. Here we used human embryonic stem cells (hESCs) to interrogate how exogenous factors modulate endogenous signaling events during specification of foregut endoderm lineages. We find that transforming growth factor β1 (TGF-β1) activates a putative human OTX2/LHX1 gene regulatory network which promotes anterior fate by antagonizing endogenous Wnt signaling. In contrast to Porcupine inhibition, TGF-β1 effects cannot be reversed by exogenous Wnt ligands, suggesting that induction of SHISA proteins and intracellular accumulation of Fzd receptors render TGF-β1-treated cells refractory to Wnt signaling. Subsequently, TGF-β1-mediated inhibition of BMP and Wnt signaling suppresses liver fate and promotes pancreas fate. Furthermore, combined TGF-β1 treatment and Wnt inhibition during pancreatic specification reproducibly and robustly enhance INSULIN+ cell yield across hESC lines. This modification of widely used differentiation protocols will enhance pancreatic β cell yield for cell-based therapeutic applications.

Uncovering the dynamics and consequences of RNA isoform changes during neuronal differentiation

Static gene expression programs have been extensively characterized in stem cells and mature human cells. However, the dynamics of RNA isoform changes upon cell-state-transitions during cell differentiation, the determinants and functional consequences have largely remained unclear. Here, we established an improved model for human neurogenesis in vitro that is amenable for systems-wide analyses of gene expression. Our multi-omics analysis reveals that the pronounced alterations in cell morphology correlate strongly with widespread changes in RNA isoform expression. Our approach identifies thousands of new RNA isoforms that are expressed at distinct differentiation stages. RNA isoforms mainly arise from exon skipping and the alternative usage of transcription start and polyadenylation sites during human neurogenesis. The transcript isoform changes can remodel the identity and functions of protein isoforms. Finally, our study identifies a set of RNA binding proteins as a potential determinant of differentiation stage-specific global isoform changes. This work supports the view of regulated isoform changes that underlie state-transitions during neurogenesis.

Quantitative and systematic NMR measurements of sequence-dependent A-T Hoogsteen dynamics uncovers unique conformational specificity in the DNA double helix

The propensities to form lowly-populated short-lived conformations of DNA could vary with sequence, providing an important source of sequence-specificity in biochemical reactions. However, comprehensively measuring how these dynamics vary with sequence is challenging. Using 1H CEST and 13C R 1 ρ NMR, we measured Watson-Crick to Hoogsteen dynamics for an A-T base pair in thirteen trinucleotide sequence contexts. The Hoogsteen population and exchange rate varied 4-fold and 16-fold, respectively, and were dependent on both the 3'- and 5'-neighbors but only weakly dependent on monovalent ion concentration (25 versus 100 mM NaCl) and pH (6.8 versus 8.0). Flexible TA and CA dinucleotide steps exhibited the highest Hoogsteen populations, and their kinetics rates strongly depended on the 3'-neighbor. In contrast, the stiffer AA and GA steps had the lowest Hoogsteen population, and their kinetics were weakly dependent on the 3'-neighbor. The Hoogsteen lifetime was especially short when G-C neighbors flanked the A-T base pair. The Hoogsteen dynamics had a distinct sequence-dependence compared to duplex stability and minor groove width. Thus, our results uncover a unique source of sequence-specificity hidden within the DNA double helix in the form of A-T Hoogsteen dynamics and establish the utility of 1H CEST to quantitively measure sequence-dependent DNA dynamics.

Possible Strategies to Reduce the Tumorigenic Risk of Reprogrammed Normal and Cancer Cells

The reprogramming of somatic cells to pluripotent stem cells has immense potential for use in regenerating or redeveloping tissues for transplantation, and the future application of this method is one of the most important research topics in regenerative medicine. These cells are generated from normal cells, adult stem cells, or neoplastic cancer cells. They express embryonic stem cell markers, such as OCT4, SOX2, and NANOG, and can differentiate into all tissue types in adults, both in vitro and in vivo. However, tumorigenicity, immunogenicity, and heterogeneity of cell populations may hamper the use of this method in medical therapeutics. The risk of cancer formation is dependent on mutations of these stemness genes during the transformation of pluripotent stem cells to cancer cells and on the alteration of the microenvironments of stem cell niches at genetic and epigenetic levels. Recent reports have shown that the generation of induced pluripotent stem cells (iPSCs) derived from human fibroblasts could be induced using chemicals, which is a safe, easy, and clinical-grade manufacturing strategy for modifying the cell fate of human cells required for regeneration therapies. This strategy is one of the future routes for the clinical application of reprogramming therapy. Therefore, this review highlights the recent progress in research focused on decreasing the tumorigenic risk of iPSCs or iPSC-derived organoids and increasing the safety of iPSC cell preparation and their application for therapeutic benefits.

NANOG controls testicular germ cell tumour stemness through regulation of MIR9-2

BACKGROUND

Testicular germ cell tumours (TGCTs) represent a clinical challenge; they are most prevalent in young individuals and are triggered by molecular mechanisms that are not fully understood. The origin of TGCTs can be traced back to primordial germ cells that fail to mature during embryonic development. These cells express high levels of pluripotency factors, including the transcription factor NANOG which is highly expressed in TGCTs. Gain or amplification of the NANOG locus is common in advanced tumours, suggesting a key role for this master regulator of pluripotency in TGCT stemness and malignancy.

METHODS

In this study, we analysed the expression of microRNAs (miRNAs) that are regulated by NANOG in TGCTs via integrated bioinformatic analyses of data from The Cancer Genome Atlas and NANOG chromatin immunoprecipitation in human embryonic stem cells. Through gain-of-function experiments, MIR9-2 was further investigated as a novel tumour suppressor regulated by NANOG. After transfection with MIR9-2 mimics, TGCT cells were analysed for cell proliferation, invasion, sensitivity to cisplatin, and gene expression signatures by RNA sequencing.

RESULTS

For the first time, we identified 86 miRNAs regulated by NANOG in TGCTs. Among these, 37 miRNAs were differentially expressed in NANOG-high tumours, and they clustered TGCTs according to their subtypes. Binding of NANOG within 2 kb upstream of the MIR9-2 locus was associated with a negative regulation. Low expression of MIR9-2 was associated with tumour progression and MIR9-2-5p was found to play a role in the control of tumour stemness. A gain of function of MIR9-2-5p was associated with reduced proliferation, invasion, and sensitivity to cisplatin in both embryonal carcinoma and seminoma tumours. MIR9-2-5p expression in TGCT cells significantly reduced the expression of genes regulating pluripotency and cell division, consistent with its functional effect on reducing cancer stemness.

CONCLUSIONS

This study provides new molecular insights into the role of NANOG as a key determinant of pluripotency in TGCTs through the regulation of MIR9-2-5p, a novel epigenetic modulator of cancer stemness. Our data also highlight the potential negative feedback mediated by MIR9-2-5p on NANOG expression, which could be exploited as a therapeutic strategy for the treatment of TGCTs.

Cis-regulatory interfaces reveal the molecular mechanisms underlying the notochord gene regulatory network of Ciona

Tissue-specific gene expression is fundamental in development and evolution, and is mediated by transcription factors and by the cis-regulatory regions (enhancers) that they control. Transcription factors and their respective tissue-specific enhancers are essential components of gene regulatory networks responsible for the development of tissues and organs. Although numerous transcription factors have been characterized from different organisms, the knowledge of the enhancers responsible for their tissue-specific expression remains fragmentary. Here we use Ciona to study the enhancers associated with ten transcription factors expressed in the notochord, an evolutionary hallmark of the chordate phylum. Our results illustrate how two evolutionarily conserved transcription factors, Brachyury and Foxa2, coordinate the deployment of other notochord transcription factors. The results of these detailed cis-regulatory analyses delineate a high-resolution view of the essential notochord gene regulatory network of Ciona, and provide a reference for studies of transcription factors, enhancers, and their roles in development, disease, and evolution.

TBXT dose sensitivity and the decoupling of nascent mesoderm specification from EMT progression in 2D human gastruloids

In the nascent mesoderm, TBXT expression must be precisely regulated to ensure that cells exit the primitive streak and pattern the anterior-posterior axis, but how varying dosage informs morphogenesis is not well understood. In this study, we define the transcriptional consequences of TBXT dosage reduction during early human gastrulation using human induced pluripotent stem cell models of gastrulation and mesoderm differentiation. Multi-omic single-nucleus RNA and single-nucleus ATAC sequencing of 2D gastruloids comprising wild-type, TBXT heterozygous or TBXT null human induced pluripotent stem cells reveal that varying TBXT dosage does not compromise the ability of a cell to differentiate into nascent mesoderm, but instead directly influences the temporal progression of the epithelial-to-mesenchymal transition with wild type transitioning first, followed by TBXT heterozygous and then TBXT null. By differentiating cells into nascent mesoderm in a monolayer format, we further illustrate that TBXT dosage directly impacts the persistence of junctional proteins and cell-cell adhesions. These results demonstrate that epithelial-to-mesenchymal transition progression can be decoupled from the acquisition of mesodermal identity in the early gastrula and shed light on the mechanisms underlying human embryogenesis.

Time-integrated BMP signaling determines fate in a stem cell model for early human development

How paracrine signals are interpreted to yield multiple cell fate decisions in a dynamic context during human development in vivo and in vitro remains poorly understood. Here we report an automated tracking method to follow signaling histories linked to cell fate in large numbers of human pluripotent stem cells (hPSCs). Using an unbiased statistical approach, we discover that measured BMP signaling history correlates strongly with fate in individual cells. We find that BMP response in hPSCs varies more strongly in the duration of signaling than the level. However, both the level and duration of signaling activity control cell fate choices only by changing the time integral. Therefore, signaling duration and level are interchangeable in this context. In a stem cell model for patterning of the human embryo, we show that signaling histories predict the fate pattern and that the integral model correctly predicts changes in cell fate domains when signaling is perturbed. Our data suggest that mechanistically, BMP signaling is integrated by SOX2.

HAND factors regulate cardiac lineage commitment and differentiation from human pluripotent stem cells

BACKGROUND

Transcription factors HAND1 and HAND2 (HAND1/2) play significant roles in cardiac organogenesis. Abnormal expression and deficiency of HAND1/2 result in severe cardiac defects. However, the function and mechanism of HAND1/2 in regulating human early cardiac lineage commitment and differentiation are still unclear.

METHODS

With NKX2.5eGFP H9 human embryonic stem cells (hESCs), we established single and double knockout cell lines for HAND1 and HAND2, respectively, whose cardiomyocyte differentiation efficiency could be monitored by assessing NKX2.5-eGFP+ cells with flow cytometry. The expression of specific markers for heart fields and cardiomyocyte subtypes was examined by quantitative PCR, western blot and immunofluorescence staining. Microelectrode array and whole-cell patch clamp were performed to determine the electrophysiological characteristics of differentiated cardiomyocytes. The transcriptomic changes of HAND knockout cells were revealed by RNA sequencing. The HAND1/2 target genes were identified and validated experimentally by integrating with HAND1/2 chromatin immunoprecipitation sequencing data.

RESULTS

Either HAND1 or HAND2 knockout did not affect the cardiomyocyte differentiation kinetics, whereas depletion of HAND1/2 resulted in delayed differentiation onset. HAND1 knockout biased cardiac mesoderm toward second heart field progenitors at the expense of first heart field progenitors, leading to increased expression of atrial and outflow tract cardiomyocyte markers, which was further confirmed by the appearance of atrial-like action potentials. By contrast, HAND2 knockout cardiomyocytes had reduced expression of atrial cardiomyocyte markers and displayed ventricular-like action potentials. HAND1/2-deficient hESCs were more inclined to second heart field lineage and its derived cardiomyocytes with atrial-like action potentials than HAND1 single knockout during differentiation. Further mechanistic investigations suggested TBX5 as one of the downstream targets of HAND1/2, whose overexpression partially restored the abnormal cardiomyocyte differentiation in HAND1/2-deficient hESCs.

CONCLUSIONS

HAND1/2 have specific and redundant roles in cardiac lineage commitment and differentiation. These findings not only reveal the essential function of HAND1/2 in cardiac organogenesis, but also provide important information on the pathogenesis of HAND1/2 deficiency-related congenital heart diseases, which could potentially lead to new therapeutic strategies.

CRISPR-assisted transcription activation by phase-separation proteins

The clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 system has been widely used for genome engineering and transcriptional regulation in many different organisms. Current CRISPR-activation (CRISPRa) platforms often require multiple components because of inefficient transcriptional activation. Here, we fused different phase-separation proteins to dCas9-VPR (dCas9-VP64-P65-RTA) and observed robust increases in transcriptional activation efficiency. Notably, human NUP98 (nucleoporin 98) and FUS (fused in sarcoma) IDR domains were best at enhancing dCas9-VPR activity, with dCas9-VPR-FUS IDR (VPRF) outperforming the other CRISPRa systems tested in this study in both activation efficiency and system simplicity. dCas9-VPRF overcomes the target strand bias and widens gRNA designing windows without affecting the off-target effect of dCas9-VPR. These findings demonstrate the feasibility of using phase-separation proteins to assist in the regulation of gene expression and support the broad appeal of the dCas9-VPRF system in basic and clinical applications.

A SNAI2/CTCF Interaction is Required for NOTCH1 Expression in Rhabdomyosarcoma

Rhabdomyosarcoma (RMS) is a pediatric malignancy of the muscle with characteristics of cells blocked in differentiation. NOTCH1 is an oncogene that promotes self-renewal and blocks differentiation in the fusion negative-RMS sub-type. However, how NOTCH1 expression is transcriptionally maintained in tumors is unknown. Analyses of SNAI2 and CTCF chromatin binding and HiC analyses revealed a conserved SNAI2/CTCF overlapping peak downstream of the NOTCH1 locus marking a sub-topologically associating domain (TAD) boundary. Deletion of the SNAI2-CTCF peak showed that it is essential for NOTCH1 expression and viability of FN-RMS cells. Reintroducing constitutively activated NOTCH1-ΔE in cells with the SNAI2-CTCF peak deleted restored cell-viability. Ablation of SNAI2 using CRISPR/Cas9 reagents resulted in the loss of majority of RD and SMS-CTR FN-RMS cells. However, the few surviving clones that repopulate cultures have recovered NOTCH1. Cells that re-establish NOTCH1 expression after SNAI2 ablation are unable to differentiate robustly as SNAI2 shRNA knockdown cells; yet, SNAI2-ablated cells continued to be exquisitely sensitive to ionizing radiation. Thus, we have uncovered a novel mechanism by which SNAI2 and CTCF maintenance of a sub-TAD boundary promotes rather than represses NOTCH1 expression. Further, we demonstrate that SNAI2 suppression of apoptosis post-radiation is independent of SNAI2/NOTCH1 effects on self-renewal and differentiation.

Statistical learning quantifies transposable element-mediated cis-regulation

BACKGROUND

Transposable elements (TEs) have colonized the genomes of most metazoans, and many TE-embedded sequences function as cis-regulatory elements (CREs) for genes involved in a wide range of biological processes from early embryogenesis to innate immune responses. Because of their repetitive nature, TEs have the potential to form CRE platforms enabling the coordinated and genome-wide regulation of protein-coding genes by only a handful of trans-acting transcription factors (TFs).

RESULTS

Here, we directly test this hypothesis through mathematical modeling and demonstrate that differences in expression at protein-coding genes alone are sufficient to estimate the magnitude and significance of TE-contributed cis-regulatory activities, even in contexts where TE-derived transcription fails to do so. We leverage hundreds of overexpression experiments and estimate that, overall, gene expression is influenced by TE-embedded CREs situated within approximately 500 kb of promoters. Focusing on the cis-regulatory potential of TEs within the gene regulatory network of human embryonic stem cells, we find that pluripotency-specific and evolutionarily young TE subfamilies can be reactivated by TFs involved in post-implantation embryogenesis. Finally, we show that TE subfamilies can be split into truly regulatorily active versus inactive fractions based on additional information such as matched epigenomic data, observing that TF binding may better predict TE cis-regulatory activity than differences in histone marks.

CONCLUSION

Our results suggest that TE-embedded CREs contribute to gene regulation during and beyond gastrulation. On a methodological level, we provide a statistical tool that infers TE-dependent cis-regulation from RNA-seq data alone, thus facilitating the study of TEs in the next-generation sequencing era.

TBXT dose sensitivity and the decoupling of nascent mesoderm specification from EMT progression in 2D human gastruloids

In the nascent mesoderm, levels of Brachyury (TBXT) expression must be precisely regulated to ensure cells exit the primitive streak and pattern the anterior-posterior axis, but how this varying dosage informs morphogenesis is not well understood. In this study, we define the transcriptional consequences of TBXT dose reduction during early human gastrulation using human induced pluripotent stem cell (hiPSC)-based models of gastrulation and mesoderm differentiation. Multiomic single-nucleus RNA and single-nucleus ATAC sequencing of 2D gastruloids comprised of WT, TBXT heterozygous (TBXT-Het), or TBXT null (TBXT-KO) hiPSCs reveal that varying TBXT dosage does not compromise a cell's ability to differentiate into nascent mesoderm, but that the loss of TBXT significantly delays the temporal progression of the epithelial to mesenchymal transition (EMT). This delay is dependent on TBXT dose, as cells heterozygous for TBXT proceed with EMT at an intermediate pace relative to WT or TBXT-KO. By differentiating iPSCs of the allelic series into nascent mesoderm in a monolayer format, we further illustrate that TBXT dose directly impacts the persistence of junctional proteins and cell-cell adhesions. These results demonstrate that EMT progression can be decoupled from the acquisition of mesodermal identity in the early gastrula and shed light on the mechanisms underlying human embryogenesis.

A body map of super-enhancers and their function in pig

Introduction

Super-enhancers (SEs) are clusters of enhancers that act synergistically to drive the high-level expression of genes involved in cell identity and function. Although SEs have been extensively investigated in humans and mice, they have not been well characterized in pigs.

Methods

Here, we identified 42,380 SEs in 14 pig tissues using chromatin immunoprecipitation sequencing, and statistics of its overall situation, studied the composition and characteristics of SE, and explored the influence of SEs characteristics on gene expression.

Results

We observed that approximately 40% of normal enhancers (NEs) form SEs. Compared to NEs, we found that SEs were more likely to be enriched with an activated enhancer and show activated functions. Interestingly, SEs showed X chromosome depletion and short interspersed nuclear element enrichment, implying that SEs play an important role in sex traits and repeat evolution. Additionally, SE-associated genes exhibited higher expression levels and stronger conservation than NE-associated genes. However, genes with the largest SEs had higher expression levels than those with the smallest SEs, indicating that SE size may influence gene expression. Moreover, we observed a negative correlation between SE gene distance and gene expression, indicating that the proximity of SEs can affect gene activity. Gene ontology enrichment and motif analysis revealed that SEs have strong tissue-specific activity. For example, the CORO2B gene with a brain-specific SE shows strong brain-specific expression, and the phenylalanine hydroxylase gene with liver-specific SEs shows strong liver-specific expression.

Discussion

In this study, we illustrated a body map of SEs and explored their functions in pigs, providing information on the composition and tissue-specific patterns of SEs. This study can serve as a valuable resource of gene regulatory and comparative analyses to the scientific community and provides a theoretical reference for genetic control mechanisms of important traits in pigs.

Engineered Nano-Bio Interfaces for Stem Cell Therapy

Engineered nanomaterials (ENMs) with different topographies provide effective nano-bio interfaces for controlling the differentiation of stem cells. The interaction of stem cells with nanoscale topographies and chemical cues in their microenvironment at the nano-bio interface can guide their fate. The use of nanotopographical cues, in particular nanorods, nanopillars, nanogrooves, nanofibers, and nanopits, as well as biochemical forces mediated factors, including growth factors, cytokines, and extracellular matrix proteins, can significantly impact stem cell differentiation. These factors were seen as very effective in determining the proliferation and spreading of stem cells. The specific outgrowth of stem cells can be decided with size variation of topographic nanomaterial along with variation in matrix stiffness and surface structure like a special arrangement. The precision chemistry enabled controlled design, synthesis, and chemical composition of ENMs can regulate stem cell behaviors. The parameters of size such as aspect ratio, diameter, and pore size of nanotopographic structures are the main factors for specific termination of stem cells. Protein corona nanoparticles (NPs) have shown a powerful facet in stem cell therapy, where combining specific proteins could facilitate a certain stem cell differentiation and cellular proliferation. Nano-bio reactions implicate the interaction between biological entities and nanoparticles, which can be used to tailor the stem cells' culmination. The ion release can also be a parameter to enhance cellular proliferation and to commit the early differentiation of stem cells. Further research is needed to fully understand the mechanisms underlying the interactions between engineered nano-bio interfaces and stem cells and to develop optimized regenerative medicine and tissue engineering designs.

Co-effects of m6A and chromatin accessibility dynamics in the regulation of cardiomyocyte differentiation

BACKGROUND

Cardiomyocyte growth and differentiation rely on precise gene expression regulation, with epigenetic modifications emerging as key players in this intricate process. Among these modifications, N6-methyladenosine (m6A) stands out as one of the most prevalent modifications on mRNA, exerting influence over mRNA metabolism and gene expression. However, the specific function of m6A in cardiomyocyte differentiation remains poorly understood.

RESULTS

We investigated the relationship between m6A modification and cardiomyocyte differentiation by conducting a comprehensive profiling of m6A dynamics during the transition from pluripotent stem cells to cardiomyocytes. Our findings reveal that while the overall m6A modification level remains relatively stable, the m6A levels of individual genes undergo significant changes throughout cardiomyocyte differentiation. We discovered the correlation between alterations in chromatin accessibility and the binding capabilities of m6A writers, erasers, and readers. The changes in chromatin accessibility influence the recruitment and activity of m6A regulatory proteins, thereby impacting the levels of m6A modification on specific mRNA transcripts.

CONCLUSION

Our data demonstrate that the coordinated dynamics of m6A modification and chromatin accessibility are prominent during the cardiomyocyte differentiation.

H3K36 methylation maintains cell identity by regulating opposing lineage programmes

The epigenetic mechanisms that maintain differentiated cell states remain incompletely understood. Here we employed histone mutants to uncover a crucial role for H3K36 methylation in the maintenance of cell identities across diverse developmental contexts. Focusing on the experimental induction of pluripotency, we show that H3K36M-mediated depletion of H3K36 methylation endows fibroblasts with a plastic state poised to acquire pluripotency in nearly all cells. At a cellular level, H3K36M facilitates epithelial plasticity by rendering fibroblasts insensitive to TGFβ signals. At a molecular level, H3K36M enables the decommissioning of mesenchymal enhancers and the parallel activation of epithelial/stem cell enhancers. This enhancer rewiring is Tet dependent and redirects Sox2 from promiscuous somatic to pluripotency targets. Our findings reveal a previously unappreciated dual role for H3K36 methylation in the maintenance of cell identity by integrating a crucial developmental pathway into sustained expression of cell-type-specific programmes, and by opposing the expression of alternative lineage programmes through enhancer methylation.

The time-resolved genomic impact of Wnt/β-catenin signaling

Wnt signaling orchestrates gene expression via its effector, β-catenin. However, it is unknown whether β-catenin binds its target genomic regions simultaneously and how this impacts chromatin dynamics to modulate cell behavior. Using a combination of time-resolved CUT&RUN against β-catenin, ATAC-seq, and perturbation assays in different cell types, we show that Wnt/β-catenin physical targets are tissue-specific, β-catenin "moves" on different loci over time, and its association to DNA accompanies changing chromatin accessibility landscapes that determine cell behavior. In particular, Wnt/β-catenin progressively shapes the chromatin of human embryonic stem cells (hESCs) as they undergo mesodermal differentiation, a behavior that we define as "plastic." In HEK293T cells, on the other hand, Wnt/β-catenin drives a transient chromatin opening, followed by re-establishment of the pre-stimulation state, a response that we define as "elastic." Future experiments shall assess whether other cell communication mechanisms, in addition to Wnt signaling, are ruled by time, cellular idiosyncrasies, and chromatin constraints. A record of this paper's transparent peer review process is included in the supplemental information.

Comprehensive tissue deconvolution of cell-free DNA by deep learning for disease diagnosis and monitoring

Plasma cell-free DNA (cfDNA) is a noninvasive biomarker for cell death of all organs. Deciphering the tissue origin of cfDNA can reveal abnormal cell death because of diseases, which has great clinical potential in disease detection and monitoring. Despite the great promise, the sensitive and accurate quantification of tissue-derived cfDNA remains challenging to existing methods due to the limited characterization of tissue methylation and the reliance on unsupervised methods. To fully exploit the clinical potential of tissue-derived cfDNA, here we present one of the largest comprehensive and high-resolution methylation atlas based on 521 noncancer tissue samples spanning 29 major types of human tissues. We systematically identified fragment-level tissue-specific methylation patterns and extensively validated them in orthogonal datasets. Based on the rich tissue methylation atlas, we develop the first supervised tissue deconvolution approach, a deep-learning-powered model, cfSort, for sensitive and accurate tissue deconvolution in cfDNA. On the benchmarking data, cfSort showed superior sensitivity and accuracy compared to the existing methods. We further demonstrated the clinical utilities of cfSort with two potential applications: aiding disease diagnosis and monitoring treatment side effects. The tissue-derived cfDNA fraction estimated from cfSort reflected the clinical outcomes of the patients. In summary, the tissue methylation atlas and cfSort enhanced the performance of tissue deconvolution in cfDNA, thus facilitating cfDNA-based disease detection and longitudinal treatment monitoring.

Discovery of Competent Chromatin Regions in Human Embryonic Stem Cells

The mechanisms underlying the ability of embryonic stem cells (ESCs) to rapidly activate lineage-specific genes during differentiation remain largely unknown. Through multiple CRISPR-activation screens, we discovered human ESCs have pre-established transcriptionally competent chromatin regions (CCRs) that support lineage-specific gene expression at levels comparable to differentiated cells. CCRs reside in the same topological domains as their target genes. They lack typical enhancer-associated histone modifications but show enriched occupancy of pluripotent transcription factors, DNA demethylation factors, and histone deacetylases. TET1 and QSER1 protect CCRs from excessive DNA methylation, while HDAC1 family members prevent premature activation. This "push and pull" feature resembles bivalent domains at developmental gene promoters but involves distinct molecular mechanisms. Our study provides new insights into pluripotency regulation and cellular plasticity in development and disease.

One sentence summary

We report a class of distal regulatory regions distinct from enhancers that confer human embryonic stem cells with the competence to rapidly activate the expression of lineage-specific genes.

Epigenetic Control of Cell Potency and Fate Determination during Mammalian Gastrulation

Pluripotent embryonic stem cells have a unique and characteristic epigenetic profile, which is critical for differentiation to all embryonic germ lineages. When stem cells exit the pluripotent state and commit to lineage-specific identities during the process of gastrulation in early embryogenesis, extensive epigenetic remodelling mediates both the switch in cellular programme and the loss of potential to adopt alternative lineage programmes. However, it remains to be understood how the stem cell epigenetic profile encodes pluripotency, or how dynamic epigenetic regulation helps to direct cell fate specification. Recent advances in stem cell culture techniques, cellular reprogramming, and single-cell technologies that can quantitatively profile epigenetic marks have led to significant insights into these questions, which are important for understanding both embryonic development and cell fate engineering. This review provides an overview of key concepts and highlights exciting new advances in the field.

Cellular population dynamics shape the route to human pluripotency

Human cellular reprogramming to induced pluripotency is still an inefficient process, which has hindered studying the role of critical intermediate stages. Here we take advantage of high efficiency reprogramming in microfluidics and temporal multi-omics to identify and resolve distinct sub-populations and their interactions. We perform secretome analysis and single-cell transcriptomics to show functional extrinsic pathways of protein communication between reprogramming sub-populations and the re-shaping of a permissive extracellular environment. We pinpoint the HGF/MET/STAT3 axis as a potent enhancer of reprogramming, which acts via HGF accumulation within the confined system of microfluidics, and in conventional dishes needs to be supplied exogenously to enhance efficiency. Our data suggest that human cellular reprogramming is a transcription factor-driven process that it is deeply dependent on extracellular context and cell population determinants.

Small-sample learning reveals propionylation in determining global protein homeostasis

Proteostasis is fundamental for maintaining organismal health. However, the mechanisms underlying its dynamic regulation and how its disruptions lead to diseases are largely unclear. Here, we conduct in-depth propionylomic profiling in Drosophila, and develop a small-sample learning framework to prioritize the propionylation at lysine 17 of H2B (H2BK17pr) to be functionally important. Mutating H2BK17 which eliminates propionylation leads to elevated total protein level in vivo. Further analyses reveal that H2BK17pr modulates the expression of 14.7-16.3% of genes in the proteostasis network, and determines global protein level by regulating the expression of genes involved in the ubiquitin-proteasome system. In addition, H2BK17pr exhibits daily oscillation, mediating the influences of feeding/fasting cycles to drive rhythmic expression of proteasomal genes. Our study not only reveals a role of lysine propionylation in regulating proteostasis, but also implements a generally applicable method which can be extended to other issues with little prior knowledge.

Multifactorial profiling of epigenetic landscapes at single-cell resolution using MulTI-Tag

Chromatin profiling at locus resolution uncovers gene regulatory features that define cell types and developmental trajectories, but it remains challenging to map and compare different chromatin-associated proteins in the same sample. Here we describe Multiple Target Identification by Tagmentation (MulTI-Tag), an antibody barcoding approach for profiling multiple chromatin features simultaneously in single cells. We optimized MulTI-Tag to retain high sensitivity and specificity, and we demonstrate detection of up to three histone modifications in the same cell: H3K27me3, H3K4me1/2 and H3K36me3. We apply MulTI-Tag to resolve distinct cell types and developmental trajectories; to distinguish unique, coordinated patterns of active and repressive element regulatory usage associated with differentiation outcomes; and to uncover associations between histone marks. Multifactorial epigenetic profiling holds promise for comprehensively characterizing cell-specific gene regulatory landscapes in development and disease.

Escape from oncogene-induced senescence is controlled by POU2F2 and memorized by chromatin scars

Although oncogene-induced senescence (OIS) is a potent tumor-suppressor mechanism, recent studies revealed that cells could escape from OIS with features of transformed cells. However, the mechanisms that promote OIS escape remain unclear, and evidence of post-senescent cells in human cancers is missing. Here, we unravel the regulatory mechanisms underlying OIS escape using dynamic multidimensional profiling. We demonstrate a critical role for AP1 and POU2F2 transcription factors in escape from OIS and identify senescence-associated chromatin scars (SACSs) as an epigenetic memory of OIS detectable during colorectal cancer progression. POU2F2 levels are already elevated in precancerous lesions and as cells escape from OIS, and its expression and binding activity to cis-regulatory elements are associated with decreased patient survival. Our results support a model in which POU2F2 exploits a precoded enhancer landscape necessary for senescence escape and reveal POU2F2 and SACS gene signatures as valuable biomarkers with diagnostic and prognostic potential.

Dynamic antagonism between key repressive pathways maintains the placental epigenome

DNA and Histone 3 Lysine 27 methylation typically function as repressive modifications and operate within distinct genomic compartments. In mammals, the majority of the genome is kept in a DNA methylated state, whereas the Polycomb repressive complexes regulate the unmethylated CpG-rich promoters of developmental genes. In contrast to this general framework, the extra-embryonic lineages display non-canonical, globally intermediate DNA methylation levels, including disruption of local Polycomb domains. Here, to better understand this unusual landscape's molecular properties, we genetically and chemically perturbed major epigenetic pathways in mouse trophoblast stem cells. We find that the extra-embryonic epigenome reflects ongoing and dynamic de novo methyltransferase recruitment, which is continuously antagonized by Polycomb to maintain intermediate, locally disordered methylation. Despite its disorganized molecular appearance, our data point to a highly controlled equilibrium between counteracting repressors within extra-embryonic cells, one that can seemingly persist indefinitely without bistable features typically seen for embryonic forms of epigenetic regulation.

Leukemia core transcriptional circuitry is a sparsely interconnected hierarchy stabilized by incoherent feed-forward loops

Lineage-defining transcription factors form densely interconnected circuits in chromatin occupancy assays, but the functional significance of these networks remains underexplored. We reconstructed the functional topology of a leukemia cell transcription network from the direct gene-regulatory programs of eight core transcriptional regulators established in pre-steady state assays coupling targeted protein degradation with nascent transcriptomics. The core regulators displayed narrow, largely non-overlapping direct transcriptional programs, forming a sparsely interconnected functional hierarchy stabilized by incoherent feed-forward loops. BET bromodomain and CDK7 inhibitors disrupted the core regulators' direct programs, acting as mixed agonists/antagonists. The network is predictive of dynamic gene expression behaviors in time-resolved assays and clinically relevant pathway activity in patient populations.

H3.3 contributes to chromatin accessibility and transcription factor binding at promoter-proximal regulatory elements in embryonic stem cells

BACKGROUND

The histone variant H3.3 is enriched at active regulatory elements such as promoters and enhancers in mammalian genomes. These regions are highly accessible, creating an environment that is permissive to transcription factor binding and the recruitment of transcriptional coactivators that establish a unique chromatin post-translational landscape. How H3.3 contributes to the establishment and function of chromatin states at these regions is poorly understood.

RESULTS

We perform genomic analyses of features associated with active promoter chromatin in mouse embryonic stem cells (ESCs) and find evidence of subtle yet widespread promoter dysregulation in the absence of H3.3. Loss of H3.3 results in reduced chromatin accessibility and transcription factor (TF) binding at promoters of expressed genes in ESCs. Likewise, enrichment of the transcriptional coactivator p300 and downstream histone H3 acetylation at lysine 27 (H3K27ac) is reduced at promoters in the absence of H3.3, along with reduced enrichment of the acetyl lysine reader BRD4. Despite the observed chromatin dysregulation, H3.3 KO ESCs maintain transcription from ESC-specific genes. However, upon undirected differentiation, H3.3 KO cells retain footprinting of ESC-specific TF motifs and fail to generate footprints of lineage-specific TF motifs, in line with their diminished capacity to differentiate.

CONCLUSIONS

H3.3 facilitates DNA accessibility, transcription factor binding, and histone post-translational modification at active promoters. While H3.3 is not required for maintaining transcription in ESCs, it does promote de novo transcription factor binding which may contribute to the dysregulation of cellular differentiation in the absence of H3.3.

Specification of human germ cell fate with enhanced progression capability supported by hindgut organoids

Human primordial germ cells (hPGCs), the precursors of sperm and eggs, are specified during weeks 2-3 after fertilization. Few studies on ex vivo and in vitro cultured human embryos reported plausible hPGCs on embryonic day (E) 12-13 and in an E16-17 gastrulating embryo. In vitro, hPGC-like cells (hPGCLCs) can be specified from the intermediary pluripotent stage or peri-gastrulation precursors. Here, we explore the broad spectrum of hPGCLC precursors and how different precursors impact hPGCLC development. We show that resetting precursors can give rise to hPGCLCs (rhPGCLCs) in response to BMP. Strikingly, rhPGCLCs co-cultured with human hindgut organoids progress at a pace reminiscent of in vivo hPGC development, unlike those derived from peri-gastrulation precursors. Moreover, rhPGCLC specification depends on both EOMES and TBXT, not just on EOMES as for peri-gastrulation hPGCLCs. Importantly, our study provides the foundation for developing efficient in vitro models of human gametogenesis.

T-REX17 is a transiently expressed non-coding RNA essential for human endoderm formation

Long non-coding RNAs (lncRNAs) have emerged as fundamental regulators in various biological processes, including embryonic development and cellular differentiation. Despite much progress over the past decade, the genome-wide annotation of lncRNAs remains incomplete and many known non-coding loci are still poorly characterized. Here, we report the discovery of a previously unannotated lncRNA that is transcribed 230 kb upstream of the SOX17 gene and located within the same topologically associating domain. We termed it T-REX17 (Transcript Regulating Endoderm and activated by soX17) and show that it is induced following SOX17 activation but its expression is more tightly restricted to early definitive endoderm. Loss of T-REX17 affects crucial functions independent of SOX17 and leads to an aberrant endodermal transcriptome, signaling pathway deregulation and epithelial to mesenchymal transition defects. Consequently, cells lacking the lncRNA cannot further differentiate into more mature endodermal cell types. Taken together, our study identified and characterized T-REX17 as a transiently expressed and essential non-coding regulator in early human endoderm differentiation.

Epigenetic and transcriptional regulations prime cell fate before division during human pluripotent stem cell differentiation

Stem cells undergo cellular division during their differentiation to produce daughter cells with a new cellular identity. However, the epigenetic events and molecular mechanisms occurring between consecutive cell divisions have been insufficiently studied due to technical limitations. Here, using the FUCCI reporter we developed a cell-cycle synchronised human pluripotent stem cell (hPSC) differentiation system for uncovering epigenome and transcriptome dynamics during the first two divisions leading to definitive endoderm. We observed that transcription of key differentiation markers occurs before cell division, while chromatin accessibility analyses revealed the early inhibition of alternative cell fates. We found that Activator protein-1 members controlled by p38/MAPK signalling are necessary for inducing endoderm while blocking cell fate shifting toward mesoderm, and that enhancers are rapidly established and decommissioned between different cell divisions. Our study has practical biomedical utility for producing hPSC-derived patient-specific cell types since p38/MAPK induction increased the differentiation efficiency of insulin-producing pancreatic beta-cells.

Temporal transcriptional control of neural induction in human induced pluripotent stem cells

Introduction

Neural induction of human induced pluripotent stem cells represents a critical switch in cell state during which pluripotency is lost and commitment to a neural lineage is initiated. Although many of the key transcription factors involved in neural induction are known, we know little of the temporal and causal relationships that are required for this state transition.

Methods

Here, we have carried out a longitudinal analysis of the transcriptome of human iPSCs undergoing neural induction. Using the temporal relationships between the changing profile of key transcription factors and subsequent changes in their target gene expression profiles, we have identified distinct functional modules operative throughout neural induction.

Results

In addition to modules that govern loss of pluripotency and gain of neural ectoderm identity, we discover other modules governing cell cycle and metabolism. Strikingly, some of these functional modules are retained throughout neural induction, even though the gene membership of the module changes. Systems analysis identifies other modules associated with cell fate commitment, genome integrity, stress response and lineage specification. We then focussed on OTX2, one of the most precociously activated transcription factors during neural induction. Our temporal analysis of OTX2 target gene expression identified several OTX2 regulated gene modules representing protein remodelling, RNA splicing and RNA processing. Further CRISPRi inhibition of OTX2 prior to neural induction promotes an accelerated loss of pluripotency and a precocious and aberrant neural induction disrupting some of the previously identified modules.

Discussion

We infer that OTX2 has a diverse role during neural induction and regulates many of the biological processes that are required for loss of pluripotency and gain of neural identity. This dynamical analysis of transcriptional changes provides a unique perspective of the widespread remodelling of the cell machinery that occurs during neural induction of human iPSCs.

Cloning by SCNT: Integrating Technical and Biology-Driven Advances

Somatic cell nuclear transfer (SCNT) into enucleated oocytes initiates nuclear reprogramming of lineage-committed cells to totipotency. Pioneer SCNT work culminated with cloned amphibians from tadpoles, while technical and biology-driven advances led to cloned mammals from adult animals. Cloning technology has been addressing fundamental questions in biology, propagating desired genomes, and contributing to the generation of transgenic animals or patient-specific stem cells. Nonetheless, SCNT remains technically complex and cloning efficiency relatively low. Genome-wide technologies revealed barriers to nuclear reprogramming, such as persistent epigenetic marks of somatic origin and reprogramming resistant regions of the genome. To decipher the rare reprogramming events that are compatible with full-term cloned development, it will likely require technical advances for large-scale production of SCNT embryos alongside extensive profiling by single-cell multi-omics. Altogether, cloning by SCNT remains a versatile technology, while further advances should continuously refresh the excitement of its applications.

ESRG is critical to maintain the cell survival and self-renewal/pluripotency of hPSCs by collaborating with MCM2 to suppress p53 pathway

The mechanisms of self-renewal and pluripotency maintenance of human pluripotent stem cells (hPSCs) have not been fully elucidated, especially for the role of those poorly characterized long noncoding RNAs (lncRNAs). ESRG is a lncRNA highly expressed in hPSCs, and its functional roles are being extensively explored in the field. Here, we identified that the transcription of ESRG can be directly regulated by OCT4, a key self-renewal factor in hPSCs. Knockdown of ESRG induces hPSC differentiation, cell cycle arrest, and apoptosis. ESRG binds to MCM2, a replication-licensing factor, to sustain its steady-state level and nuclear location, safeguarding error-free DNA replication. Further study showed that ESRG knockdown leads to MCM2 abnormalities, resulting in DNA damage and activation of the p53 pathway, ultimately impairs hPSC self-renewal and pluripotency, and induces cell apoptosis. In summary, our study suggests that ESRG, as a novel target of OCT4, plays an essential role in maintaining the cell survival and self-renewal/pluripotency of hPSCs in collaboration with MCM2 to suppress p53 signaling. These findings provide critical insights into the mechanisms underlying the maintenance of self-renewal and pluripotency in hPSCs by lncRNAs.

Protein sumoylation in normal and cancer stem cells

Stem cells with the capacity of self-renewal and differentiation play pivotal roles in normal tissues and malignant tumors. Whereas stem cells are supposed to be genetically identical to their non-stem cell counterparts, cell stemness is deliberately regulated by a dynamic network of molecular mechanisms. Reversible post-translational protein modifications (PTMs) are rapid and reversible non-genetic processes that regulate essentially all physiological and pathological process. Numerous studies have reported the involvement of post-translational protein modifications in the acquirement and maintenance of cell stemness. Recent studies underscore the importance of protein sumoylation, i.e., the covalent attachment of the small ubiquitin-like modifiers (SUMO), as a critical post-translational protein modification in the stem cell populations in development and tumorigenesis. In this review, we summarize the functions of protein sumoylation in different kinds of normal and cancer stem cells. In addition, we describe the upstream regulators and the downstream effectors of protein sumoylation associated with cell stemness. We also introduce the translational studies aiming at sumoylation to target stem cells for disease treatment. Finally, we propose future directions for sumoylation studies in stem cells.

Locus-Specific Isolation of the Nanog Chromatin Identifies Regulators Relevant to Pluripotency of Mouse Embryonic Stem Cells and Reprogramming of Somatic Cells

Pluripotency is a crucial feature of pluripotent stem cells, which are regulated by the core pluripotency network consisting of key transcription factors and signaling molecules. However, relatively less is known about the molecular mechanisms that modify the core pluripotency network. Here we used the CAPTURE (CRISPR Affinity Purification in situ of Regulatory Elements) to unbiasedly isolate proteins assembled on the Nanog promoter in mouse embryonic stem cells (mESCs), and then tested their functional relevance to the maintenance of mESCs and reprogramming of somatic cells. Gene ontology analysis revealed that the identified proteins, including many RNA-binding proteins (RBPs), are enriched in RNA-related functions and gene expression. ChIP-qPCR experiments confirmed that BCLAF1, FUBP1, MSH6, PARK7, PSIP1, and THRAP3 occupy the Nanog promoter region in mESCs. Knockdown experiments of these factors show that they play varying roles in self-renewal, pluripotency gene expression, and differentiation of mESCs as well as in the reprogramming of somatic cells. Our results show the utility of unbiased identification of chromatin-associated proteins on a pluripotency gene in mESCs and reveal the functional relevance of RBPs in ESC differentiation and somatic cell reprogramming.

Carnegie in 4D? Stem-cell-based models of human embryo development

How cells build embryos is still a major mystery. Many unresolved questions require the study of the processes that pattern and shape the embryo in live specimens, in toto, across spatial and temporal scales. In mammalian embryogenesis, this remains a major challenge as the embryo develops in utero, precluding easy accessibility. For human embryos, technical, ethical and legal limitations further hamper the in-depth investigation of embryogenesis, especially beyond gastrulation stages. This has resulted in an over-reliance on model organisms, particularly mice, to understand mammalian development. However, recent efforts show critical differences between rodent and primate embryos, including timing, architecture and transcriptional regulation. Thus, a human-centric understanding of embryogenesis is much needed. To empower this, novel in vitro approaches, which coax human pluripotent stem cells to form embryonic organoids that model embryo development, are pivotal. Here, we summarize these emergent technologies that recapitulate aspects of human development "in a dish". We show how these technologies can provide insights into the molecular, cellular and morphogenetic processes that fuel the formation of a fully formed fetus, and discuss the potential of these platforms to revolutionize our understanding of human development in health and disease. Despite their clear promise, we caution against over-interpreting the extent to which these in vitro platforms model the natural embryo. In particular, we discuss how fate, form and function - a tightly coupled trinity in vivo, can be disconnected in vitro. Finally, we propose how careful benchmarking of existing models, in combination with rational protocol design based on an increased understanding of in vivo developmental dynamics and insights from mouse in vitro models of embryo development, will help guide the establishment of better models of human embryo development.

Cardiac Development in the Presence of Cadmium: An in Vitro Study Using Human Embryonic Stem Cells and Cardiac Organoids

BACKGROUND

Exposure to cadmium (Cd) is associated with cardiovascular diseases. Maternal Cd exposure is a significant risk factor for congenital heart disease. However, mechanisms of Cd on developmental cardiotoxicity are not well defined.

OBJECTIVES

We evaluated the effects of Cd on the different stages (mesoderm, cardiac induction, cardiac function) of cardiac development using an early embryo development in vitro model and two- or three-dimensional (2- or 3D) cardiomyocyte and cardiac organoid formation models mimicking early cardiac development.

METHODS

Embryonic stem cells (ESCs) form 3D aggregates, called embryoid bodies, that recapitulate events involved with early embryogenesis (e.g., germ layer formation). This model was used for early germ layer formation and signaling pathway identification. The 2D cardiomyocyte differentiation from the NKX2-5eGFP/w human ESCs model was used to explore the effects of Cd exposure on cardiomyocyte formation and to model mesoderm differentiation and cardiac induction, allowing us to explore different developmental windows of Cd toxicity. The 3D cardiac organoid model was used in evaluating the effects of Cd exposure on contractility and cardiac development.

RESULTS

Cd (0.6μM; 110 ppb) lowered the differentiation of embryoid bodies to mesoderm via suppression of Wnt/β-catenin-signaling pathways. During early mesoderm induction, the mesoderm-associated transcription factors MESP1 and EOMES showed a transient up-regulation, which decreased later in the cardiac induction stage. Cd (0.15μM) lowered mesoderm formation and cardiac induction through suppression of the transcription factors and mesoderm marker genes HAND1, SNAI2, HOPX, and the cardiac-specific genes NKX2-5, GATA4, troponin T, and alpha-actinin. In addition, Cd-induced histone modifications for both gene activation (H3K4me3) and repression (H3K27me3), which play vital roles in regulating mesoderm commitment markers. The effects of Cd inhibition on cardiomyocyte differentiation were confirmed in 3D cardiac organoids.

DISCUSSION

In conclusion, using a human ESC-derived 2D/3D in vitro differentiation model system and cardiac organoids, we demonstrated that low-dose Cd suppressed mesoderm formation through mesoderm gene histone modification, thus inhibiting cardiomyocyte differentiation and cardiac induction. The studies provide valuable insights into cellular events and molecular mechanisms associated with Cd-induced congenital heart disease. https://doi.org/10.1289/EHP11208.

Evolution and antiviral activity of a human protein of retroviral origin

Endogenous retroviruses are abundant components of mammalian genomes descended from ancient germline infections. In several mammals, the envelope proteins encoded by these elements protect against exogenous viruses, but this activity has not been documented with endogenously expressed envelopes in humans. We report that the human genome harbors a large pool of envelope-derived sequences with the potential to restrict retroviral infection. To test this, we characterized an envelope-derived protein, Suppressyn. We found that Suppressyn is expressed in human preimplantation embryos and developing placenta using its ancestral retroviral promoter. Cell culture assays showed that Suppressyn, and its hominoid orthologs, could restrict infection by extant mammalian type D retroviruses. Our data support a generalizable model of retroviral envelope co-option for host immunity and genome defense.

Mechanisms underlying WNT-mediated priming of human embryonic stem cells

Embryogenesis is guided by a limited set of signaling pathways dynamically expressed in different places. How a context-dependent signaling response is generated has been a central question of developmental biology, which can now be addressed with in vitro models of human embryos that are derived from embryonic stem cells (hESCs). Our previous work demonstrated that during early stages of hESC differentiation, cells chronicle signaling hierarchy. Only cells that have been exposed (primed) by WNT signaling can respond to subsequent activin exposure and differentiate to mesendodermal (ME) fates. Here, we show that WNT priming does not alter SMAD2 binding nor its chromatin opening but, instead, acts by inducing the expression of the SMAD2 co-factor EOMES. Expression of EOMES is sufficient to replace WNT upstream of activin-mediated ME differentiation, thus unveiling the mechanistic basis for priming and cellular memory in early development.

Human-gained heart enhancers are associated with species-specific cardiac attributes

The heart, a vital organ which is first to develop, has adapted its size, structure and function in order to accommodate the circulatory demands for a broad range of animals. Although heart development is controlled by a relatively conserved network of transcriptional/chromatin regulators, how the human heart has evolved species-specific features to maintain adequate cardiac output and function remains to be defined. Here, we show through comparative epigenomic analysis the identification of enhancers and promoters that have gained activity in humans during cardiogenesis. These cis-regulatory elements (CREs) are associated with genes involved in heart development and function, and may account for species-specific differences between human and mouse hearts. Supporting these findings, genetic variants that are associated with human cardiac phenotypic/disease traits, particularly those differing between human and mouse, are enriched in human-gained CREs. During early stages of human cardiogenesis, these CREs are also gained within genomic loci of transcriptional regulators, potentially expanding their role in human heart development. In particular, we discovered that gained enhancers in the locus of the early human developmental regulator ZIC3 are selectively accessible within a subpopulation of mesoderm cells which exhibits cardiogenic potential, thus possibly extending the function of ZIC3 beyond its conserved left-right asymmetry role. Genetic deletion of these enhancers identified a human gained enhancer that was required for not only ZIC3 and early cardiac gene expression at the mesoderm stage but also cardiomyocyte differentiation. Overall, our results illuminate how human gained CREs may contribute to human-specific cardiac attributes, and provide insight into how transcriptional regulators may gain cardiac developmental roles through the evolutionary acquisition of enhancers.

Reorganization of 3D chromatin architecture in doxorubicin-resistant breast cancer cells

Background: Doxorubicin resistance remains a major therapeutic challenge leading to poor survival prognosis and treatment failure in breast cancer. Although doxorubicin induces massive changes in the transcriptional landscape are well known, potential diagnostic or therapeutic targets associated with the reorganization of three-dimensional (3D) chromatin architecture have not yet been systematically investigated.

Methods: Here we performed in situ high-throughput chromosome conformation capture (Hi-C) on parental and doxorubicin-resistant MCF7 (MCF7-DR) human breast cancer cells, followed by integrative analysis of HiC, ATAC-seq, RNA-seq and TCGA data.

Results: It revealed that A/B compartment switching was positively correlated to genome-wide differential gene expression. The genome of MCF7-DR cells was spatially reorganized into smaller topologically associating domains (TADs) and chromatin loops. We also revealed the contribution of increased chromatin accessibility and potential transcription factor families, including CTCF, AP-1 and bHLH, to gained TADs or loops. Intriguingly, we observed two condensed genomic regions (∼20 kb) with decreased chromatin accessibility flanking TAD boundaries, which might play a critical role in the formation or maintenance of TADs. Finally, combining data from TCGA, we identified a number of gained and lost enhancer-promoter interactions and their corresponding differentially expressed genes involved in chromatin organization and breast cancer signaling pathways, including FA2H, FOXA1 and JRKL, which might serve as potential treatment targets for breast cancer.

Conclusion: These data uncovered a close connection between 3D genome reorganization, chromatin accessibility as well as gene transcription and provide novel insights into the epigenomic mechanisms involving doxorubicin resistance in breast cancer.

Eomes function is conserved between zebrafish and mouse and controls left-right organiser progenitor gene expression via interlocking feedforward loops

The T-box family transcription factor Eomesodermin (Eomes) is present in all vertebrates, with many key roles in the developing mammalian embryo and immune system. Homozygous Eomes mutant mouse embryos exhibit early lethality due to defects in both the embryonic mesendoderm and the extraembryonic trophoblast cell lineage. In contrast, zebrafish lacking the predominant Eomes homologue A (Eomesa) do not suffer complete lethality and can be maintained. This suggests fundamental differences in either the molecular function of Eomes orthologues or the molecular configuration of processes in which they participate. To explore these hypotheses we initially analysed the expression of distinct Eomes isoforms in various mouse cell types. Next we compared the functional capabilities of these murine isoforms to zebrafish Eomesa. These experiments provided no evidence for functional divergence. Next we examined the functions of zebrafish Eomesa and other T-box family members expressed in early development, as well as its paralogue Eomesb. Though Eomes is a member of the Tbr1 subfamily we found evidence for functional redundancy with the Tbx6 subfamily member Tbx16, known to be absent from eutherians. However, Tbx16 does not appear to synergise with Eomesa cofactors Mixl1 and Gata5. Finally, we analysed the ability of Eomesa and other T-box factors to induce zebrafish left-right organiser progenitors (known as dorsal forerunner cells) known to be positively regulated by vgll4l, a gene we had previously shown to be repressed by Eomesa. Here we demonstrate that Eomesa indirectly upregulates vgll4l expression via interlocking feedforward loops, suggesting a role in establishment of left-right asymmetry. Conversely, other T-box factors could not similarly induce left-right organiser progenitors. Overall these findings demonstrate conservation of Eomes molecular function and participation in similar processes, but differential requirements across evolution due to additional co-expressed T-box factors in teleosts, albeit with markedly different molecular capabilities. Our analyses also provide insights into the role of Eomesa in left-right organiser formation in zebrafish.

Repetitive DNA symmetry elements negatively regulate gene expression in embryonic stem cells

Transcription factor (TF) binding to genomic DNA elements constitutes one of the key mechanisms that regulates gene expression program in cells. Both consensus and nonconsensus DNA sequence elements influence the recognition specificity of TFs. Based on the analysis of experimentally determined c-Myc binding preferences to genomic DNA, here we statistically predict that certain repetitive, nonconsensus DNA symmetry elements can relatively reduce TF-DNA binding preferences. This is in contrast to a different set of repetitive, nonconsensus symmetry elements that can increase the strength of TF-DNA binding. Using c-Myc enhancer reporter system containing consensus motif flanked by nonconsensus sequences in embryonic stem cells, we directly demonstrate that the enrichment in such negatively regulating repetitive symmetry elements is sufficient to reduce the gene expression level compared with native genomic sequences. Negatively regulating repetitive symmetry elements around consensus c-Myc motif and DNA sequences containing consensus c-Myc motif flanked by entirely randomized sequences show similar expression baseline. A possible explanation for this observation is that rather than complete repression, negatively regulating repetitive symmetry elements play a regulatory role in fine-tuning the reduction of gene expression, most probably by binding TFs other than c-Myc.

Epithelial-Mesenchymal Transition Suppresses AMPK and Sensitizes Cancer Cells to Pyroptosis under Energy Stress

Epithelial-mesenchymal transition (EMT) is implicated in tumor metastasis and therapeutic resistance. It remains a challenge to target cancer cells that have undergone EMT. The Snail family of key EMT-inducing transcription factors directly binds to and transcriptionally represses not only epithelial genes but also a myriad of additional genomic targets that may carry out significant biological functions. Therefore, we reasoned that EMT inherently causes various concomitant phenotypes, some of which may create targetable vulnerabilities for cancer treatment. In the present study, we found that Snail transcription factors bind to the promoters of multiple genes encoding subunits of the AMP-activated protein kinase (AMPK) complex, and expression of AMPK genes was markedly downregulated by EMT. Accordingly, high AMPK expression in tumors correlated with epithelial cell markers and low AMPK expression in tumors was strongly associated with adverse prognosis. AMPK is the principal sensor of cellular energy status. In response to energy stress, AMPK is activated and critically reprograms cellular metabolism to restore energy homeostasis and maintain cell survival. We showed that activation of AMPK by energy stress was severely impaired by EMT. Consequently, EMT cancer cells became hypersensitive to a variety of energy stress conditions and primarily underwent pyroptosis, a regulated form of necrotic cell death. Collectively, the study suggests that EMT impedes the activation of AMPK signaling induced by energy stress and sensitizes cancer cells to pyroptotic cell death under energy stress conditions. Therefore, while EMT promotes malignant progression, it concurrently induces collateral vulnerabilities that may be therapeutically exploited.

Definition of germ layer cell lineage alternative splicing programs reveals a critical role for Quaking in specifying cardiac cell fate

Alternative splicing is critical for development; however, its role in the specification of the three embryonic germ layers is poorly understood. By performing RNA-Seq on human embryonic stem cells (hESCs) and derived definitive endoderm, cardiac mesoderm, and ectoderm cell lineages, we detect distinct alternative splicing programs associated with each lineage. The most prominent splicing program differences are observed between definitive endoderm and cardiac mesoderm. Integrative multi-omics analyses link each program with lineage-enriched RNA binding protein regulators, and further suggest a widespread role for Quaking (QKI) in the specification of cardiac mesoderm. Remarkably, knockout of QKI disrupts the cardiac mesoderm-associated alternative splicing program and formation of myocytes. These changes arise in part through reduced expression of BIN1 splice variants linked to cardiac development. Mechanistically, we find that QKI represses inclusion of exon 7 in BIN1 pre-mRNA via an exonic ACUAA motif, and this is concomitant with intron removal and cleavage from chromatin. Collectively, our results uncover alternative splicing programs associated with the three germ lineages and demonstrate an important role for QKI in the formation of cardiac mesoderm.

Optimization of the TeraTox assay for preclinical teratogenicity assessment

Current animal-free methods to assess teratogenicity of drugs under development still deliver high numbers of false negatives. To improve the sensitivity of human teratogenicity prediction, we characterized the TeraTox test, a newly developed multilineage differentiation assay using 3D human-induced pluripotent stem cells. TeraTox produces primary output concentration-dependent cytotoxicity and altered gene expression induced by each test compound. These data are fed into an interpretable machine-learning model to perform prediction, which relates to the concentration-dependent human teratogenicity potential of drug candidates. We applied TeraTox to profile 33 approved pharmaceuticals and 12 proprietary drug candidates with known in vivo data. Comparing TeraTox predictions with known human or animal toxicity, we report an accuracy of 69% (specificity: 53%, sensitivity: 79%). TeraTox performed better than 2 quantitative structure-activity relationship models and had a higher sensitivity than the murine embryonic stem cell test (accuracy: 58%, specificity: 76%, and sensitivity: 46%) run in the same laboratory. The overall prediction accuracy could be further improved by combining TeraTox and mouse embryonic stem cell test results. Furthermore, patterns of altered gene expression revealed by TeraTox may help grouping toxicologically similar compounds and possibly deducing common modes of action. The TeraTox assay and the dataset described here therefore represent a new tool and a valuable resource for drug teratogenicity assessment.

Role of YAP in early ectodermal specification and a Huntington's Disease model of human neurulation

The Hippo pathway, a highly conserved signaling cascade that functions as an integrator of molecular signals and biophysical states, ultimately impinges upon the transcription coactivator Yes-associated protein 1 (YAP). Hippo-YAP signaling has been shown to play key roles both at the early embryonic stages of implantation and gastrulation, and later during neurogenesis. To explore YAP's potential role in neurulation, we used self-organizing neuruloids grown from human embryonic stem cells on micropatterned substrates. We identified YAP activation as a key lineage determinant, first between neuronal ectoderm and nonneuronal ectoderm, and later between epidermis and neural crest, indicating that YAP activity can enhance the effect of BMP4 stimulation and therefore affect ectodermal specification at this developmental stage. Because aberrant Hippo-YAP signaling has been implicated in the pathology of Huntington's Disease (HD), we used isogenic mutant neuruloids to explore the relationship between signaling and the disease. We found that HD neuruloids demonstrate ectopic activation of gene targets of YAP and that pharmacological reduction of YAP's transcriptional activity can partially rescue the HD phenotype.

Glioblastoma Embryonic-like Stem Cells Exhibit Immune-Evasive Phenotype

Simple Summary

Most glioblastoma (GBM) patients relapse after an initial response to treatment. These aggressive traits are often associated with the presence of glioma stem cells (GSCs) within the tumor bulk, which are thought to participate in GBM therapy resistance. Given GBM cellular heterogeneity, we hypothesized that GSCs might also display cellular hierarchies associated with different degrees of stemness. Based on single-cell RNAseq data from GBM patients, we identified a subpopulation of GSCs, named core-GSCs (c-GSCs), with a similar profile to embryonic stem cells and downregulation of immune-associated pathways. In addition, we developed an in vitro induced c-GSC (ic-GSC) model resembling their tumor counterpart. The characterization of immune-privileged c-GSCs provides a valuable resource to study immune evasion mechanisms in GBM and to identify potential unexplored targets to improve immunotherapy treatments.

Abstract

Background: Glioma stem cells (GSCs) have self-renewal and tumor-initiating capacities involved in drug resistance and immune evasion mechanisms in glioblastoma (GBM). Methods: Core-GSCs (c-GSCs) were identified by selecting cells co-expressing high levels of embryonic stem cell (ESC) markers from a single-cell RNA-seq patient-derived GBM dataset (n = 28). Induced c-GSCs (ic-GSCs) were generated by reprogramming GBM-derived cells (GBM-DCs) using induced pluripotent stem cell (iPSC) technology. The characterization of ic-GSCs and GBM-DCs was conducted by immunostaining, transcriptomic, and DNA methylation (DNAm) analysis. Results: We identified a GSC population (4.22% ± 0.59) exhibiting concurrent high expression of ESC markers and downregulation of immune-associated pathways, named c-GSCs. In vitro ic-GSCs presented high expression of ESC markers and downregulation of antigen presentation HLA proteins. Transcriptomic analysis revealed a strong agreement of enriched biological pathways between tumor c-GSCs and in vitro ic-GSCs (κ = 0.71). Integration of our epigenomic profiling with 833 functional ENCODE epigenetic maps identifies increased DNA methylation on HLA genes’ regulatory regions associated with polycomb repressive marks in a stem-like phenotype. Conclusions: This study unravels glioblastoma immune-evasive mechanisms involving a c-GSC population. In addition, it provides a cellular model with paired gene expression, and DNA methylation maps to explore potential therapeutic complements for GBM immunotherapy.