Genetic dissection of the glutamatergic neuron system in cerebral cortex.

Diverse types of glutamatergic pyramidal neurons mediate the myriad processing streams and output channels of the cerebral cortex1,2, yet all derive from neural progenitors of the embryonic dorsal telencephalon3,4. Here we establish genetic strategies and tools for dissecting and fate-mapping subpopulations of pyramidal neurons on the basis of their developmental and molecular programs. We leverage key transcription factors and effector genes to systematically target temporal patterning programs in progenitors and differentiation programs in postmitotic neurons. We generated over a dozen temporally inducible mouse Cre and Flp knock-in driver lines to enable the combinatorial targeting of major progenitor types and projection classes. Combinatorial strategies confer viral access to subsets of pyramidal neurons defined by developmental origin, marker expression, anatomical location and projection targets. These strategies establish an experimental framework for understanding the hierarchical organization and developmental trajectory of subpopulations of pyramidal neurons that assemble cortical processing networks and output channels.

A guide to the BRAIN Initiative Cell Census Network data ecosystem.

Characterizing cellular diversity at different levels of biological organization and across data modalities is a prerequisite to understanding the function of cell types in the brain. Classification of neurons is also essential to manipulate cell types in controlled ways and to understand their variation and vulnerability in brain disorders. The BRAIN Initiative Cell Census Network (BICCN) is an integrated network of data-generating centers, data archives, and data standards developers, with the goal of systematic multimodal brain cell type profiling and characterization. Emphasis of the BICCN is on the whole mouse brain with demonstration of prototype feasibility for human and nonhuman primate (NHP) brains. Here, we provide a guide to the cellular and spatial approaches employed by the BICCN, and to accessing and using these data and extensive resources, including the BRAIN Cell Data Center (BCDC), which serves to manage and integrate data across the ecosystem. We illustrate the power of the BICCN data ecosystem through vignettes highlighting several BICCN analysis and visualization tools. Finally, we present emerging standards that have been developed or adopted toward Findable, Accessible, Interoperable, and Reusable (FAIR) neuroscience. The combined BICCN ecosystem provides a comprehensive resource for the exploration and analysis of cell types in the brain.

A Standardized Protocol for the Safe Retrieval of Infectious Postmortem Human Brain for Studying Whole-Brain Pathology.

ABSTRACT: We describe a safe and standardized perfusion protocol for studying brain pathology in high-risk autopsies using a custom-designed low-cost infection containment chamber and high-resolution histology. The output quality was studied using the histological data from the whole cerebellum and brain stem processed using a high-resolution cryohistology pipeline at 0.5 µm per pixel, in-plane resolution with serial sections at 20-µm thickness. To understand the pathophysiology of highly infectious diseases, it is necessary to have a safe and cost-effective method of performing high-risk autopsies and a standardized perfusion protocol for preparing high-quality tissues. Using the low-cost infection containment chamber, we detail the cranial autopsy protocol and ex situ perfusion-fixation of 4 highly infectious adult human brains. The digitized high-resolution histology images of the Nissl-stained series reveal that most of the sections were free of processing artifacts, such as fixation damage, freezing artifacts, and osmotic shock, at the macrocellular and microcellular level. The quality of our protocol was also tested with the highly sensitive immunohistochemistry staining for specific protein markers. Our protocol provides a safe and effective method in high-risk autopsies that allows for the evaluation of pathogen-host interaction, the underlying pathophysiology, and the extent of the infection across the whole brain at microscopic resolutions.

An active texture-based digital atlas enables automated mapping of structures and markers across brains.

Brain atlases enable the mapping of labeled cells and projections from different brains onto a standard coordinate system. We address two issues in the construction and use of atlases. First, expert neuroanatomists ascertain the fine-scale pattern of brain tissue, the 'texture' formed by cellular organization, to define cytoarchitectural borders. We automate the processes of localizing landmark structures and alignment of brains to a reference atlas using machine learning and training data derived from expert annotations. Second, we construct an atlas that is active; that is, augmented with each use. We show that the alignment of new brains to a reference atlas can continuously refine the coordinate system and associated variance. We apply this approach to the adult murine brainstem and achieve a precise alignment of projections in cytoarchitecturally ill-defined regions across brains from different animals.

Rapid Emergence of SARS-CoV-2 in the Greater New York Metropolitan Area: Geolocation, Demographics, Positivity Rates, and Hospitalization for 46 793 Persons Tested by Northwell Health.

BACKGROUND: In March 2020, the greater New York metropolitan area became an epicenter for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. The initial evolution of case incidence has not been well characterized. METHODS: Northwell Health Laboratories tested 46 793 persons for SARS-CoV-2 from 4 March through 10 April. The primary outcome measure was a positive reverse transcription-polymerase chain reaction test for SARS-CoV-2. The secondary outcomes included patient age, sex, and race, if stated; dates the specimen was obtained and the test result; clinical practice site sources; geolocation of patient residence; and hospitalization. RESULTS: From 8 March through 10 April, a total of 26 735 of 46 793 persons (57.1%) tested positive for SARS-CoV-2. Males of each race were disproportionally more affected than females above age 25, with a progressive male predominance as age increased. Of the positive persons, 7292 were hospitalized directly upon presentation; an additional 882 persons tested positive in an ambulatory setting before subsequent hospitalization, a median of 4.8 days later. Total hospitalization rate was thus 8174 persons (30.6% of positive persons). There was a broad range (>10-fold) in the cumulative number of positive cases across individual zip codes following documented first caseincidence. Test positivity was greater for persons living in zip codes with lower annual household income. CONCLUSIONS: Our data reveal that SARS-CoV-2 incidence emerged rapidly and almost simultaneously across a broad demographic population in the region. These findings support the premise that SARS-CoV-2 infection was widely distributed prior to virus testing availability.

Traumatic microbleeds suggest vascular injury and predict disability in traumatic brain injury.

Traumatic microbleeds are small foci of hypointensity seen on T2*-weighted MRI in patients following head trauma that have previously been considered a marker of axonal injury. The linear appearance and location of some traumatic microbleeds suggests a vascular origin. The aims of this study were to: (i) identify and characterize traumatic microbleeds in patients with acute traumatic brain injury; (ii) determine whether appearance of traumatic microbleeds predict clinical outcome; and (iii) describe the pathology underlying traumatic microbleeds in an index patient. Patients presenting to the emergency department following acute head trauma who received a head CT were enrolled within 48 h of injury and received a research MRI. Disability was defined using Glasgow Outcome Scale-Extended ≤6 at follow-up. All magnetic resonance images were interpreted prospectively and were used for subsequent analysis of traumatic microbleeds. Lesions on T2* MRI were stratified based on 'linear' streak-like or 'punctate' petechial-appearing traumatic microbleeds. The brain of an enrolled subject imaged acutely was procured following death for evaluation of traumatic microbleeds using MRI targeted pathology methods. Of the 439 patients enrolled over 78 months, 31% (134/439) had evidence of punctate and/or linear traumatic microbleeds on MRI. Severity of injury, mechanism of injury, and CT findings were associated with traumatic microbleeds on MRI. The presence of traumatic microbleeds was an independent predictor of disability (P < 0.05; odds ratio = 2.5). No differences were found between patients with punctate versus linear appearing microbleeds. Post-mortem imaging and histology revealed traumatic microbleed co-localization with iron-laden macrophages, predominately seen in perivascular space. Evidence of axonal injury was not observed in co-localized histopathological sections. Traumatic microbleeds were prevalent in the population studied and predictive of worse outcome. The source of traumatic microbleed signal on MRI appeared to be iron-laden macrophages in the perivascular space tracking a network of injured vessels. While axonal injury in association with traumatic microbleeds cannot be excluded, recognizing traumatic microbleeds as a form of traumatic vascular injury may aid in identifying patients who could benefit from new therapies targeting the injured vasculature and secondary injury to parenchyma.

Relation of koniocellular layers of dorsal lateral geniculate to inferior pulvinar nuclei in common marmosets.

Traditionally, the dorsal lateral geniculate nucleus (LGN) and the inferior pulvinar (IPul) nucleus are considered as anatomically and functionally distinct thalamic nuclei. However, in several primate species it has also been established that the koniocellular (K) layers of LGN and parts of the IPul have a shared pattern of immunoreactivity for the calcium-binding protein calbindin. These calbindin-rich cells constitute a thalamic matrix system which is implicated in thalamocortical synchronisation. Further, the K layers and IPul are both involved in visual processing and have similar connections with retina and superior colliculus. Here, we confirmed the continuity between calbindin-rich cells in LGN K layers and the central lateral division of IPul (IPulCL) in marmoset monkeys. By employing a high-throughput neuronal tracing method, we found that both the K layers and IPulCL form comparable patterns of connections with striate and extrastriate cortices; these connections are largely different to those of the parvocellular and magnocellular laminae of LGN. Retrograde tracer-labelled cells and anterograde tracer-labelled axon terminals merged seamlessly from IPulCL into LGN K layers. These results support continuity between LGN K layers and IPulCL, providing an anatomical basis for functional congruity of this region of the dorsal thalamic matrix and calling into question the traditional segregation between LGN and the inferior pulvinar nucleus.

On variational solutions for whole brain serial-section histology using a Sobolev prior in the computational anatomy random orbit model.

This paper presents a variational framework for dense diffeomorphic atlas-mapping onto high-throughput histology stacks at the 20 µm meso-scale. The observed sections are modelled as Gaussian random fields conditioned on a sequence of unknown section by section rigid motions and unknown diffeomorphic transformation of a three-dimensional atlas. To regularize over the high-dimensionality of our parameter space (which is a product space of the rigid motion dimensions and the diffeomorphism dimensions), the histology stacks are modelled as arising from a first order Sobolev space smoothness prior. We show that the joint maximum a-posteriori, penalized-likelihood estimator of our high dimensional parameter space emerges as a joint optimization interleaving rigid motion estimation for histology restacking and large deformation diffeomorphic metric mapping to atlas coordinates. We show that joint optimization in this parameter space solves the classical curvature non-identifiability of the histology stacking problem. The algorithms are demonstrated on a collection of whole-brain histological image stacks from the Mouse Brain Architecture Project.

ESTIMATING DIFFEOMORPHIC MAPPINGS BETWEEN TEMPLATES AND NOISY DATA: VARIANCE BOUNDS ON THE ESTIMATED CANONICAL VOLUME FORM.

Anatomy is undergoing a renaissance driven by the availability of large digital data sets generated by light microscopy. A central computational task is to map individual data volumes to standardized templates. This is accomplished by regularized estimation of a diffeomorphic transformation between the coordinate systems of the individual data and the template, building the transformation incrementally by integrating a smooth flow field. The canonical volume form of this transformation is used to quantify local growth, atrophy, or cell density. While multiple implementations exist for this estimation, less attention has been paid to the variance of the estimated diffeomorphism for noisy data. Notably, there is an infinite dimensional unobservable space defined by those diffeomorphisms which leave the template invariant. These form the stabilizer subgroup of the diffeomorphic group acting on the template. The corresponding flat directions in the energy landscape are expected to lead to increased estimation variance. Here we show that a least-action principle used to generate geodesics in the space of diffeomor-phisms connecting the subject brain to the template removes the stabilizer. This provides reduced-variance estimates of the volume form. Using simulations we demonstrate that the asymmetric large deformation diffeomorphic mapping methods (LDDMM), which explicitly incorporate the asymmetry between idealized template images and noisy empirical images, provide lower variance estimators than their symmetrized counterparts (cf. ANTs). We derive Cramer-Rao bounds for the variances in the limit of small deformations. Analytical results are shown for the Jacobian in terms of perturbations of the vector fields and divergence of the vector field.

Cell-type-based model explaining coexpression patterns of genes in the brain.

Spatial patterns of gene expression in the vertebrate brain are not independent, as pairs of genes can exhibit complex patterns of coexpression. Two genes may be similarly expressed in one region, but differentially expressed in other regions. These correlations have been studied quantitatively, particularly for the Allen Atlas of the adult mouse brain, but their biological meaning remains obscure. We propose a simple model of the coexpression patterns in terms of spatial distributions of underlying cell types and establish its plausibility using independently measured cell-type-specific transcriptomes. The model allows us to predict the spatial distribution of cell types in the mouse brain.

Conventions and nomenclature for double diffusion encoding NMR and MRI.

Stejskal and Tanner's ingenious pulsed field gradient design from 1965 has made diffusion NMR and MRI the mainstay of most studies seeking to resolve microstructural information in porous systems in general and biological systems in particular. Methods extending beyond Stejskal and Tanner's design, such as double diffusion encoding (DDE) NMR and MRI, may provide novel quantifiable metrics that are less easily inferred from conventional diffusion acquisitions. Despite the growing interest on the topic, the terminology for the pulse sequences, their parameters, and the metrics that can be derived from them remains inconsistent and disparate among groups active in DDE. Here, we present a consensus of those groups on terminology for DDE sequences and associated concepts. Furthermore, the regimes in which DDE metrics appear to provide microstructural information that cannot be achieved using more conventional counterparts (in a model-free fashion) are elucidated. We highlight in particular DDE's potential for determining microscopic diffusion anisotropy and microscopic fractional anisotropy, which offer metrics of microscopic features independent of orientation dispersion and thus provide information complementary to the standard, macroscopic, fractional anisotropy conventionally obtained by diffusion MR. Finally, we discuss future vistas and perspectives for DDE.

De novo establishment of wild-type song culture in the zebra finch.

Culture is typically viewed as consisting of traits inherited epigenetically, through social learning. However, cultural diversity has species-typical constraints, presumably of genetic origin. A celebrated, if contentious, example is whether a universal grammar constrains syntactic diversity in human languages. Oscine songbirds exhibit song learning and provide biologically tractable models of culture: members of a species show individual variation in song and geographically separated groups have local song dialects. Different species exhibit distinct song cultures, suggestive of genetic constraints. Without such constraints, innovations and copying errors should cause unbounded variation over multiple generations or geographical distance, contrary to observations. Here we report an experiment designed to determine whether wild-type song culture might emerge over multiple generations in an isolated colony founded by isolates, and, if so, how this might happen and what type of social environment is required. Zebra finch isolates, unexposed to singing males during development, produce song with characteristics that differ from the wild-type song found in laboratory or natural colonies. In tutoring lineages starting from isolate founders, we quantified alterations in song across tutoring generations in two social environments: tutor-pupil pairs in sound-isolated chambers and an isolated semi-natural colony. In both settings, juveniles imitated the isolate tutors but changed certain characteristics of the songs. These alterations accumulated over learning generations. Consequently, songs evolved towards the wild-type in three to four generations. Thus, species-typical song culture can appear de novo. Our study has parallels with language change and evolution. In analogy to models in quantitative genetics, we model song culture as a multigenerational phenotype partly encoded genetically in an isolate founding population, influenced by environmental variables and taking multiple generations to emerge.

Co-expression profiling of autism genes in the mouse brain.

Autism spectrum disorder (ASD) is one of the most prevalent and highly heritable neurodevelopmental disorders in humans. There is significant evidence that the onset and severity of ASD is governed in part by complex genetic mechanisms affecting the normal development of the brain. To date, a number of genes have been associated with ASD. However, the temporal and spatial co-expression of these genes in the brain remain unclear. To address this issue, we examined the co-expression network of 26 autism genes from AutDB (http://mindspec.org/autdb.html), in the framework of 3,041 genes whose expression energies have the highest correlation between the coronal and sagittal images from the Allen Mouse Brain Atlas database (http://mouse.brain-map.org). These data were derived from in situ hybridization experiments conducted on male, 56-day old C57BL/6J mice co-registered to the Allen Reference Atlas, and were used to generate a normalized co-expression matrix indicating the cosine similarity between expression vectors of genes in this database. The network formed by the autism-associated genes showed a higher degree of co-expression connectivity than seen for the other genes in this dataset (Kolmogorov-Smirnov P = 5×10⁻²8). Using Monte Carlo simulations, we identified two cliques of co-expressed genes that were significantly enriched with autism genes (A Bonferroni corrected P<0.05). Genes in both these cliques were significantly over-expressed in the cerebellar cortex (P = 1×10⁻5) suggesting possible implication of this brain region in autism. In conclusion, our study provides a detailed profiling of co-expression patterns of autism genes in the mouse brain, and suggests specific brain regions and new candidate genes that could be involved in autism etiology.

NSF workshop report: discovering general principles of nervous system organization by comparing brain maps across species.

Efforts to understand nervous system structure and function have received new impetus from the federal Brain Research through Advancing Innovative Neurotechnologies (BRAIN) Initiative. Comparative analyses can contribute to this effort by leading to the discovery of general principles of neural circuit design, information processing, and gene-structure-function relationships that are not apparent from studies on single species. We here propose to extend the comparative approach to nervous system 'maps' comprising molecular, anatomical, and physiological data. This research will identify which neural features are likely to generalize across species, and which are unlikely to be broadly conserved. It will also suggest causal relationships between genes, development, adult anatomy, physiology, and, ultimately, behavior. These causal hypotheses can then be tested experimentally. Finally, insights from comparative research can inspire and guide technological development. To promote this research agenda, we recommend that teams of investigators coalesce around specific research questions and select a set of 'reference species' to anchor their comparative analyses. These reference species should be chosen not just for practical advantages, but also with regard for their phylogenetic position, behavioral repertoire, well-annotated genome, or other strategic reasons. We envision that the nervous systems of these reference species will be mapped in more detail than those of other species. The collected data may range from the molecular to the behavioral, depending on the research question. To integrate across levels of analysis and across species, standards for data collection, annotation, archiving, and distribution must be developed and respected. To that end, it will help to form networks or consortia of researchers and centers for science, technology, and education that focus on organized data collection, distribution, and training. These activities could be supported, at least in part, through existing mechanisms at NSF, NIH, and other agencies. It will also be important to develop new integrated software and database systems for cross-species data analyses. Multidisciplinary efforts to develop such analytical tools should be supported financially. Finally, training opportunities should be created to stimulate multidisciplinary, integrative research into brain structure, function, and evolution.

Establishing neuroanatomical correspondences across mouse and marmoset brain structures.

Interest in the common marmoset is growing due to evolutionarily proximity to humans compared to laboratory mice, necessitating a comparison of mouse and marmoset brain architectures, including connectivity and cell type distributions. Creating an actionable comparative platform is challenging since these brains have distinct spatial organizations and expert neuroanatomists disagree. We propose a general theoretical framework to relate named atlas compartments across taxa and use it to establish a detailed correspondence between marmoset and mice brains. Contrary to conventional wisdom that brain structures may be easier to relate at higher levels of the atlas hierarchy, we find that finer parcellations at the leaf levels offer greater reconcilability despite naming discrepancies. Utilizing existing atlases and associated literature, we created a list of leaf-level structures for both species and establish five types of correspondence between them. One-to-one relations were found between 43% of the structures in mouse and 47% in marmoset, whereas 25% of mouse and 10% of marmoset structures were not relatable. The remaining structures show a set of more complex mappings which we quantify. Implementing this correspondence with volumetric atlases of the two species, we make available a computational tool for querying and visualizing relationships between the corresponding brains. Our findings provide a foundation for computational comparative analyses of mesoscale connectivity and cell type distributions in the laboratory mouse and the common marmoset.

Histological characterization and development of mesial surface sulci in the human brain at 13-15 gestational weeks through high-resolution histology.

Cellular-level anatomical data from early fetal brain are sparse yet critical to the understanding of neurodevelopmental disorders. We characterize the organization of the human cerebral cortex between 13 and 15 gestational weeks using high-resolution whole-brain histological data sets complimented with multimodal imaging. We observed the heretofore underrecognized, reproducible presence of infolds on the mesial surface of the cerebral hemispheres. Of note at this stage, when most of the cerebrum is occupied by lateral ventricles and the corpus callosum is incompletely developed, we postulate that these mesial infolds represent the primordial stage of cingulate, callosal, and calcarine sulci, features of mesial cortical development. Our observations are based on the multimodal approach and further include histological three-dimensional reconstruction that highlights the importance of the plane of sectioning. We describe the laminar organization of the developing cortical mantle, including these infolds from the marginal to ventricular zone, with Nissl, hematoxylin and eosin, and glial fibrillary acidic protein (GFAP) immunohistochemistry. Despite the absence of major sulci on the dorsal surface, the boundaries among the orbital, frontal, parietal, and occipital cortex were very well demarcated, primarily by the cytoarchitecture differences in the organization of the subplate (SP) and intermediate zone (IZ) in these locations. The parietal region has the thickest cortical plate (CP), SP, and IZ, whereas the orbital region shows the thinnest CP and reveals an extra cell-sparse layer above the bilaminar SP. The subcortical structures show intensely GFAP-immunolabeled soma, absent in the cerebral mantle. Our findings establish a normative neurodevelopment baseline at the early stage.

Transcriptomic cytoarchitecture reveals principles of human neocortex organization.

Variation in cytoarchitecture is the basis for the histological definition of cortical areas. We used single cell transcriptomics and performed cellular characterization of the human cortex to better understand cortical areal specialization. Single-nucleus RNA-sequencing of 8 areas spanning cortical structural variation showed a highly consistent cellular makeup for 24 cell subclasses. However, proportions of excitatory neuron subclasses varied substantially, likely reflecting differences in connectivity across primary sensorimotor and association cortices. Laminar organization of astrocytes and oligodendrocytes also differed across areas. Primary visual cortex showed characteristic organization with major changes in the excitatory to inhibitory neuron ratio, expansion of layer 4 excitatory neurons, and specialized inhibitory neurons. These results lay the groundwork for a refined cellular and molecular characterization of human cortical cytoarchitecture and areal specialization.

Gamma-band synchronization in visual cortex predicts speed of change detection.

Our capacity to process and respond behaviourally to multiple incoming stimuli is very limited. To optimize the use of this limited capacity, attentional mechanisms give priority to behaviourally relevant stimuli at the expense of irrelevant distractors. In visual areas, attended stimuli induce enhanced responses and an improved synchronization of rhythmic neuronal activity in the gamma frequency band (40-70 Hz). Both effects probably improve the neuronal signalling of attended stimuli within and among brain areas. Attention also results in improved behavioural performance and shortened reaction times. However, it is not known how reaction times are related to either response strength or gamma-band synchronization in visual areas. Here we show that behavioural response times to a stimulus change can be predicted specifically by the degree of gamma-band synchronization among those neurons in monkey visual area V4 that are activated by the behaviourally relevant stimulus. When there are two visual stimuli and monkeys have to detect a change in one stimulus while ignoring the other, their reactions are fastest when the relevant stimulus induces strong gamma-band synchronization before and after the change in stimulus. This enhanced gamma-band synchronization is also followed by shorter neuronal response latencies on the fast trials. Conversely, the monkeys' reactions are slowest when gamma-band synchronization is high in response to the irrelevant distractor. Thus, enhanced neuronal gamma-band synchronization and shortened neuronal response latencies to an attended stimulus seem to have direct effects on visually triggered behaviour, reflecting an early neuronal correlate of efficient visuo-motor integration.

Alta-Cyclic: a self-optimizing base caller for next-generation sequencing.

Next-generation sequencing is limited to short read lengths and by high error rates. We systematically analyzed sources of noise in the Illumina Genome Analyzer that contribute to these high error rates and developed a base caller, Alta-Cyclic, that uses machine learning to compensate for noise factors. Alta-Cyclic substantially improved the number of accurate reads for sequencing runs up to 78 bases and reduced systematic biases, facilitating confident identification of sequence variants.

Clustering of spatial gene expression patterns in the mouse brain and comparison with classical neuroanatomy.

Spatial gene expression profiles provide a novel means of exploring the structural organization of the brain. Computational analysis of these patterns is made possible by genome-scale mapping of the C57BL/6J mouse brain in the Allen Brain Atlas. Here we describe methodology used to explore the spatial structure of gene expression patterns across a set of 3041 genes chosen on the basis of consistency across experimental observations (N=2). The analysis was performed on smoothed, co-registered 3D expression volumes for each gene obtained by aggregating cellular resolution image data. Following dimensionality and noise reduction, voxels were clustered according to similarity of expression across the gene set. We illustrate the resulting parcellations of the mouse brain for different numbers of clusters (K) and quantitatively compare these parcellations with a classically-defined anatomical reference atlas at different levels of granularity, revealing a high degree of correspondence. These observations suggest that spatial localization of gene expression offers substantial promise in connecting knowledge at the molecular level with higher-level information about brain organization.