Developmental trajectories of thalamic progenitors revealed by single-cell transcriptome profiling and Shh perturbation.

The thalamus is the principal information hub of the vertebrate brain, with essential roles in sensory and motor information processing, attention, and memory. The complex array of thalamic nuclei develops from a restricted pool of neural progenitors. We apply longitudinal single-cell RNA sequencing and regional abrogation of Sonic hedgehog (Shh) to map the developmental trajectories of thalamic progenitors, intermediate progenitors, and post-mitotic neurons as they coalesce into distinct thalamic nuclei. These data reveal that the complex architecture of the thalamus is established early during embryonic brain development through the coordinated action of four cell differentiation lineages derived from Shh-dependent and -independent progenitors. We systematically characterize the gene expression programs that define these thalamic lineages across time and demonstrate how their disruption upon Shh depletion causes pronounced locomotor impairment resembling infantile Parkinson's disease. These results reveal key principles of thalamic development and provide mechanistic insights into neurodevelopmental disorders resulting from thalamic dysfunction.

Single-cell antigen-specific landscape of CAR T infusion product identifies determinants of CD19-positive relapse in patients with ALL.

A notable number of acute lymphoblastic leukemia (ALL) patients develop CD19-positive relapse within 1 year after receiving chimeric antigen receptor (CAR) T cell therapy. It remains unclear if the long-term response is associated with the characteristics of CAR T cells in infusion products, hindering the identification of biomarkers to predict therapeutic outcomes. Here, we present 101,326 single-cell transcriptomes and surface protein landscape from the infusion products of 12 ALL patients. We observed substantial heterogeneity in the antigen-specific activation states, among which a deficiency of T helper 2 function was associated with CD19-positive relapse compared with durable responders (remission, >54 months). Proteomic data revealed that the frequency of early memory T cells, rather than activation or coinhibitory signatures, could distinguish the relapse. These findings were corroborated by independent functional profiling of 49 patients, and an integrative model was developed to predict the response. Our data unveil the molecular mechanisms that may inform strategies to boost specific T cell function to maintain long-term remission.

Single-cell multiomics dissection of basal and antigen-specific activation states of CD19-targeted CAR T cells.

BACKGROUND: Autologous T cells engineered to express a chimeric antigen receptor (CAR) specific for CD19 molecule have transformed the therapeutic landscape in patients with highly refractory leukemia and lymphoma, and the use of donor-generated allogeneic CAR T is paving the way for further breakthroughs in the treatment of cancer. However, it remains unknown how the intrinsic heterogeneities of these engineered cells mediate therapeutic efficacy and whether allogeneic products match the effectiveness of autologous therapies. METHODS: Using single-cell mRNA sequencing in conjunction with CITE-seq, we performed multiomics characterization of CAR T cells generated from healthy donor and patients with acute lymphoblastic leukemia. CAR T cells used in this study were manufactured at the University of Pennsylvania through lentiviral transduction with a CD19-4-1BB-CD3ζ construct. Besides the baseline condition, we engineered NIH-3T3 cells with human CD19 or mesothelin expression to conduct ex vivo antigen-specific or non-antigen stimulation of CAR T cells through 6-hour coculture at a 1:1 ratio. RESULTS: We delineated the global cellular and molecular CAR T landscape and identified that transcriptional CAR tonic signaling was regulated by a mixture of early activation, exhaustion signatures, and cytotoxic activities. On CD19 stimulation, we illuminated the disparities of CAR T cells derived from different origins and found that donor CAR T had more pronounced activation level in correlation with the upregulation of major histocompatibility complex class II genes compared with patient CAR T cells. This finding was independently validated in additional datasets from literature. Furthermore, GM-CSF(CSF2) expression was found to be associated with functional gene productions, but it induced little impact on the CAR T activation. CONCLUSIONS: Through integrated multiomics profiling and unbiased canonical pathway analyses, our results unveil heterogeneities in the transcriptional, phenotypic, functional, and metabolic profiles of donor and patient CAR T cells, providing mechanistic basis for ameliorating clinical outcomes and developing next-generation 'off- the-shelf' allogeneic products.

Mapping Rora expression in resting and activated CD4+ T cells.

The transcription factor Rora has been shown to be important for the development of ILC2 and the regulation of ILC3, macrophages and Treg cells. Here we investigate the role of Rora across CD4+ T cells in general, but with an emphasis on Th2 cells, both in vitro as well as in the context of several in vivo type 2 infection models. We dissect the function of Rora using overexpression and a CD4-conditional Rora-knockout mouse, as well as a RORA-reporter mouse. We establish the importance of Rora in CD4+ T cells for controlling lung inflammation induced by Nippostrongylus brasiliensis infection, and have measured the effect on downstream genes using RNA-seq. Using a systematic stimulation screen of CD4+ T cells, coupled with RNA-seq, we identify upstream regulators of Rora, most importantly IL-33 and CCL7. Our data suggest that Rora is a negative regulator of the immune system, possibly through several downstream pathways, and is under control of the local microenvironment.

Single-cell transcriptomic analysis of mIHC images via antigen mapping.

Highly multiplexed immunohistochemistry (mIHC) enables the staining and quantification of dozens of antigens in a tissue section with single-cell resolution. However, annotating cell populations that differ little in the profiled antigens or for which the antibody panel does not include specific markers is challenging. To overcome this obstacle, we have developed an approach for enriching mIHC images with single-cell RNA sequencing data, building upon recent experimental procedures for augmenting single-cell transcriptomes with concurrent antigen measurements. Spatially-resolved Transcriptomics via Epitope Anchoring (STvEA) performs transcriptome-guided annotation of highly multiplexed cytometry datasets. It increases the level of detail in histological analyses by enabling the systematic annotation of nuanced cell populations, spatial patterns of transcription, and interactions between cell types. We demonstrate the utility of STvEA by uncovering the architecture of poorly characterized cell types in the murine spleen using published cytometry and mIHC data of this organ.

SCNS: a graphical tool for reconstructing executable regulatory networks from single-cell genomic data.

BACKGROUND: Reconstruction of executable mechanistic models from single-cell gene expression data represents a powerful approach to understanding developmental and disease processes. New ambitious efforts like the Human Cell Atlas will soon lead to an explosion of data with potential for uncovering and understanding the regulatory networks which underlie the behaviour of all human cells. In order to take advantage of this data, however, there is a need for general-purpose, user-friendly and efficient computational tools that can be readily used by biologists who do not have specialist computer science knowledge. RESULTS: The Single Cell Network Synthesis toolkit (SCNS) is a general-purpose computational tool for the reconstruction and analysis of executable models from single-cell gene expression data. Through a graphical user interface, SCNS takes single-cell qPCR or RNA-sequencing data taken across a time course, and searches for logical rules that drive transitions from early cell states towards late cell states. Because the resulting reconstructed models are executable, they can be used to make predictions about the effect of specific gene perturbations on the generation of specific lineages. CONCLUSIONS: SCNS should be of broad interest to the growing number of researchers working in single-cell genomics and will help further facilitate the generation of valuable mechanistic insights into developmental, homeostatic and disease processes.

BTR: training asynchronous Boolean models using single-cell expression data.

BACKGROUND: Rapid technological innovation for the generation of single-cell genomics data presents new challenges and opportunities for bioinformatics analysis. One such area lies in the development of new ways to train gene regulatory networks. The use of single-cell expression profiling technique allows the profiling of the expression states of hundreds of cells, but these expression states are typically noisier due to the presence of technical artefacts such as drop-outs. While many algorithms exist to infer a gene regulatory network, very few of them are able to harness the extra expression states present in single-cell expression data without getting adversely affected by the substantial technical noise present. RESULTS: Here we introduce BTR, an algorithm for training asynchronous Boolean models with single-cell expression data using a novel Boolean state space scoring function. BTR is capable of refining existing Boolean models and reconstructing new Boolean models by improving the match between model prediction and expression data. We demonstrate that the Boolean scoring function performed favourably against the BIC scoring function for Bayesian networks. In addition, we show that BTR outperforms many other network inference algorithms in both bulk and single-cell synthetic expression data. Lastly, we introduce two case studies, in which we use BTR to improve published Boolean models in order to generate potentially new biological insights. CONCLUSIONS: BTR provides a novel way to refine or reconstruct Boolean models using single-cell expression data. Boolean model is particularly useful for network reconstruction using single-cell data because it is more robust to the effect of drop-outs. In addition, BTR does not assume any relationship in the expression states among cells, it is useful for reconstructing a gene regulatory network with as few assumptions as possible. Given the simplicity of Boolean models and the rapid adoption of single-cell genomics by biologists, BTR has the potential to make an impact across many fields of biomedical research.

Single-cell analysis of CD4+ T-cell differentiation reveals three major cell states and progressive acceleration of proliferation.

BACKGROUND: Differentiation of lymphocytes is frequently accompanied by cell cycle changes, interplay that is of central importance for immunity but is still incompletely understood. Here, we interrogate and quantitatively model how proliferation is linked to differentiation in CD4+ T cells. RESULTS: We perform ex vivo single-cell RNA-sequencing of CD4+ T cells during a mouse model of infection that elicits a type 2 immune response and infer that the differentiated, cytokine-producing cells cycle faster than early activated precursor cells. To dissect this phenomenon quantitatively, we determine expression profiles across consecutive generations of differentiated and undifferentiated cells during Th2 polarization in vitro. We predict three discrete cell states, which we verify by single-cell quantitative PCR. Based on these three states, we extract rates of death, division and differentiation with a branching state Markov model to describe the cell population dynamics. From this multi-scale modelling, we infer a significant acceleration in proliferation from the intermediate activated cell state to the mature cytokine-secreting effector state. We confirm this acceleration both by live imaging of single Th2 cells and in an ex vivo Th1 malaria model by single-cell RNA-sequencing. CONCLUSION: The link between cytokine secretion and proliferation rate holds both in Th1 and Th2 cells in vivo and in vitro, indicating that this is likely a general phenomenon in adaptive immunity.

Processing, visualising and reconstructing network models from single-cell data.

New single-cell technologies readily permit gene expression profiling of thousands of cells at single-cell resolution. In this review, we will discuss methods for visualisation and interpretation of single-cell gene expression data, and the computational analysis needed to go from raw data to predictive executable models of gene regulatory network function. We will focus primarily on single-cell real-time quantitative PCR and RNA-sequencing data, but much of what we cover will also be relevant to other platforms, such as the mass cytometry technology for high-dimensional single-cell proteomics.

Decoding the regulatory network of early blood development from single-cell gene expression measurements.

Reconstruction of the molecular pathways controlling organ development has been hampered by a lack of methods to resolve embryonic progenitor cells. Here we describe a strategy to address this problem that combines gene expression profiling of large numbers of single cells with data analysis based on diffusion maps for dimensionality reduction and network synthesis from state transition graphs. Applying the approach to hematopoietic development in the mouse embryo, we map the progression of mesoderm toward blood using single-cell gene expression analysis of 3,934 cells with blood-forming potential captured at four time points between E7.0 and E8.5. Transitions between individual cellular states are then used as input to develop a single-cell network synthesis toolkit to generate a computationally executable transcriptional regulatory network model of blood development. Several model predictions concerning the roles of Sox and Hox factors are validated experimentally. Our results demonstrate that single-cell analysis of a developing organ coupled with computational approaches can reveal the transcriptional programs that underpin organogenesis.

Single-cell analyses of regulatory network perturbations using enhancer-targeting TALEs suggest novel roles for PU.1 during haematopoietic specification.

Transcription factors (TFs) act within wider regulatory networks to control cell identity and fate. Numerous TFs, including Scl (Tal1) and PU.1 (Spi1), are known regulators of developmental and adult haematopoiesis, but how they act within wider TF networks is still poorly understood. Transcription activator-like effectors (TALEs) are a novel class of genetic tool based on the modular DNA-binding domains of Xanthomonas TAL proteins, which enable DNA sequence-specific targeting and the manipulation of endogenous gene expression. Here, we report TALEs engineered to target the PU.1-14kb and Scl+40kb transcriptional enhancers as efficient new tools to perturb the expression of these key haematopoietic TFs. We confirmed the efficiency of these TALEs at the single-cell level using high-throughput RT-qPCR, which also allowed us to assess the consequences of both PU.1 activation and repression on wider TF networks during developmental haematopoiesis. Combined with comprehensive cellular assays, these experiments uncovered novel roles for PU.1 during early haematopoietic specification. Finally, transgenic mouse studies confirmed that the PU.1-14kb element is active at sites of definitive haematopoiesis in vivo and PU.1 is detectable in haemogenic endothelium and early committing blood cells. We therefore establish TALEs as powerful new tools to study the functionality of transcriptional networks that control developmental processes such as early haematopoiesis.

Transcriptional hierarchies regulating early blood cell development.

Hematopoiesis represents one of the paradigmatic systems for studying stem cell biology, but our understanding of how the hematopoietic system develops during embryogenesis is still incomplete. While many lessons have been learned from studying the mouse embryo, embryonic stem cells have come to the fore as an alternative and more tractable model to recapitulate hematopoietic development. Here we review what is known about the embryonic origin of blood from these complementary systems and how transcription factor networks regulate the emergence of hematopoietic tissue from the mesoderm. Furthermore, we have performed an integrated analysis of genome-wide microarray and ChIP-seq data sets from mouse embryos and embryonic stem (ES) cell lines deficient in key regulators and demonstrate how this type of analysis can be used to reconstruct regulatory hierarchies that both confirm existing regulatory linkages and suggest additional interactions.