Hyeon-Jin Kim

Nucleoli, nuclear speckles and other compartments regulate transcription, RNA processing, and chromatin organization within the nucleus, yet the relationship of their morphology to developmental gene expression programs in vivo is poorly understood. Here, we develop a high-throughput Visual Cell Sorting (VCS) workflow for fixed cells and nuclei that combines antibody-based photoconversion; GPU-accelerated, real-time image analysis; and three-level single-cell combinatorial indexing RNA-seq (sci-RNA-seq3) to link nuclear compartment morphology to single-nucleus transcriptomes at embryo scale. We use VCS to analyze and sort over 1 million mouse embryo-derived nuclei by nucleolar, nuclear speckle, or nuclear size and construct a transcriptional atlas annotated with nuclear compartment phenotypes. Nuclear compartment size varies both between and within lineages and is shaped by proliferation and differentiation. In extracellular matrix protein-producing cell types such as fibroblasts, chondrocytes, and osteoblasts, nucleolar enlargement is uncoupled from cell cycle, and in erythroid cells exhibit a sharp nucleolar contraction preceding cell-cycle exit. We identify a 41-gene transcriptional signature whose expression tracks nucleolar size, enriched for ribosome biogenesis, mitochondrial metabolism, unfolded protein response, stress granule, and ubiquitin–proteasome pathway components. We used this nucleolar transcriptional signature to annotate mouse, zebrafish and human developmental atlases with nucleolar size information, revealing a conserved coupling between nucleolar activity and proteostasis programs. Our work establishes Visual Cell Sorting as a scalable platform for mapping image-based phenotypes to molecular programs; details the relationship between nuclear compartment phenotypes and development; and provides a transcriptional signature to estimate nucleolar size from existing single-cell datasets.

Nucleoli, nuclear speckles and other compartments regulate transcription, RNA processing, and chromatin organization within the nucleus, yet the relationship of their morphology to developmental gene expression programs in vivo is poorly understood. Here, we develop a high-throughput Visual Cell Sorting (VCS) workflow for fixed cells and nuclei that combines antibody-based photoconversion; GPU-accelerated, real-time image analysis; and three-level single-cell combinatorial indexing RNA-seq (sci-RNA-seq3) to link nuclear compartment morphology to single-nucleus transcriptomes at embryo scale. We use VCS to analyze and sort over 1 million mouse embryo-derived nuclei by nucleolar, nuclear speckle, or nuclear size and construct a transcriptional atlas annotated with nuclear compartment phenotypes. Nuclear compartment size varies both between and within lineages and is shaped by proliferation and differentiation. In extracellular matrix protein-producing cell types such as fibroblasts, chondrocytes, and osteoblasts, nucleolar enlargement is uncoupled from cell cycle, and in erythroid cells exhibit a sharp nucleolar contraction preceding cell-cycle exit. We identify a 41-gene transcriptional signature whose expression tracks nucleolar size, enriched for ribosome biogenesis, mitochondrial metabolism, unfolded protein response, stress granule, and ubiquitin–proteasome pathway components. We used this nucleolar transcriptional signature to annotate mouse, zebrafish and human developmental atlases with nucleolar size information, revealing a conserved coupling between nucleolar activity and proteostasis programs. Our work establishes Visual Cell Sorting as a scalable platform for mapping image-based phenotypes to molecular programs; details the relationship between nuclear compartment phenotypes and development; and provides a transcriptional signature to estimate nucleolar size from existing single-cell datasets.

Genetic variants often produce complex phenotypic effects that confound current assays and predictive models. We developed Variant in situ sequencing (VIS-seq), a pooled, image-based method that measures variant effects on molecular and cellular phenotypes in diverse cell types. Applying VIS-seq to ∼3,000 LMNA and PTEN variants yielded high-dimensional morphological profiles that captured variant-driven changes in protein abundance, localization, activity and cell architecture. We identified gain-of-function LMNA variants that reshape the nucleus and autism-associated PTEN variants that mislocalize. Morphological profiles predicted variant pathogenicity with near-perfect accuracy and distinguished autism-linked from tumor syndrome-linked PTEN variants. Most variants impacted a multidimensional continuum of phenotypes not recapitulated by any single functional readout. By linking protein variation to cell images at scale, we illuminate how variant effects cascade from molecular to subcellular to cell morphological phenotypes, providing a framework for resolving the complexity of variant function.

Single-cell RNA sequencing (scRNA-seq) offers a high-resolution molecular view into complex tissues, but suffers from high levels of technical noise which frustrates efforts to compare the gene expression programs of different cell types. “Spike-in” RNA standards help control for technical variation in scRNA-seq, but using them with recently developed, ultra-scalable scRNA-seq methods based on combinatorial indexing is not feasible. Here, we describe a simple and cost-effective method for normalizing transcript counts and subtracting technical variability that improves differential expression analysis in scRNA-seq. The method affixes a ladder of synthetic single-stranded DNA oligos to each cell that appears in its RNA-seq library. With improved normalization we explore chemical perturbations with broad or highly specific effects on gene regulation, including RNA pol II elongation, histone deacetylation, and activation of the glucocorticoid receptor. Our methods reveal that inhibiting histone deacetylation prevents cells from executing their canonical program of changes following glucocorticoid stimulation.

Single-cell Hi-C (scHi-C) interrogates genome-wide chromatin interaction in individual cells, allowing us to gain insights into 3D genome organization. However, the extremely sparse nature of scHi-C data poses a significant barrier to analysis, limiting our ability to tease out hidden biological information. In this work, we approach this problem by applying topic modeling to scHi-C data. Topic modeling is well-suited for discovering latent topics in a collection of discrete data. For our analysis, we generate nine different single-cell combinatorial indexed Hi-C (sci-Hi-C) libraries from five human cell lines (GM12878, H1Esc, HFF, IMR90, and HAP1), consisting over 19,000 cells. We demonstrate that topic modeling is able to successfully capture cell type differences from sci-Hi-C data in the form of “chromatin topics.” We further show enrichment of particular compartment structures associated with locus pairs in these topics. The genomes of higher organisms are intricately folded and organized in a dynamic manner that has strong implications for many biological processes. Each chromosome undergoes dramatic changes to their three dimensional conformation during the cell cycle, whereas the positioning of chromosomes within the nucleus plays an important role in controlling the activation of specific genes. Recently, it has become possible to investigate the 3D conformations of the genomes of individual cells using a high throughput sequencing assay called single cell Hi-C (scHi-C). However, data from these assays are sparse and noisy, making analysis and interpretation of scHi-C data challenging. In this work, we generated a scHi-C dataset of over 19,000 cells from five human cell lines and applied a natural language processing method called topic modeling to discover cell type-specific “chromatin” topics. We show that these topics can be used to distinguish between cells at different stages of the cell cycle and cells from different tissues based on the 3D conformation of their genomes, despite the sparsity of the data. We further show that the 3D conformations of single cells are linked to the expression of cell type-specific genes and to cell cycle-associated conformational patterns.