To provide a temporal framework for the genoarchitecture of brain development, we generated in situ hybridization data for embryonic and postnatal mouse brain at seven developmental stages for ∼2,100 genes, which were processed with an automated informatics pipeline and manually annotated. This resource comprises 434,946 images, seven reference atlases, an ontogenetic ontology, and tools to explore coexpression of genes across neurodevelopment. Gene sets coinciding with developmental phenomena were identified. A temporal shift in the principles governing the molecular organization of the brain was detected, with transient neuromeric, plate-based organization of the brain present at E11.5 and E13.5. Finally, these data provided a transcription factor code that discriminates brain structures and identifies the developmental age of a tissue, providing a foundation for eventual genetic manipulation or tracking of specific brain structures over development. The resource is available as the Allen Developing Mouse Brain Atlas (http://developingmouse.brain-map.org).

Realistic computational models of single neurons require component ion channels that reproduce experimental findings. Here, a topology-mutating genetic algorithm that searches for the best state diagram and transition-rate parameters to model macroscopic ion-channel behavior is described. Important features of the algorithm include a topology-altering strategy, automatic satisfaction of equilibrium constraints (microscopic reversibility), and multiple-protocol fitting using sequential goal programming rather than explicit weighting. Application of this genetic algorithm to design a sodium-channel model exhibiting both fast and prolonged inactivation yields a six-state model that produces realistic activity-dependent attenuation of action-potential backpropagation in current-clamp simulations of a CA1 pyramidal neuron.

Nervous systems are composed of various cell types, but the extent of cell type diversity is poorly understood. We constructed a cellular taxonomy of one cortical region, primary visual cortex, in adult mice on the basis of single-cell RNA sequencing. We identified 49 transcriptomic cell types, including 23 GABAergic, 19 glutamatergic and 7 non-neuronal types. We also analyzed cell type-specific mRNA processing and characterized genetic access to these transcriptomic types by many transgenic Cre lines. Finally, we found that some of our transcriptomic cell types displayed specific and differential electrophysiological and axon projection properties, thereby confirming that the single-cell transcriptomic signatures can be associated with specific cellular properties.

Neuronal computation involves the integration of synaptic inputs that are often distributed over expansive dendritic trees, suggesting the need for compensatory mechanisms that enable spatially disparate synapses to influence neuronal output. In hippocampal CA1 pyramidal neurons, such mechanisms have indeed been reported, which normalize either the ability of distributed synapses to drive action potential initiation in the axon or their ability to drive dendritic spiking locally. Here we report that these mechanisms can coexist, through an elegant combination of distance-dependent regulation of synapse number and synaptic expression of AMPA and NMDA receptors. Together, these complementary gradients allow individual dendrites in both the apical and basal dendritic trees of hippocampal neurons to operate as facile computational subunits capable of supporting both global integration in the soma/axon and local integration in the dendrite.

The structure and function of the human brain are highly stereotyped, implying a conserved molecular program responsible for its development, cellular structure and function. We applied a correlation-based metric called differential stability to assess reproducibility of gene expression patterning across 132 structures in six individual brains, revealing mesoscale genetic organization. The genes with the highest differential stability are highly biologically relevant, with enrichment for brain-related annotations, disease associations, drug targets and literature citations. Using genes with high differential stability, we identified 32 anatomically diverse and reproducible gene expression signatures, which represent distinct cell types, intracellular components and/or associations with neurodevelopmental and neurodegenerative disorders. Genes in neuron-associated compared to non-neuronal networks showed higher preservation between human and mouse; however, many diversely patterned genes displayed marked shifts in regulation between species. Finally, highly consistent transcriptional architecture in neocortex is correlated with resting state functional connectivity, suggesting a link between conserved gene expression and functionally relevant circuitry.

Advances in single-cell RNA-sequencing technology have resulted in a wealth of studies aiming to identify transcriptomic cell types in various biological systems. There are multiple experimental approaches to isolate and profile single cells, which provide different levels of cellular and tissue coverage. In addition, multiple computational strategies have been proposed to identify putative cell types from single-cell data. From a data generation perspective, recent single-cell studies can be classified into two groups: those that distribute reads shallowly over large numbers of cells and those that distribute reads more deeply over a smaller cell population. Although there are advantages to both approaches in terms of cellular and tissue coverage, it is unclear whether different computational cell type identification methods are better suited to one or the other experimental paradigm. This study reviews three cell type clustering algorithms, each representing one of three broad approaches, and finds that PCA-based algorithms appear most suited to low read depth data sets, whereas gene clustering-based and biclustering algorithms perform better on high read depth data sets. In addition, highly related cell classes are better distinguished by higher-depth data, given the same total number of reads; however, simultaneous discovery of distinct and similar types is better served by lower-depth, higher cell number data. Overall, this study suggests that the depth of profiling should be determined by initial assumptions about the diversity of cells in the population, and that the selection of clustering algorithm(s) is subsequently based on the depth of profiling will allow for better identification of putative transcriptomic cell types.

The brain is an organ of immense complexity. Next-generation RNA sequencing (RNA-seq) is becoming increasingly popular in the deconstruction of this complexity into distinct classes of 'cell types'. Notably, in addition to revealing the organization of this distinct cell-type landscape, the technology has also begun to illustrate that continuous variation can be found within narrowly defined cell types. Here we summarize the evidence for graded transcriptomic heterogeneity being present, widespread, and functionally relevant in the nervous system. We explain how these graded differences can map onto higher-order organizational features and how they may reframe existing interpretations of higher-order heterogeneity. Ultimately, a multimodal approach incorporating continuously variable cell types will facilitate an accurate reductionist interpretation of the nervous system.

The neocortex contains diverse populations of excitatory neurons segregated by layer and further definable by their specific cortical and subcortical projection targets. The current study describes a systematic approach to identify molecular correlates of specific projection neuron classes in mouse primary somatosensory cortex (S1), using a combination of in situ hybridization (ISH) data mining, marker gene colocalization, and combined retrograde labeling with ISH for layer-specific marker genes. First, we identified a large set of genes with specificity for each cortical layer, and that display heterogeneous patterns within those layers. Using these genes as markers, we find extensive evidence for the covariation of gene expression and projection target specificity in layer 2/3, 5, and 6, with individual genes labeling neurons projecting to specific subsets of target structures. The combination of gene expression and target specificity imply a great diversity of projection neuron classes that is similar to or greater than that of GABAergic interneurons. The covariance of these 2 phenotypic modalities suggests that these classes are both discrete and genetically specified.

The microtubule cytoskeleton has proven to be an effective target for cancer therapeutics. One class of drugs, known as microtubule stabilizing agents (MSAs), binds to microtubule polymers and stabilizes them against depolymerization. The prototype of this group of drugs, Taxol, is an effective chemotherapeutic agent used extensively in the treatment of human ovarian, breast, and lung carcinomas. Although electron crystallography and photoaffinity labeling experiments determined that the binding site for Taxol is in a hydrophobic pocket in beta-tubulin, little was known about the effects of this drug on the conformation of the entire microtubule. A recent study from our laboratory utilizing hydrogen-deuterium exchange (HDX) in concert with various mass spectrometry (MS) techniques has provided new information on the structure of microtubules upon Taxol binding. In the current study we apply this technique to determine the binding mode and the conformational effects on chicken erythrocyte tubulin (CET) of another MSA, discodermolide, whose synthetic analogues may have potential use in the clinic. We confirmed that, like Taxol, discodermolide binds to the taxane binding pocket in beta-tubulin. However, as opposed to Taxol, which has major interactions with the M-loop, discodermolide orients itself away from this loop and toward the N-terminal H1-S2 loop. Additionally, discodermolide stabilizes microtubules mainly via its effects on interdimer contacts, specifically on the alpha-tubulin side, and to a lesser extent on interprotofilament contacts between adjacent beta-tubulin subunits. Also, our results indicate complementary stabilizing effects of Taxol and discodermolide on the microtubules, which may explain the synergy observed between the two drugs in vivo.

Clusters of time series data may change location and memberships over time; in gene expression data, this occurs as groups of genes or samples respond differently to stimuli or experimental conditions at different times. In order to uncover this underlying temporal structure, we consider dynamic clusters with time-dependent parameters which split and merge over time, enabling cluster memberships to change. These interesting time-dependent structures are useful in understanding the development of organisms or complex organs, and could not be identified using traditional clustering methods. In cell cycle data, these time-dependent structure may provide links between genes and stages of the cell cycle, whilst in developmental data sets they may highlight key developmental transitions.