Cell systems最新文献

Fibrosis remains a major unmet medical need. Simplifying principles are needed to better understand fibrosis and to yield new therapeutic approaches. Fibrosis is driven by myofibroblasts that interact with macrophages. A mathematical cell-circuit model predicts two types of fibrosis: hot fibrosis driven by macrophages and myofibroblasts and cold fibrosis driven by myofibroblasts alone. Testing these concepts in cardiac fibrosis resulting from myocardial infarction (MI) and heart failure (HF), we revealed that acute MI leads to cold fibrosis whereas chronic injury (HF) leads to hot fibrosis. MI-driven cold fibrosis is conserved in pigs and humans. We computationally identified a vulnerability of cold fibrosis: the myofibroblast autocrine growth factor loop. Inhibiting this loop by targeting TIMP1 with neutralizing antibodies reduced myofibroblast proliferation and fibrosis post-MI in mice. Our study demonstrates the utility of the concepts of hot and cold fibrosis and the feasibility of a circuit-to-target approach to pinpoint a treatment strategy that reduces fibrosis. A record of this paper's transparent peer review process is included in the supplemental information.

Single-cell CRISPR screens link genetic perturbations to transcriptional states, but high-throughput methods connecting these induced changes to their regulatory foundations are limited. Here, we introduce Multiome Perturb-seq, extending single-cell CRISPR screens to simultaneously measure perturbation-induced changes in gene expression and chromatin accessibility. We apply Multiome Perturb-seq in a CRISPRi screen of 13 chromatin remodelers in human RPE-1 cells, achieving efficient assignment of sgRNA identities to single nuclei via an improved method for capturing barcode transcripts from nuclear RNA. We organize expression and accessibility measurements into coherent programs describing the integrated effects of perturbations on cell state, finding that ARID1A and SUZ12 knockdowns induce programs enriched for developmental features. Modeling of perturbation-induced heterogeneity connects accessibility changes to changes in gene expression, highlighting the value of multimodal profiling. Overall, our method provides a scalable and simply implemented system to dissect the regulatory logic underpinning cell state. A record of this paper's transparent peer review process is included in the supplemental information.

The mitogen-activated protein kinase (MAPK) pathway integrates growth factor signaling through extracellular signal-regulated kinase (ERK) to control cell proliferation. To study ERK dynamics, many researchers use an ERK activity kinase translocation reporter (KTR). Our study reveals that this ERK KTR also partially senses cyclin-dependent kinase 2 (CDK2) activity, making it appear as if ERK activity rises as cells progress through the cell cycle. Through single-cell time-lapse imaging, we identified a residual ERK KTR signal that was eliminated by selective CDK2 inhibitors, indicating crosstalk from CDK2 onto the ERK KTR. By contrast, EKAREN5, a FRET-based ERK sensor, showed no CDK2 crosstalk. A related p38 KTR is also partly affected by CDK2 activity. To address this, we developed linear and non-linear computational correction methods that subtract CDK2 signal from the ERK and p38 KTRs. These findings will allow for more accurate quantification of MAPK activities, especially for studies of actively cycling cells.

Treatment resistance poses a significant challenge in the care of cancer patients. Hirsch et al. applied computational and genomic approaches, examining gene expression dynamics from a mouse model of melanoma at single-cell resolution to reveal that semi-heritable non-genetic alterations in tumor cell populations confer adaptive resistance to treatment.

Computational prediction of the interaction of T cell receptors (TCRs) and their ligands is a grand challenge in immunology. Despite advances in high-throughput assays, specificity-labeled TCR data remain sparse. In other domains, the pre-training of language models on unlabeled data has been successfully used to address data bottlenecks. However, it is unclear how to best pre-train protein language models for TCR specificity prediction. Here, we introduce a TCR language model called SCEPTR (simple contrastive embedding of the primary sequence of T cell receptors), which is capable of data-efficient transfer learning. Through our model, we introduce a pre-training strategy combining autocontrastive learning and masked-language modeling, which enables SCEPTR to achieve its state-of-the-art performance. In contrast, existing protein language models and a variant of SCEPTR pre-trained without autocontrastive learning are outperformed by sequence alignment-based methods. We anticipate that contrastive learning will be a useful paradigm to decode the rules of TCR specificity. A record of this paper's transparent peer review process is included in the supplemental information.

While proliferating cells optimize their metabolism to produce biomass, the metabolic objectives of cells that perform non-proliferative tasks are unclear. The opposing requirements for optimizing each objective result in a trade-off that forces single cells to prioritize their metabolic needs and optimally allocate limited resources. Here, we present single-cell optimization objective and trade-off inference (SCOOTI), which infers metabolic objectives and trade-offs in biological systems by integrating bulk and single-cell omics data, using metabolic modeling and machine learning. We validated SCOOTI by identifying essential genes from CRISPR-Cas9 screens in embryonic stem cells, and by inferring the metabolic objectives of quiescent cells, during different cell-cycle phases. Applying this to embryonic cell states, we observed a decrease in metabolic entropy upon development. We further uncovered a trade-off between glutathione and biosynthetic precursors in one-cell zygote, two-cell embryo, and blastocyst cells, potentially representing a trade-off between pluripotency and proliferation. A record of this paper's transparent peer review process is included in the supplemental information.

Cancer progression is an evolutionary process driven by the selection of cells adapted to gain growth advantage. We present a formal study on the adaptation of gene expression in subclonal evolution. We model evolutionary changes in gene expression as stochastic Ornstein-Uhlenbeck processes, jointly leveraging the evolutionary history of subclones and single-cell expression data. Applying our model to sublines derived from single cells of a mouse melanoma revealed that sublines with distinct phenotypes are underlined by different patterns of gene expression adaptation, indicating non-genetic mechanisms of cancer evolution. Sublines previously observed to be resistant to anti-CTLA4 treatment showed adaptive expression of genes related to invasion and non-canonical Wnt signaling, whereas sublines that responded to treatment showed adaptive expression of genes related to proliferation and canonical Wnt signaling. Our results suggest that clonal phenotypes emerge as the result of specific adaptivity patterns of gene expression. A record of this paper's transparent peer review process is included in the supplemental information.

Deep learning is a promising strategy for modeling cis-regulatory elements. However, models trained on genomic sequences often fail to explain why the same transcription factor can activate or repress transcription in different contexts. To address this limitation, we developed an active learning approach to train models that distinguish between enhancers and silencers composed of binding sites for the photoreceptor transcription factor cone-rod homeobox (CRX). After training the model on nearly all bound CRX sites from the genome, we coupled synthetic biology with uncertainty sampling to generate additional rounds of informative training data. This allowed us to iteratively train models on data from multiple rounds of massively parallel reporter assays. The ability of the resulting models to discriminate between CRX sites with identical sequence but opposite functions establishes active learning as an effective strategy to train models of regulatory DNA. A record of this paper's transparent peer review process is included in the supplemental information.

In any given cell type, dozens of transcription factors (TFs) act in concert to control the activity of the genome by binding to specific DNA sequences in regulatory elements. Despite their considerable importance, we currently lack simple tools to directly measure the activity of many TFs in parallel. Massively parallel reporter assays (MPRAs) allow the detection of TF activities in a multiplexed fashion; however, we lack basic understanding to rationally design sensitive reporters for many TFs. Here, we use an MPRA to systematically optimize transcriptional reporters for 86 TFs and evaluate the specificity of all reporters across a wide array of TF perturbation conditions. We thus identified critical TF reporter design features and obtained highly sensitive and specific reporters for 62 TFs, many of which outperform available reporters. The resulting collection of "prime" TF reporters can be used to uncover TF regulatory networks and to illuminate signaling pathways. A record of this paper's transparent peer review process is included in the supplemental information.

T cells develop from hematopoietic progenitors in the thymus and protect against pathogens and cancer. However, the emergence of human T cell-competent blood progenitors and their subsequent specification to the T lineage have been challenging to capture in real time. Here, we leveraged a pluripotent stem cell differentiation system to understand the transcriptional dynamics and cell fate restriction events that underlie this critical developmental process. Time-resolved single-cell RNA sequencing revealed that downregulation of the multipotent hematopoietic program, upregulation of >90 lineage-associated transcription factors, and cell-cycle exit all occur within a highly coordinated developmental window. Gene-regulatory network inference uncovered a role for YBX1 in T lineage specification. We mapped the differentiation cell fate hierarchy using transcribed lineage barcoding and discovered that mast and myeloid potential bifurcate from each other early in hematopoiesis, upstream of T lineage restriction. Our systems-level analyses provide a quantitative, time-resolved model of human T cell fate specification. A record of this paper's transparent peer review process is included in the supplemental information.