Functional specialization and distributed processing across marmoset lateral prefrontal subregions.

A prominent aspect of primate lateral prefrontal cortex organization is its division into several cytoarchitecturally distinct subregions. Neurophysiological investigations in macaques have provided evidence for the functional specialization of these subregions, but an understanding of the relative representational topography of sensory, social, and cognitive processes within them remains elusive. One explanatory factor is that evidence for functional specialization has been compiled largely from a patchwork of findings across studies, in many animals, and with considerable variation in stimulus sets and tasks. Here, we addressed this by leveraging the common marmoset (Callithrix jacchus) to carry out large-scale neurophysiological mapping of the lateral prefrontal cortex using high-density microelectrode arrays, and a diverse suite of test stimuli including faces, marmoset calls, and spatial working memory task. Task-modulated units and units responsive to visual and auditory stimuli were distributed throughout the lateral prefrontal cortex, while those with saccade-related activity or face-selective responses were restricted to 8aV, 8aD, 10, 46 V, and 47. Neurons with contralateral visual receptive fields were limited to areas 8aV and 8aD. These data reveal a mixed pattern of functional specialization in the lateral prefrontal cortex, in which responses to some stimuli and tasks are distributed broadly across lateral prefrontal cortex subregions, while others are more limited in their representation.

Neural activity for complex sounds in the marmoset anterior cingulate cortex.

Vocalizations play an important role in the daily life of nonhuman primates and are likely precursors of human language. Recent functional imaging studies in the highly vocal common marmoset (Callithrix jacchus) have suggested that anterior cingulate cortex (ACC) area 32 may be a part of a vocalization-processing network but the response properties of area 32 neurons to auditory stimuli remain unknown. Here we perform electrophysiological recordings in area 32 in marmosets with high-density Neuropixels probes and characterize neuronal responses to a variety of sounds including conspecific vocalizations. Nearly half of the neurons in area 32 respond to conspecific vocalizations and other complex auditory stimuli. These responses exhibit dynamics consisting of an initially non-selective reduction in neural activity, followed by an increase in activity that immediately conveys sound selectivity. Our findings demonstrate that primate ACC area 32 processes species-specific and biologically relevant sounds.

Small data methods in omics: the power of one.

Over the last decade, biology has begun utilizing 'big data' approaches, resulting in large, comprehensive atlases in modalities ranging from transcriptomics to neural connectomics. However, these approaches must be complemented and integrated with 'small data' approaches to efficiently utilize data from individual labs. Integration of smaller datasets with major reference atlases is critical to provide context to individual experiments, and approaches toward integration of large and small data have been a major focus in many fields in recent years. Here we discuss progress in integration of small data with consortium-sized atlases across multiple modalities, and its potential applications. We then examine promising future directions for utilizing the power of small data to maximize the information garnered from small-scale experiments. We envision that, in the near future, international consortia comprising many laboratories will work together to collaboratively build reference atlases and foundation models using small data methods.

Single-cell spatial transcriptomics reveals distinct patterns of dysregulation in non-neuronal and neuronal cells induced by the Trem2R47H Alzheimer's risk gene mutation.

The R47H missense mutation of the TREM2 gene is a known risk factor for development of Alzheimer's Disease. In this study, we analyze the impact of the Trem2R47H mutation on specific cell types in multiple cortical and subcortical brain regions in the context of wild-type and 5xFAD mouse background. We profile 19 mouse brain sections consisting of wild-type, Trem2R47H, 5xFAD and Trem2R47H; 5xFAD genotypes using MERFISH spatial transcriptomics, a technique that enables subcellular profiling of spatial gene expression. Spatial transcriptomics and neuropathology data are analyzed using our custom pipeline to identify plaque and Trem2R47H-induced transcriptomic dysregulation. We initially analyze cell type-specific transcriptomic alterations induced by plaque proximity. Next, we analyze spatial distributions of disease associated microglia and astrocytes, and how they vary between 5xFAD and Trem2R47H; 5xFAD mouse models. Finally, we analyze the impact of the Trem2R47H mutation on neuronal transcriptomes. The Trem2R47H mutation induces consistent upregulation of Bdnf and Ntrk2 across many cortical excitatory neuron types, independent of amyloid pathology. Spatial investigation of genotype enriched subclusters identified spatially localized neuronal subpopulations reduced in 5xFAD and Trem2R47H; 5xFAD mice. Overall, our MERFISH spatial transcriptomics analysis identifies glial and neuronal transcriptomic alterations induced independently by 5xFAD and Trem2R47H mutations, impacting inflammatory responses in microglia and astrocytes, and activity and BDNF signaling in neurons.

ULV: A robust statistical method for clustered data, with applications to multi-subject, single-cell omics data.

Molecular and genomic technological advancements have greatly enhanced our understanding of biological processes by allowing us to quantify key biological variables such as gene expression, protein levels, and microbiome compositions. These breakthroughs have enabled us to achieve increasingly higher levels of resolution in our measurements, exemplified by our ability to comprehensively profile biological information at the single-cell level. However, the analysis of such data faces several critical challenges: limited number of individuals, non-normality, potential dropouts, outliers, and repeated measurements from the same individual. In this article, we propose a novel method, which we call U-statistic based latent variable (ULV). Our proposed method takes advantage of the robustness of rank-based statistics and exploits the statistical efficiency of parametric methods for small sample sizes. It is a computationally feasible framework that addresses all the issues mentioned above simultaneously. We show that our method controls false positives at desired significance levels. An additional advantage of ULV is its flexibility in modeling various types of single-cell data, including both RNA and protein abundance. The usefulness of our method is demonstrated in two studies: a single-cell proteomics study of acute myelogenous leukemia (AML) and a single-cell RNA study of COVID-19 symptoms. In the AML study, ULV successfully identified differentially expressed proteins that would have been missed by the pseudobulk version of the Wilcoxon rank-sum test. In the COVID-19 study, ULV identified genes associated with covariates such as age and gender, and genes that would be missed without adjusting for covariates. The differentially expressed genes identified by our method are less biased toward genes with high expression levels. Furthermore, ULV identified additional gene pathways likely contributing to the mechanisms of COVID-19 severity.

Single-cell spatial transcriptomics reveals a dystrophic trajectory following a developmental bifurcation of myoblast cell fates in facioscapulohumeral muscular dystrophy.

Facioscapulohumeral muscular dystrophy (FSHD) is linked to abnormal derepression of the transcription activator DUX4. This effect is localized to a low percentage of cells, requiring single-cell analysis. However, single-cell/nucleus RNA-seq cannot fully capture the transcriptome of multinucleated large myotubes. To circumvent these issues, we use multiplexed error-robust fluorescent in situ hybridization (MERFISH) spatial transcriptomics that allows profiling of RNA transcripts at a subcellular resolution. We simultaneously examined spatial distributions of 140 genes, including 24 direct DUX4 targets, in in vitro differentiated myotubes and unfused mononuclear cells (MNCs) of control, isogenic D4Z4 contraction mutant and FSHD patient samples, as well as the individual nuclei within them. We find myocyte nuclei segregate into two clusters defined by the expression of DUX4 target genes, which is exclusively found in patient/mutant nuclei, whereas MNCs cluster based on developmental states. Patient/mutant myotubes are found in "FSHD-hi" and "FSHD-lo" states with the former signified by high DUX4 target expression and decreased muscle gene expression. Pseudotime analyses reveal a clear bifurcation of myoblast differentiation into control and FSHD-hi myotube branches, with variable numbers of DUX4 target-expressing nuclei found in multinucleated FSHD-hi myotubes. Gene coexpression modules related to extracellular matrix and stress gene ontologies are significantly altered in patient/mutant myotubes compared with the control. We also identify distinct subpathways within the DUX4 gene network that may differentially contribute to the disease transcriptomic phenotype. Taken together, our MERFISH-based study provides effective gene network profiling of multinucleated cells and identifies FSHD-induced transcriptomic alterations during myoblast differentiation.

Laminar Dynamics of Target Selection in the Posterior Parietal Cortex of the Common Marmoset.

The lateral intraparietal area (LIP) plays a crucial role in target selection and attention in primates, but the laminar microcircuitry of this region is largely unknown. To address this, we used ultra-high density laminar electrophysiology with Neuropixels probes to record neural activity in the posterior parietal cortex (PPC) of two adult marmosets while they performed a simple visual target selection task. Our results reveal neural correlates of visual target selection in the marmoset, similar to those observed in macaques and humans, with distinct timing and profiles of activity across cell types and cortical layers. Notably, a greater proportion of neurons exhibited stimulus-related activity in superficial layers whereas a greater proportion of infragranular neurons exhibited significant postsaccadic activity. Stimulus-related activity was first observed in granular layer putative interneurons, whereas target discrimination activity emerged first in supragranular layers putative pyramidal neurons, supporting a canonical laminar circuit underlying visual target selection in marmoset PPC. These findings provide novel insights into the neural basis of visual attention and target selection in primates.

A ubiquitous spectrolaminar motif of local field potential power across the primate cortex.

The mammalian cerebral cortex is anatomically organized into a six-layer motif. It is currently unknown whether a corresponding laminar motif of neuronal activity patterns exists across the cortex. Here we report such a motif in the power of local field potentials (LFPs). Using laminar probes, we recorded LFPs from 14 cortical areas across the cortical hierarchy in five macaque monkeys. The laminar locations of recordings were histologically identified by electrolytic lesions. Across all areas, we found a ubiquitous spectrolaminar pattern characterized by an increasing deep-to-superficial layer gradient of high-frequency power peaking in layers 2/3 and an increasing superficial-to-deep gradient of alpha-beta power peaking in layers 5/6. Laminar recordings from additional species showed that the spectrolaminar pattern is highly preserved among primates-macaque, marmoset and human-but more dissimilar in mouse. Our results suggest the existence of a canonical layer-based and frequency-based mechanism for cortical computation.

Single cell spatial transcriptomics reveals distinct patterns of dysregulation in non-neuronal and neuronal cells induced by the Trem2R47H Alzheimer's risk gene mutation.

INTRODUCTION: The R47H missense mutation of the TREM2 gene is a strong risk factor for development of Alzheimer's Disease. We investigate cell-type-specific spatial transcriptomic changes induced by the Trem2R47H mutation to determine the impacts of this mutation on transcriptional dysregulation. METHODS: We profiled 15 mouse brain sections consisting of wild-type, Trem2R47H, 5xFAD and Trem2R47H; 5xFAD genotypes using MERFISH spatial transcriptomics. Single-cell spatial transcriptomics and neuropathology data were analyzed using our custom pipeline to identify plaque and Trem2R47H induced transcriptomic dysregulation. RESULTS: The Trem2R47H mutation induced consistent upregulation of Bdnf and Ntrk2 across many cortical excitatory neuron types, independent of amyloid pathology. Spatial investigation of genotype enriched subclusters identified spatially localized neuronal subpopulations reduced in 5xFAD and Trem2R47H; 5xFAD mice. CONCLUSION: Spatial transcriptomics analysis identifies glial and neuronal transcriptomic alterations induced independently by 5xFAD and Trem2R47H mutations, impacting inflammatory responses in microglia and astrocytes, and activity and BDNF signaling in neurons.

Anatomical and molecular characterization of parvalbumin-cholecystokinin co-expressing inhibitory interneurons: implications for neuropsychiatric conditions.

Inhibitory interneurons are crucial to brain function and their dysfunction is implicated in neuropsychiatric conditions. Emerging evidence indicates that cholecystokinin (CCK)-expressing interneurons (CCK+) are highly heterogenous. We find that a large subset of parvalbumin-expressing (PV+) interneurons express CCK strongly; between 40 and 56% of PV+ interneurons in mouse hippocampal CA1 express CCK. Primate interneurons also exhibit substantial PV/CCK co-expression. Mouse PV+/CCK+ and PV+/CCK- cells show distinguishable electrophysiological and molecular characteristics. Analysis of single nuclei RNA-seq and ATAC-seq data shows that PV+/CCK+ cells are a subset of PV+ cells, not of synuclein gamma positive (SNCG+) cells, and that they strongly express oxidative phosphorylation (OXPHOS) genes. We find that mitochondrial complex I and IV-associated OXPHOS gene expression is strongly correlated with CCK expression in PV+ interneurons at both the transcriptomic and protein levels. Both PV+ interneurons and dysregulation of OXPHOS processes are implicated in neuropsychiatric conditions, including autism spectrum (ASD) disorder and schizophrenia (SCZ). Analysis of human brain samples from patients with these conditions shows alterations in OXPHOS gene expression. Together these data reveal important molecular characteristics of PV-CCK co-expressing interneurons and support their implication in neuropsychiatric conditions.

Inferring neuron-neuron communications from single-cell transcriptomics through NeuronChat.

Neural communication networks form the fundamental basis for brain function. These communication networks are enabled by emitted ligands such as neurotransmitters, which activate receptor complexes to facilitate communication. Thus, neural communication is fundamentally dependent on the transcriptome. Here we develop NeuronChat, a method and package for the inference, visualization and analysis of neural-specific communication networks among pre-defined cell groups using single-cell expression data. We incorporate a manually curated molecular interaction database of neural signaling for both human and mouse, and benchmark NeuronChat on several published datasets to validate its ability in predicting neural connectivity. Then, we apply NeuronChat to three different neural tissue datasets to illustrate its functionalities in identifying interneural communication networks, revealing conserved or context-specific interactions across different biological contexts, and predicting communication pattern changes in diseased brains with autism spectrum disorder. Finally, we demonstrate NeuronChat can utilize spatial transcriptomics data to infer and visualize neural-specific cell-cell communication.

Inferring neuron-neuron communications from single-cell transcriptomics through NeuronChat.

Neural communication networks form the fundamental basis for brain function. These communication networks are enabled by emitted ligands such as neurotransmitters, which activate receptor complexes to facilitate communication. Thus, neural communication is fundamentally dependent on the transcriptome. Here we develop NeuronChat, a method and package for the inference, visualization and analysis of neural-specific communication networks among pre-defined cell groups using single-cell expression data. We incorporate a manually curated molecular interaction database of neural signaling for both human and mouse, and benchmark NeuronChat on several published datasets to validate its ability in predicting neural connectivity. Then, we apply NeuronChat to three different neural tissue datasets to illustrate its functionalities in identifying interneural communication networks, revealing conserved or context-specific interactions across different biological contexts, and predicting communication pattern changes in diseased brains with autism spectrum disorder. Finally, we demonstrate NeuronChat can utilize spatial transcriptomics data to infer and visualize neural-specific cell-cell communication.

Degenerate mapping of environmental location presages deficits in object-location encoding and memory in the 5xFAD mouse model for Alzheimer's disease.

A key challenge in developing diagnosis and treatments for Alzheimer's disease (AD) is to detect abnormal network activity at as early a stage as possible. To date, behavioral and neurophysiological investigations in AD model mice have yet to conduct a longitudinal assessment of cellular pathology, memory deficits, and neurophysiological correlates of neuronal activity. We therefore examined the temporal relationships between pathology, neuronal activities and spatial representation of environments, as well as object location memory deficits across multiple stages of development in the 5xFAD mice model and compared these results to those observed in wild-type mice. We performed longitudinal in vivo calcium imaging with miniscope on hippocampal CA1 neurons in behaving mice. We find that 5xFAD mice show amyloid plaque accumulation, depressed neuronal calcium activity during immobile states, and degenerate and unreliable hippocampal neuron spatial tuning to environmental location at early stages by 4 months of age while their object location memory (OLM) is comparable to WT mice. By 8 months of age, 5xFAD mice show deficits of OLM, which are accompanied by progressive degradation of spatial encoding and, eventually, impaired CA1 neural tuning to object-location pairings. Furthermore, depressed neuronal activity and unreliable spatial encoding at early stage are correlated with impaired performance in OLM at 8-month-old. Our results indicate the close connection between impaired hippocampal tuning to object-location and the presence of OLM deficits. The results also highlight that depressed baseline firing rates in hippocampal neurons during immobile states and unreliable spatial representation precede object memory deficits and predict memory deficits at older age, suggesting potential early opportunities for AD detecting.

Delay-related activity in marmoset prefrontal cortex.

Persistent delay-period activity in prefrontal cortex (PFC) has long been regarded as a neural signature of working memory (WM). Electrophysiological investigations in macaque PFC have provided much insight into WM mechanisms; however, a barrier to understanding is the fact that a portion of PFC lies buried within the principal sulcus in this species and is inaccessible for laminar electrophysiology or optical imaging. The relatively lissencephalic cortex of the New World common marmoset (Callithrix jacchus) circumvents such limitations. It remains unknown, however, whether marmoset PFC neurons exhibit persistent activity. Here, we addressed this gap by conducting wireless electrophysiological recordings in PFC of marmosets performing a delayed-match-to-location task on a home cage-based touchscreen system. As in macaques, marmoset PFC neurons exhibited sample-, delay-, and response-related activity that was directionally tuned and linked to correct task performance. Models constructed from population activity consistently and accurately predicted stimulus location throughout the delay period, supporting a framework of delay activity in which mnemonic representations are relatively stable in time. Taken together, our findings support the existence of common neural mechanisms underlying WM performance in PFC of macaques and marmosets and thus validate the marmoset as a suitable model animal for investigating the microcircuitry underlying WM.

Tracking longitudinal population dynamics of single neuronal calcium signal using SCOUT.

In vivo calcium imaging enables simultaneous recording of large neuronal ensembles engaged in complex operations. Many experiments require monitoring and identification of cell populations across multiple sessions. Population cell tracking across multiple sessions is complicated by non-rigid transformations induced by cell movement and imaging field shifts. We introduce SCOUT (Single-Cell spatiOtemporal longitUdinal Tracking), a fast, robust cell-tracking method utilizing multiple cell-cell similarity metrics, probabilistic inference, and an adaptive clustering methodology, to perform cell identification across multiple sessions. By comparing SCOUT with earlier cell-tracking algorithms on simulated, 1-photon, and 2-photon recordings, we show that our approach significantly improves cell-tracking quality, particularly when recordings exhibit spatial footprint movement between sessions or sub-optimal neural extraction quality.

Simultaneous functional MRI of two awake marmosets.

Social cognition is a dynamic process that requires the perception and integration of a complex set of idiosyncratic features between interacting conspecifics. Here we present a method for simultaneously measuring the whole-brain activation of two socially interacting marmoset monkeys using functional magnetic resonance imaging. MRI hardware (a radiofrequency coil and peripheral devices) and image-processing pipelines were developed to assess brain responses to socialization, both on an intra-brain and inter-brain level. Notably, the brain activation of a marmoset when viewing a second marmoset in-person versus when viewing a pre-recorded video of the same marmoset-i.e., when either capable or incapable of socially interacting with a visible conspecific-demonstrates increased activation in the face-patch network. This method enables a wide range of possibilities for potentially studying social function and dysfunction in a non-human primate model.

Interspecies activation correlations reveal functional correspondences between marmoset and human brain areas.

The common marmoset has enormous promise as a nonhuman primate model of human brain functions. While resting-state functional MRI (fMRI) has provided evidence for a similar organization of marmoset and human cortices, the technique cannot be used to map the functional correspondences of brain regions between species. This limitation can be overcome by movie-driven fMRI (md-fMRI), which has become a popular tool for noninvasively mapping the neural patterns generated by rich and naturalistic stimulation. Here, we used md-fMRI in marmosets and humans to identify whole-brain functional correspondences between the two primate species. In particular, we describe functional correlates for the well-known human face, body, and scene patches in marmosets. We find that these networks have a similar organization in both species, suggesting a largely conserved organization of higher-order visual areas between New World marmoset monkeys and humans. However, while face patches in humans and marmosets were activated by marmoset faces, only human face patches responded to the faces of other animals. Together, the results demonstrate that higher-order visual processing might be a conserved feature between humans and New World marmoset monkeys but that small, potentially important functional differences exist.

Ketamine disrupts gaze patterns during face viewing in the common marmoset.

Faces are stimuli of critical importance for primates. The common marmoset (Callithrix jacchus) is a promising model for investigations of face processing, as this species possesses oculomotor and face-processing networks resembling those of macaques and humans. Face processing is often disrupted in neuropsychiatric conditions such as schizophrenia (SZ), and thus, it is important to recapitulate underlying circuitry dysfunction preclinically. The N-methyl-d-aspartate (NMDA) noncompetitive antagonist ketamine has been used extensively to model the cognitive symptoms of SZ. Here, we investigated the effects of a subanesthetic dose of ketamine on oculomotor behavior in marmosets during face viewing. Four marmosets received systemic ketamine or saline injections while viewing phase-scrambled or intact videos of conspecifics' faces. To evaluate effects of ketamine on scan paths during face viewing, we identified regions of interest in each face video and classified locations of saccade onsets and landing positions within these areas. A preference for the snout over eye regions was observed following ketamine administration. In addition, regions in which saccades landed could be significantly predicted by saccade onset region in the saline but not the ketamine condition. Effects on saccade control were limited to an increase in saccade peak velocity in all conditions and a reduction in saccade amplitudes during viewing of scrambled videos. Thus, ketamine induced a significant disruption of scan paths during viewing of conspecific faces but limited effects on saccade motor control. These findings support the use of ketamine in marmosets for investigating changes in neural circuits underlying social cognition in neuropsychiatric disorders.NEW & NOTEWORTHY Face processing, an important social cognitive ability, is impaired in neuropsychiatric conditions such as schizophrenia. The highly social common marmoset model presents an opportunity to investigate these impairments. We administered subanesthetic doses of ketamine to marmosets to model the cognitive symptoms of schizophrenia. We observed a disruption of scan paths during viewing of conspecifics' faces. These findings support the use of ketamine in marmosets as a model for investigating social cognition in neuropsychiatric disorders.

A Novel Three-Choice Touchscreen Task to Examine Spatial Attention and Orienting Responses in Rodents.

Mammalian orienting behavior consists of coordinated movements of the eyes, head, pinnae, vibrissae, or body to attend to an external stimulus. The present study aimed to develop a novel operant task using a touch-screen system to measure spatial attention. In this task, rats were trained to nose-poke a light stimulus presented in one of three locations. The stimulus was presented more frequently in the center location to develop spatial attention bias toward the center stimulus. Changes in orienting responses were detected by measuring the animals' response accuracy and latency to stimuli at the lateral locations, following reversible unilateral chemogenetic inactivation of the superior colliculus (SC). Additionally, spontaneous turning and rotation behavior was measured using an open-field test (OFT). Our results show that right SC inactivation significantly increased the whole body turn angle in the OFT, in line with previous literature that indicated an ipsiversive orientating bias and the presence of contralateral neglect following unilateral SC lesions. In the touch screen orienting task, unilateral SC inactivation significantly increased bias toward the ipsilateral side, as measured by response frequency in various experimental conditions, and a very large left-shift of a respective psychometric function. Our results demonstrate that this novel touchscreen task is able to detect changes in spatial attention and orienting responses because of e.g. experimental manipulations or injury with very high sensitivity, while taking advantage of the touch screen technology that allows for high transferability of the task between labs and for open-source data sharing through https://www.mousebytes.ca.

Structural alterations in cortical and thalamocortical white matter tracts after recovery from prefrontal cortex lesions in macaques.

Unilateral damage to the frontoparietal network typically impairs saccade target selection within the contralesional visual hemifield. Severity of deficits and the degree of recovery have been associated with widespread network dysfunction, yet it is not clear how these behavioural and functional brain changes relate with the underlying structural white matter tracts. Here, we investigated whether recovery after unilateral prefrontal cortex (PFC) lesions was associated with changes in white matter microstructure across large-scale frontoparietal cortical and thalamocortical networks. Diffusion-weighted imaging was acquired in four male rhesus macaques at pre-lesion, week 1, and week 8-16 post-lesion when target selection deficits largely recovered. Probabilistic tractography was used to reconstruct cortical frontoparietal fiber tracts, including the superior longitudinal fasciculus (SLF) and transcallosal fibers connecting the PFC or posterior parietal cortex (PPC), as well as thalamocortical fiber tracts connecting the PFC and PPC to thalamic nuclei. We found that the two animals with small PFC lesions showed increased fractional anisotropy in both cortical and thalamocortical fiber tracts when behaviour had recovered. However, we found that fractional anisotropy decreased in cortical frontoparietal tracts after larger PFC lesions yet increased in some thalamocortical tracts at the time of behavioural recovery. These findings indicate that behavioural recovery after small PFC lesions may be supported by both cortical and subcortical areas, whereas larger PFC lesions may have induced widespread structural damage and hindered compensatory remodeling in the cortical frontoparietal network.