An integrated chromatin accessibility and transcriptome landscape of human pre-implantation embryos

DUX4 activates common and context-specific intergenic transcripts and isoforms

DUX4 regulates the expression of genic and nongenic elements and modulates chromatin accessibility during zygotic genome activation in cleavage stage embryos. Its misexpression in skeletal muscle causes facioscapulohumeral dystrophy (FSHD). By leveraging full-length RNA isoform sequencing with short-read RNA sequencing of DUX4-inducible myoblasts, we elucidate an isoform-resolved transcriptome featuring numerous unannotated isoforms from known loci and novel intergenic loci. While DUX4 activates similar programs in early embryos and FSHD muscle, the isoform usage of known DUX4 targets is notably distinct between the two contexts. DUX4 also activates hundreds of previously unannotated intergenic loci dominated by repetitive elements. The transcriptional and epigenetic profiles of these loci in myogenic and embryonic contexts indicate that the usage of DUX4-binding sites at these intergenic loci is influenced by the cellular environment. These findings demonstrate that DUX4 induces context-specific transcriptomic programs, enriching our understanding of DUX4-induced muscle pathology.

Proteogenomic reprogramming to a functional human blastomere-like stem cell state via a PARP-DUX4 regulatory axis

Here, we show that conventional human pluripotent stem cells cultured with non-specific tankyrase-PARP1-inhibited conditions underwent proteogenomic reprogramming to functional blastomere-like tankyrase/PARP inhibitor-regulated naive stem cells (TIRN-SC). TIRN-SCs concurrently expressed hundreds of pioneer factors in hybrid 2C-8C-morula-ICM programs that were augmented by induced expression of DUX4. Injection of TIRN-SCs into 8C-staged murine embryos equipotently differentiated human cells to the extra-embryonic and embryonic compartments of chimeric blastocysts and fetuses. Ectopic expression of murine-E-Cadherin in TIRN-SCs further enhanced interspecific chimeric tissue targeting. TIRN-SC-derived trophoblast stem cells efficiently generated placental chimeras. Proteome-ubiquitinome analyses revealed increased TNKS and reduced PARP1 levels and an ADP-ribosylation-deficient, hyper-ubiquitinated proteome that impacted expression of both tankyrase and PARP1 substrates. ChIP-seq of NANOG-SOX2-OCT4 and PARP1 (NSOP) revealed genome-wide NSOP co-binding at DUX4-accessible enhancers of embryonic lineage factors; suggesting a DUX4-NSOP axis regulated TIRN-SC lineage plasticity. TIRN-SCs may serve as valuable models for studying the proteogenomic regulation of pre-lineage human embryogenesis. VIDEO ABSTRACT.

Tracking single-cell evolution using clock-like chromatin accessibility loci

Single-cell chromatin accessibility sequencing (scATAC-seq) reconstructs developmental trajectory by phenotypic similarity. However, inferring the exact developmental trajectory is challenging. Previous studies showed age-associated DNA methylation (DNAm) changes in specific genomic regions, termed clock-like differential methylation loci (ClockDML). Age-associated DNAm could either result from or result in chromatin accessibility changes at ClockDML. As cells undergo mitosis, the heterogeneity of chromatin accessibility on clock-like loci is reduced, providing a measure of mitotic age. In this study, we developed a method, called EpiTrace, that counts the fraction of opened clock-like loci from scATAC-seq data to determine cell age and perform lineage tracing in various cell lineages and animal species. It shows concordance with known developmental hierarchies, correlates well with DNAm-based clocks and is complementary with mutation-based lineage tracing, RNA velocity and stemness predictions. Applying EpiTrace to scATAC-seq data reveals biological insights with clinically relevant implications, ranging from hematopoiesis, organ development, tumor biology and immunity to cortical gyrification.

Multiple Retrotransposon-mediated NF-YA Gene Duplication Events Recurred in Diverse Groups of Mammals at Different Ancestry Levels

NF-Y is a transcription factor trimer formed by Histone Fold subunits NF-YB/NF-YC and NF-YA, which confers sequence-specificity for the CCAAT box, an important cis regulatory element. The subunits are extremely conserved in all eukaryotes and in mammals they are typically encoded by single copy genes. We describe here the presence of a second NF-YA termed NF-YAr for retrogene in diverse groups of mammals (Cetacea, Ruminantia, Ursidae, Sciuridae, hippopotamus, and greater horseshoe bat). NF-YAr retrogenes are located on different chromosomes with respect to the parental gene; they are compact and intronless, or with few annotated introns. Phylogenetic and synteny analyses indicate multiple independent retrotransposition events in the different orders. Analysis of RNA-seq data of Bos taurus suggests expression confined to spermatozoa. Conservation of translation initiation signals around predicted start codons, and of 5'UTR sequences, are consistent with protein expression, suggesting that NF-YAr is a translated, retroposed NF-YA. 3D-informed structural considerations of the predicted protein sequences point at deleterious changes for CCAAT-binding and, potentially, for trimer formation. These findings indicate that multiple independent NF-YA retrotransposition events were fixed in selected orders of mammals, generating a second NF-YA with a strict tissue distribution.

Moving toward totipotency: the molecular and cellular features of totipotent and naive pluripotent stem cells

BACKGROUND

Dissecting the key molecular mechanism of embryonic development provides novel insights into embryogenesis and potential intervention strategies for clinical practices. However, the ability to study the molecular mechanisms of early embryo development in humans, such as zygotic genome activation and lineage segregation, is meaningfully constrained by methodological limitations and ethical concerns. Totipotent stem cells have an extended developmental potential to differentiate into embryonic and extraembryonic tissues, providing a suitable model for studying early embryo development. Recently, a series of ground-breaking results on stem cells have identified totipotent-like cells or induced pluripotent stem cells into totipotent-like cells.

OBJECTIVE AND RATIONALE

This review followed the PRISMA guidelines, surveys the current works of literature on totipotent, naive, and formative pluripotent stem cells, introduces the molecular and biological characteristics of those stem cells, and gives advice for future research.

SEARCH METHODS

The search method employed the terms 'totipotent' OR 'naive pluripotent stem cell' OR 'formative pluripotent stem cell' for unfiltered search on PubMed, Web of Science, and Cochrane Library. Papers included were those with information on totipotent stem cells, naive pluripotent stem cells, or formative pluripotent stem cells until June 2024 and were published in the English language. Articles that have no relevance to stem cells, or totipotent, naive pluripotent, or formative pluripotent cells were excluded.

OUTCOMES

There were 152 records included in this review. These publications were divided into four groups according to the species of the cells included in the studies: 67 human stem cell studies, 70 mouse stem cell studies, 9 porcine stem cell studies, and 6 cynomolgus stem cell studies. Naive pluripotent stem cell models have been established in other species such as porcine and cynomolgus. Human and mouse totipotent stem cells, e.g. human 8-cell-like cells, human totipotent blastomere-like cells, and mouse 2-cell-like cells, have been successfully established and exhibit high developmental potency for both embryonic and extraembryonic contributions. However, the observed discrepancies between these cells and real embryos in terms of epigenetics and transcription suggest that further research is warranted. Our results systematically reviewed the established methods, molecular characteristics, and developmental potency of different naive, formative pluripotent, and totipotent stem cells. Furthermore, we provide a parallel comparison between animal and human models, and offer recommendations for future applications to advance early embryo research and assisted reproduction technologies.

WIDER IMPLICATIONS

Totipotent cell models provide a valuable resource to understand the underlying mechanisms of embryo development and forge new paths toward future treatment of infertility and regenerative medicine. However, current in vitro cell models exhibit epigenetic and transcriptional differences from in vivo embryos, and many cell models are unstable across passages, thus imperfectly recapitulating embryonic development. In this regard, standardizing and expanding current research on totipotent stem cell models are essential to enhance our capability to resemble and decipher embryogenesis.

Maternal transcription factor OTX2 directly induces SETD1A and promotes embryonic genome activation in human pre-implantation embryos

Early embryonic development is controlled by maternal factors originating from mature oocytes. The zygotic genome is activated from a transcriptionally quiescent state through a process called embryonic genome activation (EGA), which involves the depletion and clearance of maternal factors. However, the mechanism by which maternal factors regulate EGA and embryonic development, particularly in humans, remains elusive. In this study, using tri-pronuclear (3PN) embryos and human embryonic stem cells (hESCs), we demonstrated that the maternal transcription factor Orthodenticle Homeobox 2 (OTX2), a paired-like homeobox gene, promotes EGA in human pre-implantation embryos. Knockdown of OTX2 through Trim-Away technology blocked embryonic development and minor EGA gene expression. Overexpression of OTX2 (OTX2OE) in hESCs increased transcript products, primarily at the 2-cell embryo stage genes, including genes encoding methyltransferase of histone H3K4. OTX2OE increased the level of H3K4me3 and increased the open chromatin region that co-occurs with the H3K4me3 region at the 4-cell stage in hESCs. Based on these findings in hESCs, we further verified that OTX2 directly induced the expression of SETD1A by binding to its promoter, leading to increased H3K4me3 levels in both hESCs and 3PN embryos. These findings suggest that the maternal transcription factor OTX2 regulates EGA and early embryogenesis via epigenetic mechanisms.

Comparative analysis of chromatin accessibility in porcine oocytes from follicles of different sizes

Chromatin accessibility and transcription levels during oocyte growth are important for oocyte maturation and subsequent development. However, chromatin accessibility changes in porcine oocytes during growth are unclear. The present study demonstrated that porcine oocytes derived from large follicles (LFO) exhibited higher developmental capacity than those derived from small follicles (SFO). Assay for transposase-accessible chromatin using sequencing (ATAC-seq) analysis identified 1117 and 1694 uniquely accessible chromatin peaks in LFO and SFO, respectively. Motif analysis of differential peaks revealed the top 10 significantly enriched transcription factor (TF)-binding motifs in LFO versus SFO, with only one increased peak (Spi1 binding site) and nine decreased peaks (NFYA, ATOH1, ZNF549, Foxn1, HAND2, THRB, NHLH2, FoxP1, and FoxP2 binding sites). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis identified key processes in the regulation of oocyte growth and maturation. Integration of ATAC-seq and RNA sequencing data revealed the top 10 hub genes involved in chromatin remodeling (MYSM1 and EZH2), histone modification (MYSM1, RNF2, USP1, EZH2, and MIER1), and transcription regulation (MYSM1, ASXL3, and MIER1), as well as those involved in metabolic processes and signal transduction (DOCK7, FGGY, DTL, and DNAJC6). All these genes exhibited increased expression levels in LFO versus SFO. In conclusion, the study demonstrated the dynamic nature of chromatin accessibility during porcine oocyte growth and revealed the TFs and genes closely associated with oocyte growth and maturation. These findings provide new insight into porcine oocyte growth and offer a potential strategy to enhance the in vitro developmental ability of SFO.

The role(s) of NF-Y in development and differentiation

NF-Y is a conserved sequence-specific trimeric Transcription Factor -TF- binding to the CCAAT element. We review here the role(s) in development, from pre-implantation embryo to terminally differentiated tissues, by rationalizing and commenting on genetic, genomic, epigenetic and biochemical studies. This effort brings to light the impact of NF-YA isoforms on stemness and differentiation, as well as binding to distal vs promoter proximal sites and connections with selected TFs.

The nuclear matrix stabilizes primed-specific genes in human pluripotent stem cells

The nuclear matrix, a proteinaceous gel composed of proteins and RNA, is an important nuclear structure that supports chromatin architecture, but its role in human pluripotent stem cells (hPSCs) has not been described. Here we show that by disrupting heterogeneous nuclear ribonucleoprotein U (HNRNPU) or the nuclear matrix protein, Matrin-3, primed hPSCs adopted features of the naive pluripotent state, including morphology and upregulation of naive-specific marker genes. We demonstrate that HNRNPU depletion leads to increased chromatin accessibility, reduced DNA contacts and increased nuclear size. Mechanistically, HNRNPU acts as a transcriptional co-factor that anchors promoters of primed-specific genes to the nuclear matrix with POLII to promote their expression and their RNA stability. Overall, HNRNPU promotes cell-type stability and when reduced promotes conversion to earlier embryonic states.

Quantifying Nascent Transcription in Early Embryogenesis

The early embryonic genome exists in a dormant state following fertilization, and it then subsequently undergoes broad activation of zygotic transcription at the early stages of development. A major challenge is the detection of newly made zygotic transcripts and the determination of their activation onset time due to the presence of large and predominantly maternal pool of RNAs. Here we describe a detailed method to measure the zygotic transcription during zygotic genome activation (ZGA) of Xenopus early embryos using metabolic labeling of nascent transcripts with 5-ethynyl uridine (5-EU) followed by purifying and sequencing the nascent EU-RNAs (EU-RNA-seq). This method is highly sensitive in detecting early zygotic transcripts that are not detected by total RNA-seq and determines the actual onset time of transcriptional activation for zygotic genes. The method is applicable to a wide variety of embryonic model systems and has already afforded novel insights into gene regulation in early embryogenesis.

Transposable element activity captures human pluripotent cell states

Human pluripotent stem cells (hPSCs) exist in multiple, transcriptionally distinct states and serve as powerful models for studying human development. Despite their significance, the molecular determinants and pathways governing these pluripotent states remain incompletely understood. Here, we demonstrate that transposable elements act as sensitive indicators of distinct pluripotent cell states. We engineered hPSCs with fluorescent reporters to capture the temporal expression dynamics of two state-specific transposable elements, LTR5_Hs, and MER51B. This dual reporter system enables real-time monitoring and isolation of stem cells transitioning from naïve to primed pluripotency and further towards differentiation, serving as a more accurate readout of pluripotency states compared to conventional systems. Unexpectedly, we identified a rare, metastable cell population within primed hPSCs, marked by transcripts related to preimplantation embryo development and which is associated with a DNA damage response. Moreover, our system establishes the chromatin factor NSD1 and the RNA-binding protein FUS as potent molecular safeguards of primed pluripotency. Our study introduces a novel system for investigating cellular potency and provides key insights into the regulation of embryonic development.

Chromatin accessibility reveals potential prognostic value of the peak set associated with smoking history in patients with lung adenocarcinoma

Considerable differences in molecular characteristics have been defined between non-smoker and smokers in patients with lung adenocarcinoma (LUAD), yet studies on open chromatin patterns associated with LUAD progression caused by smoking are still lacking. Here, we constructed a novel network based on correlations between each ATAC-seq peak from TCGA data using our previously developed algorithm. Subsequently, principal component analysis was performed on LUAD samples with retained peaks filtered by the correlation network, and pathway analysis was conducted to identify potential pathways involved. We identified a set of peaks that discriminated smokers in LUAD patients according to levels of exposure to tobacco quantified in pack-years. These peaks were also significantly associated with progression-free survival and overall survival of these patients. Further examination of the gene set related to those peaks revealed that the comprising genes, such as KRT19, B3GNT3, CLDN7 and CLDN3 are strongly associated with LUAD development. They are consistent with the important roles of the associated pathways in LUAD oncogenesis induced by smoking, including estrogen response, apical junction and glycolysis pathways. In summary, our study may provide valuable insights into exploring ATAC-seq peaks and understanding smoking-related LUAD carcinogenesis from a perspective of open chromatin changes.

Transcription of Endogenous Retroviruses: Broad and Precise Mechanisms of Control

Endogenous retroviruses (ERVs) are the remnants of retroviral germline infections and are highly abundant in the genomes of vertebrates. At one time considered to be nothing more than inert 'junk' within genomes, ERVs have been tolerated within host genomes over vast timescales, and their study continues to reveal complex co-evolutionary histories within their respective host species. For example, multiple instances have been characterized of ERVs having been 'borrowed' for normal physiology, from single copies to ones involved in various regulatory networks such as innate immunity and during early development. Within the cell, the accessibility of ERVs is normally tightly controlled by epigenetic mechanisms such as DNA methylation or histone modifications. However, these silencing mechanisms of ERVs are reversible, and epigenetic alterations to the chromatin landscape can thus lead to their aberrant expression, as is observed in abnormal cellular environments such as in tumors. In this review, we focus on ERV transcriptional control and draw parallels and distinctions concerning the loss of regulation in disease, as well as their precise regulation in early development.

Lactylation: the novel histone modification influence on gene expression, protein function, and disease

Lactic acid, traditionally considered as a metabolic waste product arising from glycolysis, has undergone a resurgence in scientific interest since the discovery of the Warburg effect in tumor cells. Numerous studies have proved that lactic acid could promote angiogenesis and impair the function of immune cells within tumor microenvironments. Nevertheless, the precise molecular mechanisms governing these biological functions remain inadequately understood. Recently, lactic acid has been found to induce a posttranslational modification, lactylation, that may offer insight into lactic acid's non-metabolic functions. Notably, the posttranslational modification of proteins by lactylation has emerged as a crucial mechanism by which lactate regulates cellular processes. This article provides an overview of the discovery of lactate acidification, outlines the potential "writers" and "erasers" responsible for protein lactylation, presents an overview of protein lactylation patterns across different organisms, and discusses the diverse physiological roles of lactylation. Besides, the article highlights the latest research progress concerning the regulatory functions of protein lactylation in pathological processes and underscores its scientific significance for future investigations.

Disruption of maternal vascular remodeling by a fetal endoretrovirus-derived gene in preeclampsia

BACKGROUND

Preeclampsia, one of the most lethal pregnancy-related diseases, is associated with the disruption of uterine spiral artery remodeling during placentation. However, the early molecular events leading to preeclampsia remain unknown.

RESULTS

By analyzing placentas from preeclampsia, non-preeclampsia, and twin pregnancies with selective intrauterine growth restriction, we show that the pathogenesis of preeclampsia is attributed to immature trophoblast and maldeveloped endothelial cells. Delayed epigenetic reprogramming during early extraembryonic tissue development leads to generation of excessive immature trophoblast cells. We find reduction of de novo DNA methylation in these trophoblast cells results in selective overexpression of maternally imprinted genes, including the endoretrovirus-derived gene PEG10 (paternally expressed gene 10). PEG10 forms virus-like particles, which are transferred from the trophoblast to the closely proximate endothelial cells. In normal pregnancy, only a low amount of PEG10 is transferred to maternal cells; however, in preeclampsia, excessive PEG10 disrupts maternal vascular development by inhibiting TGF-beta signaling.

CONCLUSIONS

Our study reveals the intricate epigenetic mechanisms that regulate trans-generational genetic conflict and ultimately ensure proper maternal-fetal interface formation.

The SWI/SNF ATP-dependent chromatin remodeling complex in cell lineage priming and early development

ATP dependent chromatin remodelers have pivotal roles in transcription, DNA replication and repair, and maintaining genome integrity. SWI/SNF remodelers were first discovered in yeast genetic screens for factors involved in mating type switching or for using alternative energy sources therefore termed SWI/SNF complex (short for SWItch/Sucrose NonFermentable). The SWI/SNF complexes utilize energy from ATP hydrolysis to disrupt histone-DNA interactions and shift, eject, or reposition nucleosomes making the underlying DNA more accessible to specific transcription factors and other regulatory proteins. In development, SWI/SNF orchestrates the precise activation and repression of genes at different stages, safe guards the formation of specific cell lineages and tissues. Dysregulation of SWI/SNF have been implicated in diseases such as cancer, where they can drive uncontrolled cell proliferation and tumor metastasis. Additionally, SWI/SNF defects are associated with neurodevelopmental disorders, leading to disruption of neural development and function. This review offers insights into recent developments regarding the roles of the SWI/SNF complex in pluripotency and cell lineage primining and the approaches that have helped delineate its importance. Understanding these molecular mechanisms is crucial for unraveling the intricate processes governing embryonic stem cell biology and developmental transitions and may potentially apply to human diseases linked to mutations in the SWI/SNF complex.

Uncovering uncharacterized binding of transcription factors from ATAC-seq footprinting data

Transcription factors (TFs) are crucial epigenetic regulators, which enable cells to dynamically adjust gene expression in response to environmental signals. Computational procedures like digital genomic footprinting on chromatin accessibility assays such as ATACseq can be used to identify bound TFs in a genome-wide scale. This method utilizes short regions of low accessibility signals due to steric hindrance of DNA bound proteins, called footprints (FPs), which are combined with motif databases for TF identification. However, while over 1600 TFs have been described in the human genome, only ~ 700 of these have a known binding motif. Thus, a substantial number of FPs without overlap to a known DNA motif are normally discarded from FP analysis. In addition, the FP method is restricted to organisms with a substantial number of known TF motifs. Here we present DENIS (DE Novo motIf diScovery), a framework to generate and systematically investigate the potential of de novo TF motif discovery from FPs. DENIS includes functionality (1) to isolate FPs without binding motifs, (2) to perform de novo motif generation and (3) to characterize novel motifs. Here, we show that the framework rediscovers artificially removed TF motifs, quantifies de novo motif usage during an early embryonic development example dataset, and is able to analyze and uncover TF activity in organisms lacking canonical motifs. The latter task is exemplified by an investigation of a scATAC-seq dataset in zebrafish which covers different cell types during hematopoiesis.

Endogenous retrovirus HERVH-derived lncRNA UCA1 controls human trophoblast development

Endogenous retroviruses (ERVs) are frequently reactivated in mammalian placenta. It has been proposed that ERVs contribute to shaping the gene regulatory network of mammalian trophoblasts, dominantly acting as species- and placental-specific enhancers. However, whether and how ERVs control human trophoblast development through alternative pathways remains poorly understood. Besides the well-recognized function of human endogenous retrovirus-H (HERVH) in maintaining pluripotency of early human epiblast, here we present a unique role of HERVH on trophoblast lineage development. We found that the LTR7C/HERVH subfamily exhibits an accessible chromatin state in the human trophoblast lineage. Particularly, the LTR7C/HERVH-derived Urothelial Cancer Associated 1 (UCA1), a primate-specific long non-coding RNA (lncRNA), is transcribed in human trophoblasts and promotes the proliferation of human trophoblast stem cells (hTSCs), whereas its ectopic expression compromises human trophoblast syncytialization coinciding with increased interferon signaling pathway. Importantly, UCA1 upregulation is detectable in placental samples from early-onset preeclampsia (EO-PE) patients and the transcriptome of EO-PE placenta exhibits considerable similarities to that of the syncytiotrophoblasts differentiated from UCA1-overexpressing hTSCs, supporting up-regulated UCA1 as a potential biomarker of this disease. Altogether, our data shed light on the versatile regulatory role of HERVH in early human development and provide a unique mechanism whereby ERVs exert a function in human placentation and placental syndromes.

Interaction network of human early embryonic transcription factors

Embryonic genome activation (EGA) occurs during preimplantation development and is characterized by the initiation of de novo transcription from the embryonic genome. Despite its importance, the regulation of EGA and the transcription factors involved in this process are poorly understood. Paired-like homeobox (PRDL) family proteins are implicated as potential transcriptional regulators of EGA, yet the PRDL-mediated gene regulatory networks remain uncharacterized. To investigate the function of PRDL proteins, we are identifying the molecular interactions and the functions of a subset family of the Eutherian Totipotent Cell Homeobox (ETCHbox) proteins, seven PRDL family proteins and six other transcription factors (TFs), all suggested to participate in transcriptional regulation during preimplantation. Using mass spectrometry-based interactomics methods, AP-MS and proximity-dependent biotin labeling, and chromatin immunoprecipitation sequencing we derive the comprehensive regulatory networks of these preimplantation TFs. By these interactomics tools we identify more than a thousand high-confidence interactions for the 21 studied bait proteins with more than 300 interacting proteins. We also establish that TPRX2, currently assigned as pseudogene, is a transcriptional activator.

Genome-wide ATAC-see screening identifies TFDP1 as a modulator of global chromatin accessibility

Chromatin accessibility is a hallmark of active regulatory regions and is functionally linked to transcriptional networks and cell identity. However, the molecular mechanisms and networks that govern chromatin accessibility have not been thoroughly studied. Here we conducted a genome-wide CRISPR screening combined with an optimized ATAC-see protocol to identify genes that modulate global chromatin accessibility. In addition to known chromatin regulators like CREBBP and EP400, we discovered a number of previously unrecognized proteins that modulate chromatin accessibility, including TFDP1, HNRNPU, EIF3D and THAP11 belonging to diverse biological pathways. ATAC-seq analysis upon their knockouts revealed their distinct and specific effects on chromatin accessibility. Remarkably, we found that TFDP1, a transcription factor, modulates global chromatin accessibility through transcriptional regulation of canonical histones. In addition, our findings highlight the manipulation of chromatin accessibility as an approach to enhance various cell engineering applications, including genome editing and induced pluripotent stem cell reprogramming.

Expression and function of NF-Y subunits in cancer

NF-Y is a Transcription Factor (TF) targeting the CCAAT box regulatory element. It consists of the NF-YB/NF-YC heterodimer, each containing an Histone Fold Domain (HFD), and the sequence-specific subunit NF-YA. NF-YA expression is associated with cell proliferation and absent in some post-mitotic cells. The review summarizes recent findings impacting on cancer development. The logic of the NF-Y regulome points to pro-growth, oncogenic genes in the cell-cycle, metabolism and transcriptional regulation routes. NF-YA is involved in growth/differentiation decisions upon cell-cycle re-entry after mitosis and it is widely overexpressed in tumors, the HFD subunits in some tumor types or subtypes. Overexpression of NF-Y -mostly NF-YA- is oncogenic and decreases sensitivity to anti-neoplastic drugs. The specific roles of NF-YA and NF-YC isoforms generated by alternative splicing -AS- are discussed, including the prognostic value of their levels, although the specific molecular mechanisms of activity are still to be deciphered.

Species-Specific Transcription Factors Associated with Long Terminal Repeat Promoters of Endogenous Retroviruses: A Comprehensive Review

Endogenous retroviruses (ERVs) became a part of the eukaryotic genome through endogenization millions of years ago. Moreover, they have lost their innate capability of virulence or replication. Nevertheless, in eukaryotic cells, they actively engage in various activities that may be advantageous or disadvantageous to the cells. The mechanisms by which transcription is triggered and implicated in cellular processes are complex. Owing to the diversity in the expression of transcription factors (TFs) in cells and the TF-binding motifs of viruses, the comprehensibility of ERV initiation and its impact on cellular functions are unclear. Currently, several factors are known to be related to their initiation. TFs that bind to the viral long-terminal repeat (LTR) are critical initiators. This review discusses the TFs shown to actively associate with ERV stimulation across species such as humans, mice, pigs, monkeys, zebrafish, Drosophila, and yeast. A comprehensive summary of the expression of previously reported TFs may aid in identifying similarities between animal species and endogenous viruses. Moreover, an in-depth understanding of ERV expression will assist in elucidating their physiological roles in eukaryotic cell development and in clarifying their relationship with endogenous retrovirus-associated diseases.

Single-cell multi-omics profiling of human preimplantation embryos identifies cytoskeletal defects during embryonic arrest

Human in vitro fertilized embryos exhibit low developmental capabilities, and the mechanisms that underlie embryonic arrest remain unclear. Here using a single-cell multi-omics sequencing approach, we simultaneously analysed alterations in the transcriptome, chromatin accessibility and the DNA methylome in human embryonic arrest due to unexplained reasons. Arrested embryos displayed transcriptome disorders, including a distorted microtubule cytoskeleton, increased genomic instability and impaired glycolysis, which were coordinated with multiple epigenetic reprogramming defects. We identified Aurora A kinase (AURKA) repression as a cause of embryonic arrest. Mechanistically, arrested embryos induced through AURKA inhibition resembled the reprogramming abnormalities of natural embryonic arrest in terms of the transcriptome, the DNA methylome, chromatin accessibility and H3K4me3 modifications. Mitosis-independent sequential activation of the zygotic genome in arrested embryos showed that YY1 contributed to human major zygotic genome activation. Collectively, our study decodes the reprogramming abnormalities and mechanisms of human embryonic arrest and the key regulators of zygotic genome activation.

The transcriptional activator Klf5 recruits p300-mediated H3K27ac for maintaining trophoblast stem cell pluripotency

The effective proliferation and differentiation of trophoblast stem cells (TSCs) is indispensable for the development of the placenta, which is the key to maintaining normal fetal growth during pregnancy. Kruppel-like factor 5 (Klf5) is implicated in the activation of pluripotency gene expression in embryonic stem cells (ESCs), yet its function in TSCs is poorly understood. Here, we showed that Klf5 knockdown resulted in the downregulation of core TSC-specific genes, consequently causing rapid differentiation of TSCs. Consistently, Klf5-depleted embryos lost the ability to establish TSCs in vitro. At the molecular level, Klf5 preferentially occupied the proximal promoter regions and maintained an open chromatin architecture of key TSC-specific genes. Deprivation of Klf5 impaired the enrichment of p300, a major histone acetyl transferase of H3 lysine 27 acetylation (H3K27ac), and further reduced the occupancy of H3K27ac at promoter regions, leading to decreased transcriptional activity of TSC pluripotency genes. Thus, our findings highlight a novel mechanism of Klf5 in regulating the self-renewal and differentiation of TSCs and provide a reference for understanding placental development and improving pregnancy rates.

MICA: a multi-omics method to predict gene regulatory networks in early human embryos

Recent advances in single-cell omics have transformed characterisation of cell types in challenging-to-study biological contexts. In contexts with limited single-cell samples, such as the early human embryo inference of transcription factor-gene regulatory network (GRN) interactions is especially difficult. Here, we assessed application of different linear or non-linear GRN predictions to single-cell simulated and human embryo transcriptome datasets. We also compared how expression normalisation impacts on GRN predictions, finding that transcripts per million reads outperformed alternative methods. GRN inferences were more reproducible using a non-linear method based on mutual information (MI) applied to single-cell transcriptome datasets refined with chromatin accessibility (CA) (called MICA), compared with alternative network prediction methods tested. MICA captures complex non-monotonic dependencies and feedback loops. Using MICA, we generated the first GRN inferences in early human development. MICA predicted co-localisation of the AP-1 transcription factor subunit proto-oncogene JUND and the TFAP2C transcription factor AP-2γ in early human embryos. Overall, our comparative analysis of GRN prediction methods defines a pipeline that can be applied to single-cell multi-omics datasets in especially challenging contexts to infer interactions between transcription factor expression and target gene regulation.

Epigenetic regulation of early human embryo development

Studies of mammalian development have advanced our understanding of the genetic, epigenetic, and cellular processes that orchestrate embryogenesis and have uncovered new insights into the unique aspects of human embryogenesis. Recent studies have now produced the first epigenetic maps of early human embryogenesis, stimulating new ideas about epigenetic reprogramming, cell fate control, and the potential mechanisms underpinning developmental plasticity in human embryos. In this review, we discuss these new insights into the epigenetic regulation of early human development and the importance of these processes for safeguarding development. We also highlight unanswered questions and key challenges that remain to be addressed.

Isoform-resolved transcriptome of the human preimplantation embryo

Human preimplantation development involves extensive remodeling of RNA expression and splicing. However, its transcriptome has been compiled using short-read sequencing data, which fails to capture most full-length mRNAs. Here, we generate an isoform-resolved transcriptome of early human development by performing long- and short-read RNA sequencing on 73 embryos spanning the zygote to blastocyst stages. We identify 110,212 unannotated isoforms transcribed from known genes, including highly conserved protein-coding loci and key developmental regulators. We further identify 17,964 isoforms from 5,239 unannotated genes, which are largely non-coding, primate-specific, and highly associated with transposable elements. These isoforms are widely supported by the integration of published multi-omics datasets, including single-cell 8CLC and blastoid studies. Alternative splicing and gene co-expression network analyses further reveal that embryonic genome activation is associated with splicing disruption and transient upregulation of gene modules. Together, these findings show that the human embryo transcriptome is far more complex than currently known, and will act as a valuable resource to empower future studies exploring development.

The composition dynamics of transposable elements in human blastocysts

Transposable elements (TEs) are mobile DNA sequences that can replicate themselves and play significant roles in embryo development and chromosomal structure remodeling. In this study, we investigated the variation of TEs in blastocysts with different parental genetic backgrounds. We analyzed the proportions of 1137 TEs subfamilies from six classes at the DNA level using Bowtie2 and PopoolationTE2 in 196 blastocysts with abnormal parental chromosomal diseases. Our findings revealed that the parental karyotype was the dominant factor influencing TEs frequencies. Out of the 1116 subfamilies, different frequencies were observed in blastocysts with varying parental karyotypes. The development stage of blastocysts was the second most crucial factor influencing TEs proportions. A total of 614 subfamilies exhibited different proportions at distinct blastocyst stages. Notably, subfamily members belonging to the Alu family showed a high proportion at stage 6, while those from the LINE class exhibited a high proportion at stage 3 and a low proportion at stage 6. Moreover, the proportions of some TEs subfamilies also varied depending on blastocyst karyotype, inner cell mass status, and outer trophectoderm status. We found that 48 subfamilies displayed different proportions between balanced and unbalanced blastocysts. Additionally, 19 subfamilies demonstrated varying proportions among different inner cell mass scores, and 43 subfamilies exhibited different proportions among outer trophectoderm scores. This study suggests that the composition of TEs subfamilies may be influenced by various factors and undergoes dynamic modulation during embryo development.

Genome-wide high-throughput chromosome conformation capture analysis reveals hierarchical chromatin interactions during early somatic embryogenesis

Somatic embryogenesis (SE), like zygotic embryo development, is a progressive process. Early SE is the beginning of a switch from a somatic to an embryogenic state and is an important stage for initiating chromatin reprogramming of SE. Previous studies suggest that changes in chromatin accessibility occur during early SE, although information on the 3D structure of chromatin is not yet available. Here, we present a chromosome-level genome assembly of longan (Dimocarpus longan) using PacBio combined with high-through chromosome conformation capture scaffolding, which resulted in a 446 Mb genome assembly anchored onto 15 scaffolds. During early SE, chromatin was concentrated and then decondensed, and a large number of long terminal repeat retrotransposons (LTR-RTs) were enriched in the local chromatin interaction region, suggesting LTR-RTs were involved in chromatin reorganization. Early SE was accompanied by the transformation from A to B compartments, and the interactions between B compartments were enhanced. Results from chromatin accessibility, monomethylation of histone H3 at lysine 4 (H3K4me1) modification, and transcription analyses further revealed a gene regulatory network for cell wall thickening during SE. Particularly, we found that the H3K4me1 differential peak binding motif showed abnormal activation of ethylene response factor transcription factors and participation in SE. The chromosome-level genomic and multiomics analyses revealed the 3D conformation of chromatin during early SE, providing insight into the molecular mechanisms underlying cell wall thickening and the potential regulatory networks of TFs during early SE in D. longan. These results provide additional clues for revealing the molecular mechanisms of plant SE.

Activation of lnc-ALVE1-AS1 inhibited ALV-J replication through triggering the TLR3 pathway in chicken macrophage like cell line

Endogenous retroviruses (ERVs) are remnants of the historical retroviral infections, and their derived transcripts with viral signatures are important sources of long noncoding RNAs (lncRNAs). We have previously shown that the chicken ERV-derived lncRNA lnc-ALVE1-AS1 exerts antiviral innate immunity in chicken embryo fibroblasts. However, it is not clear whether this endogenous retroviral RNA has a similar function in immune cells. Here, we found that lnc-ALVE1-AS1 was persistently inhibited in chicken macrophages after avian leukosis virus subgroup J (ALV-J) infection. Furthermore, overexpression of lnc-ALVE1-AS1 significantly inhibited the replication of exogenous ALV-J, whereas knockdown of lnc-ALVE1-AS1 promoted the replication of ALV-J in chicken macrophages. This phenomenon is attributed to the induction of antiviral innate immunity by lnc-ALVE1-AS1 in macrophages, whereas knockdown of lnc-ALVE1-AS1 had the opposite effect. Mechanistically, lnc-ALVE1-AS1 can be sensed by the cytosolic pattern recognition receptor TLR3 and trigger the type I interferons response. The present study provides novel insights into the antiviral defense of ERV-derived lncRNAs in macrophages and offers new strategies for future antiviral solutions.

Selective Translation of Maternal mRNA by eIF4E1B Controls Oocyte to Embryo Transition

Maternal messenger ribonucleic acids (mRNAs) are driven by a highly orchestrated scheme of recruitment to polysomes and translational activation. However, selecting and regulating individual mRNAs for the translation from a competitive pool of mRNAs are little-known processes. This research shows that the maternal eukaryotic translation initiation factor 4e1b (Eif4e1b) expresses during the oocyte-to-embryo transition (OET), and maternal deletion of Eif4e1b leads to multiple defects concerning oogenesis and embryonic developmental competence during OET. The linear amplification of complementary deoxyribonucleic acid (cDNA) ends, and sequencing (LACE-seq) is used to identify the distinct subset of mRNA and its CG-rich binding sites within the 5' untranslated region (UTR) targeted by eIF4E1B. The proteomics analyses indicate that eIF4E1B-specific bound genes show stronger downregulation at the protein level, which further verify a group of proteins that plays a crucial role in oocyte maturation and embryonic developmental competence is insufficiently synthesized in Eif4e1b-cKO oocytes during OET. Moreover, the biochemical results in vitro are combined to further confirm the maternal-specific translation activation model assembled by eIF4E1B and 3'UTR-associated mRNA binding proteins. The findings demonstrate the indispensability of eIF4E1B for selective translation activation in mammalian oocytes and provide a potential network regulated by eIF4E1B in OET.

Comprehensive characterization of the embryonic factor LEUTX

The paired-like homeobox transcription factor LEUTX is expressed in human preimplantation embryos between the 4- and 8-cell stages, and then silenced in somatic tissues. To characterize the function of LEUTX, we performed a multiomic characterization of LEUTX using two proteomics methods and three genome-wide sequencing approaches. Our results show that LEUTX stably interacts with the EP300 and CBP histone acetyltransferases through its 9 amino acid transactivation domain (9aaTAD), as mutation of this domain abolishes the interactions. LEUTX targets genomic cis-regulatory sequences that overlap with repetitive elements, and through these elements it is suggested to regulate the expression of its downstream genes. We find LEUTX to be a transcriptional activator, upregulating several genes linked to preimplantation development as well as 8-cell-like markers, such as DPPA3 and ZNF280A. Our results support a role for LEUTX in preimplantation development as an enhancer binding protein and as a potent transcriptional activator.

Histone remodeling reflects conserved mechanisms of bovine and human preimplantation development

How histone modifications regulate changes in gene expression during preimplantation development in any species remains poorly understood. Using CUT&Tag to overcome limiting amounts of biological material, we profiled two activating (H3K4me3 and H3K27ac) and two repressive (H3K9me3 and H3K27me3) marks in bovine oocytes, 2-, 4-, and 8-cell embryos, morula, blastocysts, inner cell mass, and trophectoderm. In oocytes, broad bivalent domains mark developmental genes, and prior to embryonic genome activation (EGA), H3K9me3 and H3K27me3 co-occupy gene bodies, suggesting a global mechanism for transcription repression. During EGA, chromatin accessibility is established before canonical H3K4me3 and H3K27ac signatures. Embryonic transcription is required for this remodeling, indicating that maternally provided products alone are insufficient for reprogramming. Last, H3K27me3 plays a major role in restriction of cellular potency, as blastocyst lineages are defined by differential polycomb repression and transcription factor activity. Notably, inferred regulators of EGA and blastocyst formation strongly resemble those described in humans, as opposed to mice. These similarities suggest that cattle are a better model than rodents to investigate the molecular basis of human preimplantation development.

Single-cell multi-omic analysis profiles defective genome activation and epigenetic reprogramming associated with human pre-implantation embryo arrest

During pre-implantation stages of mammalian development, maternally stored material promotes both the erasure of the sperm and oocyte epigenetic profiles and is responsible for concomitant genome activation. Here, we have utilized single-cell methylome and transcriptome sequencing (scM&T-seq) to quantify both mRNA expression and DNA methylation in oocytes and a developmental series of human embryos at single-cell resolution. We fully characterize embryonic genome activation and maternal transcript degradation and map key epigenetic reprogramming events in developmentally high-quality embryos. By comparing these signatures with early embryos that have undergone spontaneous cleavage-stage arrest, as determined by time-lapse imaging, we identify embryos that fail to appropriately activate their genomes or undergo epigenetic reprogramming. Our results indicate that a failure to successfully accomplish these essential milestones impedes the developmental potential of pre-implantation embryos and is likely to have important implications, similar to aneuploidy, for the success of assisted reproductive cycles.

Telomeres cooperate in zygotic genome activation by affecting DUX4/Dux transcription

Zygotic genome activation (ZGA) is initiated once the genome chromatin state is organized in the newly formed zygote. Telomeres are specialized chromatin structures at the ends of chromosomes and are reset during early embryogenesis, while the details and significance of telomere changes in preimplantation embryos remain unclear. We demonstrated that the telomere length was shortened in the minor ZGA stage and significantly elongated in the major ZGA stage of human and mouse embryos. Expression of the ZGA pioneer factor DUX4/Dux was negatively correlated with the telomere length. ATAC sequencing data revealed that the chromatin accessibility peaks on the DUX4 promoter region (i.e., the subtelomere of chromosome 4q) were transiently augmented in human minor ZGA. Reduction of telomeric heterochromatin H3K9me3 in the telomeric region also synergistically activated DUX4 expression with p53 in human embryonic stem cells. We propose herein that telomeres regulate the expression of DUX4/Dux through chromatin remodeling and are thereby involved in ZGA.

Revealing cell populations catching the early stages of human embryo development in naive pluripotent stem cell cultures

Naive human pluripotent stem cells (hPSCs) are defined as the in vitro counterpart of the human preimplantation embryo's epiblast and are used as a model system to study developmental processes. In this study, we report the discovery and characterization of distinct cell populations coexisting with epiblast-like cells in 5iLAF naive human induced PSC (hiPSC) cultures. It is noteworthy that these populations closely resemble different cell types of the human embryo at early developmental stages. While epiblast-like cells represent the main cell population, interestingly we detect a cell population with gene and transposable element expression profile closely resembling the totipotent eight-cell (8C)-stage human embryo, and three cell populations analogous to trophectoderm cells at different stages of their maturation process: transition, early, and mature stages. Moreover, we reveal the presence of cells resembling primitive endoderm. Thus, 5iLAF naive hiPSC cultures provide an excellent opportunity to model the earliest events of human embryogenesis, from the 8C stage to the peri-implantation period.

Profiling the transcriptomic signatures and identifying the patterns of zygotic genome activation - a comparative analysis between early porcine embryos and their counterparts in other three mammalian species

BACKGROUND

The transcriptional changes around zygotic genome activation (ZGA) in preimplantation embryos are critical for studying mechanisms of embryonic developmental arrest and searching for key transcription factors. However, studies on the transcription profile of porcine ZGA are limited.

RESULTS

In this study, we performed RNA sequencing in porcine in vivo developed (IVV) and somatic cell nuclear transfer (SCNT) embryo at different stages and compared the transcriptional activity of porcine embryos with mouse, bovine and human embryos. The results showed that the transcriptome map of the early porcine embryos was significantly changed at the 4-cell stage, and 5821 differentially expressed genes (DEGs) in SCNT embryos failed to be reprogrammed or activated during ZGA, which mainly enrichment to metabolic pathways. c-MYC was identified as the highest expressed transcription factor during ZGA. By treating with 10,058-F4, an inhibitor of c-MYC, the cleavage rate (38.33 ± 3.4%) and blastocyst rate (23.33 ± 4.3%) of porcine embryos were significantly lower than those of the control group (50.82 ± 2.7% and 34.43 ± 1.9%). Cross-species analysis of transcriptome during ZGA showed that pigs and bovines had the highest similarity coefficient in biological processes. KEGG pathway analysis indicated that there were 10 co-shared pathways in the four species.

CONCLUSIONS

Our results reveal that embryos with impaired developmental competence may be arrested at an early stage of development. c-MYC helps promote ZGA and preimplantation embryonic development in pigs. Pigs and bovines have the highest coefficient of similarity in biological processes during ZGA. This study provides an important reference for further studying the reprogramming regulatory mechanism of porcine embryos during ZGA.

Primate-specific transposable elements shape transcriptional networks during human development

The human genome contains more than 4.5 million inserts derived from transposable elements (TEs), the result of recurrent waves of invasion and internal propagation throughout evolution. For new TE copies to be inherited, they must become integrated in the genome of the germline or pre-implantation embryo, which requires that their source TE be expressed at these stages. Accordingly, many TEs harbor DNA binding sites for the pluripotency factors OCT4, NANOG, SOX2, and KLFs and are transiently expressed during embryonic genome activation. Here, we describe how many primate-restricted TEs have additional binding sites for lineage-specific transcription factors driving their expression during human gastrulation and later steps of fetal development. These TE integrants serve as lineage-specific enhancers fostering the transcription, amongst other targets, of KRAB-zinc finger proteins (KZFPs) of comparable evolutionary age, which in turn corral the activity of TE-embedded regulatory sequences in a similarly lineage-restricted fashion. Thus, TEs and their KZFP controllers play broad roles in shaping transcriptional networks during early human development.

The single-cell expression profile of transposable elements and transcription factors in human early biparental and uniparental embryonic development

Transposable elements (TEs) and transcription factors (TFs) are involved in the precise regulation of gene expression during the preimplantation stage. Activation of TEs is a key event for mammalian embryonic genome activation and preimplantation early embryonic development. TFs are involved in the regulation of drastic changes in gene expression patterns, but an inventory of the interplay between TEs and TFs during normal/abnormal human embryonic development is still lacking. Here we used single-cell RNA sequencing data generated from biparental and uniparental embryos to perform an integrative analysis of TE and TF expression. Our results showed that endogenous retroviruses (ERVs) are mainly expressed during the minor embryonic genome activation (EGA) process of early embryos, while Alu is gradually expressed in the middle and later stages. Some important ERVs (e.g., LTR5_Hs, MLT2A1) and Alu TEs are expressed at significantly lower levels in androgenic embryos. Integrative analysis revealed that the expression of the transcription factors CTCF and POU5F1 is correlated with the differential expression of ERV TEs. Comparative coexpression network analysis further showed distinct expression levels of important TFs (e.g., LEUTX and ZSCAN5A) in dizygotic embryos vs. parthenogenetic and androgenic embryos. This systematic investigation of TE and TF expression in human early embryonic development by single-cell RNA sequencing provides valuable insights into mammalian embryonic development.

The non-coding genome in early human development - Recent advancements

Not that long ago, the human genome was discovered to be mainly non-coding, that is comprised of DNA sequences that do not code for proteins. The initial paradigm that non-coding is also non-functional was soon overturned and today the work to uncover the functions of non-coding DNA and RNA in human early embryogenesis has commenced. Early human development is characterized by large-scale changes in genomic activity and the transcriptome that are partly driven by the coordinated activation and repression of repetitive DNA elements scattered across the genome. Here we provide examples of recent novel discoveries of non-coding DNA and RNA interactions and mechanisms that ensure accurate non-coding activity during human maternal-to-zygotic transition and lineage segregation. These include studies on small and long non-coding RNAs, transposable element regulation, and RNA tailing in human oocytes and early embryos. High-throughput approaches to dissect the non-coding regulatory networks governing early human development are a foundation for functional studies of specific genomic elements and molecules that has only begun and will provide a wider understanding of early human embryogenesis and causes of infertility.

Dynamic reprogramming of H3K9me3 at hominoid-specific retrotransposons during human preimplantation development

Reprogramming of H3K9me3-dependent heterochromatin is required for early development. How H3K9me3 is involved in early human development remains, however, largely unclear. Here, we resolve the temporal landscape of H3K9me3 during human preimplantation development and its regulation for diverse hominoid-specific retrotransposons. At the 8-cell stage, H3K9me3 reprogramming at hominoid-specific retrotransposons termed SINE-VNTR-Alu (SVA) facilitates interaction between certain promoters and SVA-derived enhancers, promoting the zygotic genome activation. In trophectoderm, de novo H3K9me3 domains prevent pluripotent transcription factors from binding to hominoid-specific retrotransposons-derived regulatory elements for inner cell mass (ICM)-specific genes. H3K9me3 re-establishment at SVA elements in the ICM is associated with higher transcription of DNA repair genes, when compared with naive human pluripotent stem cells. Our data demonstrate that species-specific reorganization of H3K9me3-dependent heterochromatin at hominoid-specific retrotransposons plays important roles during early human development, shedding light on how the epigenetic regulation for early development has evolved in mammals.

Advances in studying human gametogenesis and embryonic development in China

Reproductive medicine in China has developed rapidly since 1988 due to the support from the government and scientific exploration. However, the success rate of assisted reproduction technology (ART) is around 30-40% and many unknown "black boxes" in gametogenesis and embryo development are still present. With the development of single-cell and low-input sequencing technologies, the network of transcriptome and epigenetic regulation (DNA methylation, chromatin accessibility, and histone modifications) during the development of human primordial germ cells (PGCs), gametes and embryos has been investigated in depth. Furthermore, pre-implantation genetic testing (PGT) has also rapidly developed. In this review, we summarize and analyze China's outstanding progress in these fields.

Active endogenous retroviral elements in human pluripotent stem cells play a role in regulating host gene expression

Human endogenous retroviruses, also called LTR elements, can be bound by transcription factors and marked by different histone modifications in different biological contexts. Recently, individual LTR or certain subclasses of LTRs such as LTR7/HERVH and LTR5_Hs/HERVK families have been identified as cis-regulatory elements. However, there are still many LTR elements with unknown functions. Here, we dissected the landscape of histone modifications and regulatory map of LTRs by integrating 98 ChIP-seq data in human embryonic stem cells (ESCs), and annotated the active LTRs enriching enhancer/promoter-related histone marks. Notably, we found that MER57E3 functionally acted as proximal regulatory element to activate respective ZNF gene. Additionally, HERVK transcript could mainly function in nucleus to activate the adjacent genes. Since LTR5_Hs/LTR5 was bound by many early embryo-specific transcription factors, we further investigated the expression dynamics in different pluripotent states. LTR5_Hs/LTR5/HERVK exhibited higher expression level in naïve ESCs and extended pluripotent stem cells (EPSCs). Functionally, the LTR5_Hs/LTR5 with high activity could serve as a distal enhancer to regulate the host genes. Ultimately, our study not only provides a comprehensive regulatory map of LTRs in human ESCs, but also explores the regulatory models of MER57E3 and LTR5_Hs/LTR5 in host genome.

Sperm histone H3 lysine 4 tri-methylation serves as a metabolic sensor of paternal obesity and is associated with the inheritance of metabolic dysfunction

OBJECTIVE

Parental environmental exposures can strongly influence descendant risks for adult disease. How paternal obesity changes the sperm chromatin leading to the acquisition of metabolic disease in offspring remains controversial and ill-defined. The objective of this study was to assess (1) whether obesity induced by a high-fat diet alters sperm histone methylation; (2) whether paternal obesity can induce metabolic disturbances across generations; (3) whether there could be cumulative damage to the sperm epigenome leading to enhanced metabolic dysfunction in descendants; and (4) whether obesity-sensitive regions associate with embryonic epigenetic and transcriptomic profiles. Using a genetic mouse model of epigenetic inheritance, we investigated the role of histone H3 lysine 4 methylation (H3K4me3) in the paternal transmission of metabolic dysfunction. This transgenic mouse overexpresses the histone demethylase enzyme KDM1A in the developing germline and has an altered sperm epigenome at the level of histone H3K4 methylation. We hypothesized that challenging transgenic sires with a high-fat diet would further erode the sperm epigenome and lead to enhanced metabolic disturbances in the next generations.

METHODS

To assess whether paternal obesity can have inter- or transgenerational impacts, and if so to identify potential mechanisms of this non-genetic inheritance, we used wild-type C57BL/6NCrl and transgenic males with a pre-existing altered sperm epigenome. To induce obesity, sires were fed either a control or high-fat diet (10% or 60% kcal fat, respectively) for 10-12 weeks, then bred to wild-type C57BL/6NCrl females fed a regular diet. F1 and F2 descendants were characterized for metabolic phenotypes by examining the effects of paternal obesity by sex, on body weight, fat mass distribution, the liver transcriptome, intraperitoneal glucose, and insulin tolerance tests. To determine whether obesity altered the F0 sperm chromatin, native chromatin immunoprecipitation-sequencing targeting H3K4me3 was performed. To gain insight into mechanisms of paternal transmission, we compared our sperm H3K4me3 profiles with embryonic and placental chromatin states, histone modification, and gene expression profiles.

RESULTS

Obesity-induced alterations in H3K4me3 occurred in genes implicated in metabolic, inflammatory, and developmental processes. These processes were associated with offspring metabolic dysfunction and corresponded to genes enriched for H3K4me3 in embryos and overlapped embryonic and placenta gene expression profiles. Transgenerational susceptibility to metabolic disease was only observed when obese F0 had a pre-existing modified sperm epigenome. This coincided with increased H3K4me3 alterations in sperm and more severe phenotypes affecting their offspring.

CONCLUSIONS

Our data suggest sperm H3K4me3 might serve as a metabolic sensor that connects paternal diet with offspring phenotypes via the placenta. This non-DNA-based knowledge of inheritance has the potential to improve our understanding of how environment shapes heritability and may lead to novel routes for the prevention of disease. This study highlights the need to further study the connection between the sperm epigenome, placental development, and children's health.

SUMMARY SENTENCE

Paternal obesity impacts sperm H3K4me3 and is associated with placenta, embryonic and metabolic outcomes in descendants.

Rolling back human pluripotent stem cells to an eight-cell embryo-like stage

After fertilization, the quiescent zygote experiences a burst of genome activation that initiates a short-lived totipotent state. Understanding the process of totipotency in human cells would have broad applications. However, in contrast to in mice1,2, demonstration of the time of zygotic genome activation or the eight-cell (8C) stage in in vitro cultured human cells has not yet been reported, and the study of embryos is limited by ethical and practical considerations. Here we describe a transgene-free, rapid and controllable method for producing 8C-like cells (8CLCs) from human pluripotent stem cells. Single-cell analysis identified key molecular events and gene networks associated with this conversion. Loss-of-function experiments identified fundamental roles for DPPA3, a master regulator of DNA methylation in oocytes3, and TPRX1, a eutherian totipotent cell homeobox (ETCHbox) family transcription factor that is absent in mice4. DPPA3 induces DNA demethylation throughout the 8CLC conversion process, whereas TPRX1 is a key executor of 8CLC gene networks. We further demonstrate that 8CLCs can produce embryonic and extraembryonic lineages in vitro or in vivo in the form of blastoids5 and complex teratomas. Our approach provides a resource to uncover the molecular process of early human embryogenesis.

DUX4 is a multifunctional factor priming human embryonic genome activation

Double homeobox 4 (DUX4) is expressed at the early pre-implantation stage in human embryos. Here we show that induced human DUX4 expression substantially alters the chromatin accessibility of non-coding DNA and activates thousands of newly identified transcribed enhancer-like regions, preferentially located within ERVL-MaLR repeat elements. CRISPR activation of transcribed enhancers by C-terminal DUX4 motifs results in the increased expression of target embryonic genome activation (EGA) genes ZSCAN4 and KHDC1P1. We show that DUX4 is markedly enriched in human zygotes, followed by intense nuclear DUX4 localization preceding and coinciding with minor EGA. DUX4 knockdown in human zygotes led to changes in the EGA transcriptome but did not terminate the embryos. We also show that the DUX4 protein interacts with the Mediator complex via the C-terminal KIX binding motif. Our findings contribute to the understanding of DUX4 as a regulator of the non-coding genome.

8C-like cells capture the human zygotic genome activation program in vitro

The activation of the embryonic genome marks the first major wave of transcription in the developing organism. Zygotic genome activation (ZGA) in mouse 2-cell embryos and 8-cell embryos in humans is crucial for development. Here, we report the discovery of human 8-cell-like cells (8CLCs) among naive embryonic stem cells, which transcriptionally resemble the 8-cell human embryo. They express ZGA markers, including ZSCAN4 and LEUTX, and transposable elements, such as HERVL and MLT2A1. 8CLCs show reduced SOX2 levels and can be identified using TPRX1 and H3.Y marker proteins in vitro. Overexpression of the transcription factor DUX4 and spliceosome inhibition increase human ZGA-like transcription. Excitingly, the 8CLC markers TPRX1 and H3.Y are also expressed in ZGA-stage 8-cell human embryos and may thus be relevant in vivo. 8CLCs provide a unique opportunity to characterize human ZGA-like transcription and might provide critical insights into early events in embryogenesis in humans.

Application of ATAC-Seq for genome-wide analysis of the chromatin state at single myofiber resolution

Myofibers are the main components of skeletal muscle, which is the largest tissue in the body. Myofibers are highly adaptive and can be altered under different biological and disease conditions. Therefore, transcriptional and epigenetic studies on myofibers are crucial to discover how chromatin alterations occur in the skeletal muscle under different conditions. However, due to the heterogenous nature of skeletal muscle, studying myofibers in isolation proves to be a challenging task. Single-cell sequencing has permitted the study of the epigenome of isolated myonuclei. While this provides sequencing with high dimensionality, the sequencing depth is lacking, which makes comparisons between different biological conditions difficult. Here, we report the first implementation of single myofiber ATAC-Seq, which allows for the sequencing of an individual myofiber at a depth sufficient for peak calling and for comparative analysis of chromatin accessibility under various physiological and disease conditions. Application of this technique revealed significant differences in chromatin accessibility between resting and regenerating myofibers, as well as between myofibers from a mouse model of Duchenne Muscular Dystrophy (mdx) and wild-type (WT) counterparts. This technique can lead to a wide application in the identification of chromatin regulatory elements and epigenetic mechanisms in muscle fibers during development and in muscle-wasting diseases.

Downstream events initiated by expression of FSHD-associated DUX4: Studies of nucleocytoplasmic transport, γH2AX accumulation, and Bax/Bak-dependence

Abnormal expression in skeletal muscle of the double homeobox transcription factor DUX4 underlies pathogenesis in facioscapulohumeral muscular dystrophy (FSHD). Though multiple changes are known to be initiated by aberrant DUX4 expression, the downstream events initiated by DUX4 remain incompletely understood. In this study, we examined plausible downstream events initiated by DUX4. First, we found that nucleocytoplasmic protein export appeared to be decreased upon DUX4 expression as indicated by nuclear accumulation of a shuttle-GFP reporter. Second, building on studies from other labs, we showed that phospho(Ser139)-H2AX (γH2AX), an indicator of double-strand DNA breaks, accumulated both in human FSHD1 myotube nuclei upon endogenous DUX4 expression and in Bax-/-;Bak-/- (double knockout), SV40-immortalized mouse embryonic fibroblasts upon exogenous DUX4 expression. In contrast, DUX4-induced caspase 3/7 activation was prevented in Bax-/-;Bak-/- double knockout SV40-MEFs, but not by single knockouts of Bax, Bak, or Bid. Thus, aberrant DUX4 expression appeared to alter nucleocytoplasmic protein transport and generate double-strand DNA breaks in FSHD1 myotube nuclei, and the Bax/Bak pathway is required for DUX4-induced caspase activation but not γH2AX accumulation. These results add to our knowledge of downstream events induced by aberrant DUX4 expression and suggest possibilities for further mechanistic investigation.

Profiling of open chromatin in developing pig (Sus scrofa) muscle to identify regulatory regions

There is very little information about how the genome is regulated in domestic pigs (Sus scrofa). This lack of knowledge hinders efforts to define and predict the effects of genetic variants in pig breeding programs. To address this knowledge gap, we need to identify regulatory sequences in the pig genome starting with regions of open chromatin. We used the "Improved Protocol for the Assay for Transposase-Accessible Chromatin (Omni-ATAC-Seq)" to identify putative regulatory regions in flash-frozen semitendinosus muscle from 24 male piglets. We collected samples from the smallest-, average-, and largest-sized male piglets from each litter through five developmental time points. Of the 4661 ATAC-Seq peaks identified that represent regions of open chromatin, >50% were within 1 kb of known transcription start sites. Differential read count analysis revealed 377 ATAC-Seq defined genomic regions where chromatin accessibility differed significantly across developmental time points. We found regions of open chromatin associated with downregulation of genes involved in muscle development that were present in small-sized fetal piglets but absent in large-sized fetal piglets at day 90 of gestation. The dataset that we have generated provides a resource for studies of genome regulation in pigs and contributes valuable functional annotation information to filter genetic variants for use in genomic selection in pig breeding programs.