Xinglong Wu

Single-cell RNA-Seq profiling of human preimplantation embryos and embryonic stem cells

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.

The DNA methylation landscape of human early embryos

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.

Single-cell methylome landscapes of mouse embryonic stem cells and early embryos analyzed using reduced representation bisulfite sequencing

DNA methylation is crucial for a wide variety of biological processes, yet no technique suitable for the methylome analysis of DNA methylation at single-cell resolution is available. Here, we describe a methylome analysis technique that enables single-cell and single-base resolution DNA methylation analysis based on reduced representation bisulfite sequencing (scRRBS). The technique is highly sensitive and can detect the methylation status of up to 1.5 million CpG sites within the genome of an individual mouse embryonic stem cell (mESC). Moreover, we show that the technique can detect the methylation status of individual CpG sites in a haploid sperm cell in a digitized manner as either unmethylated or fully methylated. Furthermore, we show that the demethylation dynamics of maternal and paternal genomes after fertilization can be traced within the individual pronuclei of mouse zygotes. The demethylation process of the genic regions is faster than that of the intergenic regions in both male and female pronuclei. Our method paves the way for the exploration of the dynamic methylome landscapes of individual cells at single-base resolution during physiological processes such as embryonic development, or during pathological processes such as tumorigenesis.

Active and Passive Demethylation of Male and Female Pronuclear DNA in the Mammalian Zygote

The epigenomes of mammalian sperm and oocytes, characterized by gamete-specific 5-methylcytosine (5mC) patterns, are reprogrammed during early embryogenesis to establish full developmental potential. Previous studies have suggested that the paternal genome is actively demethylated in the zygote while the maternal genome undergoes subsequent passive demethylation via DNA replication during cleavage. Active demethylation is known to depend on 5mC oxidation by Tet dioxygenases and excision of oxidized bases by thymine DNA glycosylase (TDG). Here we show that both maternal and paternal genomes undergo widespread active and passive demethylation in zygotes before the first mitotic division. Passive demethylation was blocked by the replication inhibitor aphidicolin, and active demethylation was abrogated by deletion of Tet3 in both pronuclei. At actively demethylated loci, 5mCs were processed to unmodified cytosines. Surprisingly, the demethylation process was unaffected byxa0the deletion of TDG from the zygote, suggesting the existence of other demethylation mechanisms downstream of Tet3-mediated oxidation.

Whole-genome analysis of 5-hydroxymethylcytosine and 5-methylcytosine at base resolution in the human brain

Background5-methylcytosine (mC) can be oxidized by the tet methylcytosine dioxygenase (Tet) family of enzymes to 5-hydroxymethylcytosine (hmC), which is an intermediate of mC demethylation and may also be a stable epigenetic modification that influences chromatin structure. hmC is particularly abundant in mammalian brains but its function is currently unknown. A high-resolution hydroxymethylome map is required to fully understand the function of hmC in the human brain.ResultsWe present genome-wide and single-base resolution maps of hmC and mC in the human brain by combined application of Tet-assisted bisulfite sequencing and bisulfite sequencing. We demonstrate that hmCs increase markedly from the fetal to the adult stage, and in the adult brain, 13% of all CpGs are highly hydroxymethylated with strong enrichment at genic regions and distal regulatory elements. Notably, hmC peaks are identified at the 5′splicing sites at the exon-intron boundary, suggesting a mechanistic link between hmC and splicing. We report a surprising transcription-correlated hmC bias toward the sense strand and an mC bias toward the antisense strand of gene bodies. Furthermore, hmC is negatively correlated with H3K27me3-marked and H3K9me3-marked repressive genomic regions, and is more enriched at poised enhancers than active enhancers.ConclusionsWe provide single-base resolution hmC and mC maps in the human brain and our data imply novel roles of hmC in regulating splicing and gene expression. Hydroxymethylation is the main modification status for a large portion of CpGs situated at poised enhancers and actively transcribed regions, suggesting its roles in epigenetic tuning at these regions.

Microfluidic single-cell whole-transcriptome sequencing.

Significance RNA sequencing of single cells enables measurement of biological variation in heterogeneous cellular populations and dissection of transcriptome complexity that is masked in ensemble measurements of gene expression. The low quantity of RNA in a single cell, however, hinders efficient and consistent reverse transcription and amplification of cDNA, limiting accuracy and obscuring biological variation with high technical noise. We developed a microfluidic approach to prepare cDNA from single cells for high-throughput transcriptome sequencing. The microfluidic platform facilitates single-cell manipulation, minimizes contamination, and furthermore, provides improved detection sensitivity and measurement precision, which is necessary for differentiating biological variability from technical noise. Single-cell whole-transcriptome analysis is a powerful tool for quantifying gene expression heterogeneity in populations of cells. Many techniques have, thus, been recently developed to perform transcriptome sequencing (RNA-Seq) on individual cells. To probe subtle biological variation between samples with limiting amounts of RNA, more precise and sensitive methods are still required. We adapted a previously developed strategy for single-cell RNA-Seq that has shown promise for superior sensitivity and implemented the chemistry in a microfluidic platform for single-cell whole-transcriptome analysis. In this approach, single cells are captured and lysed in a microfluidic device, where mRNAs with poly(A) tails are reverse-transcribed into cDNA. Double-stranded cDNA is then collected and sequenced using a next generation sequencing platform. We prepared 94 libraries consisting of single mouse embryonic cells and technical replicates of extracted RNA and thoroughly characterized the performance of this technology. Microfluidic implementation increased mRNA detection sensitivity as well as improved measurement precision compared with tube-based protocols. With 0.2 M reads per cell, we were able to reconstruct a majority of the bulk transcriptome with 10 single cells. We also quantified variation between and within different types of mouse embryonic cells and found that enhanced measurement precision, detection sensitivity, and experimental throughput aided the distinction between biological variability and technical noise. With this work, we validated the advantages of an early approach to single-cell RNA-Seq and showed that the benefits of combining microfluidic technology with high-throughput sequencing will be valuable for large-scale efforts in single-cell transcriptome analysis.

Single-cell triple omics sequencing reveals genetic, epigenetic, and transcriptomic heterogeneity in hepatocellular carcinomas

Single-cell genome, DNA methylome, and transcriptome sequencing methods have been separately developed. However, to accurately analyze the mechanism by which transcriptome, genome and DNA methylome regulate each other, these omic methods need to be performed in the same single cell. Here we demonstrate a single-cell triple omics sequencing technique, scTrio-seq, that can be used to simultaneously analyze the genomic copy-number variations (CNVs), DNA methylome, and transcriptome of an individual mammalian cell. We show that large-scale CNVs cause proportional changes in RNA expression of genes within the gained or lost genomic regions, whereas these CNVs generally do not affect DNA methylation in these regions. Furthermore, we applied scTrio-seq to 25 single cancer cells derived from a human hepatocellular carcinoma tissue sample. We identified two subpopulations within these cells based on CNVs, DNA methylome, or transcriptome of individual cells. Our work offers a new avenue of dissecting the complex contribution of genomic and epigenomic heterogeneities to the transcriptomic heterogeneity within a population of cells.

Single-cell RNA-seq transcriptome analysis of linear and circular RNAs in mouse preimplantation embryos

Circular RNAs (circRNAs) are a new class of non-polyadenylated non-coding RNAs that may play important roles in many biological processes. Here we develop a single-cell universal poly(A)-independent RNA sequencing (SUPeR-seq) method to sequence both polyadenylated and non-polyadenylated RNAs from individual cells. This method exhibits robust sensitivity, precision and accuracy. We discover 2891 circRNAs and 913 novel linear transcripts in mouse preimplantation embryos and further analyze the abundance of circRNAs along development, the function of enriched genes, and sequence features of circRNAs. Our work is key to deciphering regulation mechanisms of circRNAs during mammalian early embryonic development.

Profiling DNA methylome landscapes of mammalian cells with single-cell reduced-representation bisulfite sequencing

The heterogeneity of DNA methylation within a population of cells necessitates DNA methylome profiling at single-cell resolution. Recently, we developed a single-cell reduced-representation bisulfite sequencing (scRRBS) technique in which we modified the original RRBS method by integrating all the experimental steps before PCR amplification into a single-tube reaction. These modifications enable scRRBS to provide digitized methylation information on ∼1 million CpG sites within an individual diploid mouse or human cell at single-base resolution. Compared with the single-cell bisulfite sequencing (scBS) technique, scRRBS covers fewer CpG sites, but it provides better coverage for CpG islands (CGIs), which are likely to be the most informative elements for DNA methylation. The entire procedure takes ∼3 weeks, and it requires strong molecular biology skills.

Single-cell multi-omics sequencing of mouse early embryos and embryonic stem cells

Single-cell epigenome sequencing techniques have recently been developed. However, the combination of different layers of epigenome sequencing in an individual cell has not yet been achieved. Here, we developed a single-cell multi-omics sequencing technology (single-cell COOL-seq) that can analyze the chromatin state/nucleosome positioning, DNA methylation, copy number variation and ploidy simultaneously from the same individual mammalian cell. We used this method to analyze the reprogramming of the chromatin state and DNA methylation in mouse preimplantation embryos. We found that within < 12 h of fertilization, each individual cell undergoes global genome demethylation together with the rapid and global reprogramming of both maternal and paternal genomes to a highly opened chromatin state. This was followed by decreased openness after the late zygote stage. Furthermore, from the late zygote to the 4-cell stage, the residual DNA methylation is preferentially preserved on intergenic regions of the paternal alleles and intragenic regions of maternal alleles in each individual blastomere. However, chromatin accessibility is similar between paternal and maternal alleles in each individual cell from the late zygote to the blastocyst stage. The binding motifs of several pluripotency regulators are enriched at distal nucleosome depleted regions from as early as the 2-cell stage. This indicates that the cis-regulatory elements of such target genes have been primed to an open state from the 2-cell stage onward, long before pluripotency is eventually established in the ICM of the blastocyst. Genes may be classified into homogeneously open, homogeneously closed and divergent states based on the chromatin accessibility of their promoter regions among individual cells. This can be traced to step-wise transitions during preimplantation development. Our study offers the first single-cell and parental allele-specific analysis of the genome-scale chromatin state and DNA methylation dynamics at single-base resolution in early mouse embryos and provides new insights into the heterogeneous yet highly ordered features of epigenomic reprogramming during this process.