Next-generation sequencing technology is a powerful tool for transcriptome analysis. However, under certain conditions, only a small amount of material is available, which requires more sensitive techniques that can preferably be used at the single-cell level. Here we describe a single-cell digital gene expression profiling assay. Using our mRNA-Seq assay with only a single mouse blastomere, we detected the expression of 75% (5,270) more genes than microarray techniques and identified 1,753 previously unknown splice junctions called by at least 5 reads. Moreover, 8-19% of the genes with multiple known transcript isoforms expressed at least two isoforms in the same blastomere or oocyte, which unambiguously demonstrated the complexity of the transcript variants at whole-genome scale in individual cells. Finally, for Dicer1(-/-) and Ago2(-/-) (Eif2c2(-/-)) oocytes, we found that 1,696 and 1,553 genes, respectively, were abnormally upregulated compared to wild-type controls, with 619 genes in common.

Measuring gene expression in individual cells is crucial for understanding the gene regulatory network controlling human embryonic development. Here we apply single-cell RNA sequencing (RNA-Seq) analysis to 124 individual cells from human preimplantation embryos and human embryonic stem cells (hESCs) at different passages. The number of maternally expressed genes detected in our data set is 22,687, including 8,701 long noncoding RNAs (lncRNAs), which represents a significant increase from 9,735 maternal genes detected previously by cDNA microarray. We discovered 2,733 novel lncRNAs, many of which are expressed in specific developmental stages. To address the long-standing question whether gene expression signatures of human epiblast (EPI) and in vitro hESCs are the same, we found that EPI cells and primary hESC outgrowth have dramatically different transcriptomes, with 1,498 genes showing differential expression between them. This work provides a comprehensive framework of the transcriptome landscapes of human early embryos and hESCs.

DNA methylation is a crucial element in the epigenetic regulation of mammalian embryonic development. However, its dynamic patterns have not been analysed at the genome scale in human pre-implantation embryos due to technical difficulties and the scarcity of required materials. Here we systematically profile the methylome of human early embryos from the zygotic stage through to post-implantation by reduced representation bisulphite sequencing and whole-genome bisulphite sequencing. We show that the major wave of genome-wide demethylation is complete at the 2-cell stage, contrary to previous observations in mice. Moreover, the demethylation of the paternal genome is much faster than that of the maternal genome, and by the end of the zygotic stage the genome-wide methylation level in male pronuclei is already lower than that in female pronuclei. The inverse correlation between promoter methylation and gene expression gradually strengthens during early embryonic development, reaching its peak at the post-implantation stage. Furthermore, we show that active genes, with the trimethylation of histone H3 at lysine 4 (H3K4me3) mark at the promoter regions in pluripotent human embryonic stem cells, are essentially devoid of DNA methylation in both mature gametes and throughout pre-implantation development. Finally, we also show that long interspersed nuclear elements or short interspersed nuclear elements that are evolutionarily young are demethylated to a milder extent compared to older elements in the same family and have higher abundance of transcripts, indicating that early embryos tend to retain higher residual methylation at the evolutionarily younger and more active transposable elements. Our work provides insights into the critical features of the methylome of human early embryos, as well as its functional relation to the regulation of gene expression and the repression of transposable elements.