The transcriptomic architecture of mouse Sertoli cell clone embryos reveals temporal–spatial-specific reprogramming

Fetal growth delay caused by loss of non-canonical imprinting is resolved late in pregnancy and culminates in offspring overgrowth

Germline epigenetic programming, including genomic imprinting, substantially influences offspring development. Polycomb Repressive Complex 2 (PRC2) plays an important role in Histone 3 Lysine 27 trimethylation (H3K27me3)-dependent imprinting, loss of which leads to growth and developmental changes in mouse offspring. In this study, we show that offspring from mouse oocytes lacking the PRC2 protein Embryonic Ectoderm Development (EED) were initially developmentally delayed, characterised by low blastocyst cell counts and substantial growth delay in mid-gestation embryos. This initial developmental delay was resolved as offspring underwent accelerated fetal development and growth in late gestation resulting in offspring that were similar stage and weight to controls at birth. The accelerated development and growth in offspring from Eed-null oocytes was associated with remodelling of the placenta, which involved an increase in fetal and maternal tissue size, conspicuous expansion of the glycogen-enriched cell population, and delayed parturition. Despite placental remodelling and accelerated offspring fetal growth and development, placental efficiency, and fetal blood glucose levels were low, and the fetal blood metabolome was unchanged. Moreover, while expression of the H3K27me3-imprinted gene and amino acid transporter Slc38a4 was increased, fetal blood levels of individual amino acids were similar to controls, indicating that placental amino acid transport was not enhanced. Genome-wide analyses identified extensive transcriptional dysregulation and DNA methylation changes in affected placentas, including a range of imprinted and non-imprinted genes. Together, while deletion of Eed in growing oocytes resulted in fetal growth and developmental delay and placental hyperplasia, our data indicate a remarkable capacity for offspring fetal growth to be normalised despite inefficient placental function and the loss of H3K27me3-dependent genomic imprinting.

Polar Gini Curve: A Technique to Discover Gene Expression Spatial Patterns from Single-cell RNA-seq Data

In this work, we describe the development of Polar Gini Curve, a method for characterizing cluster markers by analyzing single-cell RNA sequencing (scRNA-seq) data. Polar Gini Curve combines the gene expression and the 2D coordinates ("spatial") information to detect patterns of uniformity in any clustered cells from scRNA-seq data. We demonstrate that Polar Gini Curve can help users characterize the shape and density distribution of cells in a particular cluster, which can be generated during routine scRNA-seq data analysis. To quantify the extent to which a gene is uniformly distributed in a cell cluster space, we combine two polar Gini curves (PGCs)-one drawn upon the cell-points expressing the gene (the "foreground curve") and the other drawn upon all cell-points in the cluster (the "background curve"). We show that genes with highly dissimilar foreground and background curves tend not to uniformly distributed in the cell cluster-thus having spatially divergent gene expression patterns within the cluster. Genes with similar foreground and background curves tend to uniformly distributed in the cell cluster-thus having uniform gene expression patterns within the cluster. Such quantitative attributes of PGCs can be applied to sensitively discover biomarkers across clusters from scRNA-seq data. We demonstrate the performance of the Polar Gini Curve framework in several simulation case studies. Using this framework to analyze a real-world neonatal mouse heart cell dataset, the detected biomarkers may characterize novel subtypes of cardiac muscle cells. The source code and data for Polar Gini Curve could be found at http://discovery.informatics.uab.edu/PGC/ or https://figshare.com/projects/Polar_Gini_Curve/76749.

Loss of H3K27me3 imprinting in the Sfmbt2 miRNA cluster causes enlargement of cloned mouse placentas

Somatic cell nuclear transfer (SCNT) in mammals is an inefficient process that is frequently associated with abnormal phenotypes, especially in placentas. Recent studies demonstrated that mouse SCNT placentas completely lack histone methylation (H3K27me3)-dependent imprinting, but how it affects placental development remains unclear. Here, we provide evidence that the loss of H3K27me3 imprinting is responsible for abnormal placental enlargement and low birth rates following SCNT, through upregulation of imprinted miRNAs. When we restore the normal paternal expression of H3K27me3-dependent imprinted genes (Sfmbt2, Gab1, and Slc38a4) in SCNT placentas by maternal knockout, the placentas remain enlarged. Intriguingly, correcting the expression of clustered miRNAs within the Sfmbt2 gene ameliorates the placental phenotype. Importantly, their target genes, which are confirmed to cause SCNT-like placental histology, recover their expression level. The birth rates increase about twofold. Thus, we identify loss of H3K27me3 imprinting as an epigenetic error that compromises embryo development following SCNT.

Somatic cell nuclear transfer (SCNT) frequently results in abnormal placenta development in cloned mice. Here the authors show that loss of histone methylation (H3K27me3) imprinting in clustered Sfmbt2 miRNAs contributes to SCNT placenta defect.

Amhr2-Cre-Mediated Global Tspo Knockout

Although the role of translocator protein (TSPO) in cholesterol transport in steroid-synthesizing cells has been studied extensively, recent studies of TSPO genetic depletion have questioned its role. Amhr2-Cre mice have been used to generate Leydig cell-specific Tspo conditional knockout (cKO) mice. Using the same Cre line, we were unable to generate Tspo cKO mice possibly because of genetic linkage between Tspo and Amhr2 and coexpression of Amhr2-Cre and Tspo in early embryonic development. We found that Amhr2-Cre is expressed during preimplantation stages, resulting in global heterozygous mice (gHE; Amhr2-Cre+/-,Tspo -/+). Two gHE mice were crossed, generating Amhr2-Cre-mediated Tspo global knockout (gKO; Tspo -/-) mice. We found that 33.3% of blastocysts at E3.5 to E4.5 showed normal morphology, whereas 66.7% showed delayed development, which correlates with the expected Mendelian proportions of Tspo +/+ (25%), Tspo -/- (25%), and Tspo +/- (50%) genotypes from crossing 2 Tspo -/+ mice. Adult Tspo gKO mice exhibited disturbances in neutral lipid homeostasis and reduced intratesticular and circulating testosterone levels, but no change in circulating basal corticosterone levels. RNA-sequencing data from mouse adrenal glands and lungs revealed transcriptome changes in response to the loss of TSPO, including changes in several cholesterol-binding and transfer proteins. This study demonstrates that Amhr2-Cre can be used to produce Tspo gKO mice instead of cKO, and can serve as a new global "Cre deleter." Moreover, our results show that Tspo deletion causes delayed preimplantation embryonic development, alters neutral lipid storage and steroidogenesis, and leads to transcriptome changes that may reflect compensatory mechanisms in response to the loss of function of TSPO.

Computational analysis of single-cell transcriptomics data elucidates the stabilization of Oct4 expression in the E3.25 mouse preimplantation embryo

Our computational analysis focuses on the 32- to 64-cell mouse embryo transition, Embryonic day (E3.25), whose study in literature is concentrated mainly on the search for an early onset of the second cell-fate decision, the specification of the inner cell mass (ICM) to primitive endoderm (PE) and epiblast (EPI). We analysed single-cell (sc) microarray transcriptomics data from E3.25 using Hierarchical Optimal k-Means (HOkM) clustering, and identified two groups of ICM cells: a group of cells from embryos with less than 34 cells (E3.25-LNCs), and another group of cells from embryos with more than 33 cells (E3.25-HNCs), corresponding to two developmental stages. Although we found massive underlying heterogeneity in the ICM cells at E3.25-HNC with over 3,800 genes with transcriptomics bifurcation, many of which are PE and EPI markers, we showed that the E3.25-HNCs are neither PE nor EPI. Importantly, analysing the differently expressed genes between the E3.25-LNCs and E3.25-HNCs, we uncovered a non-autonomous mechanism, based on a minimal number of four inner-cell contacts in the ICM, which activates Oct4 in the preimplantation embryo. Oct4 is highly expressed but unstable at E3.25-LNC, and stabilizes at high level at E3.25-HNC, with Bsg highly expressed, and the chromatin remodelling program initialised to establish an early naïve pluripotent state. Our results indicate that the pluripotent state we found to exist in the ICM at E3.25-HNC is the in vivo counterpart of a new, very early pluripotent state. We compared the transcriptomics profile of this in vivo E3.25-HNC pluripotent state, together with the profiles of E3.25-LNC, E3.5 EPI and E4.5 EPI cells, with the profiles of all embryonic stem cells (ESCs) available in the GEO database from the same platform (over 600 microarrays). The shortest distance between the set of inner cells (E3.25, E3.5 and E4.5) and the ESCs is between the E3.25-HNC cells and 2i + LIF ESCs; thus, the developmental transition from 33 to 34 cells decreases dramatically the distance with the naïve ground state of the 2i + LIF ESCs. We validated the E3.25 events through analysis of scRNA-seq data from early and late 32-cell ICM cells.

Transcriptional defects and reprogramming barriers in somatic cell nuclear reprogramming as revealed by single-embryo RNA sequencing

Background

Nuclear reprogramming reinstates totipotency or pluripotency in somatic cells by changing their gene transcription profile. This technology is widely used in medicine, animal husbandry and other industries. However, certain deficiencies severely restrict the applications of this technology.

Results

Using single-embryo RNA-seq, our study provides complete transcriptome blueprints of embryos generated by cumulus cell (CC) donor nuclear transfer (NT), embryos generated by mouse embryonic fibroblast (MEF) donor NT and in vivo embryos at each stage (zygote, 2-cell, 4-cell, 8-cell, morula, and blastocyst). According to the results from further analyses, NT embryos exhibit RNA processing and translation initiation defects during the zygotic genome activation (ZGA) period, and protein kinase activity and protein phosphorylation are defective during blastocyst formation. Two thousand three constant genes are not able to be reprogrammed in CCs and MEFs. Among these constant genes, 136 genes are continuously mis-transcribed throughout all developmental stages. These 136 differential genes may be reprogramming barrier genes (RBGs) and more studies are needed to identify.

Conclusions

These embryonic transcriptome blueprints provide new data for further mechanistic studies of somatic nuclear reprogramming. These findings may improve the efficiency of somatic cell nuclear transfer.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5091-1) contains supplementary material, which is available to authorized users.

Identification of key factors conquering developmental arrest of somatic cell cloned embryos by combining embryo biopsy and single-cell sequencing

Differentiated somatic cells can be reprogrammed into totipotent embryos through somatic cell nuclear transfer. However, most cloned embryos arrest at early stages and the underlying molecular mechanism remains largely unexplored. Here, we first developed a somatic cell nuclear transfer embryo biopsy system at two- or four-cell stage, which allows us to trace the developmental fate of the biopsied embryos precisely. Then, through single-cell transcriptome sequencing of somatic cell nuclear transfer embryos with different developmental fates, we identified that inactivation of Kdm4b, a histone H3 lysine 9 trimethylation demethylase, functions as a barrier for two-cell arrest of cloned embryos. Moreover, we discovered that inactivation of another histone demethylase Kdm5b accounts for the arrest of cloned embryos at the four-cell stage through single-cell analysis. Co-injection of Kdm4b and Kdm5b can restore transcriptional profiles of somatic cell nuclear transfer embryos and greatly improve the blastocyst development (over 95%) as well as the production of cloned mice. Our study therefore provides an effective approach to identify key factors responsible for the developmental arrest of somatic cell cloned embryos.

DNA Methylation Errors in Cloned Mouse Sperm by Germ Line Barrier Evasion

The germ line reprogramming barrier resets parental epigenetic modifications according to sex, conferring totipotency to mammalian embryos upon fertilization. However, it is not known whether epigenetic errors are committed during germ line reprogramming that are then transmitted to germ cells, and consequently to offspring. We addressed this question in the present study by performing a genome-wide DNA methylation analysis using a target postbisulfite sequencing method in order to identify DNA methylation errors in cloned mouse sperm. The sperm genomes of two somatic cell-cloned mice (CL1 and CL7) contained significantly higher numbers of differentially methylated CpG sites (P = 0.0045 and P = 0.0116). As a result, they had higher numbers of differentially methylated CpG islands. However, there was no evidence that these sites were transmitted to the sperm genome of offspring. These results suggest that DNA methylation errors resulting from embryo cloning are transmitted to the sperm genome by evading the germ line reprogramming barrier.

Does mouse embryo primordial germ cell activation start before implantation as suggested by single-cell transcriptomics dynamics?

STUDY HYPOTHESIS

Does primordial germ cell (PGC) activation start before mouse embryo implantation, and does the possible regulation of the DNA (cytosine-5-)-methyltransferase 3-like (Dnmt3l) by transcription factor AP-2, gamma (TCFAP2C) have a role in this activation and in the primitive endoderm (PE)-epiblast (EPI) lineage specification?

STUDY FINDING

A burst of expression of PGC markers, such as Dppa3/Stella, Ifitm2/Fragilis, Fkbp6 and Prdm4, is observed from embryonic day (E) 3.25, and some of them, together with the late germ cell markers Zp3, Mcf2 and Morc1, become restricted to the EPI subpopulation at E4.5, while the dynamics analysis of the PE-EPI transitions in the single-cell data suggests that TCFAP2C transitorily represses Dnmt3l in EPI cells at E3.5 and such repression is withdrawn with reactivation of Dnmt3l expression in PE and EPI cells at E4.5.

WHAT IS KNOWN ALREADY

In the mouse preimplantation embryo, cells with the same phenotype take different fates based on the orchestration between topological clues (cell polarity, positional history and division orientation) and gene regulatory rules (at transcriptomics and epigenomics level), prompting the proposal of positional, stochastic and combined models explaining the specification mechanism. PGC specification starts at E6.0-6.5 post-implantation. In view of the important role of DNA methylation in developmental events, the cross-talk between some transcription factors and DNA methyltransferases is of particular relevance. TCFAP2C has a CpG DNA methylation motif that is not methylated in pluripotent cells and that could potentially bind on DNMT3L, the stimulatory DNA methyltransferase co-factor that assists in the process of de novo DNA methylation. Chromatin-immunoprecipitation analysis has demonstrated that Dnmt3l is indeed a target of TCFAP2C.

STUDY DESIGN, SAMPLES/MATERIALS, METHODS

We aimed to assess the timing of early preimplantation events and to understand better the segregation of the inner cell mass (ICM) into PE and EPI. We designed a single-cell transcriptomics dynamics computational study to identify markers of the PE-EPI bifurcation in ICM cells through searching for statistically significant (using the Student's t-test method) differently expressed genes (DEGs) between PE and EPI cells from E3.5 to E4.5. The DEGs common for E3.5 and E4.5 were used as the markers defining the steady states. We collected microarray and next-generation sequencing transcriptomics data from public databases from bulk populations and single cells from mice at E3.25, E3.5 and E4.5. The results are based on three independent single-cell transcriptomics data sets, with a fold change of 3 and P-value <0.01 for the DEG selection.

MAIN RESULTS AND THE ROLE OF CHANCE

The dynamics analysis revealed new transitory E3.5 and steady PE and EPI markers. Among the transitory E3.5 PE markers (Dnmt3l, Dusp4, Cpne8, Akap13, Dcaf12l1, Aaed1, B4galt6, BC100530, Rnpc3, Tfpi, Lgalsl, Ckap4 and Fbxl20), several (Dusp4, Akap13, Cpn8, Dcaf12l1 and Tfpi) are related to the extracellular regulated kinase pathway. We also identified new transitory E3.5 EPI markers (Sgk1, Mal, Ubxn2a, Atg16l2, Gm13102, Tcfap2c, Hexb, Slc1a1, Svip, Liph and Mier3), six new stable PE markers (Sdc4, Cpn1, Dkk1, Havcr1, F2r/Par1 and Slc7a6os) as well as three new stable EPI markers (Zp3, Mcf2 and Hexb), which are known to be late stage germ cell markers. We found that mouse PGC marker activation starts at least at E3.25 preimplantation. The transcriptomics dynamics analyses support the regulation of Dnmt3l expression by TCFAP2C.

LIMITATIONS, REASONS FOR CAUTION

Since the regulation of Dnmt3l by TCFAP2C is based on computational prediction of DNA methylation motifs, Chip-Seq and transcriptomics data, functional studies are required to validate this result.

WIDER IMPLICATIONS OF THE FINDINGS

We identified a collection of previously undescribed E3.5-specific PE and EPI markers, and new steady PE and EPI markers. Identification of these genes, many of which encode cell membrane proteins, will facilitate the isolation and characterization of early PE and EPI populations. Since it is so well established in the literature that mouse PGC specification is a post-implantation event, it was surprising for us to see activation of PGC markers as early as E3.25 preimplantation, and identify the newly found steady EPI markers as late germ cell markers. The discovery of such early activation of PGC markers has important implications in the derivation of germ cells from pluripotent cells (embryonic stem cells or induced pluripotent stem cells), since the initial stages of such derivation resemble early development. The early activation of PGC markers points out the difficulty of separating PGC cells from pluripotent populations. Collectively, our results suggest that the combining of the precision of single-cell omics data with dynamic analysis of time-series data can establish the timing of some developmental stages as earlier than previously thought.

LARGE-SCALE DATA

Not applicable.

STUDY FUNDING AND COMPETING INTERESTS

This work was supported by grants DFG15/14 and DFG15/020 from Diputación Foral de Gipuzkoa (Spain), and grant II14/00016 from I + D + I National Plan 2013-2016 (Spain) and FEDER funds. The authors declare no conflict of interest.

Reshaping the transcriptional frontier: epigenetics and somatic cell nuclear transfer

Somatic-cell nuclear transfer (SCNT) experiments have paved the way to the field of cellular reprogramming. The demonstrated ability to clone over 20 different species to date has proven that the technology is robust but very inefficient, and is prone to developmental anomalies. Yet, the offspring from cloned animals exhibit none of the abnormalities of their parents, suggesting the low efficiency and high developmental mortality are epigenetic in origin. The epigenetic barriers to reprogramming somatic cells into a totipotent embryo capable of developing into a viable offspring are significant and varied. Despite their intimate relationship, chromatin structure and transcription are often not uniformly reprogramed after nuclear transfer, and many cloned embryos develop gene expression profiles that are hybrids between the donor cell and an embryonic blastomere. Recent advances in cellular reprogramming suggest that alteration of donor-cell chromatin structure towards that found in an normal embryo is actually the rate-limiting step in successful development of SCNT embryos. Here we review the literature relevant to the transformation of a somatic-cell nucleus into an embryo capable of full-term development. Interestingly, while resetting somatic transcription and associated epigenetic marks are absolutely required for development of SCNT embryos, life does not demand perfection.

Nuclear reprogramming of sperm and somatic nuclei in eggs and oocytes

Eggs and oocytes have a prominent ability to reprogram sperm nuclei for ensuring embryonic development. The reprogramming activity that eggs/oocytes intrinsically have towards sperm is utilised to reprogram somatic nuclei injected into eggs/oocytes in nuclear transfer (NT) embryos. NT embryos of various species can give rise to cloned animals, demonstrating that eggs/oocytes can confer totipotency even to somatic nuclei. However, many studies indicate that reprogramming of somatic nuclei is not as efficient as that of sperm nuclei. In this review, we explain how and why sperm and somatic nuclei are differentially reprogrammed in eggs/oocytes. Recent studies have shown that sperm chromatin is epigenetically modified to be adequate for early embryonic development, while somatic nuclei do not have such modifications. Moreover, epigenetic memories encoded in sperm chromatin are transgenerationally inherited, implying unique roles of sperm. We also discuss whether somatic nuclei can be artificially modified to acquire sperm-like chromatin states in order to increase the efficiency of nuclear reprogramming.