Epigenomic diversity of cortical projection neurons in the mouse brain

Projection-TAGs enable multiplex projection tracing and multi-modal profiling of projection neurons

Single-cell multiomic techniques have sparked immense interest in developing a comprehensive multi-modal map of diverse neuronal cell types and their brain wide projections. However, investigating the spatial organization, transcriptional and epigenetic landscapes of brain wide projection neurons is hampered by the lack of efficient and easily adoptable tools. Here we introduce Projection-TAGs, a retrograde AAV platform that allows multiplex tagging of projection neurons using RNA barcodes. By using Projection-TAGs, we performed multiplex projection tracing of the mouse cortex and high-throughput single-cell profiling of the transcriptional and epigenetic landscapes of the cortical projection neurons. Projection-TAGs can be leveraged to obtain a snapshot of activity-dependent recruitment of distinct projection neurons and their molecular features in the context of a specific stimulus. Given its flexibility, usability, and compatibility, we envision that Projection-TAGs can be readily applied to build a comprehensive multi-modal map of brain neuronal cell types and their projections.

High-throughput mapping of single-neuron projection and molecular features by retrograde barcoded labeling

Deciphering patterns of connectivity between neurons in the brain is a critical step toward understanding brain function. Imaging-based neuroanatomical tracing identifies area-to-area or sparse neuron-to-neuron connectivity patterns, but with limited throughput. Barcode-based connectomics maps large numbers of single-neuron projections, but remains a challenge for jointly analyzing single-cell transcriptomics. Here, we established a rAAV2-retro barcode-based multiplexed tracing method that simultaneously characterizes the projectome and transcriptome at the single neuron level. We uncovered dedicated and collateral projection patterns of ventromedial prefrontal cortex (vmPFC) neurons to five downstream targets and found that projection-defined vmPFC neurons are molecularly heterogeneous. We identified transcriptional signatures of projection-specific vmPFC neurons, and verified Pou3f1 as a marker gene enriched in neurons projecting to the lateral hypothalamus, denoting a distinct subset with collateral projections to both dorsomedial striatum and lateral hypothalamus. In summary, we have developed a new multiplexed technique whose paired connectome and gene expression data can help reveal organizational principles that form neural circuits and process information.

Non-CG DNA methylation and MeCP2 stabilize repeated tuning of long genes that distinguish closely related neuron types

The extraordinary diversity of neuron types in the mammalian brain is delineated at the highest resolution by subtle gene expression differences that may require specialized molecular mechanisms to be maintained. Neurons uniquely express the longest genes in the genome and utilize neuron-enriched non-CG DNA methylation (mCA) together with the Rett syndrome protein, MeCP2, to control gene expression, but the function of these unique gene structures and machinery in regulating finely resolved neuron type-specific gene programs has not been explored. Here, we employ epigenomic and spatial transcriptomic analyses to discover a major role for mCA and MeCP2 in maintaining neuron type-specific gene programs at the finest scale of cellular resolution. We uncover differential susceptibility to MeCP2 loss in neuronal populations depending on global mCA levels and dissect methylation patterns and intragenic enhancer repression that drive overlapping and distinct gene regulation between neuron types. Strikingly, we show that mCA and MeCP2 regulate genes that are repeatedly tuned to differentiate neuron types at the highest cellular resolution, including spatially resolved, vision-dependent gene programs in the visual cortex. These repeatedly tuned genes display genomic characteristics, including long length, numerous intragenic enhancers, and enrichment for mCA, that predispose them to regulation by MeCP2. Thus, long gene regulation by the MeCP2 pathway maintains differential gene expression between closely-related neurons to facilitate the exceptional cellular diversity in the complex mammalian brain.

Enabling neighbour labelling: using synthetic biology to explore how cells influence their neighbours

Cell-cell interactions are central to development, but exploring how a change in any given cell relates to changes in the neighbour of that cell can be technically challenging. Here, we review recent developments in synthetic biology and image analysis that are helping overcome this problem. We highlight the opportunities presented by these advances and discuss opportunities and limitations in applying them to developmental model systems.

Revealing spatial multimodal heterogeneity in tissues with SpaTrio

Capturing and depicting the multimodal tissue information of tissues at the spatial scale remains a significant challenge owing to technical limitations in single-cell multi-omics and spatial transcriptomics sequencing. Here, we developed a computational method called SpaTrio that can build spatial multi-omics data by integrating these two datasets through probabilistic alignment and enabling further analysis of gene regulation and cellular interactions. We benchmarked SpaTrio using simulation datasets and demonstrated its accuracy and robustness. Next, we evaluated SpaTrio on biological datasets and showed that it could detect topological patterns of cells and modalities. SpaTrio has also been applied to multiple sets of actual data to uncover spatially multimodal heterogeneity, understand the spatiotemporal regulation of gene expression, and resolve multimodal communication among cells. Our data demonstrated that SpaTrio could accurately map single cells and reconstruct the spatial distribution of various biomolecules, providing valuable multimodal insights into spatial biology.

Brain-wide correspondence of neuronal epigenomics and distant projections

Single-cell analyses parse the brain's billions of neurons into thousands of 'cell-type' clusters residing in different brain structures1. Many cell types mediate their functions through targeted long-distance projections allowing interactions between specific cell types. Here we used epi-retro-seq2 to link single-cell epigenomes and cell types to long-distance projections for 33,034 neurons dissected from 32 different regions projecting to 24 different targets (225 source-to-target combinations) across the whole mouse brain. We highlight uses of these data for interrogating principles relating projection types to transcriptomics and epigenomics, and for addressing hypotheses about cell types and connections related to genetics. We provide an overall synthesis with 926 statistical comparisons of discriminability of neurons projecting to each target for every source. We integrate this dataset into the larger BRAIN Initiative Cell Census Network atlas, composed of millions of neurons, to link projection cell types to consensus clusters. Integration with spatial transcriptomics further assigns projection-enriched clusters to smaller source regions than the original dissections. We exemplify this by presenting in-depth analyses of projection neurons from the hypothalamus, thalamus, hindbrain, amygdala and midbrain to provide insights into properties of those cell types, including differentially expressed genes, their associated cis-regulatory elements and transcription-factor-binding motifs, and neurotransmitter use.

Single-cell DNA methylome and 3D multi-omic atlas of the adult mouse brain

Cytosine DNA methylation is essential in brain development and is implicated in various neurological disorders. Understanding DNA methylation diversity across the entire brain in a spatial context is fundamental for a complete molecular atlas of brain cell types and their gene regulatory landscapes. Here we used single-nucleus methylome sequencing (snmC-seq3) and multi-omic sequencing (snm3C-seq)1 technologies to generate 301,626 methylomes and 176,003 chromatin conformation-methylome joint profiles from 117 dissected regions throughout the adult mouse brain. Using iterative clustering and integrating with companion whole-brain transcriptome and chromatin accessibility datasets, we constructed a methylation-based cell taxonomy with 4,673 cell groups and 274 cross-modality-annotated subclasses. We identified 2.6 million differentially methylated regions across the genome that represent potential gene regulation elements. Notably, we observed spatial cytosine methylation patterns on both genes and regulatory elements in cell types within and across brain regions. Brain-wide spatial transcriptomics data validated the association of spatial epigenetic diversity with transcription and improved the anatomical mapping of our epigenetic datasets. Furthermore, chromatin conformation diversities occurred in important neuronal genes and were highly associated with DNA methylation and transcription changes. Brain-wide cell-type comparisons enabled the construction of regulatory networks that incorporate transcription factors, regulatory elements and their potential downstream gene targets. Finally, intragenic DNA methylation and chromatin conformation patterns predicted alternative gene isoform expression observed in a whole-brain SMART-seq2 dataset. Our study establishes a brain-wide, single-cell DNA methylome and 3D multi-omic atlas and provides a valuable resource for comprehending the cellular-spatial and regulatory genome diversity of the mouse brain.

Functional genomic mechanisms of opioid action and opioid use disorder: a systematic review of animal models and human studies

In the past two decades, over-prescription of opioids for pain management has driven a steep increase in opioid use disorder (OUD) and death by overdose, exerting a dramatic toll on western countries. OUD is a chronic relapsing disease associated with a lifetime struggle to control drug consumption, suggesting that opioids trigger long-lasting brain adaptations, notably through functional genomic and epigenomic mechanisms. Current understanding of these processes, however, remain scarce, and have not been previously reviewed systematically. To do so, the goal of the present work was to synthesize current knowledge on genome-wide transcriptomic and epigenetic mechanisms of opioid action, in primate and rodent species. Using a prospectively registered methodology, comprehensive literature searches were completed in PubMed, Embase, and Web of Science. Of the 2709 articles identified, 73 met our inclusion criteria and were considered for qualitative analysis. Focusing on the 5 most studied nervous system structures (nucleus accumbens, frontal cortex, whole striatum, dorsal striatum, spinal cord; 44 articles), we also conducted a quantitative analysis of differentially expressed genes, in an effort to identify a putative core transcriptional signature of opioids. Only one gene, Cdkn1a, was consistently identified in eleven studies, and globally, our results unveil surprisingly low consistency across published work, even when considering most recent single-cell approaches. Analysis of sources of variability detected significant contributions from species, brain structure, duration of opioid exposure, strain, time-point of analysis, and batch effects, but not type of opioid. To go beyond those limitations, we leveraged threshold-free methods to illustrate how genome-wide comparisons may generate new findings and hypotheses. Finally, we discuss current methodological development in the field, and their implication for future research and, ultimately, better care.

Single-nucleus genomics in outbred rats with divergent cocaine addiction-like behaviors reveals changes in amygdala GABAergic inhibition

The amygdala processes positive and negative valence and contributes to addiction, but the cell-type-specific gene regulatory programs involved are unknown. We generated an atlas of single-nucleus gene expression and chromatin accessibility in the amygdala of outbred rats with high and low cocaine addiction-like behaviors following prolonged abstinence. Differentially expressed genes between the high and low groups were enriched for energy metabolism across cell types. Rats with high addiction index (AI) showed increased relapse-like behaviors and GABAergic transmission in the amygdala. Both phenotypes were reversed by pharmacological inhibition of the glyoxalase 1 enzyme, which metabolizes methylglyoxal-a GABAA receptor agonist produced by glycolysis. Differences in chromatin accessibility between high and low AI rats implicated pioneer transcription factors in the basic helix-loop-helix, FOX, SOX and activator protein 1 families. We observed opposite regulation of chromatin accessibility across many cell types. Most notably, excitatory neurons had greater accessibility in high AI rats and inhibitory neurons had greater accessibility in low AI rats.

Functional genomics and systems biology in human neuroscience

Neuroscience research has entered a phase of key discoveries in the realm of neurogenomics owing to strong financial and intellectual support for resource building and tool development. The previous challenge of tissue heterogeneity has been met with the application of techniques that can profile individual cells at scale. Moreover, the ability to perturb genes, gene regulatory elements and neuronal activity in a cell-type-specific manner has been integrated with gene expression studies to uncover the functional underpinnings of the genome at a systems level. Although these insights have necessarily been grounded in model systems, we now have the opportunity to apply these approaches in humans and in human tissue, thanks to advances in human genetics, brain imaging and tissue collection. We acknowledge that there will probably always be limits to the extent to which we can apply the genomic tools developed in model systems to human neuroscience; however, as we describe in this Perspective, the neuroscience field is now primed with an optimal foundation for tackling this ambitious challenge. The application of systems-level network analyses to these datasets will facilitate a deeper appreciation of human neurogenomics that cannot otherwise be achieved from directly observable phenomena.

Srsf1 and Elavl1 act antagonistically on neuronal fate choice in the developing neocortex by controlling TrkC receptor isoform expression

The seat of higher-order cognitive abilities in mammals, the neocortex, is a complex structure, organized in several layers. The different subtypes of principal neurons are distributed in precise ratios and at specific positions in these layers and are generated by the same neural progenitor cells (NPCs), steered by a spatially and temporally specified combination of molecular cues that are incompletely understood. Recently, we discovered that an alternatively spliced isoform of the TrkC receptor lacking the kinase domain, TrkC-T1, is a determinant of the corticofugal projection neuron (CFuPN) fate. Here, we show that the finely tuned balance between TrkC-T1 and the better known, kinase domain-containing isoform, TrkC-TK+, is cell type-specific in the developing cortex and established through the antagonistic actions of two RNA-binding proteins, Srsf1 and Elavl1. Moreover, our data show that Srsf1 promotes the CFuPN fate and Elavl1 promotes the callosal projection neuron (CPN) fate in vivo via regulating the distinct ratios of TrkC-T1 to TrkC-TK+. Taken together, we connect spatio-temporal expression of Srsf1 and Elavl1 in the developing neocortex with the regulation of TrkC alternative splicing and transcript stability and neuronal fate choice, thus adding to the mechanistic and functional understanding of alternative splicing in vivo.

Signature morphoelectric properties of diverse GABAergic interneurons in the human neocortex

Human cortex transcriptomic studies have revealed a hierarchical organization of γ-aminobutyric acid-producing (GABAergic) neurons from subclasses to a high diversity of more granular types. Rapid GABAergic neuron viral genetic labeling plus Patch-seq (patch-clamp electrophysiology plus single-cell RNA sequencing) sampling in human brain slices was used to reliably target and analyze GABAergic neuron subclasses and individual transcriptomic types. This characterization elucidated transitions between PVALB and SST subclasses, revealed morphological heterogeneity within an abundant transcriptomic type, identified multiple spatially distinct types of the primate-specialized double bouquet cells (DBCs), and shed light on cellular differences between homologous mouse and human neocortical GABAergic neuron types. These results highlight the importance of multimodal phenotypic characterization for refinement of emerging transcriptomic cell type taxonomies and for understanding conserved and specialized cellular properties of human brain cell types.

Infernape uncovers cell type-specific and spatially resolved alternative polyadenylation in the brain

Differential polyadenylation sites (PAs) critically regulate gene expression, but their cell type-specific usage and spatial distribution in the brain have not been systematically characterized. Here, we present Infernape, which infers and quantifies PA usage from single-cell and spatial transcriptomic data and show its application in the mouse brain. Infernape uncovers alternative intronic PAs and 3'-UTR lengthening during cortical neurogenesis. Progenitor-neuron comparisons in the excitatory and inhibitory neuron lineages show overlapping PA changes in embryonic brains, suggesting that the neural proliferation-differentiation axis plays a prominent role. In the adult mouse brain, we uncover cell type-specific PAs and visualize such events using spatial transcriptomic data. Over two dozen neurodevelopmental disorder-associated genes such as Csnk2a1 and Mecp2 show differential PAs during brain development. This study presents Infernape to identify PAs from scRNA-seq and spatial data, and highlights the role of alternative PAs in neuronal gene regulation.

Perturbation of 3D nuclear architecture, epigenomic dysregulation and aging, and cannabinoid synaptopathy reconfigures conceptualization of cannabinoid pathophysiology: part 1-aging and epigenomics

Much recent attention has been directed toward the spatial organization of the cell nucleus and the manner in which three-dimensional topologically associated domains and transcription factories are epigenetically coordinated to precisely bring enhancers into close proximity with promoters to control gene expression. Twenty lines of evidence robustly implicate cannabinoid exposure with accelerated organismal and cellular aging. Aging has recently been shown to be caused by increased DNA breaks. These breaks rearrange and maldistribute the epigenomic machinery to weaken and reverse cellular differentiation, cause genome-wide DNA demethylation, reduce gene transcription, and lead to the inhibition of developmental pathways, which contribute to the progressive loss of function and chronic immune stimulation that characterize cellular aging. Both cell lineage-defining superenhancers and the superanchors that control them are weakened. Cannabis exposure phenocopies the elements of this process and reproduces DNA and chromatin breakages, reduces the DNA, RNA protein and histone synthesis, interferes with the epigenomic machinery controlling both DNA and histone modifications, induces general DNA hypomethylation, and epigenomically disrupts both the critical boundary elements and the cohesin motors that create chromatin loops. This pattern of widespread interference with developmental programs and relative cellular dedifferentiation (which is pro-oncogenic) is reinforced by cannabinoid impairment of intermediate metabolism (which locks in the stem cell-like hyper-replicative state) and cannabinoid immune stimulation (which perpetuates and increases aging and senescence programs, DNA damage, DNA hypomethylation, genomic instability, and oncogenesis), which together account for the diverse pattern of teratologic and carcinogenic outcomes reported in recent large epidemiologic studies in Europe, the USA, and elsewhere. It also accounts for the prominent aging phenotype observed clinically in long-term cannabis use disorder and the 20 characteristics of aging that it manifests. Increasing daily cannabis use, increasing use in pregnancy, and exponential dose-response effects heighten the epidemiologic and clinical urgency of these findings. Together, these findings indicate that cannabinoid genotoxicity and epigenotoxicity are prominent features of cannabis dependence and strongly indicate coordinated multiomics investigations of cannabinoid genome-epigenome-transcriptome-metabolome, chromatin conformation, and 3D nuclear architecture. Considering the well-established exponential dose-response relationships, the diversity of cannabinoids, and the multigenerational nature of the implications, great caution is warranted in community cannabinoid penetration.

The brain's dark transcriptome: Sequencing RNA in distal compartments of neurons and glia

Transcriptomic approaches are powerful strategies to map the molecular diversity of cells in the brain. Single-cell genomic atlases have now been compiled for entire mammalian brains. However, complementary techniques are only just beginning to map the subcellular transcriptomes from distal cellular compartments. We review single-cell datasets alongside subtranscriptome data from the mammalian brain to explore the development of cellular and subcellular diversity. We discuss how single-cell RNA-seq misses transcripts localized away from cell bodies, which form the 'dark transcriptome' of the brain: a collection of subtranscriptomes in dendrites, axons, growth cones, synapses, and endfeet with important roles in brain development and function. Recent advances in subcellular transcriptome sequencing are beginning to reveal these elusive pools of RNA. We outline the success stories to date in uncovering the constituent subtranscriptomes of neurons and glia, as well as present the emerging toolkit that is accelerating the pace of subtranscriptome discovery.

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.

Single-cell DNA Methylome and 3D Multi-omic Atlas of the Adult Mouse Brain

Cytosine DNA methylation is essential in brain development and has been implicated in various neurological disorders. A comprehensive understanding of DNA methylation diversity across the entire brain in the context of the brain's 3D spatial organization is essential for building a complete molecular atlas of brain cell types and understanding their gene regulatory landscapes. To this end, we employed optimized single-nucleus methylome (snmC-seq3) and multi-omic (snm3C-seq1) sequencing technologies to generate 301,626 methylomes and 176,003 chromatin conformation/methylome joint profiles from 117 dissected regions throughout the adult mouse brain. Using iterative clustering and integrating with companion whole-brain transcriptome and chromatin accessibility datasets, we constructed a methylation-based cell type taxonomy that contains 4,673 cell groups and 261 cross-modality-annotated subclasses. We identified millions of differentially methylated regions (DMRs) across the genome, representing potential gene regulation elements. Notably, we observed spatial cytosine methylation patterns on both genes and regulatory elements in cell types within and across brain regions. Brain-wide multiplexed error-robust fluorescence in situ hybridization (MERFISH2) data validated the association of this spatial epigenetic diversity with transcription and allowed the mapping of the DNA methylation and topology information into anatomical structures more precisely than our dissections. Furthermore, multi-scale chromatin conformation diversities occur in important neuronal genes, highly associated with DNA methylation and transcription changes. Brain-wide cell type comparison allowed us to build a regulatory model for each gene, linking transcription factors, DMRs, chromatin contacts, and downstream genes to establish regulatory networks. Finally, intragenic DNA methylation and chromatin conformation patterns predicted alternative gene isoform expression observed in a companion whole-brain SMART-seq3 dataset. Our study establishes the first brain-wide, single-cell resolution DNA methylome and 3D multi-omic atlas, providing an unparalleled resource for comprehending the mouse brain's cellular-spatial and regulatory genome diversity.

A molecularly defined and spatially resolved cell atlas of the whole mouse brain

In mammalian brains, tens of millions to billions of cells form complex interaction networks to enable a wide range of functions. The enormous diversity and intricate organization of cells in the brain have so far hindered our understanding of the molecular and cellular basis of its functions. Recent advances in spatially resolved single-cell transcriptomics have allowed systematic mapping of the spatial organization of molecularly defined cell types in complex tissues1-3. However, these approaches have only been applied to a few brain regions1-11 and a comprehensive cell atlas of the whole brain is still missing. Here, we imaged a panel of >1,100 genes in ~8 million cells across the entire adult mouse brain using multiplexed error-robust fluorescence in situ hybridization (MERFISH)12 and performed spatially resolved, single-cell expression profiling at the whole-transcriptome scale by integrating MERFISH and single-cell RNA-sequencing (scRNA-seq) data. Using this approach, we generated a comprehensive cell atlas of >5,000 transcriptionally distinct cell clusters, belonging to ~300 major cell types, in the whole mouse brain with high molecular and spatial resolution. Registration of the MERFISH images to the common coordinate framework (CCF) of the mouse brain further allowed systematic quantifications of the cell composition and organization in individual brain regions defined in the CCF. We further identified spatial modules characterized by distinct cell-type compositions and spatial gradients featuring gradual changes in the gene-expression profiles of cells. Finally, this high-resolution spatial map of cells, with a transcriptome-wide expression profile associated with each cell, allowed us to infer cell-type-specific interactions between several hundred pairs of molecularly defined cell types and predict potential molecular (ligand-receptor) basis and functional implications of these cell-cell interactions. These results provide rich insights into the molecular and cellular architecture of the brain and a valuable resource for future functional investigations of neural circuits and their dysfunction in diseases.

Diversity in spatial frequency, temporal frequency, and speed tuning across mouse visual cortical areas and layers

The mouse visual system consists of several visual cortical areas thought to be specialized for different visual features and/or tasks. Previous studies have revealed differences between primary visual cortex (V1) and other higher visual areas, namely, anterolateral (AL) and posteromedial (PM), and their tuning preferences for spatial and temporal frequency. However, these differences have primarily been characterized using methods that are biased toward superficial layers of cortex, such as two-photon calcium imaging. Fewer studies have investigated cell types in deeper layers of these areas and their tuning preferences. Because superficial versus deep-layer neurons and different types of deep-layer neurons are known to have different feedforward and feedback inputs and outputs, comparing the tuning preferences of these groups is important for understanding cortical visual information processing. In this study, we used extracellular electrophysiology and two-photon calcium imaging targeted toward two different layer 5 cell classes to characterize their tuning properties in V1, AL, and PM. We find that deep-layer neurons, similar to superficial layer neurons, are also specialized for different spatial and temporal frequencies, with the strongest differences between AL and V1, and AL and PM, but not V1 and PM. However, we note that the deep-layer neuron populations preferred a larger range of SFs and TFs compared to previous studies. We also find that extratelencephalically projecting layer 5 neurons are more direction selective than intratelencephalically projecting layer 5 neurons.

Homeorhesis: envisaging the logic of life trajectories in molecular research on trauma and its effects

What sets someone on a life trajectory? This question is at the heart of studies of 21st-century neurosciences that build on scientific models developed over the last 150 years that attempt to link psychopathology risk and human development. Historically, this research has documented persistent effects of singular, negative life experiences on people's subsequent development. More recently, studies have documented neuromolecular effects of early life adversity on life trajectories, resulting in models that frame lives as disproportionately affected by early negative experiences. This view is dominant, despite little evidence of the stability of the presumably early-developed molecular traits and their potential effects on phenotypes. We argue that in the context of gaps in knowledge and the need for scientists to reason across molecular and phenotypic scales, as well as time spans that can extend beyond an individual's life, specific interpretative frameworks shape the ways in which individual scientific findings are assessed. In the process, scientific reasoning oscillates between understandings of cellular homeostasis and organisms' homeorhesis, or life trajectory. Biologist and historian François Jacob described this framework as the "attitude" that researchers bring to bear on their "objects" of study. Through an analysis of, first, historical and contemporary scientific literature and then ethnographic research with neuroscientists, we consider how early life trauma came to be associated with specific psychological and neurobiological effects grounded in understandings of life trajectories. We conclude with a consideration of the conceptual, ontological, and ethical implications of interpreting life trajectories as the result of the persistence of long-embodied biological traits, persistent life environments, or both.

Identification of visual cortex cell types and species differences using single-cell RNA sequencing

The primate neocortex exerts high cognitive ability and strong information processing capacity. Here, we establish a single-cell RNA sequencing dataset of 133,454 macaque visual cortical cells. It covers major cortical cell classes including 25 excitatory neuron types, 37 inhibitory neuron types and all glial cell types. We identified layer-specific markers including HPCAL1 and NXPH4, and also identified two cell types, an NPY-expressing excitatory neuron type that expresses the dopamine receptor D3 gene; and a primate specific activity-dependent OSTN + sensory neuron type. Comparisons of our dataset with humans and mice show that the gene expression profiles differ between species in relation to genes that are implicated in the synaptic plasticity and neuromodulation of excitatory neurons. The comparisons also revealed that glutamatergic neurons may be more diverse across species than GABAergic neurons and non-neuronal cells. These findings pave the way for understanding how the primary cortex fulfills the high-cognitive functions.

Temporally divergent regulatory mechanisms govern neuronal diversification and maturation in the mouse and marmoset neocortex

Mammalian neocortical neurons span one of the most diverse cell type spectra of any tissue. Cortical neurons are born during embryonic development, and their maturation extends into postnatal life. The regulatory strategies underlying progressive neuronal development and maturation remain unclear. Here we present an integrated single-cell epigenomic and transcriptional analysis of individual mouse and marmoset cortical neuron classes, spanning both early postmitotic stages of identity acquisition and later stages of neuronal plasticity and circuit integration. We found that, in both species, the regulatory strategies controlling early and late stages of pan-neuronal development diverge. Early postmitotic neurons use more widely shared and evolutionarily conserved molecular regulatory programs. In contrast, programs active during later neuronal maturation are more brain- and neuron-specific and more evolutionarily divergent. Our work uncovers a temporal shift in regulatory choices during neuronal diversification and maturation in both mice and marmosets, which likely reflects unique evolutionary constraints on distinct events of neuronal development in the neocortex.

What is a cell type and how to define it?

Cell types are the basic functional units of an organism. Cell types exhibit diverse phenotypic properties at multiple levels, making them challenging to define, categorize, and understand. This review provides an overview of the basic principles of cell types rooted in evolution and development and discusses approaches to characterize and classify cell types and investigate how they contribute to the organism's function, using the mammalian brain as a primary example. I propose a roadmap toward a conceptual framework and knowledge base of cell types that will enable a deeper understanding of the dynamic changes of cellular function under healthy and diseased conditions.

Exploring epigenomic mechanisms of neural tube defects using multi-omics methods and data

Neural tube defects (NTDs) are a heterogeneous set of malformations attributed to disruption in normal neural tube closure during early embryogenesis. An in-depth understanding of NTD etiology and mechanisms remains elusive, however. Among the proposed mechanisms, epigenetic changes are thought to play an important role in the formation of NTDs. Epigenomics covers a wide spectrum of genomic DNA sequence modifications that can be investigated via high-throughput techniques. Recent advances in epigenomic technologies have enabled epigenetic studies of congenital malformations and facilitated the integration of big data into the understanding of NTDs. Herein, we review clinical epigenomic data that focuses on DNA methylation, histone modification, and miRNA alterations in human neural tissues, placental tissues, and leukocytes to explore potential mechanisms by which candidate genes affect human NTD pathogenesis. We discuss the links between epigenomics and gene regulatory mechanisms, and the effects of epigenetic alterations in human tissues on neural tube closure.

A novel DNA methylation-driver gene signature for long-term survival prediction of hepatitis-positive hepatocellular carcinoma patients

BACKGROUND

Abnormal DNA methylation is one of the most general epigenetic modifications in hepatocellular carcinoma (HCC). Recent research showed that DNA methylation was a prognostic indicator of all-cause HCC and nonviral HCC. However, whether DNA methylation-driver genes could be used for predicting survival, the probability of hepatitis-positive HCC remains unclear.

METHODS

In this study, DNA methylation-driver genes (MDGs) were screened by a joint analysis of methylome and transcriptome data of 142 hepatitis-positive HCC patients. Subsequently, a prognostic risk score and nomogram were constructed. Finally, correlation analyses between the risk score and signaling pathways and immunity were conducted by GSVA and CIBERSORT.

RESULTS

Through random forest screening and Cox progression analysis, 10 prognostic methylation-driver genes (AC008271.1, C11orf53, CASP8, F2RL2, GBP5, LUCAT1, RP11-114B7.6, RP11-149I23.3, RP11-383 J24.1, and SLC35G2) were screened out. As a result, a prognostic risk score signature was constructed. The independent value of the risk score for prognosis prediction were addressed in the TCGA-HCC and the China-HCC cohorts. Next, clinicopathological features were analyzed and HBV status and histological grade were screened to construct a nomogram together with the risk score. The prognostic efficiency of the nomogram was validated by the calibration curves and the concordance index (C index: 0.829, 95% confidence interval: 0.794-0.864), while its clinical application ability was confirmed by decision curve analysis (DCA). At last, the relationship between the risk score and signaling pathways, as well as the correlations between immune cells were elucidated preliminary.

CONCLUSIONS

Taken together, our study explored a novel DNA methylation-driver gene risk score signature and an efficient nomogram for long-term survival prediction of hepatitis-positive HCC patients.

Glucocorticoid Receptor-Regulated Enhancers Play a Central Role in the Gene Regulatory Networks Underlying Drug Addiction

Substance abuse and addiction represent a significant public health problem that impacts multiple dimensions of society, including healthcare, the economy, and the workforce. In 2021, over 100,000 drug overdose deaths were reported in the US, with an alarming increase in fatalities related to opioids and psychostimulants. Understanding the fundamental gene regulatory mechanisms underlying addiction and related behaviors could facilitate more effective treatments. To explore how repeated drug exposure alters gene regulatory networks in the brain, we combined capped small (cs)RNA-seq, which accurately captures nascent-like initiating transcripts from total RNA, with Hi-C and single nuclei (sn)ATAC-seq. We profiled initiating transcripts in two addiction-related brain regions, the prefrontal cortex (PFC) and the nucleus accumbens (NAc), from rats that were never exposed to drugs or were subjected to prolonged abstinence after oxycodone or cocaine intravenous self-administration (IVSA). Interrogating over 100,000 active transcription start regions (TSRs) revealed that most TSRs had hallmarks of bonafide enhancers and highlighted the KLF/SP1, RFX, and AP1 transcription factors families as central to establishing brain-specific gene regulatory programs. Analysis of rats with addiction-like behaviors versus controls identified addiction-associated repression of transcription at regulatory enhancers recognized by nuclear receptor subfamily 3 group C (NR3C) factors, including glucocorticoid receptors. Cell-type deconvolution analysis using snATAC-seq uncovered a potential role of glial cells in driving the gene regulatory programs associated with addiction-related phenotypes. These findings highlight the power of advanced transcriptomics methods to provide insight into how addiction perturbs gene regulatory programs in the brain.

Brain Networks of Connectionally Unique Basolateral Amygdala Cell Types

Different brain regions structurally interconnected through networks regulate behavior output. Therefore, understanding the functional organization of the brain in health and disease necessitates a foundational anatomic roadmap to its network organization. To provide this to the research community, our lab has systematically traced thousands of pathways in the mouse brain and has applied computational measures to determine the network architecture of major brain systems. Toward this effort, the brain-wide networks of the basolateral amygdalar complex (BLA) were recently generated. The data revealed uniquely connected cell types within the same BLA nucleus that were constituents of distinct neural networks. Here, we elaborate on how these connectionally unique BLA cell types fit within the larger cortico-basal ganglia and limbic networks that were previously described by our team. The significance and utility of high quality, detailed anatomic data is also discussed.

Precise Mapping of Otp Expressing Cells Across Different Pallial Regions Throughout Ontogenesis Using Otp-Specific Reporter Transgenic Mice

Taking advantage of two Otp-specific reporter lines of transgenic mice (Otp-eGFP and Otp-Cre; Rpl22-HA), we identify and describe different Otp cell populations across various pallial regions, including the pallial amygdala, the piriform cortex, the mesocortex, the neocortex, and the hippocampal complex. Some of these populations can be followed throughout development, suggesting migration from external sources (for example, those of the pallial amygdala and at least some of the cingulate cortex). Other cells become visible during postnatal development (some of those in the neocortex and hippocampal formation) or in adulthood (those of the parahippocampal lobe), and seem to be produced locally. We discuss the possible role of Otp in these different populations during different moments of ontogenesis. We also analyze the connectivity patterns of some of these cells and discuss their functional implications. For example, our data suggest that Otp cells of the pallial amygdala might be engaged in networks with other Otp cells of the medial amygdala with the same embryonic origin, and may regulate specific aspects of social behavior. Regarding Otp cells in the parahippocampal lobe, they seem to be projection neurons and may regulate hippocampal function during spatial navigation and memory formation. The two reporter transgenic mice employed here provide very powerful tools for high precision studies on these different Otp cells of the pallium, but careful attention should be paid to the age and to differences between lines.

Congenital anomaly epidemiological correlates of Δ8THC across USA 2003-16: panel regression and causal inferential study

Δ8-Tetrahydrocannabinol (Δ8THC) is marketed in many US states as 'legal weed'. Concerns exist relating to class-wide genotoxic cannabinoid effects. We conducted an epidemiological investigation of Δ8THC-related genotoxicity expressed as 57 congenital anomaly (CA) rates (CARs) in the USA. CARs were taken from the Centers for Disease Control, Atlanta, Georgia. Drug exposure data were taken from the National Survey of Drug Use and Health, with a response rate of 74.1%. Ethnicity and income data were taken from the US Census Bureau. National cannabinoid exposure was taken from Drug Enforcement Agency publications and multiplied by state cannabis use data to derive state-based estimates of Δ8THC exposure. At bivariate continuous analysis, Δ8THC was associated with 23 CAs on raw CA rates, 33 CARs after correction for early termination for anomaly estimates and 41 on a categorical analysis comparing the highest and lowest exposure quintiles. At inverse probability weighted multivariable additive and interactive models lagged to 0, 2 and 4 years, Δ8THC was linked with 39, 8, 4 and 9 CAs. Chromosomal, cardiovascular, gastrointestinal, genitourinary, limb, central nervous system (CNS) and face systems were particularly affected. The minimum E-values ranged to infinity. Both the number of anomalies implicated and the effect sizes demonstrated were much greater for Δ8THC than for tobacco and alcohol combined. Δ8THC appears epidemiologically to be more strongly associated with many CAs than for tobacco and alcohol and is consistent with a cannabinoid class genotoxic/epigenotoxic effect. Quantitative causality criteria were fulfilled, and causal relationships either for Δ8THC or for cannabinoid/s, for which it is a surrogate marker, may be in operation.

Cortical Cartography: Mapping Arealization Using Single-Cell Omics Technology

The cerebral cortex derives its cognitive power from a modular network of specialized areas processing a multitude of information. The assembly and organization of these regions is vital for human behavior and perception, as evidenced by the prevalence of area-specific phenotypes that manifest in neurodevelopmental and psychiatric disorders. Generations of scientists have examined the architecture of the human cortex, but efforts to capture the gene networks which drive arealization have been hampered by the lack of tractable models of human neurodevelopment. Advancements in "omics" technologies, imaging, and computational power have enabled exciting breakthroughs into the molecular and structural characteristics of cortical areas, including transcriptomic, epigenomic, metabolomic, and proteomic profiles of mammalian models. Here we review the single-omics atlases that have shaped our current understanding of cortical areas, and their potential to fuel a new era of multi-omic single-cell endeavors to interrogate both the developing and adult human cortex.

Morphological diversity of single neurons in molecularly defined cell types

Dendritic and axonal morphology reflects the input and output of neurons and is a defining feature of neuronal types1,2, yet our knowledge of its diversity remains limited. Here, to systematically examine complete single-neuron morphologies on a brain-wide scale, we established a pipeline encompassing sparse labelling, whole-brain imaging, reconstruction, registration and analysis. We fully reconstructed 1,741 neurons from cortex, claustrum, thalamus, striatum and other brain regions in mice. We identified 11 major projection neuron types with distinct morphological features and corresponding transcriptomic identities. Extensive projectional diversity was found within each of these major types, on the basis of which some types were clustered into more refined subtypes. This diversity follows a set of generalizable principles that govern long-range axonal projections at different levels, including molecular correspondence, divergent or convergent projection, axon termination pattern, regional specificity, topography, and individual cell variability. Although clear concordance with transcriptomic profiles is evident at the level of major projection type, fine-grained morphological diversity often does not readily correlate with transcriptomic subtypes derived from unsupervised clustering, highlighting the need for single-cell cross-modality studies. Overall, our study demonstrates the crucial need for quantitative description of complete single-cell anatomy in cell-type classification, as single-cell morphological diversity reveals a plethora of ways in which different cell types and their individual members may contribute to the configuration and function of their respective circuits.

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.

DNA methylation atlas of the mouse brain at single-cell resolution

Mammalian brain cells show remarkable diversity in gene expression, anatomy and function, yet the regulatory DNA landscape underlying this extensive heterogeneity is poorly understood. Here we carry out a comprehensive assessment of the epigenomes of mouse brain cell types by applying single-nucleus DNA methylation sequencing1,2 to profile 103,982 nuclei (including 95,815 neurons and 8,167 non-neuronal cells) from 45 regions of the mouse cortex, hippocampus, striatum, pallidum and olfactory areas. We identified 161 cell clusters with distinct spatial locations and projection targets. We constructed taxonomies of these epigenetic types, annotated with signature genes, regulatory elements and transcription factors. These features indicate the potential regulatory landscape supporting the assignment of putative cell types and reveal repetitive usage of regulators in excitatory and inhibitory cells for determining subtypes. The DNA methylation landscape of excitatory neurons in the cortex and hippocampus varied continuously along spatial gradients. Using this deep dataset, we constructed an artificial neural network model that precisely predicts single neuron cell-type identity and brain area spatial location. Integration of high-resolution DNA methylomes with single-nucleus chromatin accessibility data3 enabled prediction of high-confidence enhancer-gene interactions for all identified cell types, which were subsequently validated by cell-type-specific chromatin conformation capture experiments4. By combining multi-omic datasets (DNA methylation, chromatin contacts, and open chromatin) from single nuclei and annotating the regulatory genome of hundreds of cell types in the mouse brain, our DNA methylation atlas establishes the epigenetic basis for neuronal diversity and spatial organization throughout the mouse cerebrum.