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.

A novel approach to single-cell analysis reveals intrinsic differences in immune marker expression in unstimulated BALB/c and C57BL/6 macrophages.

The origin of functional heterogeneity among macrophages, key innate immune system components, is still debated. While mouse strains differ in their immune responses, the range of gene expression variation among their pre-stimulation macrophages is unknown. With a novel approach to scRNA-seq analysis, we reveal the gene expression variation in unstimulated macrophage populations from BALB/c and C57BL/6 mice. We show that intrinsic strain-to-strain differences are detectable before stimulation and we place the unstimulated single cells within the gene expression landscape of stimulated macrophages. C57BL/6 mice show stronger evidence of macrophage polarization than BALB/c mice, which may contribute to their relative resistance to pathogens. Our computational methods can be generally adopted to uncover biological variation between cell populations.

Bayesian inference of gene expression states from single-cell RNA-seq data.

Despite substantial progress in single-cell RNA-seq (scRNA-seq) data analysis methods, there is still little agreement on how to best normalize such data. Starting from the basic requirements that inferred expression states should correct for both biological and measurement sampling noise and that changes in expression should be measured in terms of fold changes, we here derive a Bayesian normalization procedure called Sanity (SAmpling-Noise-corrected Inference of Transcription activitY) from first principles. Sanity estimates expression values and associated error bars directly from raw unique molecular identifier (UMI) counts without any tunable parameters. Using simulated and real scRNA-seq datasets, we show that Sanity outperforms other normalization methods on downstream tasks, such as finding nearest-neighbor cells and clustering cells into subtypes. Moreover, we show that by systematically overestimating the expression variability of genes with low expression and by introducing spurious correlations through mapping the data to a lower-dimensional representation, other methods yield severely distorted pictures of the data.

Tead transcription factors differentially regulate cortical development.

Neural stem cells (NSCs) generate neurons of the cerebral cortex with distinct morphologies and functions. How specific neuron production, differentiation and migration are orchestrated is unclear. Hippo signaling regulates gene expression through Tead transcription factors (TFs). We show that Hippo transcriptional coactivators Yap1/Taz and the Teads have distinct functions during cortical development. Yap1/Taz promote NSC maintenance and Satb2+ neuron production at the expense of Tbr1+ neuron generation. However, Teads have moderate effects on NSC maintenance and do not affect Satb2+ neuron differentiation. Conversely, whereas Tead2 blocks Tbr1+ neuron formation, Tead1 and Tead3 promote this early fate. In addition, we found that Hippo effectors regulate neuronal migration to the cortical plate (CP) in a reciprocal fashion, that ApoE, Dab2 and Cyr61 are Tead targets, and these contribute to neuronal fate determination and migration. Our results indicate that multifaceted Hippo signaling is pivotal in different aspects of cortical development.

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.

Automated incorporation of pairwise dependency in transcription factor binding site prediction using dinucleotide weight tensors.

Gene regulatory networks are ultimately encoded by the sequence-specific binding of (TFs) to short DNA segments. Although it is customary to represent the binding specificity of a TF by a position-specific weight matrix (PSWM), which assumes each position within a site contributes independently to the overall binding affinity, evidence has been accumulating that there can be significant dependencies between positions. Unfortunately, methodological challenges have so far hindered the development of a practical and generally-accepted extension of the PSWM model. On the one hand, simple models that only consider dependencies between nearest-neighbor positions are easy to use in practice, but fail to account for the distal dependencies that are observed in the data. On the other hand, models that allow for arbitrary dependencies are prone to overfitting, requiring regularization schemes that are difficult to use in practice for non-experts. Here we present a new regulatory motif model, called dinucleotide weight tensor (DWT), that incorporates arbitrary pairwise dependencies between positions in binding sites, rigorously from first principles, and free from tunable parameters. We demonstrate the power of the method on a large set of ChIP-seq data-sets, showing that DWTs outperform both PSWMs and motif models that only incorporate nearest-neighbor dependencies. We also demonstrate that DWTs outperform two previously proposed methods. Finally, we show that DWTs inferred from ChIP-seq data also outperform PSWMs on HT-SELEX data for the same TF, suggesting that DWTs capture inherent biophysical properties of the interactions between the DNA binding domains of TFs and their binding sites. We make a suite of DWT tools available at dwt.unibas.ch, that allow users to automatically perform 'motif finding', i.e. the inference of DWT motifs from a set of sequences, binding site prediction with DWTs, and visualization of DWT 'dilogo' motifs.

Quantifying the strength of miRNA-target interactions.

We quantify the strength of miRNA-target interactions with MIRZA, a recently introduced biophysical model. We show that computationally predicted energies of interaction correlate strongly with the energies of interaction estimated from biochemical measurements of Michaelis-Menten constants. We further show that the accuracy of the MIRZA model can be improved taking into account recently emerged experimental data types. In particular, we use chimeric miRNA-mRNA sequences to infer a MIRZA-CHIMERA model and we provide a framework for inferring a similar model from measurements of rate constants of miRNA-mRNA interaction in the context of Argonaute proteins. Finally, based on a simple model of miRNA-based regulation, we discuss the importance of interaction energy and its variability between targets for the modulation of miRNA target expression in vivo.