A single-nucleus transcriptomic atlas of medium spiny neurons in the rat nucleus accumbens.

Neural processing of rewarding stimuli involves several distinct regions, including the nucleus accumbens (NAc). The majority of NAc neurons are GABAergic projection neurons known as medium spiny neurons (MSNs). MSNs are broadly defined by dopamine receptor expression, but evidence suggests that a wider array of subtypes exist. To study MSN heterogeneity, we analyzed single-nucleus RNA sequencing data from the largest available rat NAc dataset. Analysis of 48,040 NAc MSN nuclei identified major populations belonging to the striosome and matrix compartments. Integration with mouse and human data indicated consistency across species and disease-relevance scoring using genome-wide association study results revealed potentially differential roles for MSN populations in substance use disorders. Additional high-resolution clustering identified 34 transcriptomically distinct subtypes of MSNs definable by a limited number of marker genes. Together, these data demonstrate the diversity of MSNs in the NAc and provide a basis for more targeted genetic manipulation of specific populations.

Differential contributions of circadian clock genes to cell survival in bipolar disorder patient derived neuronal progenitor cells distinguishes lithium responders and non-responders.

Bipolar disorder (BD) is characterized by disrupted circadian rhythms and neuronal loss. Lithium is neuroprotective and used to treat BD, but outcomes are variable. Past research identified that circadian rhythms in BD patient neurons are associated with lithium response (Li-R) or non-response (Li-NR). However, the underlying cellular mechanisms remain unknown. To study interactions among circadian clock genes and cell survival, and their role in BD and predicting lithium response, we tested selected genes (PER1, BMAL1 and REV-ERBα) and small molecule modulators of ROR/REV-ERB nuclear receptors in models of cell survival using mouse neurons and stem-cell derived neuronal progenitor cells (NPC) from BD patients and controls. In apoptosis assays using staurosporine (STS), lithium was neuroprotective. Knockdown of PER1, BMAL1 and REV-ERBα modified cell survival across models. In NPCs, reduced expression of PER1 and BMAL1 led to more extensive cell death in Li-NR vs. Li-R. Reduced REV-ERBα expression caused more extensive cell death in BD vs. control NPCs, without distinguishing Li-R and Li-NR. In IMHN, The REV-ERB agonist GSK4112 had strong effects on circadian rhythm amplitude, and was neuroprotective in mouse neurons and control NPCs, but not in BD NPCs. Expression of cell survival genes following STS and GSK4112 treatments revealed BD-associated, and Li-R associated differences in expression profiles. We conclude that the neuroprotective response to lithium is similar in NPCs from Li-R and Li-NR. However, knockdown of circadian clock genes or stimulation of REV-ERBs reveal distinct contributions to cell death in BD patient NPCs, some of which distinguish Li-R and Li-NR.

GWAS Meta-Analysis of Suicide Attempt: Identification of 12 Genome-Wide Significant Loci and Implication of Genetic Risks for Specific Health Factors.

OBJECTIVE: Suicidal behavior is heritable and is a major cause of death worldwide. Two large-scale genome-wide association studies (GWASs) recently discovered and cross-validated genome-wide significant (GWS) loci for suicide attempt (SA). The present study leveraged the genetic cohorts from both studies to conduct the largest GWAS meta-analysis of SA to date. Multi-ancestry and admixture-specific meta-analyses were conducted within groups of significant African, East Asian, and European ancestry admixtures. METHODS: This study comprised 22 cohorts, including 43,871 SA cases and 915,025 ancestry-matched controls. Analytical methods across multi-ancestry and individual ancestry admixtures included inverse variance-weighted fixed-effects meta-analyses, followed by gene, gene-set, tissue-set, and drug-target enrichment, as well as summary-data-based Mendelian randomization with brain expression quantitative trait loci data, phenome-wide genetic correlation, and genetic causal proportion analyses. RESULTS: Multi-ancestry and European ancestry admixture GWAS meta-analyses identified 12 risk loci at p values <5×10-8. These loci were mostly intergenic and implicated DRD2, SLC6A9, FURIN, NLGN1, SOX5, PDE4B, and CACNG2. The multi-ancestry SNP-based heritability estimate of SA was 5.7% on the liability scale (SE=0.003, p=5.7×10-80). Significant brain tissue gene expression and drug set enrichment were observed. There was shared genetic variation of SA with attention deficit hyperactivity disorder, smoking, and risk tolerance after conditioning SA on both major depressive disorder and posttraumatic stress disorder. Genetic causal proportion analyses implicated shared genetic risk for specific health factors. CONCLUSIONS: This multi-ancestry analysis of suicide attempt identified several loci contributing to risk and establishes significant shared genetic covariation with clinical phenotypes. These findings provide insight into genetic factors associated with suicide attempt across ancestry admixture populations, in veteran and civilian populations, and in attempt versus death.

Placental microRNA methylome signatures may serve as biomarkers and therapeutic targets for prenatally opioid-exposed infants with neonatal opioid withdrawal syndrome.

Introduction: The neonate exposed to opioids in utero faces a constellation of withdrawal symptoms postpartum commonly called neonatal opioid withdrawal syndrome (NOWS). The incidence of NOWS has increased in recent years due to the opioid epidemic. MicroRNAs (miRNAs) are small non-coding RNA molecules that play a crucial role in gene regulation. Epigenetic variations in microRNAs (miRNAs) and their impact on addiction-related processes is a rapidly evolving area of research. Methods: The Illumina Infinium Methylation EPIC BeadChip was used to analyze DNA methylation levels of miRNA-encoding genes in 96 human placental tissues to identify miRNA gene methylation profiles as-sociated with NOWS: 32 from mothers whose prenatally opioid-exposed infants required pharmacologic management for NOWS, 32 from mothers whose prenatally opioid-exposed infants did not require treat-ment for NOWS, and 32 unexposed controls. Results: The study identified 46 significantly differentially methylated (FDR p-value ≤ 0.05) CpGs associated with 47 unique miRNAs, with a receiver operating characteristic (ROC) area under the curve (AUC) ≥0.75 including 28 hypomethylated and 18 hypermethylated CpGs as potentially associated with NOWS. These dysregulated microRNA methylation patterns may be a contributing factor to NOWS pathogenesis. Conclusion: This is the first study to analyze miRNA methylation profiles in NOWS infants and illustrates the unique role miRNAs might have in diagnosing and treating the disease. Furthermore, these data may provide a step toward feasible precision medicine for NOWS babies as well.

Contributions of circadian clock genes to cell survival in fibroblast models of lithium-responsive bipolar disorder.

Bipolar disorder (BD) is characterized by mood episodes, disrupted circadian rhythms and gray matter reduction in the brain. Lithium is an effective pharmacotherapy for BD, but not all patients respond to treatment. Lithium has neuroprotective properties and beneficial effects on circadian rhythms that may distinguish lithium responders (Li-R) from non-responders (Li-NR). The circadian clock regulates molecular pathways involved in apoptosis and cell survival, but how this overlap impacts BD and/or lithium responsiveness is unknown. In primary fibroblasts from Li-R/Li-NR BD patients and controls, we found patterns of co-expression among circadian clock and cell survival genes that distinguished BD vs. control, and Li-R vs. Li-NR cells. In cellular models of apoptosis using staurosporine (STS), lithium preferentially protected fibroblasts against apoptosis in BD vs. control samples, regardless of Li-R/Li-NR status. When examining the effects of lithium treatment of cells in vitro, caspase activation by lithium correlated with period alteration, but the relationship differed in control, Li-R and Li-NR samples. Knockdown of Per1 and Per3 in mouse fibroblasts altered caspase activity, cell death and circadian rhythms in an opposite manner. In BD cells, genetic variation in PER1 and PER3 predicted sensitivity to apoptosis in a manner consistent with knockdown studies. We conclude that distinct patterns of coordination between circadian clock and cell survival genes in BD may help predict lithium response.

Focal adhesion is associated with lithium response in bipolar disorder: evidence from a network-based multi-omics analysis.

Lithium (Li) is one of the most effective drugs for treating bipolar disorder (BD), however, there is presently no way to predict response to guide treatment. The aim of this study is to identify functional genes and pathways that distinguish BD Li responders (LR) from BD Li non-responders (NR). An initial Pharmacogenomics of Bipolar Disorder study (PGBD) GWAS of lithium response did not provide any significant results. As a result, we then employed network-based integrative analysis of transcriptomic and genomic data. In transcriptomic study of iPSC-derived neurons, 41 significantly differentially expressed (DE) genes were identified in LR vs NR regardless of lithium exposure. In the PGBD, post-GWAS gene prioritization using the GWA-boosting (GWAB) approach identified 1119 candidate genes. Following DE-derived network propagation, there was a highly significant overlap of genes between the top 500- and top 2000-proximal gene networks and the GWAB gene list (Phypergeometric = 1.28E-09 and 4.10E-18, respectively). Functional enrichment analyses of the top 500 proximal network genes identified focal adhesion and the extracellular matrix (ECM) as the most significant functions. Our findings suggest that the difference between LR and NR was a much greater effect than that of lithium. The direct impact of dysregulation of focal adhesion on axon guidance and neuronal circuits could underpin mechanisms of response to lithium, as well as underlying BD. It also highlights the power of integrative multi-omics analysis of transcriptomic and genomic profiling to gain molecular insights into lithium response in BD.

Common Genetic Variation and Age of Onset of Anorexia Nervosa.

Background: Genetics and biology may influence the age of onset of anorexia nervosa (AN). The aims of this study were to determine whether common genetic variation contributes to age of onset of AN and to investigate the genetic associations between age of onset of AN and age at menarche. Methods: A secondary analysis of the Psychiatric Genomics Consortium genome-wide association study (GWAS) of AN was performed, which included 9335 cases and 31,981 screened controls, all from European ancestries. We conducted GWASs of age of onset, early-onset AN (<13 years), and typical-onset AN, and genetic correlation, genetic risk score, and Mendelian randomization analyses. Results: Two loci were genome-wide significant in the typical-onset AN GWAS. Heritability estimates (single nucleotide polymorphism-h 2) were 0.01-0.04 for age of onset, 0.16-0.25 for early-onset AN, and 0.17-0.25 for typical-onset AN. Early- and typical-onset AN showed distinct genetic correlation patterns with putative risk factors for AN. Specifically, early-onset AN was significantly genetically correlated with younger age at menarche, and typical-onset AN was significantly negatively genetically correlated with anthropometric traits. Genetic risk scores for age of onset and early-onset AN estimated from independent GWASs significantly predicted age of onset. Mendelian randomization analysis suggested a causal link between younger age at menarche and early-onset AN. Conclusions: Our results provide evidence consistent with a common variant genetic basis for age of onset and implicate biological pathways regulating menarche and reproduction.

Single Nucleus Transcriptome Data from Alzheimer's Disease Mouse Models Yield New Insight into Pathophysiology.

BACKGROUND: 5XFAD humanized mutant mice and Trem2 knockout (T2KO) mice are two mouse models relevant to the study of Alzheimer's disease (AD)-related pathology. OBJECTIVE: To determine hippocampal transcriptomic and polyadenylation site usage alterations caused by genetic mutations engineered in 5XFAD and T2KO mice. METHODS: Employing a publicly available single-nucleus RNA sequencing dataset, we used Seurat and Sierra analytic programs to identify differentially expressed genes (DEGs) and differential transcript usage (DTU), respectively, in hippocampal cell types from each of the two mouse models. We analyzed cell type-specific DEGs further using Ingenuity Pathway Analysis (IPA). RESULTS: We identified several DEGs in both neuronal and glial cell subtypes in comparisons of wild type (WT) versus 5XFAD and WT versus T2KO mice, including Ttr, Fth1, Pcsk1n, Malat1, Rpl37, Rtn1, Sepw1, Uba52, Mbp, Arl6ip5, Gm26917, Vwa1, and Pgrmc1. We also observed DTU in common between the two comparisons in neuronal and glial subtypes, specifically in the genes Prnp, Rbm4b, Pnisr, Opcml, Cpne7, Adgrb1, Gabarapl2, Ubb, Ndfip1, Car11, and Stmn4. IPA identified three statistically significant canonical pathways that appeared in multiple cell types and that overlapped between 5XFAD and T2KO comparisons to WT, including 'FXR/RXR Activation', 'LXR/RXR Activation', and 'Acute Phase Response Signaling'. CONCLUSION: DEG, DTU, and IPA findings, derived from two different mouse models of AD, highlight the importance of energy imbalance and inflammatory processes in specific hippocampal cell types, including subtypes of neurons and glial cells, in the development of AD-related pathology. Additional studies are needed to further characterize these findings.

Multi-trait genome-wide association study of opioid addiction: OPRM1 and beyond.

Opioid addiction (OA) is moderately heritable, yet only rs1799971, the A118G variant in OPRM1, has been identified as a genome-wide significant association with OA and independently replicated. We applied genomic structural equation modeling to conduct a GWAS of the new Genetics of Opioid Addiction Consortium (GENOA) data together with published studies (Psychiatric Genomics Consortium, Million Veteran Program, and Partners Health), comprising 23,367 cases and effective sample size of 88,114 individuals of European ancestry. Genetic correlations among the various OA phenotypes were uniformly high (rg > 0.9). We observed the strongest evidence to date for OPRM1: lead SNP rs9478500 (p = 2.56 × 10-9). Gene-based analyses identified novel genome-wide significant associations with PPP6C and FURIN. Variants within these loci appear to be pleiotropic for addiction and related traits.

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.

Artificial Intelligence and Circulating Cell-Free DNA Methylation Profiling: Mechanism and Detection of Alzheimer's Disease.

Background: Despite extensive efforts, significant gaps remain in our understanding of Alzheimer's disease (AD) pathophysiology. Novel approaches using circulating cell-free DNA (cfDNA) have the potential to revolutionize our understanding of neurodegenerative disorders. Methods: We performed DNA methylation profiling of cfDNA from AD patients and compared them to cognitively normal controls. Six Artificial Intelligence (AI) platforms were utilized for the diagnosis of AD while enrichment analysis was used to elucidate the pathogenesis of AD. Results: A total of 3684 CpGs were significantly (adj. p-value < 0.05) differentially methylated in AD versus controls. All six AI algorithms achieved high predictive accuracy (AUC = 0.949−0.998) in an independent test group. As an example, Deep Learning (DL) achieved an AUC (95% CI) = 0.99 (0.95−1.0), with 94.5% sensitivity and specificity. Conclusion: We describe numerous epigenetically altered genes which were previously reported to be differentially expressed in the brain of AD sufferers. Genes identified by AI to be the best predictors of AD were either known to be expressed in the brain or have been previously linked to AD. We highlight enrichment in the Calcium signaling pathway, Glutamatergic synapse, Hedgehog signaling pathway, Axon guidance and Olfactory transduction in AD sufferers. To the best of our knowledge, this is the first reported genome-wide DNA methylation study using cfDNA to detect AD.

Correction of depression-associated circadian rhythm abnormalities is associated with lithium response in bipolar disorder.

BACKGROUND: Bipolar disorder (BD) is characterized by episodes of depression and mania and disrupted circadian rhythms. Lithium is an effective therapy for BD, but only 30%-40% of patients are fully responsive. Preclinical models show that lithium alters circadian rhythms. However, it is unknown if the circadian rhythm effects of lithium are essential to its therapeutic properties. METHODS: In secondary analyses of a multi-center, prospective, trial of lithium for BD, we examined the relationship between circadian rhythms and therapeutic response to lithium. Using standardized instruments, we measured morningness, diurnal changes in mood, sleep, and energy (circadian rhythm disturbances) in a cross-sectional study of 386 BD subjects with varying lithium exposure histories. Next, we tracked symptoms of depression and mania prospectively over 12 weeks in a subset of 88 BD patients initiating treatment with lithium. Total, circadian, and affective mood symptoms were scored separately and analyzed. RESULTS: Subjects with no prior lithium exposure had the most circadian disruption, while patients stable on lithium monotherapy had the least. Patients who were stable on lithium with another drug or unstable on lithium showed intermediate levels of disruption. Treatment with lithium for 12 weeks yielded significant reductions in total and affective depression symptoms. Lithium responders (Li-Rs) showed improvement in circadian symptoms of depression, but non-responders did not. There was no difference between Li-Rs and nonresponders in affective, circadian, or total symptoms of mania. CONCLUSIONS: Exposure to lithium is associated with reduced circadian disruption. Lithium response at 12 weeks was selectively associated with the reduction of circadian depressive symptoms. We conclude that stabilization of circadian rhythms may be an important feature of lithium's therapeutic effects. CLINICAL TRIALS REGISTRY: NCT0127253.

Analysis of differential gene expression and transcript usage in hippocampus of Apoe null mutant mice: Implications for Alzheimer's disease.

A dataset of single-nucleus RNA sequencing (snRNAseq) data was analyzed using Seurat, Sierra, and Ingenuity Pathway Analysis (IPA) programs to assess differentially expressed genes (DEGs) and differential transcript usage (DTU) in mouse hippocampal cell types. Seurat identified DEGs between the wild type (WT) and Apoe knockout (EKO) mice. IPA identified 11 statistically significant canonical pathways in >1 cell type. Sierra identified Sipa1l1 with DTU between WT and EKO samples. Analysis of the Sipa1l1 peak region identified an alternative non-canonical polyadenylation signal and a putative cytoplasmic polyadenylation element. APOE regulation of gene transcription and co-transcriptional RNA processing may underlie Alzheimer's disease.

Dissecting the Shared Genetic Architecture of Suicide Attempt, Psychiatric Disorders, and Known Risk Factors.

BACKGROUND: Suicide is a leading cause of death worldwide, and nonfatal suicide attempts, which occur far more frequently, are a major source of disability and social and economic burden. Both have substantial genetic etiology, which is partially shared and partially distinct from that of related psychiatric disorders. METHODS: We conducted a genome-wide association study (GWAS) of 29,782 suicide attempt (SA) cases and 519,961 controls in the International Suicide Genetics Consortium (ISGC). The GWAS of SA was conditioned on psychiatric disorders using GWAS summary statistics via multitrait-based conditional and joint analysis, to remove genetic effects on SA mediated by psychiatric disorders. We investigated the shared and divergent genetic architectures of SA, psychiatric disorders, and other known risk factors. RESULTS: Two loci reached genome-wide significance for SA: the major histocompatibility complex and an intergenic locus on chromosome 7, the latter of which remained associated with SA after conditioning on psychiatric disorders and replicated in an independent cohort from the Million Veteran Program. This locus has been implicated in risk-taking behavior, smoking, and insomnia. SA showed strong genetic correlation with psychiatric disorders, particularly major depression, and also with smoking, pain, risk-taking behavior, sleep disturbances, lower educational attainment, reproductive traits, lower socioeconomic status, and poorer general health. After conditioning on psychiatric disorders, the genetic correlations between SA and psychiatric disorders decreased, whereas those with nonpsychiatric traits remained largely unchanged. CONCLUSIONS: Our results identify a risk locus that contributes more strongly to SA than other phenotypes and suggest a shared underlying biology between SA and known risk factors that is not mediated by psychiatric disorders.

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.

Targeting the coronavirus nucleocapsid protein through GSK-3 inhibition.

The coronaviruses responsible for severe acute respiratory syndrome (SARS-CoV), COVID-19 (SARS-CoV-2), Middle East respiratory syndrome-CoV, and other coronavirus infections express a nucleocapsid protein (N) that is essential for viral replication, transcription, and virion assembly. Phosphorylation of N from SARS-CoV by glycogen synthase kinase 3 (GSK-3) is required for its function and inhibition of GSK-3 with lithium impairs N phosphorylation, viral transcription, and replication. Here we report that the SARS-CoV-2 N protein contains GSK-3 consensus sequences and that this motif is conserved in diverse coronaviruses, raising the possibility that SARS-CoV-2 may be sensitive to GSK-3 inhibitors, including lithium. We conducted a retrospective analysis of lithium use in patients from three major health systems who were PCR-tested for SARS-CoV-2. We found that patients taking lithium have a significantly reduced risk of COVID-19 (odds ratio = 0.51 [0.35-0.74], P = 0.005). We also show that the SARS-CoV-2 N protein is phosphorylated by GSK-3. Knockout of GSK3A and GSK3B demonstrates that GSK-3 is essential for N phosphorylation. Alternative GSK-3 inhibitors block N phosphorylation and impair replication in SARS-CoV-2 infected lung epithelial cells in a cell-type-dependent manner. Targeting GSK-3 may therefore provide an approach to treat COVID-19 and future coronavirus outbreaks.

Genetics and prescription opioid use (GaPO): study design for consenting a cohort from an existing biobank to identify clinical and genetic factors influencing prescription opioid use and abuse.

BACKGROUND: Prescription opioids (POs) are commonly used to treat moderate to severe chronic pain in the health system setting. Although they improve quality of life for many patients, more work is needed to identify both the clinical and genetic factors that put certain individuals at high risk for developing opioid use disorder (OUD) following use of POs for pain relief. With a greater understanding of important risk factors, physicians will be better able to identify patients at highest risk for developing OUD for whom non-opioid alternative therapies and treatments should be considered. METHODS: We are conducting a prospective observational study that aims to identify the clinical and genetic factors most stongly associated with OUD. The study design leverages an existing biobank that includes whole exome sequencing and array genotyping. The biobank is maintained within an integrated health system, allowing for the large-scale capture and integration of genetic and non-genetic data. Participants are enrolled into the health system biobank via informed consent and then into a second study that focuses on opioid medication use. Data capture includes validated self-report surveys measuring addiction severity, depression, anxiety, and nicotine use, as well as additional clinical, prescription, and brain imaging data extracted from electronic health records. DISCUSSION: We will harness this multimodal data capture to establish meaningful patient phenotypes in order to understand the genetic and non-genetic contributions to OUD.

Maternal opioid use disorder: Placental transcriptome analysis for neonatal opioid withdrawal syndrome.

Excessive prenatal opioid exposure may lead to the development of Neonatal Opioid Withdrawal Syndrome (NOWS). RNA-seq was done on 64 formalin-fixed paraffin-embedded placental tissue samples from 32 mothers with opioid use disorder, with newborns with NOWS that required treatment, and 32 prenatally unexposed controls. We identified 93 differentially expressed genes in the placentas of infants with NOWS compared to unexposed controls. There were 4 up- and 89 downregulated genes. Among these, 7 genes CYP1A1, APOB, RPH3A, NRXN1, LINC01206, AL157396.1, UNC80 achieved an FDR p-value of <0.01. The remaining 87 genes were significant with FDR p-value <0.05. The 4 upregulated, CYP1A1, FP671120.3, RAD1, RN7SL856P, and the 10 most significantly downregulated genes were RNA5SP364, GRIN2A, UNC5D, DMBT1P1, MIR3976HG, LINC02199, LINC02822, PANTR1, AC012178.1, CTNNA2. Ingenuity Pathway Analysis identified the 7 most likely to play an important role in the etiology of NOWS. Our study expands insights into the genetic mechanisms of NOWS development.

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.