Search

Summary: Key mechanisms underlying chronic pain occur within the dorsal horn. Genome-wide association studies (GWASs) have identified genetic variants predisposed to chronic pain. However, most of these variants lie within regulatory non-coding regions that have not been linked to spinal cord biology. Here, we take a multi-species approach to determine whether chronic pain variants impact the regulatory genomics of dorsal horn neurons. First, we generate a large rhesus macaque single-nucleus RNA sequencing (snRNA-seq) atlas and integrate it with available human and mouse datasets to produce a single unified, species-conserved atlas of neuron subtypes. Cellular-resolution spatial transcriptomics in mouse shows the precise laminar location of these neuron subtypes, consistent with our analysis of neuron-subtype-selective markers in macaque. Using this cross-species framework, we generate a mouse single-nucleus open chromatin atlas of regulatory elements that shows strong and selective relationships between the neuron-subtype-specific chromatin regions and variants from major chronic pain GWASs.

Abstract In brain, the striatum is a heterogenous region involved in reward and goal-directed behaviors. Striatal dysfunction is linked to psychiatric disorders, including opioid use disorder (OUD). Striatal subregions are divided based on neuroanatomy, each with unique roles in OUD. In OUD, the dorsal striatum is involved in altered reward processing, formation of habits, and development of negative affect during withdrawal. Using single nuclei RNA-sequencing, we identified both canonical (e.g., dopamine receptor subtype) and less abundant cell populations (e.g., interneurons) in human dorsal striatum. Pathways related to neurodegeneration, interferon response, and DNA damage were significantly enriched in striatal neurons of individuals with OUD. DNA damage markers were also elevated in striatal neurons of opioid-exposed rhesus macaques. Sex-specific molecular differences in glial cell subtypes associated with chronic stress were found in OUD, particularly female individuals. Together, we describe different cell types in human dorsal striatum and identify cell type-specific alterations in OUD.

Abstract Background Evolutionary conservation is an invaluable tool for inferring functional significance in the genome, including regions that are crucial across many species and those that have undergone convergent evolution. Computational methods to test for sequence conservation are dominated by algorithms that examine the ability of one or more nucleotides to align across large evolutionary distances. While these nucleotide alignment-based approaches have proven powerful for protein-coding genes and some non-coding elements, they fail to capture conservation of many enhancers, distal regulatory elements that control spatial and temporal patterns of gene expression. The function of enhancers is governed by a complex, often tissue- and cell type-specific code that links combinations of transcription factor binding sites and other regulation-related sequence patterns to regulatory activity. Thus, function of orthologous enhancer regions can be conserved across large evolutionary distances, even when nucleotide turnover is high. Results We present a new machine learning-based approach for evaluating enhancer conservation that leverages the combinatorial sequence code of enhancer activity rather than relying on the alignment of individual nucleotides. We first train a convolutional neural network model that can predict tissue-specific open chromatin, a proxy for enhancer activity, across mammals. Next, we apply that model to distinguish instances where the genome sequence would predict conserved function versus a loss of regulatory activity in that tissue. We present criteria for systematically evaluating model performance for this task and use them to demonstrate that our models accurately predict tissue-specific conservation and divergence in open chromatin between primate and rodent species, vastly out-performing leading nucleotide alignment-based approaches. We then apply our models to predict open chromatin at orthologs of brain and liver open chromatin regions across hundreds of mammals and find that brain enhancers associated with neuron activity have a stronger tendency than the general population to have predicted lineage-specific open chromatin. Conclusion The framework presented here provides a mechanism to annotate tissue-specific regulatory function across hundreds of genomes and to study enhancer evolution using predicted regulatory differences rather than nucleotide-level conservation measurements.

Abstract Severe and persistent disruptions to sleep and circadian rhythms are common in people with opioid use disorder (OUD). Preclinical evidence suggests altered molecular rhythms in the brain modulate opioid reward and relapse. However, whether molecular rhythms are disrupted in the brains of people with OUD remained an open question, critical to understanding the role of circadian rhythms in opioid addiction. Using subjects’ times of death as a marker of time of day, we investigated transcriptional rhythms in the brains of subjects with OUD compared to unaffected comparison subjects. We discovered rhythmic transcripts in both the dorsolateral prefrontal cortex (DLPFC) and nucleus accumbens (NAc), key brain areas involved in OUD, that were largely distinct between OUD and unaffected subjects. Fewer rhythmic transcripts were identified in DLPFC of subjects with OUD compared to unaffected subjects, whereas in the NAc, nearly double the number of rhythmic transcripts was identified in subjects with OUD. In NAc of subjects with OUD, rhythmic transcripts peaked either in the evening or near sunrise, and were associated with an opioid, dopamine, and GABAergic neurotransmission. Associations with altered neurotransmission in NAc were further supported by co-expression network analysis which identified OUD-specific modules enriched for transcripts involved in dopamine, GABA, and glutamatergic synaptic functions. Additionally, rhythmic transcripts in DLPFC and NAc of subjects with OUD were enriched for genomic loci associated with sleep-related GWAS traits, including sleep duration and insomnia. Collectively, our findings connect transcriptional rhythm changes in opioidergic, dopaminergic, GABAergic signaling in the human brain to sleep-related traits in opioid addiction.

We profile genome-wide histone 3 lysine 27 acetylation (H3K27ac) of 3 major brain cell types from hippocampus and dorsolateral prefrontal cortex (dlPFC) of subjects with and without Alzheimer’s Disease (AD). We confirm that single nucleotide polymorphisms (SNPs) associated with late onset AD (LOAD) show a strong tendency to reside in microglia-specific gene regulatory elements. Despite this significant colocalization, we find that microglia harbor more acetylation changes associated with age than with amyloid-β (Aβ) load. In contrast, we detect that an oligodendrocyte-enriched glial (OEG) population contains the majority of differentially acetylated peaks associated with Aβ load. These differential peaks reside near both early onset risk genes (APP, PSEN1, PSEN2) and late onset AD risk loci (including BIN1, PICALM, CLU, ADAM10, ADAMTS4, SORL1, FERMT2), Aβ processing genes (BACE1), as well as genes involved in myelinating and oligodendrocyte development processes. Interestingly, a number of LOAD risk loci associated with differentially acetylated risk genes contain H3K27ac peaks that are specifically enriched in OEG. These findings implicate oligodendrocyte gene regulation as a potential mechanism by which early onset and late onset risk genes mediate their effects, and highlight the deregulation of myelinating processes in AD. More broadly, our dataset serves as a resource for the study of functional effects of genetic variants and cell type specific gene regulation in AD.

Recent discoveries of extreme cellular diversity in the brain warrant rapid development of technologies to access specific cell populations within heterogeneous tissue. Available approaches for engineering-targeted technologies for new neuron subtypes are low yield, involving intensive transgenic strain or virus screening. Here, we present Specific Nuclear-Anchored Independent Labeling (SNAIL), an improved virus-based strategy for cell labeling and nuclear isolation from heterogeneous tissue. SNAIL works by leveraging machine learning and other computational approaches to identify DNA sequence features that confer cell type-specific gene activation and then make a probe that drives an affinity purification-compatible reporter gene. As a proof of concept, we designed and validated two novel SNAIL probes that target parvalbumin-expressing (PV+) neurons. Nuclear isolation using SNAIL in wild-type mice is sufficient to capture characteristic open chromatin features of PV+ neurons in the cortex, striatum, and external globus pallidus. The SNAIL framework also has high utility for multispecies cell probe engineering; expression from a mouse PV+ SNAIL enhancer sequence was enriched in PV+ neurons of the macaque cortex. Expansion of this technology has broad applications in cell type-specific observation, manipulation, and therapeutics across species and disease models.

Background: Adeno-associated virus (AAV)-mediated gene therapies are rapidly advancing to the clinic, and AAV engineering has resulted in vectors with increased ability to deliver therapeutic genes. Although the choice of vector is critical, quantitative comparison of AAVs, especially in large animals, remains challenging. Methods: Here, we developed an efficient single-cell AAV engineering pipeline (scAAVengr) to simultaneously quantify and rank efficiency of competing AAV vectors across all cell types in the same animal. Results: To demonstrate proof-of-concept for the scAAVengr workflow, we quantified – with cell-type resolution – the abilities of naturally occurring and newly engineered AAVs to mediate gene expression in primate retina following intravitreal injection. A top performing variant identified using this pipeline, K912, was used to deliver SaCas9 and edit the rhodopsin gene in macaque retina, resulting in editing efficiency similar to infection rates detected by the scAAVengr workflow. scAAVengr was then used to identify top-performing AAV variants in mouse brain, heart, and liver following systemic injection. Conclusions: These results validate scAAVengr as a powerful method for development of AAV vectors. Funding: This work was supported by funding from the Ford Foundation, NEI/NIH, Research to Prevent Blindness, Foundation Fighting Blindness, UPMC Immune Transplant and Therapy Center, and the Van Sloun fund for canine genetic research.

Calcium-Response Factor (CaRF) was first identified as a transcription factor based on its affinity for a neuronal-selective calcium-response element (CaRE1) in the gene encoding Brain-Derived Neurotrophic Factor (BDNF). However, because CaRF shares no homology with other transcription factors, its properties and gene targets have remained unknown. Here we show that the DNA binding domain of CaRF has been highly conserved across evolution and that CaRF binds DNA directly in a sequence-specific manner in the absence of other eukaryotic cofactors. Using a binding site selection screen we identify a high-affinity consensus CaRF response element (cCaRE) that shares significant homology with the CaRE1 element of Bdnf. In a genome-wide chromatin immunoprecipitation analysis (ChIP-Seq), we identified 176 sites of CaRF-specific binding (peaks) in neuronal genomic DNA. 128 of these peaks are within 10kB of an annotated gene, and 60 are within 1kB of an annotated transcriptional start site. At least 138 of the CaRF peaks contain a common 10-bp motif with strong statistical similarity to the cCaRE, and we provide evidence predicting that CaRF can bind independently to at least 64.5% of these motifs in vitro. Analysis of this set of putative CaRF targets suggests the enrichment of genes that regulate intracellular signaling cascades. Finally we demonstrate that expression of a subset of these target genes is altered in the cortex of Carf knockout (KO) mice. Together these data strongly support the characterization of CaRF as a unique transcription factor and provide the first insight into the program of CaRF-regulated transcription in neurons.