Temporal and sequential transcriptional dynamics define lineage shifts in corticogenesis

The cerebral cortex contains billions of neurons, and their disorganization or misspecification leads to neurodevelopmental disorders. Understanding how the plethora of projection neuron subtypes are generated by cortical neural stem cells (NSCs) is a major challenge. Here, we focused on elucidating the transcriptional landscape of murine embryonic NSCs, basal progenitors (BPs), and newborn neurons (NBNs) throughout cortical development. We uncover dynamic shifts in transcriptional space over time and heterogeneity within each progenitor population. We identified signature hallmarks of NSC, BP, and NBN clusters and predict active transcriptional nodes and networks that contribute to neural fate specification. We find that the expression of receptors, ligands, and downstream pathway components is highly dynamic over time and throughout the lineage implying differential responsiveness to signals. Thus, we provide an expansive compendium of gene expression during cortical development that will be an invaluable resource for studying neural developmental processes and neurodevelopmental disorders

TGF-β induces oncofetal fibronectin, which in turn modulates TGF-β superfamily signaling in endothelial cells

Gene splicing profiles are frequently altered in cancer, and the splice variants of fibronectin (FN) that contain the extra-domains A (EDA) or B (EDB), referred to as EDA+FN or EDB+FN, are highly upregulated in tumor vasculature. Transforming growth factor β (TGF-β) signaling has been attributed a pivotal role in glioblastoma, with TGF-β promoting angiogenesis and vessel remodeling. By using immunohistochemistry staining, we observed that the oncofetal FN isoforms EDA+FN and EDB+FN are expressed in glioblastoma vasculature. single-cell gene expression profiling of tumors by using CD31 and α-smooth muscle actin (αSMA) as markers for endothelial cells, and pericytes and vascular smooth muscle cells (VSMCs), respectively, confirmed the predominant expression of FN, EDA+FN and EDB+FN in the vascular compartment of glioblastoma. Specifically, within the CD31-positive cell population, we identified a positive correlation between the expression of EDA+FN and EDB+FN, and of molecules associated with TGF-β signaling. Further, TGF-β induced EDA+FN and EDB+FN in human cerebral microvascular endothelial cells and glioblastoma-derived endothelial cells in a SMAD3- and SMAD4-dependent manner. In turn, we found that FN modulated TGF-β superfamily signaling in endothelial cells via the EDA and EDB, pointing towards a bidirectional influence of oncofetal FN and TGF-β superfamily signaling

Single-cell mRNA profiling reveals the hierarchical response of miRNA targets to miRNA induction

miRNAs are small RNAs that regulate gene expression post‐transcriptionally. By repressing the translation and promoting the degradation of target mRNAs, miRNAs may reduce the cell‐to‐cell variability in protein expression, induce correlations between target expression levels, and provide a layer through which targets can influence each other's expression as “competing RNAs” (ceRNAs). However, experimental evidence for these behaviors is limited. Combining mathematical modeling with RNA sequencing of individual human embryonic kidney cells in which the expression of two distinct miRNAs was induced over a wide range, we have inferred parameters describing the response of hundreds of miRNA targets to miRNA induction. Individual targets have widely different response dynamics, and only a small proportion of predicted targets exhibit high sensitivity to miRNA induction. Our data reveal for the first time the response parameters of the entire network of endogenous miRNA targets to miRNA induction, demonstrating that miRNAs correlate target expression and at the same time increase the variability in expression of individual targets across cells. The approach is generalizable to other miRNAs and post‐transcriptional regulators to improve the understanding of gene expression dynamics in individual cell types.ISSN:1744-429

Single-cell Analysis Reveals Inter- and Intratumour Heterogeneity in Metastatic Breast Cancer

Metastasis is the leading cause of cancer-related deaths of breast cancer patients. Some cancer cells in a tumour go through successive steps, referred to as the metastatic cascade, and give rise to metastases at a distant site. We know that the plasticity and heterogeneity of cancer cells play critical roles in metastasis but the precise underlying molecular mechanisms remain elusive. Here we aimed to identify molecular mechanisms of metastasis during colonization, one of the most important yet poorly understood steps of the cascade. We performed single-cell RNA-Seq (scRNA-Seq) on tumours and matched lung macrometastases of patient-derived xenografts of breast cancer. After correcting for confounding factors such as the cell cycle and the percentage of detected genes (PDG), we identified cells in three states in both tumours and metastases. Gene-set enrichment analysis revealed biological processes specific to proliferation and invasion in two states. Our findings suggest that these states are a balance between epithelial-to-mesenchymal (EMT) and mesenchymal-to-epithelial transitions (MET) traits that results in so-called partial EMT phenotypes. Analysis of the top differentially expressed genes (DEGs) between these cell states revealed a common set of partial EMT transcription factors (TFs) controlling gene expression, including ZNF750, OVOL2, TP63, TFAP2C and HEY2. Our data suggest that the TFs related to EMT delineate different cell states in tumours and metastases. The results highlight the marked interpatient heterogeneity of breast cancer but identify common features of single cells from five models of metastatic breast cancer.ISSN:1083-3021ISSN:1573-703

Temporal and sequential transcriptional dynamics define lineage shifts in corticogenesis

The cerebral cortex contains billions of neurons, and their disorganization or misspecification leads to neurodevelopmental disorders. Understanding how the plethora of projection neuron subtypes are generated by cortical neural stem cells (NSCs) is a major challenge. Here, we focused on elucidating the transcriptional landscape of murine embryonic NSCs, basal progenitors (BPs), and newborn neurons (NBNs) throughout cortical development. We uncover dynamic shifts in transcriptional space over time and heterogeneity within each progenitor population. We identified signature hallmarks of NSC, BP, and NBN clusters and predict active transcriptional nodes and networks that contribute to neural fate specification. We find that the expression of receptors, ligands, and downstream pathway components is highly dynamic over time and throughout the lineage implying differential responsiveness to signals. Thus, we provide an expansive compendium of gene expression during cortical development that will be an invaluable resource for studying neural developmental processes and neurodevelopmental disorders.ISSN:0261-4189ISSN:1460-207

SCIM: Universal Single-Cell Matching with Unpaired Feature Sets

Motivation Recent technological advances have led to an increase in the production and availability of single-cell data. The ability to integrate a set of multi-technology measurements would allow the identification of biologically or clinically meaningful observations through the unification of the perspectives afforded by each technology. In most cases, however, profiling technologies consume the used cells and thus pairwise correspondences between datasets are lost. Due to the sheer size single-cell datasets can acquire, scalable algorithms that are able to universally match single-cell measurements carried out in one cell to its corresponding sibling in another technology are needed. Results We propose Single-Cell data Integration via Matching (SCIM), a scalable approach to recover such correspondences in two or more technologies. SCIM assumes that cells share a common (low-dimensional) underlying structure and that the underlying cell distribution is approximately constant across technologies. It constructs a technology-invariant latent space using an auto-encoder framework with an adversarial objective. Multi-modal datasets are integrated by pairing cells across technologies using a bipartite matching scheme that operates on the low-dimensional latent representations. We evaluate SCIM on a simulated cellular branching process and show that the cell-to-cell matches derived by SCIM reflect the same pseudotime on the simulated dataset. Moreover, we apply our method to two real-world scenarios, a melanoma tumor sample and a human bone marrow sample, where we pair cells from a scRNA dataset to their sibling cells in a CyTOF dataset achieving 93% and 84% cell-matching accuracy for each one of the samples respectively. Availability https://github.com/ratschlab/sci

Establishing standardized immune phenotyping of metastatic melanoma by digital pathology

CD8+ tumor-infiltrating T cells can be regarded as one of the most relevant predictive biomarkers in immune-oncology. Highly infiltrated tumors, referred to as inflamed (clinically “hot”), show the most favorable response to immune checkpoint inhibitors in contrast to tumors with a scarce immune infiltrate called immune desert or excluded (clinically “cold”). Nevertheless, quantitative and reproducible methods examining their prevalence within tumors are lacking. We therefore established a computational diagnostic algorithm to quantitatively measure spatial densities of tumor-infiltrating CD8+ T cells by digital pathology within the three known tumor compartments as recommended by the International Immuno-Oncology Biomarker Working Group in 116 prospective metastatic melanomas of the Swiss Tumor Profiler cohort. Workflow robustness was confirmed in 33 samples of an independent retrospective validation cohort. The introduction of the intratumoral tumor center compartment proved to be most relevant for establishing an immune diagnosis in metastatic disease, independent of metastatic site. Cut-off values for reproducible classification were defined and successfully assigned densities into the respective immune diagnostic category in the validation cohort with high sensitivity, specificity, and precision. We provide a robust diagnostic algorithm based on intratumoral and stromal CD8+ T-cell densities in the tumor center compartment that translates spatial densities of tumor-infiltrating CD8+ T cells into the clinically relevant immune diagnostic categories “inflamed”, “excluded”, and “desert”. The consideration of the intratumoral tumor center compartment allows immune phenotyping in the clinically highly relevant setting of metastatic lesions, even if the invasive margin compartment is not captured in biopsy material