Search

Supplemental Legends for Movies

Supplementary Movie 1

Supplementary Movie 2

Supplementary Movie 3

Supplementary Figures

Supplementary Movie 5

Supplementary Movie 8

Supplementary Movie 4

Source Code for analyses

Supplementary Movie 6

Tables of proteomics and CD4/8 staining

Supplementary Figure Legends

Supplementary Movie 9

Supplementary Movie 7

SummarySARS-CoV-2 Omicron BA.2.75 has diversified into multiple subvariants with additional spike mutations, and several are expanding in prevalence, particularly CH.1.1 and BN.1. Here, we investigated the viral receptor affinities and neutralization evasion properties of major BA.2.75 subvariants actively circulating in different regions worldwide. We found two distinct evolutionary pathways and three newly identified mutations that shaped the virological features of these subvariants. One phenotypic group exhibited a discernible decrease in viral receptor affinities, but a noteworthy increase in resistance to antibody neutralization, as exemplified by CH.1.1, which is apparently as resistant as XBB.1.5. In contrast, a second group demonstrated a substantial increase in viral receptor affinity but only a moderate increase in antibody evasion, as exemplified by BN.1. We also observed that all prevalent SARS-CoV-2 variants in the circulation presently, except for BN.1, exhibit profound levels of antibody evasion, suggesting this is the dominant determinant of virus transmissibility today.

Pluripotency can be induced in somatic cells by the expression of the four “Yamanaka” factors OCT4, KLF4, SOX2, and MYC. However, even in homogeneous conditions, usually only a rare subset of cells admit reprogramming, and the molecular characteristics of this subset remain unknown. Here, we apply retrospective clone tracing to identify and characterize the individual human fibroblast cells that are primed for reprogramming. These fibroblasts showed markers of increased cell cycle speed and decreased fibroblast activation. Knockdown of a fibroblast activation factor identified by our analysis led to increased reprogramming efficiency, identifying it as a barrier to reprogramming. Changing the frequency of reprogramming by inhibiting the activity of LSD1 led to an enlarging of the pool of cells that were primed for reprogramming. Our results show that even homogeneous cell populations can exhibit heritable molecular variability that can dictate whether individual rare cells will reprogram or not.

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.

Even amongst genetically identical cancer cells, therapy resistance often only emerges from a very small subset of those cells. Much effort has gone into uncovering the molecular differences in rare individual cells in the initial population that may allow certain cells to become therapy resistant; however, comparatively little is known about variability in the resistant outcomes themselves. 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 fate types. These different resistant types are largely predetermined by molecular differences between cells before addition of drug and not by extrinsic cell-specific microenvironmental factors. Changes in dose and kind of drug can, however, switch the resistant fate type of an initial cell, even resulting in the generation and elimination of certain fate types. Diversity in resistant fates was observed across several single-cell-derived cancer cell lines and types treated with a variety of drugs. Cell fate diversity as a result of variability in intrinsic cell states may be a generic feature of response to external cues.

Background Cardiac differentiation of human-induced pluripotent stem (hiPS) cells consistently produces a mixed population of cardiomyocytes and non-cardiac cell types, even when using well-characterized protocols. We sought to determine whether different cell types might result from intrinsic differences in hiPS cells prior to the onset of differentiation. Results By associating individual differentiated cells that share a common hiPS cell precursor, we tested whether expression variability is predetermined from the hiPS cell state. In a single experiment, cells that shared a progenitor were more transcriptionally similar to each other than to other cells in the differentiated population. However, when the same hiPS cells were differentiated in parallel, we did not observe high transcriptional similarity across differentiations. Additionally, we found that substantial cell death occurs during differentiation in a manner that suggested all cells were equally likely to survive or die, suggesting that there is no intrinsic selection bias for cells descended from particular hiPS cell progenitors. We thus wondered how cells grow spatially during differentiation, so we labeled cells by expression of marker genes and found that cells expressing the same marker tended to occur in patches. Our results suggest that cell type determination across multiple cell types, once initiated, is maintained in a cell-autonomous manner for multiple divisions. Conclusions Altogether, our results show that while substantial heterogeneity exists in the initial hiPS cell population, it is not responsible for the variability observed in differentiated outcomes; instead, factors specifying the various cell types likely act during a window that begins shortly after the seeding of hiPS cells for differentiation.

Cardiac directed differentiation of human induced pluripotent stem cells consistently produces a mixed population of cardiomyocytes and non-cardiac cell types even when using very well-characterized protocols. We wondered whether differentiated cell types might result from intrinsic differences in hiPS cells prior to the onset of differentiation. By associating individual differentiated cells that share a common hiPS cell precursor, we were able to test whether expression variability in differentiated cells was pre-determined from the hiPS cell state. Although within a single experiment, differentiated cells that share an hiPS cell progenitor were more transcriptionally similar to each other than to other cells in the differentiated population, when the same hiPS cells were differentiated in parallel, we did not observe high transcriptional similarity across differentiations. Additionally, we found that substantial cell death occurred during differentiation in a manner that suggested that all cells were equally likely to survive or die, suggesting that there was no intrinsic selection bias for cells descended from particular hiPS cell progenitors. These results led us to wonder about how cells grow out spatially during the directed differentiation process. Labeling cells by their expression of a few canonical cell type marker genes, we showed that cells expressing the same marker tended to occur in patches observable by visual inspection, suggesting that cell type determination across multiple cell types, once initiated, is maintained in a cell-autonomous manner for multiple divisions. Altogether, our results show that while there is substantial heterogeneity in the initial hiPS cell population, that heterogeneity is not responsible for heterogeneous outcomes, and that the window during which cell type specification occurs is likely to begin shortly after the seeding of hiPS cells for differentiation.

Our ability to identify the particular transcription factors that maintain cell type is limited. Identification of factors by their cell type-specific expression or their participation in developmental regulation has been only modestly successful. We hypothesized that because cell type is often resilient to perturbations, the transcriptional response to perturbations would identify identity-maintaining factors. We developed Perturbation Panel Profiling (P3) as a framework for perturbing cells in dozens of conditions and measuring gene expression responsiveness transcriptome-wide. Applying P3 to human iPSC-derived cardiac myocytes showed that transcription factors known to function in cardiac differentiation and maintenance were among the most frequently up-regulated (most responsive). We reasoned that one potential function of responsive genes may be to maintain cellular identity. We identified responsive transcription factors in fibroblasts using P3 and found that suppressing their expression led to enhanced reprogramming efficiency. We propose that responsiveness to perturbations is a property of factors that help maintain cellular identity.

Physical changes in skin are among the most visible signs of aging. We found that young dermal fibroblasts secrete high levels of extracellular matrix (ECM) constituents, including proteoglycans, glycoproteins, and cartilage-linking proteins. The most abundantly secreted was HAPLN1, a hyaluronic and proteoglycan link protein. HAPLN1 was lost in aged fibroblasts, resulting in a more aligned ECM that promoted metastasis of melanoma cells. Reconstituting HAPLN1 inhibited metastasis in an aged microenvironment, in 3-D skin reconstruction models, and in vivo. Intriguingly, aged fibroblast-derived matrices had the opposite effect on the migration of T cells, inhibiting their motility. HAPLN1 treatment of aged fibroblasts restored motility of mononuclear immune cells, while impeding that of polymorphonuclear immune cells, which in turn affected regulatory T-cell recruitment. These data suggest that although age-related physical changes in the ECM can promote tumor cell motility, they may adversely affect the motility of some immune cells, resulting in an overall change in the immune microenvironment. Understanding the physical changes in aging skin may provide avenues for more effective therapy for older patients with melanoma. Significance: These data shed light on the mechanochemical interactions that occur between aged skin, tumor, and immune cell populations, which may affect tumor metastasis and immune cell infiltration, with implications for the efficacy of current therapies for melanoma. See related commentary by Marie and Merlino, p. 19. This article is highlighted in the In This Issue feature, p. 1

Summary Identifying the particular transcription factors that maintain cell type in vitro is important for manipulating cell type. Identifying such transcription factors by their cell-type-specific expression or their involvement in developmental regulation has had limited success. We hypothesized that because cell type is often resilient to perturbations, the transcriptional response to perturbations would reveal identity-maintaining transcription factors. We developed perturbation panel profiling (P3) as a framework for perturbing cells across many conditions and measuring gene expression responsiveness transcriptome-wide. In human iPSC-derived cardiac myocytes, P3 showed that transcription factors important for cardiac myocyte differentiation and maintenance were among the most frequently upregulated (most responsive). We reasoned that one function of responsive genes may be to maintain cellular identity. We identified responsive transcription factors in fibroblasts using P3 and found that suppressing their expression led to enhanced reprogramming. We propose that responsiveness to perturbations is a property of transcription factors that help maintain cellular identity in vitro. A record of this paper’s transparent peer review process is included in the supplemental information.

Conversion of adenosine to inosine is a frequent type of RNA editing, but important details about the biology of this conversion remain unknown because of a lack of imaging tools. We developed inoFISH to directly visualize and quantify adenosine-to-inosine-edited transcripts in situ. We found that editing of the GRIA2, EIF2AK2, and NUP43 transcripts is uncorrelated with nuclear localization and paraspeckle association. Further, NUP43 exhibits constant editing levels between single cells, while GRIA2 editing levels vary.

Cellular plasticity describes the ability of cells to transition from one set of phenotypes to another. In the context of cancer therapeutics, plasticity refers to transient fluctuations in the molecular state of tumor cells, driving the formation of rare cells primed to survive drug treatment and ultimately reprogram into a stably resistant fate. However, the biological processes governing this cellular plasticity remain unknown. We used CRISPR/Cas9 genetic screens to reveal genes that affect cell fate decisions by altering cellular plasticity across a range of functional categories. We found that cellular plasticity and cell fate decision making can be decoupled in that factors can affect cell fate decisions in both plasticity-dependent and independent manners. We discovered a novel mode of altering resistance based on cellular plasticity that, contrary to known mechanisms, pushes cells towards a more differentiated state. We further confirmed our prediction that manipulating cellular plasticity before the addition of the main therapy would result in changes in therapy resistance more than concurrent administration. Together, our results indicate that identifying pathways modulating cellular plasticity has the potential to alter cell fate decisions and may provide a new avenue for treating drug resistance.

Nat. Biotechnol.; 10.1038/nbt.4286; corrected online 18 December 2018 In the supplementary information originally posted, there were errors in the primer sequences in Supplementary Table 1. The following sequences have been corrected online in all sheets: backbone sequences for MM2, MM2B and MM2C; adapter sequences for MM2, MM2B and MM2C; and MM2B_right arm sequence.

Methods for detecting single nucleic acids in cell and tissues, such as fluorescence in situ hybridization (FISH), are limited by relatively low signal intensity and nonspecific probe binding. Here we present click-amplifying FISH (clampFISH), a method for fluorescence detection of nucleic acids that achieves high specificity and high-gain (>400-fold) signal amplification. ClampFISH probes form a 'C' configuration upon hybridization to the sequence of interest in a double helical manner. The ends of the probes are ligated together using bio-orthogonal click chemistry, effectively locking the probes around the target. Iterative rounds of hybridization and click amplify the fluorescence intensity. We show that clampFISH enables the detection of RNA species with low-magnification microscopy and in RNA-based flow cytometry. Additionally, we show that the modular design of clampFISH probes allows multiplexing of RNA and DNA detection, that the locking mechanism prevents probe detachment in expansion microscopy, and that clampFISH can be applied in tissue samples.

Extensive cell-to-cell variation exists even among putatively identical cells, and there is great interest in understanding how the properties of transcription relate to this heterogeneity. Differential expression from the two gene copies in diploid cells could potentially contribute, yet our ability to measure from which gene copy individual RNAs originated remains limited, particularly in the context of tissues. Here, we demonstrate quantitative, single molecule allele-specific RNA FISH adapted for use on tissue sections, allowing us to determine the chromosome of origin of individual RNA molecules in formaldehyde-fixed tissues. We used this method to visualize the allele-specific expression of Xist and multiple autosomal genes in mouse kidney. By combining these data with mathematical modeling, we evaluated models for allele-specific heterogeneity, in particular demonstrating that apparent expression from only one of the alleles in single cells can arise as a consequence of low-level mRNA abundance and transcriptional bursting.

Biomarkers from multiple modalities have been shown to correlate with cognition in Parkinson's disease (PD) and in Alzheimer's disease (AD). However, the relationships of these markers with each other, and the use of multiple markers in concert to predict an outcome of interest, are areas that are much less explored. Our objectives in this study were (1) to evaluate relationships among 17 biomarkers previously reported to associate with cognition in PD or AD and (2) to test performance of a five-biomarker classifier trained to recognize AD in identifying PD with dementia (PDD). To do this, we evaluated a cross-sectional cohort of PD patients (n = 75) across a spectrum of cognitive abilities. All PD participants had 17 baseline biomarkers from clinical, genetic, biochemical, and imaging modalities measured, and correlations among biomarkers were assessed by Spearman's rho and by hierarchical clustering. We found that internal correlation among all 17 candidate biomarkers was modest, showing a maximum pairwise correlation coefficient of 0.51. However, a five-marker subset panel derived from AD (CSF total tau, CSF phosphorylated tau, CSF amyloid beta 42, APOE genotype, and SPARE-AD imaging score) discriminated cognitively normal PD patients vs. PDD patients with 80% accuracy, when employed in a classifier originally trained to recognize AD. Thus, an AD-derived biomarker signature may identify PDD patients with moderately high accuracy, suggesting mechanisms shared with AD in some PDD patients. Based on five measures readily obtained during life, this AD-derived signature may prove useful in identifying PDD patients most likely to respond to AD-based crossover therapies.

Small-scale molecular systems biology, by which we mean the understanding of a how a few parts work together to control a particular biological process, is predicated on the assumption that cellular regulation is arranged in a circuit-like structure. Results from the omics revolution have upset this vision to varying degrees by revealing a high degree of interconnectivity, making it difficult to develop a simple, circuit-like understanding of regulatory processes. We here outline the limitations of the small-scale systems biology approach with examples from research into genetic algorithms, genetics, transcriptional network analysis, and genomics. We also discuss the difficulties associated with deriving true understanding from the analysis of large data sets and propose that the development of new, intelligent, computational tools may point to a way forward. Throughout, we intentionally oversimplify and talk about things in which we have little expertise, and it is likely that many of our arguments are wrong on one level or another. We do believe, however, that developing a true understanding via molecular systems biology will require a fundamental rethinking of our approach, and our goal is to provoke thought along these lines.

Conditional mutagenesis is becoming a method of choice for studying gene function, but constructing conditional alleles is often laborious, limited by target gene structure, and at times, prone to incomplete conditional ablation. To address these issues, we developed a technology termed conditionals by inversion (COIN). Before activation, COINs contain an inverted module (COIN module) that lies inertly within the antisense strand of a resident gene. When inverted into the sense strand by a site-specific recombinase, the COIN module causes termination of the target gene’s transcription and simultaneously provides a reporter for tracking this event. COIN modules can be inserted into natural introns (intronic COINs) or directly into coding exons as part of an artificial intron (exonic COINs), greatly simplifying allele design and increasing flexibility over previous conditional KO approaches. Detailed analysis of over 20 COIN alleles establishes the reliability of the method and its broad applicability to any gene, regardless of exon–intron structure. Our extensive testing provides rules that help ensure success of this approach and also explains why other currently available conditional approaches often fail to function optimally. Finally, the ability to split exons using the COIN’s artificial intron opens up engineering modalities for the generation of multifunctional alleles.