Memory Sequencing Reveals Heritable Single-Cell Gene Expression Programs Associated with Distinct Cellular Behaviors.

Non-genetic factors can cause individual cells to fluctuate substantially in gene expression levels over time. It remains unclear whether these fluctuations can persist for much longer than the time of one cell division. Current methods for measuring gene expression in single cells mostly rely on single time point measurements, making the duration of gene expression fluctuations or cellular memory difficult to measure. Here, we combined Luria and Delbrück's fluctuation analysis with population-based RNA sequencing (MemorySeq) for identifying genes transcriptome-wide whose fluctuations persist for several divisions. MemorySeq revealed multiple gene modules that expressed together in rare cells within otherwise homogeneous clonal populations. These rare cell subpopulations were associated with biologically distinct behaviors like proliferation in the face of anti-cancer therapeutics. The identification of non-genetic, multigenerational fluctuations can reveal new forms of biological memory in single cells and suggests that non-genetic heritability of cellular state may be a quantitative property.

ClampFISH 2.0 enables rapid, scalable amplified RNA detection in situ.

RNA labeling in situ has enormous potential to visualize transcripts and quantify their levels in single cells, but it remains challenging to produce high levels of signal while also enabling multiplexed detection of multiple RNA species simultaneously. Here, we describe clampFISH 2.0, a method that uses an inverted padlock design to efficiently detect many RNA species and exponentially amplify their signals at once, while also reducing the time and cost compared with the prior clampFISH method. We leverage the increased throughput afforded by multiplexed signal amplification and sequential detection to detect 10 different RNA species in more than 1 million cells. We also show that clampFISH 2.0 works in tissue sections. We expect that the advantages offered by clampFISH 2.0 will enable many applications in spatial transcriptomics.

Variability within rare cell states enables multiple paths toward drug resistance.

Molecular differences between individual cells can lead to dramatic differences in cell fate, such as death versus survival of cancer cells upon drug treatment. These originating differences remain largely hidden due to difficulties in determining precisely what variable molecular features lead to which cellular fates. Thus, we developed Rewind, a methodology that combines genetic barcoding with RNA fluorescence in situ hybridization to directly capture rare cells that give rise to cellular behaviors of interest. Applying Rewind to BRAFV600E melanoma, we trace drug-resistant cell fates back to single-cell gene expression differences in their drug-naive precursors (initial frequency of ~1:1,000-1:10,000 cells) and relative persistence of MAP kinase signaling soon after drug treatment. Within this rare subpopulation, we uncover a rich substructure in which molecular differences among several distinct subpopulations predict future differences in phenotypic behavior, such as proliferative capacity of distinct resistant clones after drug treatment. Our results reveal hidden, rare-cell variability that underlies a range of latent phenotypic outcomes upon drug exposure.

Genetic screening for single-cell variability modulators driving therapy resistance.

Cellular plasticity describes the ability of cells to transition from one set of phenotypes to another. In melanoma, transient fluctuations in the molecular state of tumor cells mark the formation of rare cells primed to survive BRAF inhibition and reprogram into a stably drug-resistant fate. However, the biological processes governing cellular priming remain unknown. We used CRISPR-Cas9 genetic screens to identify genes that affect cell fate decisions by altering cellular plasticity. We found that many factors can independently affect cellular priming and fate decisions. We discovered a new plasticity-based mode of increasing resistance to BRAF inhibition that pushes cells towards a more differentiated state. Manipulating cellular plasticity through inhibition of DOT1L before the addition of the BRAF inhibitor resulted in more therapy resistance than concurrent administration. Our results indicate that modulating cellular plasticity can alter cell fate decisions and may prove useful for treating drug resistance in other cancers.

Subcellular Detection of SARS-CoV-2 RNA in Human Tissue Reveals Distinct Localization in Alveolar Type 2 Pneumocytes and Alveolar Macrophages.

The widespread coronavirus disease 2019 (COVID-19) is caused by infection with the novel coronavirus SARS-CoV-2. Currently, we have limited understanding of which cells become infected with SARS-CoV-2 in human tissues and where viral RNA localizes on the subcellular level. Here, we present a platform for preparing autopsy tissue for visualizing SARS-CoV-2 RNA using RNA fluorescence in situ hybridization (FISH) with amplification by hybridization chain reaction. We developed probe sets that target different regions of SARS-CoV-2 (including ORF1a and N), as well as probe sets that specifically target SARS-CoV-2 subgenomic mRNAs. We validated these probe sets in cell culture and tissues (lung, lymph node, and placenta) from infected patients. Using this technology, we observe distinct subcellular localization patterns of the ORF1a and N regions. In human lung tissue, we performed multiplexed RNA FISH HCR for SARS-CoV-2 and cell-type-specific marker genes. We found viral RNA in cells containing the alveolar type 2 (AT2) cell marker gene (SFTPC) and the alveolar macrophage marker gene (MARCO) but did not identify viral RNA in cells containing the alveolar type 1 (AT1) cell marker gene (AGER). Moreover, we observed distinct subcellular localization patterns of viral RNA in AT2 cells and alveolar macrophages. In sum, we demonstrate the use of RNA FISH HCR for visualizing different RNA species from SARS-CoV-2 in cell lines and FFPE (formalin fixation and paraffin embedding) autopsy specimens. We anticipate that this platform could be broadly useful for studying SARS-CoV-2 pathology in tissues, as well as extended for other applications, including investigating the viral life cycle, viral diagnostics, and drug screening. IMPORTANCE Here, we developed an in situ RNA detection assay for RNA generated by the SARS-CoV-2 virus. We found viral RNA in lung, lymph node, and placenta samples from pathology specimens from COVID patients. Using high-magnification microscopy, we can visualize the subcellular distribution of these RNA in single cells.

Multiplexed detection of SARS-CoV-2 genomic and subgenomic RNA using in situ hybridization.

The widespread Coronavirus Disease 2019 (COVID-19) is caused by infection with the novel coronavirus SARS-CoV-2. Currently, we have a limited toolset available for visualizing SARS-CoV-2 in cells and tissues, particularly in tissues from patients who died from COVID-19. Generally, single-molecule RNA FISH techniques have shown mixed results in formalin fixed paraffin embedded tissues such as those preserved from human autopsies. Here, we present a platform for preparing autopsy tissue for visualizing SARS-CoV-2 RNA using RNA FISH with amplification by hybridization chain reaction (HCR). We developed probe sets that target different regions of SARS-CoV-2 (including ORF1a and N) as well as probe sets that specifically target SARS-CoV-2 subgenomic mRNAs. We validated these probe sets in cell culture and tissues (lung, lymph node, and placenta) from infected patients. Using this technology, we observe distinct subcellular localization patterns of the ORF1a and N regions, with the ORF1a concentrated around the nucleus and the N showing a diffuse distribution across the cytoplasm. In human lung tissue, we performed multiplexed RNA FISH HCR for SARS-CoV-2 and cell-type specific marker genes. We found viral RNA in cells containing the alveolar type 2 (AT2) cell marker gene (SFTPC) and the alveolar macrophage marker gene (MARCO), but did not identify viral RNA in cells containing the alveolar type 1 (AT1) cell marker gene (AGER). Moreover, we observed distinct subcellular localization patterns of viral RNA in AT2 cells and alveolar macrophages, consistent with phagocytosis of infected cells. In sum, we demonstrate the use of RNA FISH HCR for visualizing different RNA species from SARS-CoV-2 in cell lines and FFPE autopsy specimens. Furthermore, we multiplex this assay with probes for cellular genes to determine what cell-types are infected within the lung. We anticipate that this platform could be broadly useful for studying SARS-CoV-2 pathology in tissues as well as extended for other applications including investigating the viral life cycle, viral diagnostics, and drug screening.

Gene regulation gravitates toward either addition or multiplication when combining the effects of two signals.

Two different cell signals often affect transcription of the same gene. In such cases, it is natural to ask how the combined transcriptional response compares to the individual responses. The most commonly used mechanistic models predict additive or multiplicative combined responses, but a systematic genome-wide evaluation of these predictions is not available. Here, we analyzed the transcriptional response of human MCF-7 cells to retinoic acid and TGF-ß, applied individually and in combination. The combined transcriptional responses of induced genes exhibited a range of behaviors, but clearly favored both additive and multiplicative outcomes. We performed paired chromatin accessibility measurements and found that increases in accessibility were largely additive. There was some association between super-additivity of accessibility and multiplicative or super-multiplicative combined transcriptional responses, while sub-additivity of accessibility associated with additive transcriptional responses. Our findings suggest that mechanistic models of combined transcriptional regulation must be able to reproduce a range of behaviors.

Gene Networks with Transcriptional Bursting Recapitulate Rare Transient Coordinated High Expression States in Cancer.

Non-genetic transcriptional variability is a potential mechanism for therapy resistance in melanoma. Specifically, rare subpopulations of cells occupy a transient pre-resistant state characterized by coordinated high expression of several genes and survive therapy. How might these rare states arise and disappear within the population? It is unclear whether the canonical models of probabilistic transcriptional pulsing can explain this behavior, or if it requires special, hitherto unidentified mechanisms. We show that a minimal model of transcriptional bursting and gene interactions can give rise to rare coordinated high expression states. These states occur more frequently in networks with low connectivity and depend on three parameters. While entry into these states is initiated by a long transcriptional burst that also triggers entry of other genes, the exit occurs through independent inactivation of individual genes. Together, we demonstrate that established principles of gene regulation are sufficient to describe this behavior and argue for its more general existence. A record of this paper's transparent peer review process is included in the Supplemental Information.