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.

Genetic Variation in Type 1 Diabetes Reconfigures the 3D Chromatin Organization of T Cells and Alters Gene Expression.

Genetics is a major determinant of susceptibility to autoimmune disorders. Here, we examined whether genome organization provides resilience or susceptibility to sequence variations, and how this would contribute to the molecular etiology of an autoimmune disease. We generated high-resolution maps of linear and 3D genome organization in thymocytes of NOD mice, a model of type 1 diabetes (T1D), and the diabetes-resistant C57BL/6 mice. Multi-enhancer interactions formed at genomic regions harboring genes with prominent roles in T cell development in both strains. However, diabetes risk-conferring loci coalesced enhancers and promoters in NOD, but not C57BL/6 thymocytes. 3D genome mapping of NODxC57BL/6 F1 thymocytes revealed that genomic misfolding in NOD mice is mediated in cis. Moreover, immune cells infiltrating the pancreas of humans with T1D exhibited increased expression of genes located on misfolded loci in mice. Thus, genetic variation leads to altered 3D chromatin architecture and associated changes in gene expression that may underlie autoimmune pathology.

Diverse clonal fates emerge upon drug treatment of homogeneous cancer cells.

Even among genetically identical cancer cells, resistance to therapy frequently emerges from a small subset of those cells1-7. Molecular differences in rare individual cells in the initial population enable certain cells to become resistant to therapy7-9; however, comparatively little is known about the variability in the resistance outcomes. Here we develop and apply FateMap, a framework that combines DNA barcoding with single-cell RNA sequencing, to reveal the fates of hundreds of thousands of clones exposed to anti-cancer therapies. We show that resistant clones emerging from single-cell-derived cancer cells adopt molecularly, morphologically and functionally distinct resistant types. These resistant types are largely predetermined by molecular differences between cells before drug addition and not by extrinsic factors. Changes in the dose and type of drug can switch the resistant type of an initial cell, resulting in the generation and elimination of certain resistant types. Samples from patients show evidence for the existence of these resistant types in a clinical context. We observed diversity in resistant types across several single-cell-derived cancer cell lines and cell types treated with a variety of drugs. The diversity of resistant types as a result of the variability in intrinsic cell states may be a generic feature of responses to external cues.

p53 mediates target gene association with nuclear speckles for amplified RNA expression.

Nuclear speckles are prominent nuclear bodies that contain proteins and RNA involved in gene expression. Although links between nuclear speckles and gene activation are emerging, the mechanisms regulating association of genes with speckles are unclear. We find that speckle association of p53 target genes is driven by the p53 transcription factor. Focusing on p21, a key p53 target, we demonstrate that speckle association boosts expression by elevating nascent RNA amounts. p53-regulated speckle association did not depend on p53 transactivation functions but required an intact proline-rich domain and direct DNA binding, providing mechanisms within p53 for regulating gene-speckle association. Beyond p21, a substantial subset of p53 targets have p53-regulated speckle association. Strikingly, speckle-associating p53 targets are more robustly activated and occupy a distinct niche of p53 biology compared with non-speckle-associating p53 targets. Together, our findings illuminate regulated speckle association as a mechanism used by a transcription factor to boost gene expression.

Transcriptional Burst Initiation and Polymerase Pause Release Are Key Control Points of Transcriptional Regulation.

Transcriptional regulation occurs via changes to rates of different biochemical steps of transcription, but it remains unclear which rates are subject to change upon biological perturbation. Biochemical studies have suggested that stimuli predominantly affect the rates of RNA polymerase II (Pol II) recruitment and polymerase release from promoter-proximal pausing. Single-cell studies revealed that transcription occurs in discontinuous bursts, suggesting that features of such bursts like frequency and intensity could also be regulated. We combined Pol II chromatin immunoprecipitation sequencing (ChIP-seq) and single-cell transcriptional measurements to show that an independently regulated burst initiation step is required before polymerase recruitment can occur. Using a number of global and targeted transcriptional regulatory perturbations, we showed that biological perturbations regulated both burst initiation and polymerase pause release rates but seemed not to regulate polymerase recruitment rate. Our results suggest that transcriptional regulation primarily acts by changing the rates of burst initiation and polymerase pause release.

Group 1 Innate Lymphoid Cell Lineage Identity Is Determined by a cis-Regulatory Element Marked by a Long Non-coding RNA.

Commitment to the innate lymphoid cell (ILC) lineage is determined by Id2, a transcriptional regulator that antagonizes T and B cell-specific gene expression programs. Yet how Id2 expression is regulated in each ILC subset remains poorly understood. We identified a cis-regulatory element demarcated by a long non-coding RNA (lncRNA) that controls the function and lineage identity of group 1 ILCs, while being dispensable for early ILC development and homeostasis of ILC2s and ILC3s. The locus encoding this lncRNA, which we termed Rroid, directly interacted with the promoter of its neighboring gene, Id2, in group 1 ILCs. Moreover, the Rroid locus, but not the lncRNA itself, controlled the identity and function of ILC1s by promoting chromatin accessibility and deposition of STAT5 at the promoter of Id2 in response to interleukin (IL)-15. Thus, non-coding elements responsive to extracellular cues unique to each ILC subset represent a key regulatory layer for controlling the identity and function of ILCs.

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.

The BET Protein BRD2 Cooperates with CTCF to Enforce Transcriptional and Architectural Boundaries.

Bromodomain and extraterminal motif (BET) proteins are pharmacologic targets for the treatment of diverse diseases, yet the roles of individual BET family members remain unclear. We find that BRD2, but not BRD4, co-localizes with the architectural/insulator protein CCCTC-binding factor (CTCF) genome-wide. CTCF recruits BRD2 to co-bound sites whereas BRD2 is dispensable for CTCF occupancy. Disruption of a CTCF/BRD2-occupied element positioned between two unrelated genes enables regulatory influence to spread from one gene to another, suggesting that CTCF and BRD2 form a transcriptional boundary. Accordingly, single-molecule mRNA fluorescence in situ hybridization (FISH) reveals that, upon site-specific CTCF disruption or BRD2 depletion, expression of the two genes becomes increasingly correlated. HiC shows that BRD2 depletion weakens boundaries co-occupied by CTCF and BRD2, but not those that lack BRD2. These findings indicate that BRD2 supports boundary activity, and they raise the possibility that pharmacologic BET inhibitors can influence gene expression in part by perturbing domain boundary function.

Catheter ablation versus medical therapy in atrial fibrillation: an umbrella review of meta-analyses of randomized clinical trials.

This umbrella review synthesizes data from 17 meta-analyses investigating the comparative outcomes of catheter ablation (CA) and medical treatment (MT) for atrial fibrillation (AF). Outcomes assessed were mortality, risk of hospitalization, AF recurrence, cardiovascular events, pulmonary vein stenosis, major bleeding, and changes in left ventricular ejection fraction (LVEF) and MLHFQ score. The findings indicate that CA significantly reduces overall mortality and cardiovascular hospitalization with high strength of evidence. The risk of AF recurrence was notably lower with CA, with moderate strength of evidence. Two associations reported an increased risk of pulmonary vein stenosis and major bleeding with CA, supported by high strength of evidence. Improved LVEF and a positive change in MLHFQ were also associated with CA. Among patients with AF and heart failure, CA appears superior to MT for reducing mortality, improving LVEF, and reducing cardiovascular rehospitalizations. In nonspecific populations, CA reduced mortality and improved LVEF but had higher complication rates. Our findings suggest that CA might offer significant benefits in managing AF, particularly in patients with heart failure. However, the risk of complications, including pulmonary vein stenosis and major bleeding, is notable. Further research in understudied populations may help refine these conclusions.

Corrigendum: Rare cell variability and drug-induced reprogramming as a mode of cancer drug resistance.

This corrects the article DOI: 10.1038/nature22794.

What's Luck Got to Do with It: Single Cells, Multiple Fates, and Biological Nondeterminism.

The field of single-cell biology has morphed from a philosophical digression at its inception, to a playground for quantitative biologists, to a major area of biomedical research. The last several years have witnessed an explosion of new technologies, allowing us to apply even more of the modern molecular biology toolkit to single cells. Conceptual progress, however, has been comparatively slow. Here, we provide a framework for classifying both the origins of the differences between individual cells and the consequences of those differences. We discuss how the concept of "random" differences is context dependent, and propose that rigorous definitions of inputs and outputs may bring clarity to the discussion. We also categorize ways in which probabilistic behavior may influence cellular function, highlighting studies that point to exciting future directions in the field.

Enhancer Regulation of Transcriptional Bursting Parameters Revealed by Forced Chromatin Looping.

Mammalian genes transcribe RNA not continuously, but in bursts. Transcriptional output can be modulated by altering burst fraction or burst size, but how regulatory elements control bursting parameters remains unclear. Single-molecule RNA FISH experiments revealed that the ß-globin enhancer (LCR) predominantly augments transcriptional burst fraction of the ß-globin gene with modest stimulation of burst size. To specifically measure the impact of long-range chromatin contacts on transcriptional bursting, we forced an LCR-ß-globin promoter chromatin loop. We observed that raising contact frequencies increases burst fraction but not burst size. In cells in which two developmentally distinct LCR-regulated globin genes are cotranscribed in cis, burst sizes of both genes are comparable. However, allelic co-transcription of both genes is statistically disfavored, suggesting mutually exclusive LCR-gene contacts. These results are consistent with competition between the ß-type globin genes for LCR contacts and suggest that LCR-promoter loops are formed and released with rapid kinetics.

Visualizing allele-specific expression in single cells reveals epigenetic mosaicism in an H19 loss-of-imprinting mutant.

Imprinting is a classic mammalian epigenetic phenomenon that results in expression from a single parental allele. Imprinting defects can lead to inappropriate expression from the normally silenced allele, but it remains unclear whether every cell in a mutant organism follows the population average, which would have profound implications for human imprinting disorders. Here, we apply a new fluorescence in situ hybridization method that measures allele-specific expression in single cells to address this question in mutants exhibiting aberrant H19/Igf2 (insulin-like growth factor 2) imprinting. We show that mutant primary embryonic mouse fibroblasts are comprised of two subpopulations: one expressing both H19 alleles and another expressing only the maternal copy. Only in the latter cell population is Igf2 expression detected. Furthermore, the two subpopulations are stable in that cells do not interconvert between the two expression patterns. Combined small input methylation analysis and transcriptional imaging revealed that these two mutant subpopulations exhibit distinct methylation patterns at their imprinting control regions. Consistently, pharmacological inhibition of DNA methylation reduced the proportion of monoallelic cells. Importantly, we observed that the same two subpopulations are also present in vivo within murine cardiac tissue. Our results establish that imprinting disorders can display striking single-cell heterogeneity in their molecular phenotypes and suggest that such heterogeneity may underlie epigenetic mosaicism in human imprinting disorders.

A hyperactive transcriptional state marks genome reactivation at the mitosis-G1 transition.

During mitosis, RNA polymerase II (Pol II) and many transcription factors dissociate from chromatin, and transcription ceases globally. Transcription is known to restart in bulk by telophase, but whether de novo transcription at the mitosis-G1 transition is in any way distinct from later in interphase remains unknown. We tracked Pol II occupancy genome-wide in mammalian cells progressing from mitosis through late G1. Unexpectedly, during the earliest rounds of transcription at the mitosis-G1 transition, ∼50% of active genes and distal enhancers exhibit a spike in transcription, exceeding levels observed later in G1 phase. Enhancer-promoter chromatin contacts are depleted during mitosis and restored rapidly upon G1 entry but do not spike. Of the chromatin-associated features examined, histone H3 Lys27 acetylation levels at individual loci in mitosis best predict the mitosis-G1 transcriptional spike. Single-molecule RNA imaging supports that the mitosis-G1 transcriptional spike can constitute the maximum transcriptional activity per DNA copy throughout the cell division cycle. The transcriptional spike occurs heterogeneously and propagates to cell-to-cell differences in mature mRNA expression. Our results raise the possibility that passage through the mitosis-G1 transition might predispose cells to diverge in gene expression states.

Influence of short bouts of stair climbing on young adults during prolonged sitting on posture, discomfort, and musculoskeletal performance outcomes: a counterbalanced pilot randomised crossover trial.

Prolonged sitting is postulated to influence musculoskeletal performance (cervical flexor endurance, balance, and agility), discomfort and alter cervical spine angles during work-based computer use. Stair climbing breaks may be a great addition at typical and home offices however remain unexplored for its impact on musculoskeletal performance. In our counterbalanced pilot crossover trial, 24 adults were randomised to three interventions: (1) prolonged sitting, (2) interrupted by 2 min of self-paced, and (3) externally paced stair climbing for 2 h. Cervical spine angles were measured every 30 min while balance, agility, endurance, and discomfort were assessed before and after 120 min. Stair climbing interruptions have favourable effects on agility (F = 8.12, p = 0.009, ηp2 = 0.26) and musculoskeletal discomfort, but failed to improve other musculoskeletal outcomes associated with prolonged sitting. Brief stair climbing interruptions are effective in improving discomfort and agility while pragmatic trials are warranted for translated effects.

Stair climbing interventions are known for their cardiovascular benefits, however their impact on musculoskeletal performance remains unclear. In our randomised crossover pilot study, we explored the effects of stair climbing interruptions on cervical posture, endurance, agility, and balance during 120 min of prolonged sitting. Stair climbing interruptions were favourable in reducing time to complete agility tasks and neck/shoulder discomfort score compared to prolonged sitting. However other performance measures remained unaltered.

Changes in chromatin accessibility are not concordant with transcriptional changes for single-factor perturbations.

A major goal in the field of transcriptional regulation is the mapping of changes in the binding of transcription factors to the resultant changes in gene expression. Recently, methods for measuring chromatin accessibility have enabled us to measure changes in accessibility across the genome, which are thought to correspond to transcription factor-binding events. In concert with RNA-sequencing, these data in principle enable such mappings; however, few studies have looked at their concordance over short-duration treatments with specific perturbations. Here, we used tandem, bulk ATAC-seq, and RNA-seq measurements from MCF-7 breast carcinoma cells to systematically evaluate the concordance between changes in accessibility and changes in expression in response to retinoic acid and TGF-ß. We found two classes of genes whose expression showed a significant change: those that showed some changes in the accessibility of nearby chromatin, and those that showed virtually no change despite strong changes in expression. The peaks associated with genes in the former group had lower baseline accessibility prior to exposure to signal. Focusing the analysis specifically on peaks with motifs for transcription factors associated with retinoic acid and TGF-ß signaling did not reduce the lack of correspondence. Analysis of paired chromatin accessibility and gene expression data from distinct paths along the hematopoietic differentiation trajectory showed a much stronger correspondence, suggesting that the multifactorial biological processes associated with differentiation may lead to changes in chromatin accessibility that reflect rather than driving altered transcriptional status. Together, these results show many gene expression changes can happen independently of changes in the accessibility of local chromatin in the context of a single-factor perturbation.

Nature, nurture, or chance: stochastic gene expression and its consequences.

Gene expression is a fundamentally stochastic process, with randomness in transcription and translation leading to cell-to-cell variations in mRNA and protein levels. This variation appears in organisms ranging from microbes to metazoans, and its characteristics depend both on the biophysical parameters governing gene expression and on gene network structure. Stochastic gene expression has important consequences for cellular function, being beneficial in some contexts and harmful in others. These situations include the stress response, metabolism, development, the cell cycle, circadian rhythms, and aging.

Rare cell variability and drug-induced reprogramming as a mode of cancer drug resistance.

Therapies that target signalling molecules that are mutated in cancers can often have substantial short-term effects, but the emergence of resistant cancer cells is a major barrier to full cures. Resistance can result from secondary mutations, but in other cases there is no clear genetic cause, raising the possibility of non-genetic rare cell variability. Here we show that human melanoma cells can display profound transcriptional variability at the single-cell level that predicts which cells will ultimately resist drug treatment. This variability involves infrequent, semi-coordinated transcription of a number of resistance markers at high levels in a very small percentage of cells. The addition of drug then induces epigenetic reprogramming in these cells, converting the transient transcriptional state to a stably resistant state. This reprogramming begins with a loss of SOX10-mediated differentiation followed by activation of new signalling pathways, partially mediated by the activity of the transcription factors JUN and/or AP-1 and TEAD. Our work reveals the multistage nature of the acquisition of drug resistance and provides a framework for understanding resistance dynamics in single cells. We find that other cell types also exhibit sporadic expression of many of these same marker genes, suggesting the existence of a general program in which expression is displayed in rare subpopulations of cells.

Single mammalian cells compensate for differences in cellular volume and DNA copy number through independent global transcriptional mechanisms.

Individual mammalian cells exhibit large variability in cellular volume, even with the same absolute DNA content, and so must compensate for differences in DNA concentration in order to maintain constant concentration of gene expression products. Using single-molecule counting and computational image analysis, we show that transcript abundance correlates with cellular volume at the single-cell level due to increased global transcription in larger cells. Cell fusion experiments establish that increased cellular content itself can directly increase transcription. Quantitative analysis shows that this mechanism measures the ratio of cellular volume to DNA content, most likely through sequestration of a transcriptional factor to DNA. Analysis of transcriptional bursts reveals a separate mechanism for gene dosage compensation after DNA replication that enables proper transcriptional output during early and late S phase. Our results provide a framework for quantitatively understanding the relationships among DNA content, cell size, and gene expression variability in single cells.

Dynamic enhancer-gene body contacts during transcription elongation.

Enhancers govern transcription through multiple mechanisms, including the regulation of elongation by RNA polymerase II (RNAPII). We characterized the dynamics of looped enhancer contacts during synchronous transcription elongation. We found that many distal enhancers form stable contacts with their target promoters during the entire interval of elongation. Notably, we detected additional dynamic enhancer contacts throughout the gene bodies that track with elongating RNAPII and the leading edge of RNA synthesis. These results support a model in which the gene body changes its position relative to a stable enhancer-promoter complex, which has broad ramifications for enhancer function and architectural models of transcriptional elongation.