Coordination of olfactory receptor choice with guidance receptor expression and function in olfactory sensory neurons.

Olfactory sensory neurons choose to express a single odorant receptor (OR) from a large gene repertoire and extend axons to reproducible, OR-specific locations within the olfactory bulb. This developmental process produces a topographically organized map of odorant experience in the brain. The axon guidance mechanisms that generate this pattern of connectivity, as well as those that coordinate OR choice and axonal guidance receptor expression, are incompletely understood. We applied the powerful approach of single-cell RNA-seq on newly born olfactory sensory neurons (OSNs) in young zebrafish larvae to address these issues. Expression profiles were generated for 56 individual Olfactory Marker Protein (OMP) positive sensory neurons by single-cell (SC) RNA-seq. We show that just as in mouse OSNs, mature zebrafish OSNs typically express a single predominant OR transcript. Our previous work suggests that OSN targeting is related to the OR clade from which a sensory neuron chooses to express its odorant receptor. We categorized each of the mature cells based on the clade of their predominantly expressed OR. Transcripts expressed at higher levels in each of three clade-related categories were identified using Penalized Linear Discriminant Analysis (PLDA). A genome-wide approach was used to identify membrane-associated proteins that are most likely to have guidance-related activity. We found that OSNs that choose to express an OR from a particular clade also express specific subsets of potential axon guidance genes and transcription factors. We validated our identification of candidate axon guidance genes for one clade of OSNs using bulk RNA-seq from a subset of transgene-labeled neurons that project to a single protoglomerulus. The differential expression patterns of selected candidate guidance genes were confirmed using fluorescent in situ hybridization. Most importantly, we observed axonal mistargeting in knockouts of three candidate axonal guidance genes identified in this analysis: nrp1a, nrp1b, and robo2. In each case, targeting errors were detected in the subset of axons that normally express these transcripts at high levels, and not in the axons that express them at low levels. Our findings demonstrate that specific, functional, axonal guidance related genes are expressed in subsets of OSNs that that can be categorized by their patterns of OR expression.

PIVOT: platform for interactive analysis and visualization of transcriptomics data.

BACKGROUND: Many R packages have been developed for transcriptome analysis but their use often requires familiarity with R and integrating results of different packages requires scripts to wrangle the datatypes. Furthermore, exploratory data analyses often generate multiple derived datasets such as data subsets or data transformations, which can be difficult to track. RESULTS: Here we present PIVOT, an R-based platform that wraps open source transcriptome analysis packages with a uniform user interface and graphical data management that allows non-programmers to interactively explore transcriptomics data. PIVOT supports more than 40 popular open source packages for transcriptome analysis and provides an extensive set of tools for statistical data manipulations. A graph-based visual interface is used to represent the links between derived datasets, allowing easy tracking of data versions. PIVOT further supports automatic report generation, publication-quality plots, and program/data state saving, such that all analysis can be saved, shared and reproduced. CONCLUSIONS: PIVOT will allow researchers with broad background to easily access sophisticated transcriptome analysis tools and interactively explore transcriptome datasets.

Single-cell transcriptomics and functional target validation of brown adipocytes show their complex roles in metabolic homeostasis.

Brown adipocytes (BAs) are specialized for adaptive thermogenesis and, upon sympathetic stimulation, activate mitochondrial uncoupling protein (UCP)-1 and oxidize fatty acids to generate heat. The capacity for brown adipose tissue (BAT) to protect against obesity and metabolic disease is recognized, yet information about which signals activate BA, besides ß3-adrenergic receptor stimulation, is limited. Using single-cell transcriptomics, we confirmed the presence of mRNAs encoding traditional BAT markers (i.e., UCP1, expressed in 100% of BAs Adrb3, expressed in <50% of BAs) in mouse and have shown single-cell variability (>1000-fold) in their expression at both the mRNA and protein levels. We further identified mRNAs encoding novel markers, orphan GPCRs, and many receptors that bind the classic neurotransmitters, neuropeptides, chemokines, cytokines, and hormones. The transcriptome variability between BAs suggests a much larger range of responsiveness of BAT than previously recognized and that not all BAs function identically. We examined the in vivo functional expression of 12 selected receptors by microinjecting agonists into live mouse BAT and analyzing the metabolic response. In this manner, we expanded the number of known receptors on BAs at least 25-fold, while showing that the expression of classic BA markers is more complex and variable than previously thought.

Assessing characteristics of RNA amplification methods for single cell RNA sequencing.

BACKGROUND: Recently, measurement of RNA at single cell resolution has yielded surprising insights. Methods for single-cell RNA sequencing (scRNA-seq) have received considerable attention, but the broad reliability of single cell methods and the factors governing their performance are still poorly known. RESULTS: Here, we conducted a large-scale control experiment to assess the transfer function of three scRNA-seq methods and factors modulating the function. All three methods detected greater than 70% of the expected number of genes and had a 50% probability of detecting genes with abundance greater than 2 to 4 molecules. Despite the small number of molecules, sequencing depth significantly affected gene detection. While biases in detection and quantification were qualitatively similar across methods, the degree of bias differed, consistent with differences in molecular protocol. Measurement reliability increased with expression level for all methods and we conservatively estimate measurements to be quantitative at an expression level greater than ~5-10 molecules. CONCLUSIONS: Based on these extensive control studies, we propose that RNA-seq of single cells has come of age, yielding quantitative biological information.

Serotonergic neuron regulation informed by in vivo single-cell transcriptomics.

Despite the recognized importance of the dorsal raphe (DR) serotonergic (5-HT) nuclei in the pathophysiology of depression and anxiety, the molecular components/putative drug targets expressed by these neurons are poorly characterized. Utilizing the promoter of an ETS domain transcription factor that is a stable marker of 5-HT neurons (Pet-1) to drive 5-HT neuronal expression of YFP, we identified 5-HT neurons in live acute slices. We isolated RNA from single 5-HT neurons in the ventromedial and lateral wings of the DR and performed single-cell RNA-Seq analysis identifying >500 G-protein coupled receptors (GPCRs) including receptors for classical transmitters, lipid signals, and peptides as well as dozens of orphan-GPCRs. Using these data to inform our selection of receptors to assess, we found that oxytocin and lysophosphatidic acid 1 receptors are translated and active in costimulating, with the α1-adrenergic receptor, the firing of DR 5-HT neurons, while the effects of histamine are inhibitory and exerted at H3 histamine receptors. The inhibitory histamine response provides evidence for tonic in vivo histamine inhibition of 5-HT neurons. This study illustrates that unbiased single-cell transcriptomics coupled with functional analyses provides novel insights into how neurons and neuronal systems are regulated.

Systematic Sampling of the Female Reproductive System for Molecular Characterization.

As part of the National Institutes of Health Human BioMolecular Atlas Program to develop a global platform to map the 37 trillion cells in the adult human body, we are generating a comprehensive molecular characterization of the female reproductive system. Data gathered from multiple single-cell/single-nucleus and spatial molecular assays will be used to build a 3D molecular atlas. Herein, we describe our multistep protocol, beginning with an optimized organ procurement workflow that maintains functional characteristics of the uterus, ovaries, and fallopian tubes by perfusing these organs with preservation solution. We have also developed a structured tissue sampling procedure that retains information on individual-level anatomic, physiologic, and individual diversity of the female reproductive system, toward full exploration of the function and structure of female reproductive cells. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Preparation and preservation of the female reproductive system (ovaries, fallopian tubes, and uterus) prior to procurement Basic Protocol 2: Removal of the female reproductive system en bloc Basic Protocol 3: Postsurgical dissection of ovaries Basic Protocol 4: Postsurgical dissection of fallopian tubes Basic Protocol 5: Postsurgical dissection of cervix Basic Protocol 6: Postsurgical dissection of uterine body Support Protocol 1: OCT-embedded tissue protocol Support Protocol 2: Tissue fixation protocol Support Protocol 3: Snap-frozen tissue protocol Basic Protocol 7: Tissue slice preparation for Visium analysis Support Protocol 4: Hematoxylin and eosin staining for 10X Visium imaging Basic Protocol 8: Manual tissue dissociation for Multiome analysis Basic Protocol 9: Tissue dissociation for Multiome analysis using S2 Singulator.

CHEX-seq detects single-cell genomic single-stranded DNA with catalytical potential.

Genomic DNA (gDNA) undergoes structural interconversion between single- and double-stranded states during transcription, DNA repair and replication, which is critical for cellular homeostasis. We describe "CHEX-seq" which identifies the single-stranded DNA (ssDNA) in situ in individual cells. CHEX-seq uses 3'-terminal blocked, light-activatable probes to prime the copying of ssDNA into complementary DNA that is sequenced, thereby reporting the genome-wide single-stranded chromatin landscape. CHEX-seq is benchmarked in human K562 cells, and its utilities are demonstrated in cultures of mouse and human brain cells as well as immunostained spatially localized neurons in brain sections. The amount of ssDNA is dynamically regulated in response to perturbation. CHEX-seq also identifies single-stranded regions of mitochondrial DNA in single cells. Surprisingly, CHEX-seq identifies single-stranded loci in mouse and human gDNA that catalyze porphyrin metalation in vitro, suggesting a catalytic activity for genomic ssDNA. We posit that endogenous DNA enzymatic activity is a function of genomic ssDNA.

Primary Cell Culture of Live Neurosurgically Resected Aged Adult Human Brain Cells and Single Cell Transcriptomics.

Investigation of human CNS disease and drug effects has been hampered by the lack of a system that enables single-cell analysis of live adult patient brain cells. We developed a culturing system, based on a papain-aided procedure, for resected adult human brain tissue removed during neurosurgery. We performed single-cell transcriptomics on over 300 cells, permitting identification of oligodendrocytes, microglia, neurons, endothelial cells, and astrocytes after 3 weeks in culture. Using deep sequencing, we detected over 12,000 expressed genes, including hundreds of cell-type-enriched mRNAs, lncRNAs and pri-miRNAs. We describe cell-type- and patient-specific transcriptional hierarchies. Single-cell transcriptomics on cultured live adult patient derived cells is a prime example of the promise of personalized precision medicine. Because these cells derive from subjects ranging in age into their sixties, this system permits human aging studies previously possible only in rodent systems.

Deep sequencing reveals cell-type-specific patterns of single-cell transcriptome variation.

BACKGROUND: Differentiation of metazoan cells requires execution of different gene expression programs but recent single-cell transcriptome profiling has revealed considerable variation within cells of seeming identical phenotype. This brings into question the relationship between transcriptome states and cell phenotypes. Additionally, single-cell transcriptomics presents unique analysis challenges that need to be addressed to answer this question. RESULTS: We present high quality deep read-depth single-cell RNA sequencing for 91 cells from five mouse tissues and 18 cells from two rat tissues, along with 30 control samples of bulk RNA diluted to single-cell levels. We find that transcriptomes differ globally across tissues with regard to the number of genes expressed, the average expression patterns, and within-cell-type variation patterns. We develop methods to filter genes for reliable quantification and to calibrate biological variation. All cell types include genes with high variability in expression, in a tissue-specific manner. We also find evidence that single-cell variability of neuronal genes in mice is correlated with that in rats consistent with the hypothesis that levels of variation may be conserved. CONCLUSIONS: Single-cell RNA-sequencing data provide a unique view of transcriptome function; however, careful analysis is required in order to use single-cell RNA-sequencing measurements for this purpose. Technical variation must be considered in single-cell RNA-sequencing studies of expression variation. For a subset of genes, biological variability within each cell type appears to be regulated in order to perform dynamic functions, rather than solely molecular noise.

Lifestyle Modification in Blood Pressure Study II (LIMBS): study protocol of a randomized controlled trial assessing the efficacy of a 24 week structured yoga program versus lifestyle modification on blood pressure reduction.

Hypertension is a major public health issue affecting 68 million adults in the United States. Lifestyle modifications including complementary therapies such as the movement based mind body practice of yoga have become increasingly popular in the United States and have been considered as a potential alternative to medication in blood pressure reduction. We completed a pilot study in 2009 which showed meaningful decreases in 24-hour ambulatory blood pressure readings after a 12 week period of yoga participation. Based on data from our pilot study we are now completing The Lifestyle Modification and Blood Pressure Study (LIMBS II) which is a phase 2 randomized controlled trial designed to determine the effects of yoga therapy and enhanced lifestyle modification on lowering blood pressure in pre-hypertensive and stage 1 hypertensive subjects. Using 24-hour ambulatory blood pressure monitoring, LIMBS II aims to compare the effects on blood pressure reduction in subjects randomized for 24 weeks to one of the three following groups: yoga therapy versus blood pressure education program (sodium restriction and walking program) versus a combination program that involves components of both groups. LIMBS II will also examine the impact that changes in blood pressure have on cerebral blood flow. If successful, the LIMBS study will determine if yoga therapy combined with enhanced lifestyle modification will result in clinically meaningful decreases in blood pressure and thus can be implemented as an alternative to drug therapy for patients with prehypertension and stage 1 hypertension.

Cytoplasmic intron sequence-retaining transcripts can be dendritically targeted via ID element retrotransposons.

RNA precursors give rise to mRNA after splicing of intronic sequences traditionally thought to occur in the nucleus. Here, we show that intron sequences are retained in a number of dendritically-targeted mRNAs, by using microarray and Illumina sequencing of isolated dendritic mRNA as well as in situ hybridization. Many of the retained introns contain ID elements, a class of SINE retrotransposon. A portion of these SINEs confers dendritic targeting to exogenous and endogenous transcripts showing the necessity of ID-mediated mechanisms for the targeting of different transcripts to dendrites. ID elements are capable of selectively altering the distribution of endogenous proteins, providing a link between intronic SINEs and protein function. As such, the ID element represents a common dendritic targeting element found across multiple RNAs. Retention of intronic sequence is a more general phenomenon than previously thought and plays a functional role in the biology of the neuron, partly mediated by co-opted repetitive sequences.