Phenotypic variation of transcriptomic cell types in mouse motor cortex.

Cortical neurons exhibit extreme diversity in gene expression as well as in morphological and electrophysiological properties1,2. Most existing neural taxonomies are based on either transcriptomic3,4 or morpho-electric5,6 criteria, as it has been technically challenging to study both aspects of neuronal diversity in the same set of cells7. Here we used Patch-seq8 to combine patch-clamp recording, biocytin staining, and single-cell RNA sequencing of more than 1,300 neurons in adult mouse primary motor cortex, providing a morpho-electric annotation of almost all transcriptomically defined neural cell types. We found that, although broad families of transcriptomic types (those expressing Vip, Pvalb, Sst and so on) had distinct and essentially non-overlapping morpho-electric phenotypes, individual transcriptomic types within the same family were not well separated in the morpho-electric space. Instead, there was a continuum of variability in morphology and electrophysiology, with neighbouring transcriptomic cell types showing similar morpho-electric features, often without clear boundaries between them. Our results suggest that neuronal types in the neocortex do not always form discrete entities. Instead, neurons form a hierarchy that consists of distinct non-overlapping branches at the level of families, but can form continuous and correlated transcriptomic and morpho-electrical landscapes within families.

Genomic encoding of transcriptional burst kinetics.

Mammalian gene expression is inherently stochastic1,2, and results in discrete bursts of RNA molecules that are synthesized from each allele3-7. Although transcription is known to be regulated by promoters and enhancers, it is unclear how cis-regulatory sequences encode transcriptional burst kinetics. Characterization of transcriptional bursting, including the burst size and frequency, has mainly relied on live-cell4,6,8 or single-molecule RNA fluorescence in situ hybridization3,5,8,9 recordings of selected loci. Here we determine transcriptome-wide burst frequencies and sizes for endogenous mouse and human genes using allele-sensitive single-cell RNA sequencing. We show that core promoter elements affect burst size and uncover synergistic effects between TATA and initiator elements, which were masked at mean expression levels. Notably, we provide transcriptome-wide evidence that enhancers control burst frequencies, and demonstrate that cell-type-specific gene expression is primarily shaped by changes in burst frequencies. Together, our data show that burst frequency is primarily encoded in enhancers and burst size in core promoters, and that allelic single-cell RNA sequencing is a powerful model for investigating transcriptional kinetics.

Compound models and Pearson residuals for normalization of single-cell RNA-seq data without UMIs.

Before downstream analysis can reveal biological signals in single-cell RNA sequencing data, normalization and variance stabilization are required to remove technical noise. Recently, Pearson residuals based on negative binomial models have been suggested as an efficient normalization approach. These methods were developed for UMI-based sequencing protocols, where unique molecular identifiers (UMIs) help to remove PCR amplification noise by keeping track of the original molecules. In contrast, full-length protocols such as Smart-seq2 lack UMIs and retain amplification noise, making negative binomial models inapplicable. Here, we extend Pearson residuals to such read count data by modeling them as a compound process: we assume that the captured RNA molecules follow the negative binomial distribution, but are replicated according to an amplification distribution. Based on this model, we introduce compound Pearson residuals and show that they can be analytically obtained without explicit knowledge of the amplification distribution. Further, we demonstrate that compound Pearson residuals lead to a biologically meaningful gene selection and low-dimensional embeddings of complex Smart-seq2 datasets. Finally, we empirically study amplification distributions across several sequencing protocols, and suggest that they can be described by a broken power law. We show that the resulting compound distribution captures overdispersion and zero-inflation patterns characteristic of read count data. In summary, compound Pearson residuals provide an efficient and effective way to normalize read count data based on simple mechanistic assumptions.

Transcriptional kinetics and molecular functions of long noncoding RNAs.

An increasing number of long noncoding RNAs (lncRNAs) have experimentally confirmed functions, yet little is known about their transcriptional dynamics and it is challenging to determine their regulatory effects. Here, we used allele-sensitive single-cell RNA sequencing to demonstrate that, compared to messenger RNAs, lncRNAs have twice as long duration between two transcriptional bursts. Additionally, we observed increased cell-to-cell variability in lncRNA expression due to lower frequency bursting producing larger numbers of RNA molecules. Exploiting heterogeneity in asynchronously growing cells, we identified and experimentally validated lncRNAs with cell state-specific functions involved in cell cycle progression and apoptosis. Finally, we identified cis-functioning lncRNAs and showed that knockdown of these lncRNAs modulated the nearby protein-coding gene's transcriptional burst frequency or size. In summary, we identified distinct transcriptional regulation of lncRNAs and demonstrated a role for lncRNAs in the regulation of mRNA transcriptional bursting.

A DNA-nanoassembly-based approach to map membrane protein nanoenvironments.

Most proteins at the plasma membrane are not uniformly distributed but localize to dynamic domains of nanoscale dimensions. To investigate their functional relevance, there is a need for methods that enable comprehensive analysis of the compositions and spatial organizations of membrane protein nanodomains in cell populations. Here we describe the development of a non-microscopy-based method for ensemble analysis of membrane protein nanodomains. The method, termed nanoscale deciphering of membrane protein nanodomains (NanoDeep), is based on the use of DNA nanoassemblies to translate membrane protein organization information into a DNA sequencing readout. Using NanoDeep, we characterized the nanoenvironments of Her2, a membrane receptor of critical relevance in cancer. Importantly, we were able to modulate by design the inventory of proteins analysed by NanoDeep. NanoDeep has the potential to provide new insights into the roles of the composition and spatial organization of protein nanoenvironments in the regulation of membrane protein function.

Cell type composition and circuit organization of clonally related excitatory neurons in the juvenile mouse neocortex.

Clones of excitatory neurons derived from a common progenitor have been proposed to serve as elementary information processing modules in the neocortex. To characterize the cell types and circuit diagram of clonally related excitatory neurons, we performed multi-cell patch clamp recordings and Patch-seq on neurons derived from Nestin-positive progenitors labeled by tamoxifen induction at embryonic day 10.5. The resulting clones are derived from two radial glia on average, span cortical layers 2-6, and are composed of a random sampling of transcriptomic cell types. We find an interaction between shared lineage and connection type: related neurons are more likely to be connected vertically across cortical layers, but not laterally within the same layer. These findings challenge the view that related neurons show uniformly increased connectivity and suggest that integration of vertical intra-clonal input with lateral inter-clonal input may represent a developmentally programmed connectivity motif supporting the emergence of functional circuits.

Layer 4 of mouse neocortex differs in cell types and circuit organization between sensory areas.

Layer 4 (L4) of mammalian neocortex plays a crucial role in cortical information processing, yet a complete census of its cell types and connectivity remains elusive. Using whole-cell recordings with morphological recovery, we identified one major excitatory and seven inhibitory types of neurons in L4 of adult mouse visual cortex (V1). Nearly all excitatory neurons were pyramidal and all somatostatin-positive (SOM+) non-fast-spiking interneurons were Martinotti cells. In contrast, in somatosensory cortex (S1), excitatory neurons were mostly stellate and SOM+ interneurons were non-Martinotti. These morphologically distinct SOM+ interneurons corresponded to different transcriptomic cell types and were differentially integrated into the local circuit with only S1 neurons receiving local excitatory input. We propose that cell type specific circuit motifs, such as the Martinotti/pyramidal and non-Martinotti/stellate pairs, are used across the cortex as building blocks to assemble cortical circuits.

Author Correction: Layer 4 of mouse neocortex differs in cell types and circuit organization between sensory areas.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Microscale Bioreactors for in situ characterization of GI epithelial cell physiology.

The development of in vitro artificial small intestines that realistically mimic in vivo systems will enable vast improvement of our understanding of the human gut and its impact on human health. Synthetic in vitro models can control specific parameters, including (but not limited to) cell types, fluid flow, nutrient profiles and gaseous exchange. They are also "open" systems, enabling access to chemical and physiological information. In this work, we demonstrate the importance of gut surface topography and fluid flow dynamics which are shown to impact epithelial cell growth, proliferation and intestinal cell function. We have constructed a small intestinal bioreactor using 3-D printing and polymeric scaffolds that mimic the 3-D topography of the intestine and its fluid flow. Our results indicate that TEER measurements, which are typically high in static 2-D Transwell apparatuses, is lower in the presence of liquid sheer and 3-D topography compared to a flat scaffold and static conditions. There was also increased cell proliferation and discovered localized regions of elevated apoptosis, specifically at the tips of the villi, where there is highest sheer. Similarly, glucose was actively transported (as opposed to passive) and at higher rates under flow.