Single nucleus transcriptomic analysis of rat nucleus accumbens reveals cell type-specific patterns of gene expression associated with volitional morphine intake.

Opioid exposure is known to cause transcriptomic changes in the nucleus accumbens (NAc). However, no studies to date have investigated cell type-specific transcriptomic changes associated with volitional opioid taking. Here, we use single nucleus RNA sequencing (snRNAseq) to comprehensively characterize cell type-specific alterations of the NAc transcriptome in rats self-administering morphine. One cohort of male Brown Norway rats was injected with acute morphine (10 mg/kg, i.p.) or saline. A second cohort of rats was allowed to self-administer intravenous morphine (1.0 mg/kg/infusion) for 10 consecutive days. Each morphine-experienced rat was paired with a yoked saline control rat. snRNAseq libraries were generated from NAc punches and used to identify cell type-specific gene expression changes associated with volitional morphine taking. We identified 1106 differentially expressed genes (DEGs) in the acute morphine group, compared to 2453 DEGs in the morphine self-administration group, across 27 distinct cell clusters. Importantly, we identified 1329 DEGs that were specific to morphine self-administration. DEGs were identified in novel clusters of astrocytes, oligodendrocytes, and D1R- and D2R-expressing medium spiny neurons in the NAc. Cell type-specific DEGs included Rgs9, Celf5, Oprm1, and Pde10a. Upregulation of Rgs9 and Celf5 in D2R-expressing neurons was validated by RNAscope. Approximately 85% of all oligodendrocyte DEGs, nearly all of which were associated with morphine taking, were identified in two subtypes. Bioinformatic analyses identified cell type-specific upstream regulatory mechanisms of the observed transcriptome alterations and downstream signaling pathways, including both novel and previously identified molecular pathways. These findings show that volitional morphine taking is associated with distinct cell type-specific transcriptomic changes in the rat NAc and highlight specific striatal cell populations and novel molecular substrates that could be targeted to reduce compulsive opioid taking.

Inherited L1 Retrotransposon Insertions Associated With Risk for Schizophrenia and Bipolar Disorder.

Studies of the genetic heritability of schizophrenia and bipolar disorder examining single nucleotide polymorphisms (SNPs) and copy number variations have failed to explain a large portion of the genetic liability, resulting in substantial missing heritability. Long interspersed element 1 (L1) retrotransposons are a type of inherited polymorphic variant that may be associated with risk for schizophrenia and bipolar disorder. We performed REBELseq, a genome wide assay for L1 sequences, on DNA from male and female persons with schizophrenia and controls (n = 63 each) to identify inherited L1 insertions and validated priority insertions. L1 insertions of interest were genotyped in DNA from a replication cohort of persons with schizophrenia, bipolar disorder, and controls (n = 2268 each) to examine differences in carrier frequencies. We identified an inherited L1 insertion in ARHGAP24 and a quadallelic SNP (rs74169643) inside an L1 insertion in SNTG2 that are associated with risk for developing schizophrenia and bipolar disorder (all odds ratios ~1.2). Pathway analysis identified 15 gene ontologies that were differentially affected by L1 burden, including multiple ontologies related to glutamatergic signaling and immune function, which have been previously associated with schizophrenia. These findings provide further evidence supporting the role of inherited repetitive genetic elements in the heritability of psychiatric disorders.

Differential expression and transcription factor binding associated with genotype at a pharmacogenetic variant in OPRD1.

BACKGROUND: The functional mechanism is unknown for many genetic variants associated with substance use disorder phenotypes. Rs678849, an intronic variant in the delta-opioid receptor gene (OPRD1), has been found to predict regional brain volume, addiction risk, and the efficacy of buprenorphine/naloxone in treating opioid use disorder. The variant has also been implicated as an expression quantitative trait locus (eQTL) for several genes. OBJECTIVES: The objective of this study was to identify functional differences between the two alleles of rs678849 in vitro. We hypothesized that the two alleles of rs678849 would have different effects on transcriptional activity due to differential interactions with transcription factors. METHODS: 15bp regions containing the C or T alleles of rs678849 were cloned into luciferase constructs and transfected into BE(2)C neuroblastoma cells to test the effect on transcription. Electrophoretic mobility shift assays (EMSA) using nuclear lysates from BE(2)C cell or human postmortem medial prefrontal cortex were used to identify proteins that differentially bound the two alleles. RESULTS: At 24 hours post-transfection, the C allele construct had significantly lower luciferase expression than the T allele construct and empty vector control (ANOVA p < .001). Proteomic analysis and supershift assays identified XRCC6 as a transcription factor specifically binding the C allele, whereas hnRNP D0 was found to specifically bind the T allele. CONCLUSION: These functional differences between the C and T alleles may help explain the psychiatric and neurological phenotype differences predicted by rs678849 genotype and the potential role of the variant as an eQTL.

Investigation of long interspersed element-1 retrotransposons as potential risk factors for idiopathic temporal lobe epilepsy.

OBJECTIVE: To determine if long interspersed element-1 (L1) retrotransposons convey risk for idiopathic temporal lobe epilepsy (TLE). METHODS: Surgically resected temporal cortex from individuals with TLE (N = 33) and postmortem temporal cortex from individuals with no known neurological disease (N = 33) were analyzed for L1 content by Restriction Enzyme Based Enriched L1Hs sequencing (REBELseq). Expression of three KCNIP4 splice variants was assessed by droplet digital PCR (ddPCR). Protein ANalysis THrough Evolutionary Relationships (PANTHER) was used to determine ontologies and pathways for lists of genes harboring L1 insertions. RESULTS: We identified novel L1 insertions specific to individuals with TLE, and others specific to controls. Although there were no statistically significant differences between cases and controls in the numbers of known and novel L1 insertions, PANTHER analyses of intragenic L1 insertions showed statistically significant enrichments for epilepsy-relevant gene ontologies in both cases and controls. Gene ontologies "neuron projection development" and "calcium ion transmembrane transport" were among those found only in individuals with TLE. We confirmed novel L1 insertions in several genes associated with seizures/epilepsy, including a de novo somatic L1 retrotransposition in KCNIP4 that occurred after neural crest formation in one patient. However, ddPCR results suggest this de novo L1 did not alter KCNIP4 mRNA expression. SIGNIFICANCE: Given current data from this small cohort, we conclude that L1 elements, either rare heritable germline insertions or de novo somatic retrotranspositions, may contribute only minimally to overall genetic risk for idiopathic TLE. We suggest that further studies in additional patients and additional brain regions are warranted.

The role of demographic history and selection in shaping genetic diversity of the Galápagos penguin (Spheniscus mendiculus).

Although many studies have documented the effects of demographic bottlenecks on the genetic diversity of natural populations, there is conflicting evidence of the roles that genetic drift and selection may play in driving changes in genetic variation at adaptive loci. We analyzed genetic variation at microsatellite and mitochondrial loci in conjunction with an adaptive MHC class II locus in the Galápagos penguin (Spheniscus mendiculus), a species that has undergone serial demographic bottlenecks associated with El Niño events through its evolutionary history. We compared levels of variation in the Galápagos penguin to those of its congener, the Magellanic penguin (Spheniscus magellanicus), which has consistently maintained a large population size and thus was used as a non-bottlenecked control. The comparison of neutral and adaptive markers in these two demographically distinct species allowed assessment of the potential role of balancing selection in maintaining levels of MHC variation during bottleneck events. Our analysis suggests that the lack of genetic diversity at both neutral and adaptive loci in the Galápagos penguin likely resulted from its restricted range, relatively low abundance, and history of demographic bottlenecks. The Galápagos penguin revealed two MHC alleles, one mitochondrial haplotype, and six alleles across five microsatellite loci, which represents only a small fraction of the diversity detected in Magellanic penguins. Despite the decreased genetic diversity in the Galápagos penguin, results revealed signals of balancing selection at the MHC, which suggest that selection can mitigate some of the effects of genetic drift during bottleneck events. Although Galápagos penguin populations have persisted for a long time, increased frequency of El Niño events due to global climate change, as well as the low diversity exhibited at immunological loci, may put this species at further risk of extinction.