A data-driven single-cell and spatial transcriptomic map of the human prefrontal cortex.

The molecular organization of the human neocortex historically has been studied in the context of its histological layers. However, emerging spatial transcriptomic technologies have enabled unbiased identification of transcriptionally defined spatial domains that move beyond classic cytoarchitecture. We used the Visium spatial gene expression platform to generate a data-driven molecular neuroanatomical atlas across the anterior-posterior axis of the human dorsolateral prefrontal cortex. Integration with paired single-nucleus RNA-sequencing data revealed distinct cell type compositions and cell-cell interactions across spatial domains. Using PsychENCODE and publicly available data, we mapped the enrichment of cell types and genes associated with neuropsychiatric disorders to discrete spatial domains.

The gene expression landscape of the human locus coeruleus revealed by single-nucleus and spatially-resolved transcriptomics.

Norepinephrine (NE) neurons in the locus coeruleus (LC) make long-range projections throughout the central nervous system, playing critical roles in arousal and mood, as well as various components of cognition including attention, learning, and memory. The LC-NE system is also implicated in multiple neurological and neuropsychiatric disorders. Importantly, LC-NE neurons are highly sensitive to degeneration in both Alzheimer's and Parkinson's disease. Despite the clinical importance of the brain region and the prominent role of LC-NE neurons in a variety of brain and behavioral functions, a detailed molecular characterization of the LC is lacking. Here, we used a combination of spatially-resolved transcriptomics and single-nucleus RNA-sequencing to characterize the molecular landscape of the LC region and the transcriptomic profile of LC-NE neurons in the human brain. We provide a freely accessible resource of these data in web-accessible and downloadable formats.

High frequency of HIV precursor-target-specific B cells in sub-Saharan populations.

HIV gp120 engineered outer domain germline-targeting version 8 (eOD-GT8) was designed specifically to engage naive B cell precursors of VRC01-class antibodies. However, the frequency and affinity of naive B cell precursors able to recognize eOD-GT8 have been evaluated only in U.S. populations. HIV infection is disproportionally concentrated in sub-Saharan Africa, so we seek to characterize naive B cells able to recognize eOD-GT8 in sub-Saharan cohorts. We demonstrate that people from sub-Saharan Africa have a higher or equivalent frequency of naive B cells able to engage eOD-GT8 compared with people from the U.S. Genetically, the higher frequency of eOD-GT8-positive cells is accompanied by a higher level of naive B cells with gene signatures characteristic of the VRC01 class, as well as other CD4bs-directed antibodies. Our study demonstrates that vaccination with eOD-GT8 in sub-Saharan Africa could be successful at expanding and establishing a pool of CD4bs-directed memory B cells from naive precursors.

Integrated single cell and unsupervised spatial transcriptomic analysis defines molecular anatomy of the human dorsolateral prefrontal cortex.

Generation of a molecular neuroanatomical map of the human prefrontal cortex reveals novel spatial domains and cell-cell interactions relevant for psychiatric disease. The molecular organization of the human neocortex has been historically studied in the context of its histological layers. However, emerging spatial transcriptomic technologies have enabled unbiased identification of transcriptionally-defined spatial domains that move beyond classic cytoarchitecture. Here we used the Visium spatial gene expression platform to generate a data-driven molecular neuroanatomical atlas across the anterior-posterior axis of the human dorsolateral prefrontal cortex (DLPFC). Integration with paired single nucleus RNA-sequencing data revealed distinct cell type compositions and cell-cell interactions across spatial domains. Using PsychENCODE and publicly available data, we map the enrichment of cell types and genes associated with neuropsychiatric disorders to discrete spatial domains. Finally, we provide resources for the scientific community to explore these integrated spatial and single cell datasets at research.libd.org/spatialDLPFC/.

VistoSeg: Processing utilities for high-resolution images for spatially resolved transcriptomics data.

Spatially resolved transcriptomics (SRT) is a growing field that links gene expression to anatomical context. SRT approaches that use next-generation sequencing (NGS) combine RNA sequencing with histological or fluorescent imaging to generate spatial maps of gene expression in intact tissue sections. These technologies directly couple gene expression measurements with high-resolution histological or immunofluorescent images that contain rich morphological information about the tissue under study. While broad access to NGS-based spatial transcriptomic technology is now commercially available through the Visium platform from the vendor 10× Genomics, computational tools for extracting image-derived metrics for integration with gene expression data remain limited. We developed VistoSeg as a MATLAB pipeline to process, analyze and interactively visualize the high-resolution images generated in the Visium platform. VistoSeg outputs can be easily integrated with accompanying transcriptomic data to facilitate downstream analyses in common programing languages including R and Python. VistoSeg provides user-friendly tools for integrating image-derived metrics from histological and immunofluorescent images with spatially resolved gene expression data. Integration of this data enhances the ability to understand the transcriptional landscape within tissue architecture. VistoSeg is freely available at http://research.libd.org/VistoSeg/.

spatialLIBD: an R/Bioconductor package to visualize spatially-resolved transcriptomics data.

BACKGROUND: Spatially-resolved transcriptomics has now enabled the quantification of high-throughput and transcriptome-wide gene expression in intact tissue while also retaining the spatial coordinates. Incorporating the precise spatial mapping of gene activity advances our understanding of intact tissue-specific biological processes. In order to interpret these novel spatial data types, interactive visualization tools are necessary. RESULTS: We describe spatialLIBD, an R/Bioconductor package to interactively explore spatially-resolved transcriptomics data generated with the 10x Genomics Visium platform. The package contains functions to interactively access, visualize, and inspect the observed spatial gene expression data and data-driven clusters identified with supervised or unsupervised analyses, either on the user's computer or through a web application. CONCLUSIONS: spatialLIBD is available at https://bioconductor.org/packages/spatialLIBD . It is fully compatible with SpatialExperiment and the Bioconductor ecosystem. Its functionality facilitates analyzing and interactively exploring spatially-resolved data from the Visium platform.

Reconstruction of dynamic regulatory networks reveals signaling-induced topology changes associated with germ layer specification.

Elucidating regulatory relationships between transcription factors (TFs) and target genes is fundamental to understanding how cells control their identity and behavior. Unfortunately, existing computational gene regulatory network (GRN) reconstruction methods are imprecise, computationally burdensome, and fail to reveal dynamic regulatory topologies. Here, we present Epoch, a reconstruction tool that uses single-cell transcriptomics to accurately infer dynamic networks. We apply Epoch to identify the dynamic networks underpinning directed differentiation of mouse embryonic stem cells (ESCs) guided by multiple signaling pathways, and we demonstrate that modulating these pathways drives topological changes that bias cell fate potential. We also find that Peg3 rewires the pluripotency network to favor mesoderm specification. By integrating signaling pathways with GRNs, we trace how Wnt activation and PI3K suppression govern mesoderm and endoderm specification, respectively. Finally, we identify regulatory circuits of patterning and axis formation that distinguish in vitro and in vivo mesoderm specification.

Single-nucleus transcriptome analysis reveals cell-type-specific molecular signatures across reward circuitry in the human brain.

Single-cell gene expression technologies are powerful tools to study cell types in the human brain, but efforts have largely focused on cortical brain regions. We therefore created a single-nucleus RNA-sequencing resource of 70,615 high-quality nuclei to generate a molecular taxonomy of cell types across five human brain regions that serve as key nodes of the human brain reward circuitry: nucleus accumbens, amygdala, subgenual anterior cingulate cortex, hippocampus, and dorsolateral prefrontal cortex. We first identified novel subpopulations of interneurons and medium spiny neurons (MSNs) in the nucleus accumbens and further characterized robust GABAergic inhibitory cell populations in the amygdala. Joint analyses across the 107 reported cell classes revealed cell-type substructure and unique patterns of transcriptomic dynamics. We identified discrete subpopulations of D1- and D2-expressing MSNs in the nucleus accumbens to which we mapped cell-type-specific enrichment for genetic risk associated with both psychiatric disease and addiction.

A single cell transcriptional portrait of embryoid body differentiation and comparison to progenitors of the developing embryo.

Directed differentiation of pluripotent stem cells provides an accessible system to model development. However, the distinct cell types that emerge, their dynamics, and their relationship to progenitors in the early embryo has been difficult to decipher because of the cellular heterogeneity inherent to differentiation. Here, we used a combination of bulk RNA-Seq, single cell RNA-Seq, and bioinformatics analyses to dissect the cell types that emerge during directed differentiation of mouse embryonic stem cells as embryoid bodies and we compared them to spatially and temporally resolved transcriptional profiles of early embryos. Our single cell analyses of the day 4 embryoid bodies revealed three populations which had retained related yet distinct pluripotent signatures that resemble the pre- or post-implantation epiblast, one population of presumptive neuroectoderm, one population of mesendoderm, and four populations of neural progenitors. By day 6, the neural progenitors predominated the embryoid bodies, but both a small population of pluripotent-like cells and an anterior mesoderm-like Brachyury-expressing population were present. By comparing the day 4 and day 6 populations, we identified candidate differentiation paths, transcription factors, and signaling pathways that mark the in vitro correlate of the transition from the mid-to-late primitive streak stage.

Impacts of adelphophagic development on variation in offspring size, duration of development, and temperature-mediated plasticity.

Adelphophagic development, where embryos consume sibling embryos or nurse eggs, is particularly common in marine caenogastropods and some families of polychaetes. When exogenous nutrition is provided before hatching, egg size and hatching size can be uncoupled, but advantages and constraints of adelphophagic development compared to development from large eggs are unknown. Here we examine temperature-mediated plasticity in offspring size, brooding duration, and fecundity in the adelphophagic marine gastropod Crepidula cf. onyx. We use these data combined with previously published data on two planktotrophic Crepidula and two Crepidula species that develop from large eggs to test hypotheses about the consequences of adelphophagic development and patterns of variation in offspring size. In Crepidula cf. onyx, egg size shows no significant effect of temperature. Hatching size is significantly larger at 28 °C than at 23 °C but proceeds from fewer eggs per capsule at 28 °C. Hatching size is therefore decoupled from both egg size and the number of eggs per capsule. Although development is faster at the higher temperature, broods are produced roughly every 26-27 days at both temperatures. Increased rate of development has been cited as a potential advantage of adelphophagic development in muricids, but the adelphophagic C. cf. onyx did not develop more quickly than C. atrasolea or C. ustulatulina, species that produce similarly sized hatchlings from large eggs. Comparisons across calyptraeid species support the role of adelphophagy in increasing variance in offspring size. This increased variability is primarily expressed within broods or among broods from the same female, not among females.