NANOG initiates epiblast fate through the coordination of pluripotency genes expression

A single-cell transcriptome atlas of human euploid and aneuploid blastocysts

Aneuploidy is frequently detected in early human embryos as a major cause of early pregnancy failure. However, how aneuploidy affects cellular function remains elusive. Here, we profiled the transcriptomes of 14,908 single cells from 203 human euploid and aneuploid blastocysts involving autosomal and sex chromosomes. Nearly all of the blastocysts contained four lineages. In aneuploid chromosomes, 19.5% ± 1.2% of the expressed genes showed a dosage effect, and 90 dosage-sensitive domains were identified. Aneuploidy leads to prevalent genome-wide transcriptome alterations. Common effects, including apoptosis, were identified, especially in monosomies, partially explaining the lower cell numbers in autosomal monosomies. We further identified lineage-specific effects causing unstable epiblast development in aneuploidies, which was accompanied by the downregulation of TGF-β and FGF signaling, which resulted in insufficient trophectoderm maturation. Our work provides crucial insights into the molecular basis of human aneuploid blastocysts and may shed light on the cellular interaction during blastocyst development.

Longitudinal profiling of human androgenotes through single-cell analysis unveils paternal gene expression dynamics in early embryo development

STUDY QUESTION

How do transcriptomics vary in haploid human androgenote embryos at single cell level in the first four cell cycles of embryo development?

SUMMARY ANSWER

Gene expression peaks at the fourth cell cycle, however some androcytes exhibit unique transcriptional behaviors.

WHAT IS KNOWN ALREADY

The developmental potential of an embryo is determined by the competence of the oocyte and the sperm. However, studies of the contribution of the paternal genome using pure haploid androgenotes are very scarce.

STUDY DESIGN, SIZE, DURATION

This study was performed analyzing the single-cell transcriptomic sequencing of 38 androcytes obtained from 10 androgenote bioconstructs previously produced in vitro (de Castro et al., 2023). These results were analyzed through different bioinformatics software such as g: Profiler, GSEA, Cytoscape, and Reactome.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Single cell sequencing was used to obtain the transcriptomic profiles of the different androcytes. The results obtained were compared between the different cycles studied using the DESeq2 program and functional enrichment pathways using g: Profiler, Cytoscape, and Reactome.

MAIN RESULTS AND THE ROLE OF CHANCE

A wave of paternally driven transcriptomic activation was found during the third-cell cycle, with 1128 upregulated and 225 downregulated genes and the fourth-cell cycle, with 1373 upregulated and 286 downregulated genes, compared to first-cell cycle androcytes. Differentially expressed routes related to cell differentiation, DNA-binding transcription, RNA biosynthesis and RNA polymerase II transcription regulatory complex, and cell death were found in the third and fourth with respect to the first-cell cycle. Conversely, in the fourth cell cycle, 153 downregulated and 332 upregulated genes were found compared with third cell cycle, associated with differentially expressed processes related to E-box binding and zinc finger protein 652 (ZNF652) transcription factor. Further, significant overexpression of LEUTX, PRAMEF1, DUXA, RFPL4A, TRIM43, and ZNF675 found in androgenotes, compared to biparental embryos, highlights the paternal contributions to zygote genome activation.

LARGE SCALE DATA

All raw sequencing data are available through the Gene Expression Omnibus (GEO) under accessions number: GSE216501.

LIMITATIONS, REASONS FOR CAUTION

Extrapolation of biological events from uniparental constructs to biparental embryos should be done with caution. Maternal and paternal genomes do not act independently of each other in a natural condition. The absence of one genome may affect gene transcription of the other. In this sense, the haploid condition of the bioconstructs could mask the transcriptomic patterns of the single cells.

WIDER IMPLICATIONS OF THE FINDINGS

The results obtained demonstrated the level of involvement of the human paternal haploid genome in the early stages of embryo development as well as its evolution at the transcriptomic level, laying the groundwork for the use of these bioconstructs as reliable models to dispel doubts about the genetic role played by the paternal genome in the early cycles of embryo development.

STUDY FUNDING/COMPETING INTEREST(S)

This study was funded by Instituto de Salud Carlos III (ISCIII) through the project 'PI22/00924', co-funded by European Regional Development Fund (ERDF); 'A way to make Europe'. F.D. was supported by the Spanish Ministry of Economy and Competitiveness through the Miguel Servet program (CPII018/00002). M.J.E. was supported by Instituto de Salud Carlos III (PI19/00577 [M.J.E.]) and FI20/00086. P.dC. was supported by a predoctoral grant for training in research into health (PFIS PI19/00577) from the Instituto de Salud Carlos III. All authors declare having no conflict of interest with regard to this trial.

Branching topology of the human embryo transcriptome revealed by Entropy Sort Feature Weighting

Analysis of single cell transcriptomics (scRNA-seq) data is typically performed after subsetting to highly variable genes (HVGs). Here, we show that Entropy Sorting provides an alternative mathematical framework for feature selection. On synthetic datasets, continuous Entropy Sort Feature Weighting (cESFW) outperforms HVG selection in distinguishing cell-state-specific genes. We apply cESFW to six merged scRNA-seq datasets spanning human early embryo development. Without smoothing or augmenting the raw counts matrices, cESFW generates a high-resolution embedding displaying coherent developmental progression from eight-cell to post-implantation stages and delineating 15 distinct cell states. The embedding highlights sequential lineage decisions during blastocyst development, while unsupervised clustering identifies branch point populations obscured in previous analyses. The first branching region, where morula cells become specified for inner cell mass or trophectoderm, includes cells previously asserted to lack a developmental trajectory. We quantify the relatedness of different pluripotent stem cell cultures to distinct embryo cell types and identify marker genes of naïve and primed pluripotency. Finally, by revealing genes with dynamic lineage-specific expression, we provide markers for staging progression from morula to blastocyst.

Unlocking trophectoderm mysteries: In vivo and in vitro perspectives on human and mouse trophoblast fate induction

In recent years, the pursuit of inducing the trophoblast stem cell (TSC) state has gained prominence as a compelling research objective, illuminating the establishment of the trophoblast lineage and unlocking insights into early embryogenesis. In this review, we examine how advancements in diverse technologies, including in vivo time course transcriptomics, cellular reprogramming to TSC state, chemical induction of totipotent stem-cell-like state, and stem-cell-based embryo-like structures, have enriched our insights into the intricate molecular mechanisms and signaling pathways that define the mouse and human trophectoderm/TSC states. We delve into disparities between mouse and human trophectoderm/TSC fate establishment, with a special emphasis on the intriguing role of pluripotency in this context. Additionally, we re-evaluate recent findings concerning the potential of totipotent-stem-like cells and embryo-like structures to fully manifest the trophectoderm/trophoblast lineage's capabilities. Lastly, we briefly discuss the potential applications of induced TSCs in pregnancy-related disease modeling.

The Dynamics of Histone Modifications during Mammalian Zygotic Genome Activation

Mammalian fertilization initiates the reprogramming of oocytes and sperm, forming a totipotent zygote. During this intricate process, the zygotic genome undergoes a maternal-to-zygotic transition (MZT) and subsequent zygotic genome activation (ZGA), marking the initiation of transcriptional control and gene expression post-fertilization. Histone modifications are pivotal in shaping cellular identity and gene expression in many mammals. Recent advances in chromatin analysis have enabled detailed explorations of histone modifications during ZGA. This review delves into conserved and unique regulatory strategies, providing essential insights into the dynamic changes in histone modifications and their variants during ZGA in mammals. The objective is to explore recent advancements in leading mechanisms related to histone modifications governing this embryonic development phase in depth. These considerations will be useful for informing future therapeutic approaches that target epigenetic regulation in diverse biological contexts. It will also contribute to the extensive areas of evolutionary and developmental biology and possibly lay the foundation for future research and discussion on this seminal topic.

NANOG expression in parthenogenetic porcine blastocysts is required for intact lineage specification and pluripotency

OBJECTIVE

Nanog homeobox (NANOG) is a core transcription factor that contributes to pluripotency along with octamer binding transcription factor-4 (OCT4) and sex determining region-Y box-2 (SOX2). It is an epiblast lineage marker in mammalian pre-implantation embryos and exhibits a species-specific expression pattern. Therefore, it is important to understand the lineage of NANOG, the trophectoderm, and the primitive endoderm in the pig embryo.

METHODS

A loss- and gain-of-function analysis was done to determine the role of NANOG in lineage specification in parthenogenetic porcine blastocysts. We analyzed the relationship between NANOG and pluripotent core transcription factors and other lineage makers.

RESULTS

In NANOG-null late blastocysts, OCT4-, SOX2-, and SOX17-positive cells were decreased, whereas GATA binding protein 6 (GATA6)-positive cells were increased. Quantitative real-time polymerase chain reaction revealed that the expression of SOX2 was decreased in NANOG-null blastocysts, whereas that of primitive endoderm makers, except SOX17, was increased. In NANOG-overexpressing blastocysts, caudal type homeobox 2 (CDX2-), SOX17-, and GATA6-positive cells were decreased. The results indicated that the expression of primitive endoderm markers and trophectoderm-related genes was decreased.

CONCLUSION

Taken together, the results demonstrate that NANOG is involved in the epiblast and primitive endoderm differentiation and is essential for maintaining pluripotency within the epiblast.

Integrative analysis of single-cell embryo data reveals transcriptome signatures for the human pre-implantation inner cell mass

As the source of embryonic stem cells (ESCs), inner cell mass (ICM) can form all tissues of the embryo proper, however, its role in early human lineage specification remains controversial. Although a stepwise differentiation model has been proposed suggesting the existence of ICM as a distinct developmental stage, the underlying molecular mechanism remains unclear. In the present study, we perform an integrated analysis on the public human preimplantation embryonic single-cell transcriptomic data and apply a trajectory inference algorithm to measure the cell plasticity. In our results, ICM population can be clearly discriminated on the dimension-reduced graph and confirmed by compelling evidences, thus validating the two-step hypothesis of lineage commitment. According to the branch probabilities and differentiation potential, we determine the precise time points for two lineage segregations. Further analysis on gene expression dynamics and regulatory network indicates that transcription factors including GSC, PRDM1, and SPIC may underlie the decisions of ICM fate. In addition, new human ICM marker genes, such as EPHA4 and CCR8 are discovered and validated by immunofluorescence. Given the potential clinical applications of ESCs, our analysis provides a further understanding of human ICM cells and facilitates the exploration of more unique characteristics in early human development.

Spatial positioning of preimplantation mouse embryo cells is regulated by mTORC1 and m7G-cap-dependent translation at the 8- to 16-cell transition

Preimplantation mouse embryo development involves temporal-spatial specification and segregation of three blastocyst cell lineages: trophectoderm, primitive endoderm and epiblast. Spatial separation of the outer-trophectoderm lineage from the two other inner-cell-mass (ICM) lineages starts with the 8- to 16-cell transition and concludes at the 32-cell stages. Accordingly, the ICM is derived from primary and secondary contributed cells; with debated relative EPI versus PrE potencies. We report generation of primary but not secondary ICM populations is highly dependent on temporal activation of mammalian target of Rapamycin (mTOR) during 8-cell stage M-phase entry, mediated via regulation of the 7-methylguanosine-cap (m7G-cap)-binding initiation complex (EIF4F) and linked to translation of mRNAs containing 5' UTR terminal oligopyrimidine (TOP-) sequence motifs, as knockdown of identified TOP-like motif transcripts impairs generation of primary ICM founders. However, mTOR inhibition-induced ICM cell number deficits in early blastocysts can be compensated by the late blastocyst stage, after inhibitor withdrawal; compensation likely initiated at the 32-cell stage when supernumerary outer cells exhibit molecular characteristics of inner cells. These data identify a novel mechanism specifically governing initial spatial segregation of mouse embryo blastomeres, that is distinct from those directing subsequent inner cell formation, contributing to germane segregation of late blastocyst lineages.

The Human Early Maternal–Embryonic Interactome

Background: Single cell transcriptomics offers an avenue for predicting, with improved accuracy, the gene networks that are involved in the establishment of the first direct cell–cell interactions between the blastocyst and the maternal luminal epithelium. We hypothesised that in silico modelling of the maternal–embryonic interface may provide a causal model of these interactions, leading to the identification of genes associated with a successful initiation of implantation. Methods: Bulk and single cell RNA-sequencing of endometrial epithelium and scRNAseq of day 6 and 7 trophectoderm (TE) were used to model the initial encounter between the blastocyst and the maternal uterine lining epithelium in silico. In silico modelling of the maternal–embryonic interface was performed using hypernetwork (HN) analysis of genes mediating endometrial–TE interactions and the wider endometrial epithelial transcriptome. A hypernetwork analysis identifies genes that co-ordinate the expression of many other genes to derive a higher order interaction likely to be causally linked to the function. Potential interactions of TE with non-ciliated luminal cells, ciliated cells, and glandular cells were examined. Results: Prominent epithelial activities include secretion, endocytosis, ion transport, adhesion, and immune modulation. Three highly correlated clusters of 25, 22 and 26 TE-interacting epithelial surface genes were identified, each with distinct properties. Genes in both ciliated and non-ciliated luminal epithelial cells and glandular cells exhibit significant functional associations. Ciliated cells are predicted to bind to TE via galectin–glycan interaction. Day 6 and day 7 embryonic–epithelial interactomes are largely similar. The removal of aneuploid TE-derived mRNA invoked only subtle differences. No direct interaction with the maternal gland epithelial cell surface is predicted. These functional differences validate the in silico segregation of phenotypes. Single cell analysis of the epithelium revealed significant change with the cycle phase, but differences in the cell phenotype between individual donors were also present. Conclusions: A hypernetwork analysis can identify epithelial gene clusters that show correlated change during the menstrual cycle and can be interfaced with TE genes to predict pathways and processes occurring during the initiation of embryo–epithelial interaction in the mid-secretory phase. The data are on a scale that is realistic for functional dissection using current ex vivo human implantation models. A focus on luminal epithelial cells may allow a resolution to the current bottleneck of endometrial receptivity testing based on tissue lysates, which is confounded by noise from multiple diverse cell populations.

The role of Wnt signaling in the development of the epiblast and axial progenitors

Understanding how the body plan is established during embryogenesis remains a fundamental biological question. The Wnt/β-catenin signaling pathway plays a crucial and highly conserved role in body plan formation, functioning to polarize the primary anterior-posterior (AP) or head-to-tail body axis in most metazoans. In this chapter, we focus on the roles that the mammalian Wnt/β-catenin pathway plays to prepare the pluripotent epiblast for gastrulation, and to elicit the emergence of multipotent axial progenitors from the caudal epiblast. Interactions between Wnt and retinoic acid (RA), another powerful family of developmental signaling molecules, in axial progenitors will also be discussed. Gastrulation movements and somitogenesis result in the anterior displacement of the RA source (the rostral somites and lateral plate mesoderm (LPM)), from the posterior Wnt source (the primitive streak (PS)), leading to the establishment of antiparallel gradients of RA and Wnt that control the self-renewal and successive differentiation of neck, trunk and tail progenitors.

Initial source of heterogeneity in a model for cell fate decision in the early mammalian embryo

During development, cells from a population of common progenitors evolve towards different fates characterized by distinct levels of specific transcription factors, a process known as cell differentiation. This evolution is governed by gene regulatory networks modulated by intercellular signalling. In order to evolve towards distinct fates, cells forming the population of common progenitors must display some heterogeneity. We applied a modelling approach to obtain insights into the possible sources of cell-to-cell variability initiating the specification of cells of the inner cell mass into epiblast or primitive endoderm cells in early mammalian embryo. At the single-cell level, these cell fates correspond to three possible steady states of the model. A combination of numerical simulations and bifurcation analyses predicts that the behaviour of the model is preserved with respect to the source of variability and that cell-cell coupling induces the emergence of multiple steady states associated with various cell fate configurations, and to a distribution of the levels of expression of key transcription factors. Statistical analysis of these time-dependent distributions reveals differences in the evolutions of the variance-to-mean ratios of key variables of the system, depending on the simulated source of variability, and, by comparison with experimental data, points to the rate of synthesis of the key transcription factor NANOG as a likely initial source of heterogeneity.

Capturing Transitional Pluripotency through Proline Metabolism

In this paper, we summarize the current knowledge of the role of proline metabolism in the control of the identity of Embryonic Stem Cells (ESCs). An imbalance in proline metabolism shifts mouse ESCs toward a stable naïve-to-primed intermediate state of pluripotency. Proline-induced cells (PiCs), also named primitive ectoderm-like cells (EPLs), are phenotypically metastable, a trait linked to a rapid and reversible relocalization of E-cadherin from the plasma membrane to intracellular membrane compartments. The ESC-to-PiC transition relies on the activation of Erk and Tgfβ/Activin signaling pathways and is associated with extensive remodeling of the transcriptome, metabolome and epigenome. PiCs maintain several properties of naïve pluripotency (teratoma formation, blastocyst colonization and 3D gastruloid development) and acquire a few traits of primed cells (flat-shaped colony morphology, aerobic glycolysis metabolism and competence for primordial germ cell fate). Overall, the molecular and phenotypic features of PiCs resemble those of an early-primed state of pluripotency, providing a robust model to study the role of metabolic perturbations in pluripotency and cell fate decisions.