Addiction-Associated Genetic Variants Implicate Brain Cell Type- and Region-Specific Cis-Regulatory Elements in Addiction Neurobiology.

Recent large genome-wide association studies have identified multiple confident risk loci linked to addiction-associated behavioral traits. Most genetic variants linked to addiction-associated traits lie in noncoding regions of the genome, likely disrupting cis-regulatory element (CRE) function. CREs tend to be highly cell type-specific and may contribute to the functional development of the neural circuits underlying addiction. Yet, a systematic approach for predicting the impact of risk variants on the CREs of specific cell populations is lacking. To dissect the cell types and brain regions underlying addiction-associated traits, we applied stratified linkage disequilibrium score regression to compare genome-wide association studies to genomic regions collected from human and mouse assays for open chromatin, which is associated with CRE activity. We found enrichment of addiction-associated variants in putative CREs marked by open chromatin in neuronal (NeuN+) nuclei collected from multiple prefrontal cortical areas and striatal regions known to play major roles in reward and addiction. To further dissect the cell type-specific basis of addiction-associated traits, we also identified enrichments in human orthologs of open chromatin regions of female and male mouse neuronal subtypes: cortical excitatory, D1, D2, and PV. Last, we developed machine learning models to predict mouse cell type-specific open chromatin, enabling us to further categorize human NeuN+ open chromatin regions into cortical excitatory or striatal D1 and D2 neurons and predict the functional impact of addiction-associated genetic variants. Our results suggest that different neuronal subtypes within the reward system play distinct roles in the variety of traits that contribute to addiction.SIGNIFICANCE STATEMENT We combine statistical genetic and machine learning techniques to find that the predisposition to for nicotine, alcohol, and cannabis use behaviors can be partially explained by genetic variants in conserved regulatory elements within specific brain regions and neuronal subtypes of the reward system. Our computational framework can flexibly integrate open chromatin data across species to screen for putative causal variants in a cell type- and tissue-specific manner for numerous complex traits.

Inferring mammalian tissue-specific regulatory conservation by predicting tissue-specific differences in open chromatin.

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.

Cell Type-Specific Oxidative Stress Genomic Signatures in the Globus Pallidus of Dopamine-Depleted Mice.

Neuron subtype dysfunction is a key contributor to neurologic disease circuits, but identifying associated gene regulatory pathways is complicated by the molecular complexity of the brain. For example, parvalbumin-expressing (PV+) neurons in the external globus pallidus (GPe) are critically involved in the motor deficits of dopamine-depleted mouse models of Parkinson's disease, where cell type-specific optogenetic stimulation of PV+ neurons over other neuron populations rescues locomotion. Despite the distinct roles these cell types play in the neural circuit, the molecular correlates remain unknown because of the difficulty of isolating rare neuron subtypes. To address this issue, we developed a new viral affinity purification strategy, Cre-Specific Nuclear Anchored Independent Labeling, to isolate Cre recombinase-expressing (Cre+) nuclei from the adult mouse brain. Applying this technology, we performed targeted assessments of the cell type-specific transcriptomic and epigenetic effects of dopamine depletion on PV+ and PV- cells within three brain regions of male and female mice: GPe, striatum, and cortex. We found GPe PV+ neuron-specific gene expression changes that suggested increased hypoxia-inducible factor 2α signaling. Consistent with transcriptomic data, regions of open chromatin affected by dopamine depletion within GPe PV+ neurons were enriched for hypoxia-inducible factor family binding motifs. The gene expression and epigenomic experiments performed on PV+ neurons isolated by Cre-Specific Nuclear Anchored Independent Labeling identified a transcriptional regulatory network mediated by the neuroprotective factor Hif2a as underlying neural circuit differences in response to dopamine depletion.SIGNIFICANCE STATEMENT Cre-Specific Nuclear Anchored Independent Labeling is an enhanced, virus-based approach to isolate nuclei of a specific cell type for transcriptome and epigenome interrogation that decreases dependency on transgenic animals. Applying this technology to GPe parvalbumin-expressing neurons in a mouse model of Parkinson's disease, we discovered evidence for an upregulation of the oxygen homeostasis maintaining pathway involving Hypoxia-inducible factor 2α. These results provide new insight into how neuron subtypes outside the substantia nigra pars compacta may be compensating at a molecular level for differences in the motor production neural circuit during the progression of Parkinson's disease. Furthermore, they emphasize the utility of cell type-specific technologies, such as Cre-Specific Nuclear Anchored Independent Labeling, for isolated assessment of specific neuron subtypes in complex systems.

An mRNA-based platform for the delivery of pathogen-specific IgA into mucosal secretions.

Colonization of the gut and airways by pathogenic bacteria can lead to local tissue destruction and life-threatening systemic infections, especially in immunologically compromised individuals. Here, we describe an mRNA-based platform enabling delivery of pathogen-specific immunoglobulin A (IgA) monoclonal antibodies into mucosal secretions. The platform consists of synthetic mRNA encoding IgA heavy, light, and joining (J) chains, packaged in lipid nanoparticles (LNPs) that express glycosylated, dimeric IgA with functional activity in vitro and in vivo. Importantly, mRNA-derived IgA had a significantly greater serum half-life and a more native glycosylation profile in mice than did a recombinantly produced IgA. Expression of an mRNA encoded Salmonella-specific IgA in mice resulted in intestinal localization and limited Peyer's patch invasion. The same mRNA-LNP technology was used to express a Pseudomonas-specific IgA that protected from a lung challenge. Leveraging the mRNA antibody technology as a means to intercept bacterial pathogens at mucosal surfaces opens up avenues for prophylactic and therapeutic interventions.

Transcriptional and anatomical diversity of medium spiny neurons in the primate striatum.

Medium spiny neurons (MSNs) constitute the vast majority of striatal neurons and the principal interface between dopamine reward signals and functionally diverse cortico-basal ganglia circuits. Information processing in these circuits is dependent on distinct MSN types: cell types that are traditionally defined according to their projection targets or dopamine receptor expression. Single-cell transcriptional studies have revealed greater MSN heterogeneity than predicted by traditional circuit models, but the transcriptional landscape in the primate striatum remains unknown. Here, we set out to establish molecular definitions for MSN subtypes in Rhesus monkeys and to explore the relationships between transcriptionally defined subtypes and anatomical subdivisions of the striatum. Our results suggest at least nine MSN subtypes, including dorsal striatum subtypes associated with striosome and matrix compartments, ventral striatum subtypes associated with the nucleus accumbens shell and olfactory tubercle, and an MSN-like cell type restricted to µ-opioid receptor rich islands in the ventral striatum. Although each subtype was demarcated by discontinuities in gene expression, continuous variation within subtypes defined gradients corresponding to anatomical locations and, potentially, functional specializations. These results lay the foundation for achieving cell-type-specific transgenesis in the primate striatum and provide a blueprint for investigating circuit-specific information processing.

scAAVengr, a transcriptome-based pipeline for quantitative ranking of engineered AAVs with single-cell resolution.

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.

Gene therapy is an experimental approach to treating disease that involves altering faulty genes or replacing them with new, working copies. Most often, the new genetic material is delivered into cells using a modified virus that no longer causes disease, called a viral vector. Virus-mediated gene therapies are currently being explored for degenerative eye diseases, such as retinitis pigmentosa, and neurological disorders, like Alzheimer's and Parkinson's disease. A number of gene therapies have also been approved for treating some rare cancers, blood disorders and a childhood form of motor neuron disease. Despite the promise of virus-mediated gene therapy, there are significant hurdles to its widespread success. Viral vectors need to deliver enough genetic material to the right cells without triggering an immune response or causing serious side effects. Selecting an optimal vector is key to achieving this. A type of viruses called adeno-associated viruses (AAV) are prime candidates, partly because they can be easily engineered. However, accurately comparing the safety and efficacy of newly engineered AAVs is difficult, due to variation between test subjects and the labor and cost involved in careful testing. Öztürk et al. addressed this issue by developing an experimental pipeline called scAAVengr for comparing gene therapy vectors head-to-head. The process involves tagging potential AAV vectors with unique genetic barcodes, which can then be detected and quantified in individual cells using a technique called single-cell RNA sequencing. This means that when several vectors are used to infect lab-grown cells or a test animal at the same time, they can be tracked. The vectors can then be ranked on their ability to infect specific cell types and deliver useful genetic material. Using scAAVengr, Öztürk et al. compared viral vectors designed to target the light-sensitive cells of the retina, which allow animals to see. First, a set of promising viral vectors were evaluated using the scAAVengr pipeline in the eyes of marmosets and macaques, two small primates. Precise levels and locations of gene delivery were quantified. The top-performing vector was then identified and used to deliver Cas9, a genome editing tool, to primate retinas. Öztürk et al. also used scAAVengr to compare viral vectors in mice, analysing the vectors' ability to deliver their genetic cargo to the brain, heart, and liver. These experiments demonstrated that scAAVengr can be used to evaluate vectors in multiple tissues and in different organisms. In summary, this work outlines a method for identifying and precisely quantifying the performance of top-performing viral vectors for gene therapy. By aiding the selection of optimal viral vectors, the scAAVengr pipeline could help to improve the success of preclinical studies and early clinical trials testing gene therapies.

Selecting the most appropriate time points to profile in high-throughput studies.

Biological systems are increasingly being studied by high throughput profiling of molecular data over time. Determining the set of time points to sample in studies that profile several different types of molecular data is still challenging. Here we present the Time Point Selection (TPS) method that solves this combinatorial problem in a principled and practical way. TPS utilizes expression data from a small set of genes sampled at a high rate. As we show by applying TPS to study mouse lung development, the points selected by TPS can be used to reconstruct an accurate representation for the expression values of the non selected points. Further, even though the selection is only based on gene expression, these points are also appropriate for representing a much larger set of protein, miRNA and DNA methylation changes over time. TPS can thus serve as a key design strategy for high throughput time series experiments. Supporting Website: www.sb.cs.cmu.edu/TPS.

Tradeoffs between Dense and Replicate Sampling Strategies for High-Throughput Time Series Experiments.

An important experimental design question for high-throughput time series studies is the number of replicates required for accurate reconstruction of the profiles. Due to budget and sample availability constraints, more replicates imply fewer time points and vice versa. We analyze the performance of dense and replicate sampling by developing a theoretical framework that focuses on a restricted yet expressive set of possible curves over a wide range of noise levels and by analyzing real expression data. For both the theoretical analysis and experimental data, we observe that, under reasonable noise levels, autocorrelations in the time series data allow dense sampling to better determine the correct levels of non-sampled points when compared to replicate sampling. A Java implementation of our framework can be used to determine the best replicate strategy given the expected noise. These results provide theoretical support to the large number of high-throughput time series experiments that do not use replicates.