Dynamic transcriptional symmetry-breaking in pre-implantation mammalian embryo development revealed by single-cell RNA-seq

Living in a noisy world-origins of gene expression noise and its impact on cellular decision-making

The expression level of a gene can vary between genetically identical cells under the same environmental condition-a phenomenon referred to as gene expression noise. Several studies have now elucidated a central role of transcription factors in the generation of expression noise. Transcription factors, as the key components of gene regulatory networks, drive many important cellular decisions in response to cellular and environmental signals. Therefore, a very relevant question is how expression noise impacts gene regulation and influences cellular decision-making. In this Review, we summarize the current understanding of the molecular origins of expression noise, highlighting the role of transcription factors in this process, and discuss the ways in which noise can influence cellular decision-making. As advances in single-cell technologies open new avenues for studying expression noise as well as gene regulatory circuits, a better understanding of the influence of noise on cellular decisions will have important implications for many biological processes.

New insights into the epitranscriptomic control of pluripotent stem cell fate

Each cell in the human body has a distinguishable fate. Pluripotent stem cells are challenged with a myriad of lineage differentiation options. Defects are more likely to be fatal to stem cells than to somatic cells due to the broad impact of the former on early development. Hence, a detailed understanding of the mechanisms that determine the fate of stem cells is needed. The mechanisms by which human pluripotent stem cells, although not fully equipped with complex chromatin structures or epigenetic regulatory mechanisms, accurately control gene expression and are important to the stem cell field. In this review, we examine the events driving pluripotent stem cell fate and the underlying changes in gene expression during early development. In addition, we highlight the role played by the epitranscriptome in the regulation of gene expression that is necessary for each fate-related event.

NANOG initiates epiblast fate through the coordination of pluripotency genes expression

The epiblast is the source of all mammalian embryonic tissues and of pluripotent embryonic stem cells. It differentiates alongside the primitive endoderm in a “salt and pepper” pattern from inner cell mass (ICM) progenitors during the preimplantation stages through the activity of NANOG, GATA6 and the FGF pathway. When and how epiblast lineage specification is initiated is still unclear. Here, we show that the coordinated expression of pluripotency markers defines epiblast identity. Conversely, ICM progenitor cells display random cell-to-cell variability in expression of various pluripotency markers, remarkably dissimilar from the epiblast signature and independently from NANOG, GATA6 and FGF activities. Coordination of pluripotency markers expression fails in Nanog and Gata6 double KO (DKO) embryos. Collectively, our data suggest that NANOG triggers epiblast specification by ensuring the coordinated expression of pluripotency markers in a subset of cells, implying a stochastic mechanism. These features are likely conserved, as suggested by analysis of human embryos.

Pluripotent epiblast cells segregate from primitive endoderm in the blastocyst inner cell mass (ICM). Here the authors show that mosaic epiblast differentiation during mouse and human preimplantation development initiates stochastically in ICM progenitors, independently of the FGF pathway, and requires NANOG activity

Developmental capacity is unevenly distributed among single blastomeres of 2-cell and 4-cell stage mouse embryos

During preimplantation development, mammalian embryo cells (blastomeres) cleave, gradually losing their potencies and differentiating into three primary cell lineages: epiblast (EPI), trophectoderm (TE), and primitive endoderm (PE). The exact moment at which cells begin to vary in their potency for multilineage differentiation still remains unknown. We sought to answer the question of whether single cells isolated from 2- and 4-cell embryos differ in their ability to generate the progenitors and cells of blastocyst lineages. We revealed that twins were often able to develop into blastocysts containing inner cell masses (ICMs) with PE and EPI cells. Despite their capacity to create a blastocyst, the twins differed in their ability to produce EPI, PE, and TE cell lineages. In contrast, quadruplets rarely formed normal blastocysts, but instead developed into blastocysts with ICMs composed of only one cell lineage or completely devoid of an ICM altogether. We also showed that quadruplets have unequal capacities to differentiate into TE, PE, and EPI lineages. These findings could explain the difficulty of creating monozygotic twins and quadruplets from 2- and 4-cell stage mouse embryos.

Asymmetric expression of maternal mRNA governs first cell-fate decision

The goal of preimplantation development is to establish the fates of the embryonic and extra-embryonic cells. However, when and how cell fates are determined during early mammalian embryonic development remains unclear. We report that the high mobility group (HMG) protein family member HMGA1 was distributed differentially in mouse two-cell blastomeres. Knockdown of Hmga1 expression in one of the two cells reduced the number of cells contributing to the inner cell mass (ICM), suggesting that differential distribution of HMGA1 in the blastomeres in two-cell mouse embryos affected the selection of embryonic cell lineages. Mechanistically, HMGA1 promotes the expression of the ICM-specific gene Sox2. The results of this study show that mouse embryos demonstrate heterogeneity as early as the two-cell stage, and that these differences are related to cell-fate differentiation in early mouse embryos.

Totipotency of mouse zygotes extends to single blastomeres of embryos at the four-cell stage

In multicellular organisms, oocytes and sperm undergo fusion during fertilization and the resulting zygote gives rise to a new individual. The ability of zygotes to produce a fully formed individual from a single cell when placed in a supportive environment is known as totipotency. Given that totipotent cells are the source of all multicellular organisms, a better understanding of totipotency may have a wide-ranging impact on biology. The precise delineation of totipotent cells in mammals has remained elusive, however, although zygotes and single blastomeres of embryos at the two-cell stage have been thought to be the only totipotent cells in mice. We now show that a single blastomere of two- or four-cell mouse embryos can give rise to a fertile adult when placed in a uterus, even though blastomere isolation disturbs the transcriptome of derived embryos. Single blastomeres isolated from embryos at the eight-cell or morula stages and cultured in vitro manifested pronounced defects in the formation of epiblast and primitive endoderm by the inner cell mass and in the development of blastocysts, respectively. Our results thus indicate that totipotency of mouse zygotes extends to single blastomeres of embryos at the four-cell stage.

CARM1/PRMT4: Making Its Mark beyond Its Function as a Transcriptional Coactivator

Coactivator-associated arginine methyltransferase 1 (CARM1), identified 20 years ago as a coregulator of transcription, is an enzyme that catalyzes arginine methylation of proteins. Beyond its well-established involvement in the regulation of transcription, the physiological functions of CARM1 are still poorly understood. However, recent studies have revealed novel roles of CARM1 in autophagy, metabolism, paraspeckles, and early development. In addition, CARM1 is emerging as an attractive therapeutic target and a drug response biomarker for certain types of cancer. Here, we provide a comprehensive overview of the structure of CARM1 and its post-translational modifications, its various functions, apart from transcriptional coactivation, and its involvement in cancer.

Distinct molecular profile and restricted stem cell potential defines the prospective human cranial neural crest from embryonic stem cell state

Neural crest cells are an embryonic multipotent stem cell population. Recent studies in model organisms have suggested that neural crest cells are specified earlier than previously thought, at blastula stages. However, the molecular dynamics of early neural crest specification, and functional changes from pluripotent precursors to early specified NC, remain to be elucidated. In this report, we utilized a robust human model of cranial neural crest formation to address the distinct molecular character of the earliest stages of neural crest specification and assess the functional differences from its embryonic stem cell precursor. Our human neural crest model reveals a rapid change in the epigenetic state of neural crest and pluripotency genes, accompanied by changes in gene expression upon Wnt-based induction from embryonic stem cells. These changes in gene expression are directly regulated by the transcriptional activity of β-catenin. Furthermore, prospective cranial neural crest cells are characterized by restricted stem cell potential compared to embryonic stem cells. Our results suggest that human neural crest induced by Wnt/β-catenin signaling from human embryonic stem cells rapidly acquire a prospective neural crest cell state defined by a unique molecular signature and endowed with limited potential compared to pluripotent stem cells.

Measuring and Modeling Single-Cell Heterogeneity and Fate Decision in Mouse Embryos

Cellular heterogeneity is a property of any living system; however, its relationship with cellular fate decision remains an open question. Recent technological advances have enabled valuable insights, especially in complex systems such as the mouse embryo. In this review, we discuss recent studies that characterize cellular heterogeneity at different levels during mouse development, from the two-cell stage up to gastrulation. In addition to key experimental findings, we review mathematical modeling approaches that help researchers interpret these findings. Disentangling the role of heterogeneity in cell fate decision will likely rely on the refined integration of experiments, large-scale omics data, and mathematical modeling, complemented by the use of synthetic embryos and gastruloids as promising in vitro models.

A DNA methylation state transition model reveals the programmed epigenetic heterogeneity in human pre-implantation embryos

BACKGROUND

During mammalian early embryogenesis, expression and epigenetic heterogeneity emerge before the first cell fate determination, but the programs causing such determinate heterogeneity are largely unexplored.

RESULTS

Here, we present MethylTransition, a novel DNA methylation state transition model, for characterizing methylation changes during one or a few cell cycles at single-cell resolution. MethylTransition involves the creation of a transition matrix comprising three parameters that represent the probabilities of DNA methylation-modifying activities in order to link the methylation states before and after a cell cycle. We apply MethylTransition to single-cell DNA methylome data from human pre-implantation embryogenesis and elucidate that the DNA methylation heterogeneity that emerges at promoters during this process is largely an intrinsic output of a program with unique probabilities of DNA methylation-modifying activities. Moreover, we experimentally validate the effect of the initial DNA methylation on expression heterogeneity in pre-implantation mouse embryos.

CONCLUSIONS

Our study reveals the programmed DNA methylation heterogeneity during human pre-implantation embryogenesis through a novel mathematical model and provides valuable clues for identifying the driving factors of the first cell fate determination during this process.

Identification of the Lineage Markers and Inhibition of DAB2 in In Vitro Fertilized Porcine Embryos

Specification of embryonic lineages is an important question in the field of early development. Numerous studies analyzed the expression patterns of the candidate transcripts and proteins in humans and mice and clearly determined the markers of each lineage. To overcome the limitations of human and mouse embryos, the expression of the marker transcripts in each cell has been investigated using in vivo embryos in pigs. In vitro produced embryos are more accessible, can be rapidly processed with low cost. Therefore, we analyzed the characteristics of lineage markers and the effects of the DAB2 gene (trophectoderm marker) in in vitro fertilized porcine embryos. We investigated the expression levels of the marker genes during embryonic stages and distribution of the marker proteins was assayed in day 7 blastocysts. Then, the shRNA vectors were injected into the fertilized embryos and the differences in the marker transcripts were analyzed. Marker transcripts showed diverse patterns of expression, and each embryonic lineage could be identified with localization of marker proteins. In DAB2-shRNA vectors injected embryos, HNF4A and PDGFRA were upregulated. DAB2 protein level was lower in shRNA-injected embryos without significant differences. Our results will contribute to understanding of the mechanisms of embryonic lineage specification in pigs.

Keratins are asymmetrically inherited fate determinants in the mammalian embryo

To implant in the uterus, the mammalian embryo first specifies two cell lineages: the pluripotent inner cell mass that forms the fetus, and the outer trophectoderm layer that forms the placenta1. In many organisms, asymmetrically inherited fate determinants drive lineage specification2, but this is not thought to be the case during early mammalian development. Here we show that intermediate filaments assembled by keratins function as asymmetrically inherited fate determinants in the mammalian embryo. Unlike F-actin or microtubules, keratins are the first major components of the cytoskeleton that display prominent cell-to-cell variability, triggered by heterogeneities in the BAF chromatin-remodelling complex. Live-embryo imaging shows that keratins become asymmetrically inherited by outer daughter cells during cell division, where they stabilize the cortex to promote apical polarization and YAP-dependent expression of CDX2, thereby specifying the first trophectoderm cells of the embryo. Together, our data reveal a mechanism by which cell-to-cell heterogeneities that appear before the segregation of the trophectoderm and the inner cell mass influence lineage fate, via differential keratin regulation, and identify an early function for intermediate filaments in development.

Trophoblast attachment to the endometrial epithelium elicits compartment-specific transcriptional waves in an in-vitro model

RESEARCH QUESTION

Which are the early compartment-specific transcriptional responses of the trophoblast and the endometrial epithelium throughout early attachment during implantation?

DESIGN

An endometrial epithelium proxy (cell line Ishikawa) was co-cultured with spheroids of a green fluorescent protein (GFP) expressing trophoblast cell line (JEG-3). After 0, 8 and 24 h of co-culture, the compartments were sorted by fluorescence-activated cell sorting; GFP+ (trophoblast), GFP- (epithelium) and non-co-cultured control populations were analysed (in triplicate) by RNA-seq and gene set enrichment analysis (GSEA).

RESULTS

Trophoblast challenge induced a wave of transcriptional changes in the epithelium that resulted in 295 differentially regulated genes involving epithelial to mesenchymal transition (EMT), cell movement, apoptosis, hypoxia, inflammation, allograft rejection, myogenesis and cell signalling at 8 h. Interestingly, many of the enriched pathways were subsequently de-enriched by 24 h (i.e. EMT, cell movement, allograft rejection, myogenesis and cell signalling). In the trophoblast, the co-culture induced more transcriptional changes and regulation of a variety of pathways. A total of 1247 and 481 genes were differentially expressed after 8 h and from 8 to 24 h, respectively. Angiogenesis and hypoxia were over-represented at both stages, while EMT and cell signalling only were at 8 h; from 8 to 24 h, inflammation and oestrogen response were enriched, while proliferation was under-represented.

CONCLUSIONS

Successful attachment produced a series of dynamic changes in gene expression, characterized by an overall early and transient transcriptional up-regulation in the receptive epithelium, in contrast to a more dynamic transcriptional response in the trophoblast.

Common principles of early mammalian embryo self-organisation

Pre-implantation mammalian development unites extreme plasticity with a robust outcome: the formation of a blastocyst, an organised multi-layered structure ready for implantation. The process of blastocyst formation is one of the best-known examples of self-organisation. The first three cell lineages in mammalian development specify and arrange themselves during the morphogenic process based on cell-cell interactions. Despite decades of research, the unifying principles driving early mammalian development are still not fully defined. Here, we discuss the role of physical forces, and molecular and cellular mechanisms, in driving self-organisation and lineage formation that are shared between eutherian mammals.

Insights into epigenetic patterns in mammalian early embryos

Mammalian fertilization begins with the fusion of two specialized gametes, followed by major epigenetic remodeling leading to the formation of a totipotent embryo. During the development of the pre-implantation embryo, precise reprogramming progress is a prerequisite for avoiding developmental defects or embryonic lethality, but the underlying molecular mechanisms remain elusive. For the past few years, unprecedented breakthroughs have been made in mapping the regulatory network of dynamic epigenomes during mammalian early embryo development, taking advantage of multiple advances and innovations in low-input genome-wide chromatin analysis technologies. The aim of this review is to highlight the most recent progress in understanding the mechanisms of epigenetic remodeling during early embryogenesis in mammals, including DNA methylation, histone modifications, chromatin accessibility and 3D chromatin organization.

Specification of the First Mammalian Cell Lineages In Vivo and In Vitro

Our understanding of how the first mammalian cell lineages arise has been shaped largely by studies of the preimplantation mouse embryo. Painstaking work over many decades has begun to reveal how a single totipotent cell is transformed into a multilayered structure representing the foundations of the body plan. Here, we review how the first lineage decision is initiated by epigenetic regulation but consolidated by the integration of morphological features and transcription factor activity. The establishment of pluripotent and multipotent stem cell lines has enabled deeper analysis of molecular and epigenetic regulation of cell fate decisions. The capability to assemble these stem cells into artificial embryos is an exciting new avenue of research that offers a long-awaited window into cell fate specification in the human embryo. Together, these approaches are poised to profoundly increase our understanding of how the first lineage decisions are made during mammalian embryonic development.

Lessons from expanded potential of embryonic stem cells: Moving toward totipotency

Embryonic stem cells possess fascinating capacity of self-renewal and developmental potential, leading to significant progress in understanding the molecular basis of pluripotency, disease modeling, and reprogramming technology. Recently, 2-cell-like embryonic stem cells (ESCs) and expanded potential stem cells or extended pluripotent stem cells (EPSCs) generated from early-cleavage embryos display some features of totipotent embryos. These cell lines provide valuable in vitro models to study underlying principles of totipotency, cell plasticity, and lineage segregation. In this review, we summarize the current progress in this filed and highlight the application potentials of these cells in the future.

Blastula stage specification of avian neural crest

Cell fate specification defines the earliest steps towards a distinct cell lineage. Neural crest, a multipotent stem cell population, is thought to be specified from the ectoderm, but its varied contributions defy canons of segregation potential and challenges its embryonic origin. Aiming to resolve this conflict, we have assayed the earliest specification of neural crest using blastula stage chick embryos. Specification assays on isolated chick epiblast explants identify an intermediate region specified towards the neural crest cell fate. Furthermore, low density culture suggests that the specification of intermediate cells towards the neural crest lineage is independent of contact mediated induction and Wnt-ligand induced signaling, but is, however, dependent on transcriptional activity of β-catenin. Finally, we have validated the regional identity of the intermediate region towards the neural crest cell fate using fate map studies. Our results suggest a model of neural crest specification within a restricted epiblast region in blastula stage chick embryos.

Transitions in the cell-fate induction induced by colored noise associated with the inductive stimulus

The cell-fate induction based on the saddle-node bifurcation is undoubtedly a very important concept in developmental biology, which provides a possible mechanism to explain the intrinsic irreversibility in the developmental process. In this paper, the effect of a colored noise, which is associated with the inductive stimulus, on the saddle-node landscape of cell-fate induction is investigated, especially, the effect of the change of correlation time of colored noise on cell-fate induction. The main results show clearly that the change of correlation time of colored noise could induce the transitions of the system. This implies that the colored noise associated with inductive stimulus may have a profound effect on the saddle-node bifurcation landscape of cell-fate induction. This will also help us to understand more deeply the role of cell-fate induction in developmental biology.

Light-induced injury in mouse embryos revealed by single-cell RNA sequencing

Background

Light exposure is a common stress factor in in vitro manipulation of embryos in the reproductive center. Many studies have shown the deleterious effects of high-intensity light exposure in different animal embryos. However, no transcriptomic studies have explored the light-induced injury and response in preimplantation embryos.

Results

Here, we adopt different time-courses and illumination intensities to treat mouse embryos at the 2-cell stage and evaluate their effects on blastulation. Meanwhile, single-cell transcriptomes from the 2-cell to blastocyst stage were analyzed after high-intensity light exposure. These data show that cells at each embryonic stage can be categorized into different light conditions. Further analyses of differentially expressed genes and GO terms revealed the light-induced injury as well as the potential repair response after high-intensity lighting. Maternal-to-zygote transition is also affected by the failure to remove maternal RNAs and deactivate zygotic genome expression.

Conclusion

Our work revealed an integrated response to high-intensity lighting, involving morphological changes, long-lasting injury effects, and intracellular damage repair mechanisms.

CARM1 and Paraspeckles Regulate Pre-implantation Mouse Embryo Development

Nuclear architecture has never been carefully examined during early mammalian development at the stages leading to establishment of the embryonic and extra-embryonic lineages. Heterogeneous activity of the methyltransferase CARM1 during these stages results in differential methylation of histone H3R26 to modulate establishment of these two lineages. Here we show that CARM1 accumulates in nuclear granules at the 2- to 4-cell stage transition in the mouse embryo, with the majority corresponding to paraspeckles. The paraspeckle component Neat1 and its partner p54nrb are required for CARM1's association with paraspeckles and for H3R26 methylation. Conversely, CARM1 also influences paraspeckle organization. Depletion of Neat1 or p54nrb results in arrest at the 16- to 32-cell stage, with elevated expression of transcription factor Cdx2, promoting differentiation into the extra-embryonic lineage. This developmental arrest occurs at an earlier stage than following CARM1 depletion, indicating that paraspeckles act upstream of CARM1 but also have additional earlier roles in fate choice.

Sperm RNA code programmes the metabolic health of offspring

Mammalian sperm RNA is increasingly recognized as an additional source of paternal hereditary information beyond DNA. Environmental inputs, including an unhealthy diet, mental stresses and toxin exposure, can reshape the sperm RNA signature and induce offspring phenotypes that relate to paternal environmental stressors. Our understanding of the categories of sperm RNAs (such as tRNA-derived small RNAs, microRNAs, ribosomal RNA-derived small RNAs and long non-coding RNAs) and associated RNA modifications is expanding and has begun to reveal the functional diversity and information capacity of these molecules. However, the coding mechanism endowed by sperm RNA structures and by RNA interactions with DNA and other epigenetic factors remains unknown. How sperm RNA-encoded information is decoded in early embryos to control offspring phenotypes also remains unclear. Complete deciphering of the 'sperm RNA code' with regard to metabolic control could move the field towards translational applications and precision medicine, and this may lead to prevention of intergenerational transmission of obesity and type 2 diabetes mellitus susceptibility.

Retrospective analysis: reproducibility of interblastomere differences of mRNA expression in 2-cell stage mouse embryos is remarkably poor due to combinatorial mechanisms of blastomere diversification

STUDY QUESTION

What is the prevalence, reproducibility and biological significance of transcriptomic differences between sister blastomeres of the mouse 2-cell embryo?

SUMMARY ANSWER

Sister 2-cell stage blastomeres are distinguishable from each other by mRNA analysis, attesting to the fact that differentiation starts mostly early in the mouse embryo; however, the interblastomere differences are poorly reproducible and invoke the combinatorial effects of known and new mechanisms of blastomere diversification.

WHAT IS KNOWN ALREADY

Transcriptomic datasets for single blastomeres in mice have been available for years but have never been systematically analysed together, although such an analysis may shed light onto some unclarified topics of early mammalian development. Two unknowns that remain are at which stage embryonic blastomeres start to diversify from each other and what is the molecular origin of that difference. At the earliest postzygotic stage, the 2-cell stage, opinions differ regarding the answer to these questions; one group claims that the first zygotic division yields two equal blastomeres capable of forming a full organism (totipotency) and another group claims evidence for interblastomere differences reminiscent of the prepatterning found in embryos of lower taxa. Regarding the molecular origin of interblastomere differences, there are four prevalent models which invoke (1) oocyte anisotropy, (2) sperm entry point, (3) partition errors of the transcript pool and (4) asynchronous embryonic genome activation in the two blastomeres.

STUDY DESIGN, SIZE, DURATION

Seven transcriptomic studies published between 2011 and 2017 were eligible for retrospective analysis, since both blastomeres of the mouse 2-cell embryo had been analysed individually regarding the original pair associations and since the datasets were made available in public repositories. Five of these studies, encompassing a total of 43 pairs of sister blastomeres, were selected for further analyses based on high interblastomere correlations of mRNA levels. A double cut-off was used to select mRNAs that had robust interblastomere differences both within and between embryos (hits). The hits of each study were compared and contrasted with the hits of the other studies using Venn diagrams. The hits shared by at least four of five studies were analysed further by bioinformatics.

PARTICIPANTS/MATERIALS, SETTING, METHODS

PubMed was systematically examined for mRNA expression profiles of single 2-cell stage blastomeres in addition to publicly available microarray datasets (GEO, ArrayExpress). Based on the original normalizations, data from seven studies were screened for pairwise sample correlation at the gene level (Spearman), and the top five datasets with the highest correlation were subjected to hierarchical cluster analysis. Interblastomere differences of gene expression were expressed as a ratio of the higher to the lower mRNA level for each pair of blastomeres. A double cut-off was used to make the call of interblastomere difference, accepting genes with mRNA ratios above 2 when observed in at least 50% of the pairs, and discarding the other genes. The proportion of interblastomere differences common to at least four of the five datasets was calculated. Finally, the corresponding gene, pathway and enrichment analyses were performed utilizing PANTHER and GORILLA platforms.

MAIN RESULTS AND THE ROLE OF CHANCE

An average of 17% of genes within the datasets are differently expressed between sister blastomeres, a proportion which falls to 1% when considering the differences that are common to at least four of the five studies. Housekeeping mRNAs were not included in the 17% and 1% gene lists, suggesting that the interblastomere differences do not occur simply by chance. The 1% of shared interblastomere differences comprise 100 genes, of which 35 are consistent with at least one of the four prevalent models of sister blastomere diversification. Bioinformatics analysis of the remaining 65 genes that are not consistent with the four models suggests that at least one more mechanism is at play, potentially related to the endomembrane system. Although there are many dimensions to the issue of reproducibility (biological, experimental, analytical), we consider that the sister blastomeres are poised to escape high interblastomere correlations of mRNA levels, because at least five sources of diversity superimpose on each other, accounting for at least 25 = 32 different states. As a result, interblastomere mRNA differences of a given 2-cell embryo are necessarily difficult to reproduce in another 2-cell embryo.

LARGE SCALE DATA

Data were as provided by the original studies (GSE21688, GSE22182, GSE27396, GSE45719, GSE57249, E-MTAB-3321, GSE94050).

LIMITATIONS, REASONS FOR CAUTION

The original studies present similarities (e.g. fertilization in vivo after ovarian stimulation) as well as differences (e.g. mouse strains, method and timing of blastomere separation). We identified robust mRNA differences between the sister blastomeres, but these differences are underestimated because our double cut-off method works with thresholds and affords more protection against false positives than false negatives. Regarding the false negatives, transcriptome analysis may have captured only part of the interblastomere differences due to: (1) the 2-fold cut-off not being sensitive enough to detect the remaining part of the interblastomere differences, (2) the detection limit of the transcriptomic methods not being sufficient, or (3) interblastomere differences being oblivious to transcriptomic identification because transcriptional changes are oscillatory or because differences are mediated non-transcriptionally or post-transcriptionally. Regarding the false positives, it seems unlikely that a difference was found just by chance for the same group of transcripts due to the same technical error, given that different laboratories produced the data.

WIDER IMPLICATIONS OF THE FINDINGS

It is clear that the sister blastomeres are distinguishable from each other by mRNA analysis even at the 2-cell stage; however, efforts to identify large stable patterns may be in vain. This elicits thoughts about the wisdom of adding new transcriptomic datasets to the ones that already exist; if all transcriptomic datasets produced so far show a reproducibility of 1%, then any future study would probably face the same issue again. Possibly, a solid identification of the 'large stable pattern that should be there but was not found' requires an even larger dataset than the sum of the seven datasets considered here. Conversely, small stable patterns may be easier to identify, but their biological relevance is less obvious. Alternatively, interblastomere differences may not be mediated by nucleic acids but by other cellular components.

STUDY FUNDING/COMPETING INTEREST(S)

This study was supported by the Deutsche Forschungsgemeinschaft (grant DFG BO 2540-4-3 to M.B. and grant NO 413/3-3 to V.N.). The authors declare that they have no competing financial interests.

Asymmetric Inheritance of Cell Fate Determinants: Focus on RNA

During the last decade, and mainly primed by major developments in high-throughput sequencing technologies, the catalogue of RNA molecules harbouring regulatory functions has increased at a steady pace. Current evidence indicates that hundreds of mammalian RNAs have regulatory roles at several levels, including transcription, translation/post-translation, chromatin structure, and nuclear architecture, thus suggesting that RNA molecules are indeed mighty controllers in the flow of biological information. Therefore, it is logical to suggest that there must exist a series of molecular systems that safeguard the faithful inheritance of RNA content throughout cell division and that those mechanisms must be tightly controlled to ensure the successful segregation of key molecules to the progeny. Interestingly, whilst a handful of integral components of mammalian cells seem to follow a general pattern of asymmetric inheritance throughout division, the fate of RNA molecules largely remains a mystery. Herein, we will discuss current concepts of asymmetric inheritance in a wide range of systems, including prions, proteins, and finally RNA molecules, to assess overall the biological impact of RNA inheritance in cellular plasticity and evolutionary fitness.

Systems based analysis of human embryos and gene networks involved in cell lineage allocation

BACKGROUND

Little is understood of the molecular mechanisms involved in the earliest cell fate decision in human development, leading to the establishment of the trophectoderm (TE) and inner cell mass (ICM) stem cell population. Notably, there is a lack of understanding of how transcriptional networks arise during reorganisation of the embryonic genome post-fertilisation.

RESULTS

We identified a hierarchical structure of preimplantation gene network modules around the time of embryonic genome activation (EGA). Using network models along with eukaryotic initiation factor (EIF) and epigenetic-associated gene expression we defined two sets of blastomeres that exhibited diverging tendencies towards ICM or TE. Analysis of the developmental networks demonstrated stage specific EIF expression and revealed that histone modifications may be an important epigenetic regulatory mechanism in preimplantation human embryos. Comparison to published RNAseq data confirmed that during EGA the individual 8-cell blastomeres are transcriptionally primed for the first lineage decision in development towards ICM or TE.

CONCLUSIONS

Using multiple systems biology approaches to compare developmental stages in the early human embryo with single cell transcript data from blastomeres, we have shown that blastomeres considered to be totipotent are not transcriptionally equivalent. Furthermore we have linked the developmental interactome to individual blastomeres and to later cell lineage. This has clinical implications for understanding the impact of fertility treatments and developmental programming of long term health.

Faster, higher, stronger: timely and robust cell fate/identity commitment in stem cell lineages

Precise specification of cell fate or identity within stem cell lineages is critical for ensuring correct stem cell lineage progression and tissue homeostasis. Failure to specify cell fate or identity in a timely and robust manner can result in developmental abnormalities and diseases such as cancer. However, the molecular basis of timely cell fate/identity specification is only beginning to be understood. In this review, we discuss key regulatory strategies employed in cell fate specification and highlight recent results revealing how timely and robust cell fate/identity commitment is achieved through transcriptional control.

Transcriptome profiling of human oocytes experiencing recurrent total fertilization failure

There exist some patients who face recurrent total fertilization failure during assisted reproduction treatment, but the pathological mechanism underlying is elusive. Here, by using sc-RNA-seq method, the transcriptome profiles of ten abnormally fertilized zygotes were assessed, including five zygotes from one patient with recurrent Poly-PN zygotes, and five zygotes from a patient with pronuclear fusion failure. Four zygotes with three pronuclear (Tri-PN) were collected from four different patients as controls. After that, we identified 951 and 1697 significantly differentially expressed genes (SDEGs) in Poly-PN and PN arrest zygotes, respectively as compared with the control group. KEGG analyses indicated down regulated genes in the Poly-PN group included oocyte meiosis related genes, such as PPP2R1B, YWHAZ, MAD2L1, SPDYC, SKP1 and CDC27, together with genes associated with RNA processing, such as SF3B1, LOC645691, MAGOHB, PHF5A, PRPF18, DDX5, THOC1 and BAT1. In contrast, down regulated genes in the PN arrest group, included cell cycle genes, such as E2F4, DBF4, YWHAB, SKP2, CDC23, SMC3, CDC25A, CCND3, BUB1B, MDM2, CCNA2 and CDC7, together with homologous recombination related genes, such as NBN, XRCC3, SHFM1, RAD54B and RAD51. Thus, our work provides a better understanding of transcriptome profiles underlying RTFF, although it based on a limited number of patients.

Integrated Analysis of Quantitative Proteome and Transcriptional Profiles Reveals the Dynamic Function of Maternally Expressed Proteins After Parthenogenetic Activation of Buffalo Oocyte

Maternal-effect genes are especially critical for early embryonic development after fertilization and until massive activation of the embryonic genome occurs. By applying a tandem mass tag (TMT)-labeled quantitative proteomics combined with RNA sequencing approach, the proteome of the buffalo was quantitatively analyzed during parthenogenesis of mature oocytes and the two-cell stage embryo. Of 1908 quantified proteins, 123 differed significantly. The transcriptome was analyzed eight stages (GV, MII, 2-cell, 4-cell, 8-cell, 16-cell, morula, blastocyst) of Buffalo using the RNA sequencing approach, and a total of 3567 unique genes were identified to be differently expressed between all consecutive stages of pre-implantation development. Validation of proteomics results (TUBB3, CTNNA1, CDH3, MAP2K1), which are involved in tight junction and gap junction, revealing that the maternal expression of the proteins possibly plays a role in the formation of cellular junctions firstly after parthenogenetic activation. Correlation and hierarchical analyses of transcriptional profiles and the expression of NPM2 and NLRP5 mRNA of buffalo in vitro developed oocytes and parthenogenetic embryos indicated that the "maternal-to-zygotic transition" (MZT) process might exist in the model of parthenogenesis, which is similar to a normally fertilized embryo, and may occur between the 8-cell to 16-cell stage. These data provide a rich resource for further studies on maternal proteins and genes and are conducive to improving nuclear transfer technology.

Rainbow-Seq: Combining Cell Lineage Tracing with Single-Cell RNA Sequencing in Preimplantation Embryos

We developed the Rainbow-seq technology to trace cell division history and reveal single-cell transcriptomes. With distinct fluorescent protein genes as lineage markers, Rainbow-seq enables each single-cell RNA sequencing (RNA-seq) experiment to simultaneously decode the lineage marker genes and read single-cell transcriptomes. We triggered lineage tracking in each blastomere at the 2-cell stage, observed microscopically inequivalent contributions of the progeny to the two embryonic poles at the blastocyst stage, and analyzed every single cell at either 4- or 8-cell stage with deep paired-end sequencing of full-length transcripts. Although lineage difference was not marked unequivocally at a single-gene level, it became clear when the transcriptome was analyzed as a whole. Moreover, several groups of novel transcript isoforms with embedded repeat sequences exhibited lineage difference, suggesting a possible link between DNA demethylation and cell fate decision. Rainbow-seq bridged a critical gap between division history and single-cell RNA-seq assays.

Deconstructing and reconstructing the mouse and human early embryo

The emergence of form and function during mammalian embryogenesis is a complex process that involves multiple regulatory levels. The foundations of the body plan are laid throughout the first days of post-implantation development as embryonic stem cells undergo symmetry breaking and initiate lineage specification, in a process that coincides with a global morphological reorganization of the embryo. Here, we review experimental models and how they have shaped our current understanding of the post-implantation mammalian embryo.

Building Principles for Constructing a Mammalian Blastocyst Embryo

The self-organisation of a fertilised egg to form a blastocyst structure, which consists of three distinct cell lineages (trophoblast, epiblast and hypoblast) arranged around an off-centre cavity, is unique to mammals. While the starting point (the zygote) and endpoint (the blastocyst) are similar in all mammals, the intervening events have diverged. This review examines and compares the descriptive and functional data surrounding embryonic gene activation, symmetry-breaking, first and second lineage establishment, and fate commitment in a wide range of mammalian orders. The exquisite detail known from mouse embryogenesis, embryonic stem cell studies and the wealth of recent single cell transcriptomic experiments are used to highlight the building principles underlying early mammalian embryonic development.

Tracing the origin of heterogeneity and symmetry breaking in the early mammalian embryo

A fundamental question in developmental and stem cell biology concerns the origin and nature of signals that initiate asymmetry leading to pattern formation and self-organization. Instead of having prominent pre-patterning determinants as present in model organisms (worms, sea urchin, frog), we propose that the mammalian embryo takes advantage of more subtle cues such as compartmentalized intracellular reactions that generate micro-scale inhomogeneity, which is gradually amplified over several cellular generations to drive pattern formation while keeping developmental plasticity. It is therefore possible that by making use of compartmentalized information followed by its amplification, mammalian embryos would follow general principle of development found in other organisms in which the spatial cue is more robustly presented.

Making lineage decisions with biological noise: Lessons from the early mouse embryo

Understanding how individual cells make fate decisions that lead to the faithful formation and homeostatic maintenance of tissues is a fundamental goal of contemporary developmental and stem cell biology. Seemingly uniform populations of stem cells and multipotent progenitors display a surprising degree of heterogeneity, primarily originating from the inherent stochastic nature of molecular processes underlying gene expression. Despite this heterogeneity, lineage decisions result in tissues of a defined size and with consistent proportions of differentiated cell types. Using the early mouse embryo as a model we review recent developments that have allowed the quantification of molecular intercellular heterogeneity during cell differentiation. We first discuss the relationship between these heterogeneities and developmental cellular potential. We then review recent theoretical approaches that formalize the mechanisms underlying fate decisions in the inner cell mass of the blastocyst stage embryo. These models build on our extensive knowledge of the genetic control of fate decisions in this system and will become essential tools for a rigorous understanding of the connection between noisy molecular processes and reproducible outcomes at the multicellular level. We conclude by suggesting that cell-to-cell communication provides a mechanism to exploit and buffer intercellular variability in a self-organized process that culminates in the reproducible formation of the mature mammalian blastocyst stage embryo that is ready for implantation into the maternal uterus. This article is categorized under: Gene Expression and Transcriptional Hierarchies > Cellular Differentiation Establishment of Spatial and Temporal Patterns > Regulation of Size, Proportion, and Timing Gene Expression and Transcriptional Hierarchies > Gene Networks and Genomics Gene Expression and Transcriptional Hierarchies > Quantitative Methods and Models.

Long non-coding RNAs involved in the regulatory network during porcine pre-implantation embryonic development and iPSC induction

Long non-coding RNAs (lncRNA) play a key role in the orchestration of transcriptional regulation during development and many other cellular processes. The importance of the regulatory co-expression network was highlighted in the identification of the mechanism of these processes in humans and mice. However, elucidation of the properties of porcine lncRNAs involved in the regulatory network during pre-implantation embryonic development and fibroblast reprogramming to induced pluripotent stem cell (iPSC) has been limited to date. Using a weighted gene co-expression network analysis, we constructed the regulatory network and determined that the novel lncRNAs were functionally involved in key events of embryonic development during the pre-implantation period; moreover, reprogramming could be delineated by a small number of potentially functional modules of co-expressed genes. These findings indicate that lncRNAs may be involved in the transcriptional regulation of zygotic genome activation, first lineage segregation and somatic reprogramming to pluripotency. Furthermore, we performed a conservation and synteny analysis with the significant lncRNAs involved in these vital events and validated the results via experimental assays. In summary, the current findings provide a valuable resource to dissect the protein coding gene and lncRNA regulatory networks that underlie the progressive development of embryos and somatic reprogramming.

Using single-cell genomics to understand developmental processes and cell fate decisions

High-throughput -omics techniques have revolutionised biology, allowing for thorough and unbiased characterisation of the molecular states of biological systems. However, cellular decision-making is inherently a unicellular process to which "bulk" -omics techniques are poorly suited, as they capture ensemble averages of cell states. Recently developed single-cell methods bridge this gap, allowing high-throughput molecular surveys of individual cells. In this review, we cover core concepts of analysis of single-cell gene expression data and highlight areas of developmental biology where single-cell techniques have made important contributions. These include understanding of cell-to-cell heterogeneity, the tracing of differentiation pathways, quantification of gene expression from specific alleles, and the future directions of cell lineage tracing and spatial gene expression analysis.

Transcriptional Regulation and Genes Involved in First Lineage Specification During Preimplantation Development

The successful development from a single-cell zygote into a complex multicellular organism requires precise coordination of multiple cell-fate decisions. The very first of these is lineage specification into the inner cell mass (ICM) and trophectoderm (TE) during mammalian preimplantation development. In mouse embryos, transcription factors (TFs) such as Oct4, Sox2, and Nanog are enriched in cells of ICM, which gives rise to the fetus and yolk sac. Conversely, TFs such as Cdx2 and Eomes become highly upregulated in TE, which contribute to the placenta. Here, we review the current understanding of key transcriptional control mechanisms and genes responsible for these distinct differences during the first cell lineage specification. In particular, we highlight recent insights gained through advances in genome manipulation, live imaging, single-cell transcriptomics, and loss-of-function studies.

Cell Fate Decisions During Preimplantation Mammalian Development

The early mouse embryo offers a phenomenal system to dissect how changes in the mechanisms controlling cell fate are integrated with morphogenetic events at the single-cell level. New technologies based on live imaging have enabled the discovery of dynamic changes in the regulation of single genes, transcription factors, and epigenetic mechanisms directing early cell fate decision in the early embryo. Here, we review recent progress in linking molecular dynamic events occurring at the level of the single cell in vivo, to some of the key morphogenetic changes regulating early mouse development.

Understanding development and stem cells using single cell-based analyses of gene expression

In recent years, genome-wide profiling approaches have begun to uncover the molecular programs that drive developmental processes. In particular, technical advances that enable genome-wide profiling of thousands of individual cells have provided the tantalizing prospect of cataloging cell type diversity and developmental dynamics in a quantitative and comprehensive manner. Here, we review how single-cell RNA sequencing has provided key insights into mammalian developmental and stem cell biology, emphasizing the analytical approaches that are specific to studying gene expression in single cells.

Totipotency segregates between the sister blastomeres of two-cell stage mouse embryos

Following fertilization in mammals, it is generally accepted that totipotent cells are exclusive to the zygote and to each of the two blastomeres originating from the first mitotic division. This model of totipotency was inferred from a minority of cases in which blastomeres produced monozygotic twins in mice. Was this due to experimental limitation or biological constraint? Here we removed experimental obstacles and achieved reliable quantification of the prevalence of dual totipotency among mouse two-cell stage blastomeres. We separated the blastomeres of 1,252 two-cell embryos, preserving 1,210 of the pairs. Two classes of monozygotic twins became apparent at the blastocyst stage: 27% formed a functional epiblast in both members (concordant), and 73% did so in only one member of the pair (discordant) - a partition that proved insensitive to oocyte quality, sperm-entry point, culture environment and pattern of cleavage. In intact two-cell embryos, the ability of sister blastomeres to generate epiblast was also skewed. Class discovery clustering of the individual blastomeres' and blastocysts' transcriptomes points to an innate origin of concordance and discordance rather than developmental acquisition. Our data place constraints on the commonly accepted idea that totipotency is allocated equally between the two-cell stage blastomeres in mice.

Single-Cell Landscape of Transcriptional Heterogeneity and Cell Fate Decisions during Mouse Early Gastrulation

The mouse inner cell mass (ICM) segregates into the epiblast and primitive endoderm (PrE) lineages coincident with implantation of the embryo. The epiblast subsequently undergoes considerable expansion of cell numbers prior to gastrulation. To investigate underlying regulatory principles, we performed systematic single-cell RNA sequencing (seq) of conceptuses from E3.5 to E6.5. The epiblast shows reactivation and subsequent inactivation of the X chromosome, with Zfp57 expression associated with reactivation and inactivation together with other candidate regulators. At E6.5, the transition from epiblast to primitive streak is linked with decreased expression of polycomb subunits, suggesting a key regulatory role. Notably, our analyses suggest elevated transcriptional noise at E3.5 and within the non-committed epiblast at E6.5, coinciding with exit from pluripotency. By contrast, E6.5 primitive streak cells became highly synchronized and exhibit a shortened G1 cell-cycle phase, consistent with accelerated proliferation. Our study systematically charts transcriptional noise and uncovers molecular processes associated with early lineage decisions.

New insights into human primordial germ cells and early embryonic development from single-cell analysis

Human preimplantation developmental studies are difficult to accomplish due to associated ethical and moral issues. Preimplantation cells are rare and exist only in transient cell states. From a single cell, it is very challenging to analyse the origination of the heterogeneity and complexity inherent to the human body. However, recent advances in single-cell technology and data analysis have provided new insights into the process of early human development and germ cell specification. In this Review, we examine the latest single-cell datasets of human preimplantation embryos and germ cell development, compare them to bulk cell analyses, and interpret their biological implications.

Embryoids, organoids and gastruloids: new approaches to understanding embryogenesis

Cells have an intrinsic ability to self-assemble and self-organize into complex and functional tissues and organs. By taking advantage of this ability, embryoids, organoids and gastruloids have recently been generated in vitro, providing a unique opportunity to explore complex embryological events in a detailed and highly quantitative manner. Here, we examine how such approaches are being used to answer fundamental questions in embryology, such as how cells self-organize and assemble, how the embryo breaks symmetry, and what controls timing and size in development. We also highlight how further improvements to these exciting technologies, based on the development of quantitative platforms to precisely follow and measure subcellular and molecular events, are paving the way for a more complete understanding of the complex events that help build the human embryo.

Exploratory bioinformatics investigation reveals importance of "junk" DNA in early embryo development

BACKGROUND

Instead of testing predefined hypotheses, the goal of exploratory data analysis (EDA) is to find what data can tell us. Following this strategy, we re-analyzed a large body of genomic data to study the complex gene regulation in mouse pre-implantation development (PD).

RESULTS

Starting with a single-cell RNA-seq dataset consisting of 259 mouse embryonic cells derived from zygote to blastocyst stages, we reconstructed the temporal and spatial gene expression pattern during PD. The dynamics of gene expression can be partially explained by the enrichment of transposable elements in gene promoters and the similarity of expression profiles with those of corresponding transposons. Long Terminal Repeats (LTRs) are associated with transient, strong induction of many nearby genes at the 2-4 cell stages, probably by providing binding sites for Obox and other homeobox factors. B1 and B2 SINEs (Short Interspersed Nuclear Elements) are correlated with the upregulation of thousands of nearby genes during zygotic genome activation. Such enhancer-like effects are also found for human Alu and bovine tRNA SINEs. SINEs also seem to be predictive of gene expression in embryonic stem cells (ESCs), raising the possibility that they may also be involved in regulating pluripotency. We also identified many potential transcription factors underlying PD and discussed the evolutionary necessity of transposons in enhancing genetic diversity, especially for species with longer generation time.

CONCLUSIONS

Together with other recent studies, our results provide further evidence that many transposable elements may play a role in establishing the expression landscape in early embryos. It also demonstrates that exploratory bioinformatics investigation can pinpoint developmental pathways for further study, and serve as a strategy to generate novel insights from big genomic data.