Volatile working memory representations crystallize with practice.

Working memory, the process through which information is transiently maintained and manipulated over a brief period, is essential for most cognitive functions1-4. However, the mechanisms underlying the generation and evolution of working-memory neuronal representations at the population level over long timescales remain unclear. Here, to identify these mechanisms, we trained head-fixed mice to perform an olfactory delayed-association task in which the mice made decisions depending on the sequential identity of two odours separated by a 5 s delay. Optogenetic inhibition of secondary motor neurons during the late-delay and choice epochs strongly impaired the task performance of the mice. Mesoscopic calcium imaging of large neuronal populations of the secondary motor cortex (M2), retrosplenial cortex (RSA) and primary motor cortex (M1) showed that many late-delay-epoch-selective neurons emerged in M2 as the mice learned the task. Working-memory late-delay decoding accuracy substantially improved in the M2, but not in the M1 or RSA, as the mice became experts. During the early expert phase, working-memory representations during the late-delay epoch drifted across days, while the stimulus and choice representations stabilized. In contrast to single-plane layer 2/3 (L2/3) imaging, simultaneous volumetric calcium imaging of up to 73,307 M2 neurons, which included superficial L5 neurons, also revealed stabilization of late-delay working-memory representations with continued practice. Thus, delay- and choice-related activities that are essential for working-memory performance drift during learning and stabilize only after several days of expert performance.

A method to find temporal structure of neuronal coactivity patterns with across-trial correlations.

BACKGROUND: The across-trial correlation of neurons' coactivity patterns emerges to be important for information coding, but methods for finding their temporal structures remain largely unexplored. NEW METHOD: In the present study, we propose a method to find time clusters in which coactivity patterns of neurons are correlated across trials. We transform the multidimensional neural activity at each timing into a coactivity pattern of binary states, and predict the coactivity patterns at different timings. We devise a method suitable for these coactivity pattern predictions, call general event prediction. Cross-temporal prediction accuracy is then used to estimate across-trial correlations between coactivity patterns at two timings. We extract time clusters from the cross-temporal prediction accuracy by a modified k-means algorithm. RESULTS: The feasibility of the proposed method is verified through simulations based on ground truth. We apply the proposed method to a calcium imaging dataset recorded from the motor cortex of mice, and demonstrate time clusters of motor cortical coactivity patterns during a motor task. COMPARISON WITH EXISTING METHODS: While the existing cosine similarity method, which does not account for across-trial correlation, shows temporal structures only for contralateral neural responses, the proposed method reveals those for both contralateral and ipsilateral neural responses, demonstrating the effect of across-trial correlations. CONCLUSIONS: This study introduces a novel method for measuring the temporal structure of neuronal ensemble activity.

Selective plasticity of layer 2/3 inputs onto distal forelimb controlling layer 5 corticospinal neurons with skilled grasp motor training.

Layer 5 neurons of the neocortex receive their principal inputs from layer 2/3 neurons. We seek to identify the nature and extent of the plasticity of these projections with motor learning. Using optogenetic and viral intersectional tools to selectively stimulate distinct neuronal subsets in rat primary motor cortex, we simultaneously record from pairs of corticospinal neurons associated with distinct features of motor output control: distal forelimb vs. proximal forelimb. Activation of Channelrhodopsin2-expressing layer 2/3 afferents onto layer 5 in untrained animals produces greater monosynaptic excitation of neurons controlling the proximal forelimb. Following skilled grasp training, layer 2/3 inputs onto corticospinal neurons controlling the distal forelimb associated with skilled grasping become significantly stronger. Moreover, peak excitatory response amplitude nearly doubles while latency shortens, and excitatory-to-inhibitory latencies become significantly prolonged. These findings demonstrate distinct, highly segregated, and cell-specific plasticity of layer 2/3 projections during skilled grasp motor learning.

Planar amorphous silicon carbide microelectrode arrays for chronic recording in rat motor cortex.

Chronic implantation of intracortical microelectrode arrays (MEAs) capable of recording from individual neurons can be used for the development of brain-machine interfaces. However, these devices show reduced recording capabilities under chronic conditions due, at least in part, to the brain's foreign body response (FBR). This creates a need for MEAs that can minimize the FBR to possibly enable long-term recording. A potential approach to reduce the FBR is the use of MEAs with reduced cross-sectional geometries. Here, we fabricated 4-shank amorphous silicon carbide (a-SiC) MEAs and implanted them into the motor cortex of seven female Sprague-Dawley rats. Each a-SiC MEA shank was 8 µm thick by 20 µm wide and had sixteen sputtered iridium oxide film (SIROF) electrodes (4 per shank). A-SiC was chosen as the fabrication base for its high chemical stability, good electrical insulation properties, and amenability to thin film fabrication. Electrochemical analysis and neural recordings were performed weekly for 4 months. MEAs were characterized pre-implantation in buffered saline and in vivo using electrochemical impedance spectroscopy and cyclic voltammetry at 50 mV/s and 50,000 mV/s. Neural recordings were analyzed for single unit activity. At the end of the study, animals were sacrificed for immunohistochemical analysis. We observed statistically significant, but small, increases in 1 and 30 kHz impedance values and 50,000 mV/s charge storage capacity over the 16-week implantation period. Slow sweep 50 mV/s CV and 1 Hz impedance did not significantly change over time. Impedance values increased from 11.6 MΩ to 13.5 MΩ at 1 Hz, 1.2 MΩ-2.9 MΩ at 1 kHz, and 0.11 MΩ-0.13 MΩ at 30 kHz over 16 weeks. The median charge storage capacity of the implanted electrodes at 50 mV/s was 58.1 mC/cm2 on week 1 and 55.9 mC/cm2 on week 16, and at 50,000 mV/s, 4.27 mC/cm2 on week 1 and 5.93 mC/cm2 on week 16. Devices were able to record neural activity from 92% of all active channels at the beginning of the study, At the study endpoint, a-SiC devices were still recording single-unit activity on 51% of electrochemically active electrode channels. In addition, we observed that the signal-to-noise ratio experienced a small decline of -0.19 per week. We also classified observed units as fast and slow repolarizing based on the trough-to-peak time. Although the overall presence of single units declined, fast and slow repolarizing units declined at a similar rate. At recording electrode depth, immunohistochemistry showed minimal tissue response to the a-SiC devices, as indicated by statistically insignificant differences in activated glial cell response between implanted brains slices and contralateral sham slices at 150 µm away from the implant location, as evidenced by GFAP staining. NeuN staining revealed the presence of neuronal cell bodies close to the implantation site, again statistically not different from a contralateral sham slice. These results warrant further investigation of a-SiC MEAs for future long-term implantation neural recording studies.

Vocal learning-associated convergent evolution in mammalian proteins and regulatory elements.

Vocal production learning ("vocal learning") is a convergently evolved trait in vertebrates. To identify brain genomic elements associated with mammalian vocal learning, we integrated genomic, anatomical, and neurophysiological data from the Egyptian fruit bat (Rousettus aegyptiacus) with analyses of the genomes of 215 placental mammals. First, we identified a set of proteins evolving more slowly in vocal learners. Then, we discovered a vocal motor cortical region in the Egyptian fruit bat, an emergent vocal learner, and leveraged that knowledge to identify active cis-regulatory elements in the motor cortex of vocal learners. Machine learning methods applied to motor cortex open chromatin revealed 50 enhancers robustly associated with vocal learning whose activity tended to be lower in vocal learners. Our research implicates convergent losses of motor cortex regulatory elements in mammalian vocal learning evolution.

Conserved and divergent gene regulatory programs of the mammalian neocortex.

Divergence of cis-regulatory elements drives species-specific traits1, but how this manifests in the evolution of the neocortex at the molecular and cellular level remains unclear. Here we investigated the gene regulatory programs in the primary motor cortex of human, macaque, marmoset and mouse using single-cell multiomics assays, generating gene expression, chromatin accessibility, DNA methylome and chromosomal conformation profiles from a total of over 200,000 cells. From these data, we show evidence that divergence of transcription factor expression corresponds to species-specific epigenome landscapes. We find that conserved and divergent gene regulatory features are reflected in the evolution of the three-dimensional genome. Transposable elements contribute to nearly 80% of the human-specific candidate cis-regulatory elements in cortical cells. Through machine learning, we develop sequence-based predictors of candidate cis-regulatory elements in different species and demonstrate that the genomic regulatory syntax is highly preserved from rodents to primates. Finally, we show that epigenetic conservation combined with sequence similarity helps to uncover functional cis-regulatory elements and enhances our ability to interpret genetic variants contributing to neurological disease and traits.

Preserved neural dynamics across animals performing similar behaviour.

Animals of the same species exhibit similar behaviours that are advantageously adapted to their body and environment. These behaviours are shaped at the species level by selection pressures over evolutionary timescales. Yet, it remains unclear how these common behavioural adaptations emerge from the idiosyncratic neural circuitry of each individual. The overall organization of neural circuits is preserved across individuals1 because of their common evolutionarily specified developmental programme2-4. Such organization at the circuit level may constrain neural activity5-8, leading to low-dimensional latent dynamics across the neural population9-11. Accordingly, here we suggested that the shared circuit-level constraints within a species would lead to suitably preserved latent dynamics across individuals. We analysed recordings of neural populations from monkey and mouse motor cortex to demonstrate that neural dynamics in individuals from the same species are surprisingly preserved when they perform similar behaviour. Neural population dynamics were also preserved when animals consciously planned future movements without overt behaviour12 and enabled the decoding of planned and ongoing movement across different individuals. Furthermore, we found that preserved neural dynamics extend beyond cortical regions to the dorsal striatum, an evolutionarily older structure13,14. Finally, we used neural network models to demonstrate that behavioural similarity is necessary but not sufficient for this preservation. We posit that these emergent dynamics result from evolutionary constraints on brain development and thus reflect fundamental properties of the neural basis of behaviour.

Thalamus drives vocal onsets in the zebra finch courtship song.

While motor cortical circuits contain information related to specific movement parameters1, long-range inputs also have a critical role in action execution2,3. Thalamic projections can shape premotor activity2-6 and have been suggested7 to mediate the selection of short, stereotyped actions comprising more complex behaviours8. However, the mechanisms by which thalamus interacts with motor cortical circuits to execute such movement sequences remain unknown. Here we find that thalamic drive engages a specific subpopulation of premotor neurons within the zebra finch song nucleus HVC (proper name) and that these inputs are critical for the progression between vocal motor elements (that is, 'syllables'). In vivo two-photon imaging of thalamic axons in HVC showed robust song-related activity, and online perturbations of thalamic function caused song to be truncated at syllable boundaries. We used thalamic stimulation to identify a sparse set of thalamically driven neurons within HVC, representing ~15% of the premotor neurons within that network. Unexpectedly, this population of putative thalamorecipient neurons is robustly active immediately preceding syllable onset, leading to the possibility that thalamic input can initiate individual song components through selectively targeting these 'starter cells'. Our findings highlight the motor thalamus as a director of cortical dynamics in the context of an ethologically relevant behavioural sequence.

Layer-specific pain relief pathways originating from primary motor cortex.

The primary motor cortex (M1) is involved in the control of voluntary movements and is extensively mapped in this capacity. Although the M1 is implicated in modulation of pain, the underlying circuitry and causal underpinnings remain elusive. We unexpectedly unraveled a connection from the M1 to the nucleus accumbens reward circuitry through a M1 layer 6-mediodorsal thalamus pathway, which specifically suppresses negative emotional valence and associated coping behaviors in neuropathic pain. By contrast, layer 5 M1 neurons connect with specific cell populations in zona incerta and periaqueductal gray to suppress sensory hypersensitivity without altering pain affect. Thus, the M1 employs distinct, layer-specific pathways to attune sensory and aversive-emotional components of neuropathic pain, which can be exploited for purposes of pain relief.

Neuronal population dynamics during motor plan cancellation in nonhuman primates.

To understand the cortical neuronal dynamics behind movement generation and control, most studies have focused on tasks where actions were planned and then executed using different instances of visuomotor transformations. However, to fully understand the dynamics related to movement control, one must also study how movements are actively inhibited. Inhibition, indeed, represents the first level of control both when different alternatives are available and only one solution could be adopted and when it is necessary to maintain the current position. We recorded neuronal activity from a multielectrode array in the dorsal premotor cortex (PMd) of monkeys performing a countermanding reaching task that requires, in a subset of trials, them to cancel a planned movement before its onset. In the analysis of the neuronal state space of PMd, we found a subspace in which activities conveying temporal information were confined during active inhibition and position holding. Movement execution required activities to escape from this subspace toward an orthogonal subspace and, furthermore, surpass a threshold associated with the maturation of the motor plan. These results revealed further details in the neuronal dynamics underlying movement control, extending the hypothesis that neuronal computation confined in an "output-null" subspace does not produce movements.

Cell-type-specific responses to associative learning in the primary motor cortex.

The primary motor cortex (M1) is known to be a critical site for movement initiation and motor learning. Surprisingly, it has also been shown to possess reward-related activity, presumably to facilitate reward-based learning of new movements. However, whether reward-related signals are represented among different cell types in M1, and whether their response properties change after cue-reward conditioning remains unclear. Here, we performed longitudinal in vivo two-photon Ca2+ imaging to monitor the activity of different neuronal cell types in M1 while mice engaged in a classical conditioning task. Our results demonstrate that most of the major neuronal cell types in M1 showed robust but differential responses to both the conditioned cue stimulus (CS) and reward, and their response properties undergo cell-type-specific modifications after associative learning. PV-INs' responses became more reliable to the CS, while VIP-INs' responses became more reliable to reward. Pyramidal neurons only showed robust responses to novel reward, and they habituated to it after associative learning. Lastly, SOM-INs' responses emerged and became more reliable to both the CS and reward after conditioning. These observations suggest that cue- and reward-related signals are preferentially represented among different neuronal cell types in M1, and the distinct modifications they undergo during associative learning could be essential in triggering different aspects of local circuit reorganization in M1 during reward-based motor skill learning.

Global organization of neuronal activity only requires unstructured local connectivity.

Modern electrophysiological recordings simultaneously capture single-unit spiking activities of hundreds of neurons spread across large cortical distances. Yet, this parallel activity is often confined to relatively low-dimensional manifolds. This implies strong coordination also among neurons that are most likely not even connected. Here, we combine in vivo recordings with network models and theory to characterize the nature of mesoscopic coordination patterns in macaque motor cortex and to expose their origin: We find that heterogeneity in local connectivity supports network states with complex long-range cooperation between neurons that arises from multi-synaptic, short-range connections. Our theory explains the experimentally observed spatial organization of covariances in resting state recordings as well as the behaviorally related modulation of covariance patterns during a reach-to-grasp task. The ubiquity of heterogeneity in local cortical circuits suggests that the brain uses the described mechanism to flexibly adapt neuronal coordination to momentary demands.

Temporal controls over inter-areal cortical projection neuron fate diversity.

Interconnectivity between neocortical areas is critical for sensory integration and sensorimotor transformations1-6. These functions are mediated by heterogeneous inter-areal cortical projection neurons (ICPN), which send axon branches across cortical areas as well as to subcortical targets7-9. Although ICPN are anatomically diverse10-14, they are molecularly homogeneous15, and how the diversity of their anatomical and functional features emerge during development remains largely unknown. Here we address this question by linking the connectome and transcriptome in developing single ICPN of the mouse neocortex using a combination of multiplexed analysis of projections by sequencing16,17 (MAPseq, to identify single-neuron axonal projections) and single-cell RNA sequencing (to identify corresponding gene expression). Focusing on neurons of the primary somatosensory cortex (S1), we reveal a protracted unfolding of the molecular and functional differentiation of motor cortex-projecting ([Formula: see text]) ICPN compared with secondary somatosensory cortex-projecting ([Formula: see text]) ICPN. We identify SOX11 as a temporally differentially expressed transcription factor in [Formula: see text] versus [Formula: see text] ICPN. Postnatal manipulation of SOX11 expression in S1 impaired sensorimotor connectivity and disrupted selective exploratory behaviours in mice. Together, our results reveal that within a single cortical area, different subtypes of ICPN have distinct postnatal paces of molecular differentiation, which are subsequently reflected in distinct circuit connectivities and functions. Dynamic differences in the expression levels of a largely generic set of genes, rather than fundamental differences in the identity of developmental genetic programs, may thus account for the emergence of intra-type diversity in cortical neurons.

Resurgent Na+ currents promote ultrafast spiking in projection neurons that drive fine motor control.

The underlying mechanisms that promote precise spiking in upper motor neurons controlling fine motor skills are not well understood. Here we report that projection neurons in the adult zebra finch song nucleus RA display robust high-frequency firing, ultra-narrow spike waveforms, superfast Na+ current inactivation kinetics, and large resurgent Na+ currents (INaR). These properties of songbird pallial motor neurons closely resemble those of specialized large pyramidal neurons in mammalian primary motor cortex. They emerge during the early phases of song development in males, but not females, coinciding with a complete switch of Na+ channel subunit expression from Navß3 to Navß4. Dynamic clamping and dialysis of Navß4's C-terminal peptide into juvenile RA neurons provide evidence that Navß4, and its associated INaR, promote neuronal excitability. We thus propose that INaR modulates the excitability of upper motor neurons that are required for the execution of fine motor skills.

Spatially resolved cell atlas of the mouse primary motor cortex by MERFISH.

A mammalian brain is composed of numerous cell types organized in an intricate manner to form functional neural circuits. Single-cell RNA sequencing allows systematic identification of cell types based on their gene expression profiles and has revealed many distinct cell populations in the brain1,2. Single-cell epigenomic profiling3,4 further provides information on gene-regulatory signatures of different cell types. Understanding how different cell types contribute to brain function, however, requires knowledge of their spatial organization and connectivity, which is not preserved in sequencing-based methods that involve cell dissociation. Here we used a single-cell transcriptome-imaging method, multiplexed error-robust fluorescence in situ hybridization (MERFISH)5, to generate a molecularly defined and spatially resolved cell atlas of the mouse primary motor cortex. We profiled approximately 300,000 cells in the mouse primary motor cortex and its adjacent areas, identified 95 neuronal and non-neuronal cell clusters, and revealed a complex spatial map in which not only excitatory but also most inhibitory neuronal clusters adopted laminar organizations. Intratelencephalic neurons formed a largely continuous gradient along the cortical depth axis, in which the gene expression of individual cells correlated with their cortical depths. Furthermore, we integrated MERFISH with retrograde labelling to probe projection targets of neurons of the mouse primary motor cortex and found that their cortical projections formed a complex network in which individual neuronal clusters project to multiple target regions and individual target regions receive inputs from multiple neuronal clusters.

Comparative cellular analysis of motor cortex in human, marmoset and mouse.

The primary motor cortex (M1) is essential for voluntary fine-motor control and is functionally conserved across mammals1. Here, using high-throughput transcriptomic and epigenomic profiling of more than 450,000 single nuclei in humans, marmoset monkeys and mice, we demonstrate a broadly conserved cellular makeup of this region, with similarities that mirror evolutionary distance and are consistent between the transcriptome and epigenome. The core conserved molecular identities of neuronal and non-neuronal cell types allow us to generate a cross-species consensus classification of cell types, and to infer conserved properties of cell types across species. Despite the overall conservation, however, many species-dependent specializations are apparent, including differences in cell-type proportions, gene expression, DNA methylation and chromatin state. Few cell-type marker genes are conserved across species, revealing a short list of candidate genes and regulatory mechanisms that are responsible for conserved features of homologous cell types, such as the GABAergic chandelier cells. This consensus transcriptomic classification allows us to use patch-seq (a combination of whole-cell patch-clamp recordings, RNA sequencing and morphological characterization) to identify corticospinal Betz cells from layer 5 in non-human primates and humans, and to characterize their highly specialized physiology and anatomy. These findings highlight the robust molecular underpinnings of cell-type diversity in M1 across mammals, and point to the genes and regulatory pathways responsible for the functional identity of cell types and their species-specific adaptations.

A transcriptomic and epigenomic cell atlas of the mouse primary motor cortex.

Single-cell transcriptomics can provide quantitative molecular signatures for large, unbiased samples of the diverse cell types in the brain1-3. With the proliferation of multi-omics datasets, a major challenge is to validate and integrate results into a biological understanding of cell-type organization. Here we generated transcriptomes and epigenomes from more than 500,000 individual cells in the mouse primary motor cortex, a structure that has an evolutionarily conserved role in locomotion. We developed computational and statistical methods to integrate multimodal data and quantitatively validate cell-type reproducibility. The resulting reference atlas-containing over 56 neuronal cell types that are highly replicable across analysis methods, sequencing technologies and modalities-is a comprehensive molecular and genomic account of the diverse neuronal and non-neuronal cell types in the mouse primary motor cortex. The atlas includes a population of excitatory neurons that resemble pyramidal cells in layer 4 in other cortical regions4. We further discovered thousands of concordant marker genes and gene regulatory elements for these cell types. Our results highlight the complex molecular regulation of cell types in the brain and will directly enable the design of reagents to target specific cell types in the mouse primary motor cortex for functional analysis.

Isoform cell-type specificity in the mouse primary motor cortex.

Full-length SMART-seq1 single-cell RNA sequencing can be used to measure gene expression at isoform resolution, making possible the identification of specific isoform markers for different cell types. Used in conjunction with spatial RNA capture and gene-tagging methods, this enables the inference of spatially resolved isoform expression for different cell types. Here, in a comprehensive analysis of 6,160 mouse primary motor cortex cells assayed with SMART-seq, 280,327 cells assayed with MERFISH2 and 94,162 cells assayed with 10x Genomics sequencing3, we find examples of isoform specificity in cell types-including isoform shifts between cell types that are masked in gene-level analysis-as well as examples of transcriptional regulation. Additionally, we show that isoform specificity helps to refine cell types, and that a multi-platform analysis of single-cell transcriptomic data leveraging multiple measurements provides a comprehensive atlas of transcription in the mouse primary motor cortex that improves on the possibilities offered by any single technology.

Cellular anatomy of the mouse primary motor cortex.

An essential step toward understanding brain function is to establish a structural framework with cellular resolution on which multi-scale datasets spanning molecules, cells, circuits and systems can be integrated and interpreted1. Here, as part of the collaborative Brain Initiative Cell Census Network (BICCN), we derive a comprehensive cell type-based anatomical description of one exemplar brain structure, the mouse primary motor cortex, upper limb area (MOp-ul). Using genetic and viral labelling, barcoded anatomy resolved by sequencing, single-neuron reconstruction, whole-brain imaging and cloud-based neuroinformatics tools, we delineated the MOp-ul in 3D and refined its sublaminar organization. We defined around two dozen projection neuron types in the MOp-ul and derived an input-output wiring diagram, which will facilitate future analyses of motor control circuitry across molecular, cellular and system levels. This work provides a roadmap towards a comprehensive cellular-resolution description of mammalian brain architecture.

A multimodal cell census and atlas of the mammalian primary motor cortex.

Here we report the generation of a multimodal cell census and atlas of the mammalian primary motor cortex as the initial product of the BRAIN Initiative Cell Census Network (BICCN). This was achieved by coordinated large-scale analyses of single-cell transcriptomes, chromatin accessibility, DNA methylomes, spatially resolved single-cell transcriptomes, morphological and electrophysiological properties and cellular resolution input-output mapping, integrated through cross-modal computational analysis. Our results advance the collective knowledge and understanding of brain cell-type organization1-5. First, our study reveals a unified molecular genetic landscape of cortical cell types that integrates their transcriptome, open chromatin and DNA methylation maps. Second, cross-species analysis achieves a consensus taxonomy of transcriptomic types and their hierarchical organization that is conserved from mouse to marmoset and human. Third, in situ single-cell transcriptomics provides a spatially resolved cell-type atlas of the motor cortex. Fourth, cross-modal analysis provides compelling evidence for the transcriptomic, epigenomic and gene regulatory basis of neuronal phenotypes such as their physiological and anatomical properties, demonstrating the biological validity and genomic underpinning of neuron types. We further present an extensive genetic toolset for targeting glutamatergic neuron types towards linking their molecular and developmental identity to their circuit function. Together, our results establish a unifying and mechanistic framework of neuronal cell-type organization that integrates multi-layered molecular genetic and spatial information with multi-faceted phenotypic properties.