We developed ASTAR-Seq (Assay for Single-cell Transcriptome and Accessibility Regions) integrated with automated microfluidic chips, which allows for parallel sequencing of transcriptome and chromatin accessibility within the same single-cell. Using ASTAR-Seq, we profiled 192 mESCs cultured in serum+LIF and 2i medium, 424 human cell lines including BJ, K562, JK1, and Jurkat, and 480 primary cells undergoing erythroblast differentiation. Integrative analysis using Coupled NMF identified the distinct sub-populations and uncovered sets of regulatory regions and the respective target genes determining their distinctions. Analysis of epigenetic regulomes further unravelled the key transcription factors responsible for the heterogeneity observed.

Cis Regulatory Elements (CREs) regulate the expression of the genes in their genomic neighborhoods and influence cellular processes such as cell-fate maintenance and differentiation. To date, there remain major gaps in the functional characterization of CREs and the identification of its target genes in the cellular native environment. In this study, we performed a Features Oriented CRISPR Utilized Systematic (FOCUS) screen of OCT4-bound CREs using CRISPR/Cas9 to identify functional enhancers important for pluripotency maintenance in mouse ES cells. From the initial 235 candidates tested, 16 CREs were identified to be essential stem cell enhancers. Using RNA-seq and genomic 4C-seq, we further uncovered a complex network of candidate CREs and their downstream target genes, which supports the growth and self-renewal of mESCs. Notably, an essential enhancer, CRE111, and its target, Lrrc31, form the important switch to modulate the LIF-JAK1-STAT3 signaling pathway.

To unravel the mechanism of human cellular reprogramming process at single-cell resolution, we performed parallel scRNA-Seq and scATAC-Seq analysis. Our analysis reveals that the cells undergoing reprogramming proceed in an asynchronous trajectory and diversify into heterogeneous sub-populations. BDD2-C8 fluorescent probe staining and negative staining for CD13, CD44 and CD201 markers, could enrich for the GDF3+ early reprogrammed cells. Combinatory usage of the surface markers enables the fine segregation of the early-intermediate cells with diverse reprogramming propensities. scATAC-Seq analysis further uncovered the genomic partitions and transcription factors responsible for the regulatory phasing of reprogramming process. Binary choice between a FOSL1 or a TEAD4-centric regulatory network determines the outcome of a successful reprogramming. Altogether, our study illuminates the multitude of diverse routes transversed by individual reprogramming cells and presents an integrative roadmap for identifying the mechanistic part-list of the reprogramming machinery.