Mimicking opioid analgesia in cortical pain circuits

The anterior cingulate cortex plays a pivotal role in the cognitive and affective aspects of pain perception. Both endogenous and exogenous opioid signaling within the cingulate mitigate cortical nociception, reducing pain unpleasantness. However, the specific functional and molecular identities of cells mediating opioid analgesia in the cingulate remain elusive. Given the complexity of pain as a sensory and emotional experience, and the richness of ethological pain-related behaviors, we developed a standardized, deep-learning platform for deconstructing the behavior dynamics associated with the affective component of pain in mice-LUPE (Light aUtomated Pain Evaluator). LUPE removes human bias in behavior quantification and accelerated analysis from weeks to hours, which we leveraged to discover that morphine altered attentional and motivational pain behaviors akin to affective analgesia in humans. Through activity-dependent genetics and single-nuclei RNA sequencing, we identified specific ensembles of nociceptive cingulate neuron-types expressing mu-opioid receptors. Tuning receptor expression in these cells bidirectionally modulated morphine analgesia. Moreover, we employed a synthetic opioid receptor promoter-driven approach for cell-type specific optical and chemical genetic viral therapies to mimic morphine's pain-relieving effects in the cingulate, without reinforcement. This approach offers a novel strategy for precision pain management by targeting a key nociceptive cortical circuit with on-demand, non-addictive, and effective analgesia.

A nociceptive amygdala-striatal pathway for chronic pain aversion

The basolateral amygdala (BLA) is essential for assigning positive or negative valence to sensory stimuli. Noxious stimuli that cause pain are encoded by an ensemble of nociceptive BLA projection neurons (BLAnoci ensemble). However, the role of the BLAnoci ensemble in mediating behavior changes and the molecular signatures and downstream targets distinguishing this ensemble remain poorly understood. Here, we show that the same BLAnoci ensemble neurons are required for both acute and chronic neuropathic pain behavior. Using single nucleus RNA-sequencing, we characterized the effect of acute and chronic pain on the BLA and identified enrichment for genes with known functions in axonal and synaptic organization and pain perception. We thus examined the brain-wide targets of the BLAnoci ensemble and uncovered a previously undescribed nociceptive hotspot of the nucleus accumbens shell (NAcSh) that mirrors the stability and specificity of the BLAnoci ensemble and is recruited in chronic pain. Notably, BLAnoci ensemble axons transmit acute and neuropathic nociceptive information to the NAcSh, highlighting this nociceptive amygdala-striatal circuit as a unique pathway for affective-motivational responses across pain states.

Calcitonin receptor signaling in nucleus accumbens D1R- and D2R-expressing medium spiny neurons bidirectionally alters opioid taking in male rats

The high rates of relapse associated with current medications used to treat opioid use disorder (OUD) necessitate research that expands our understanding of the neural mechanisms regulating opioid taking to identify molecular substrates that could be targeted by novel pharmacotherapies to treat OUD. Recent studies show that activation of calcitonin receptors (CTRs) is sufficient to reduce the rewarding effects of addictive drugs in rodents. However, the role of central CTR signaling in opioid-mediated behaviors has not been studied. Here, we used single nuclei RNA sequencing (snRNA-seq), fluorescent in situ hybridization (FISH), and immunohistochemistry (IHC) to characterize cell type-specific patterns of CTR expression in the nucleus accumbens (NAc), a brain region that plays a critical role in voluntary drug taking. Using these approaches, we identified CTRs expressed on D1R- and D2R-expressing medium spiny neurons (MSNs) in the medial shell subregion of the NAc. Interestingly, Calcr transcripts were expressed at higher levels in D2R- versus D1R-expressing MSNs. Cre-dependent viral-mediated miRNA knockdown of CTRs in transgenic male rats was then used to determine the functional significance of endogenous CTR signaling in opioid taking. We discovered that reduced CTR expression specifically in D1R-expressing MSNs potentiated/augmented opioid self-administration. In contrast, reduced CTR expression specifically in D2R-expressing MSNs attenuated opioid self-administration. These findings highlight a novel cell type-specific mechanism by which CTR signaling in the ventral striatum bidirectionally modulates voluntary opioid taking and support future studies aimed at targeting central CTR-expressing circuits to treat OUD.

Ventral hippocampus neurons encode meal-related memory

The ability to encode and retrieve meal-related information is critical to efficiently guide energy acquisition and consumption, yet the underlying neural processes remain elusive. Here we reveal that ventral hippocampus (HPCv) neuronal activity dynamically elevates during meal consumption and this response is highly predictive of subsequent performance in a foraging-related spatial memory task. Targeted recombination-mediated ablation of HPCv meal-responsive neurons impairs foraging-related spatial memory without influencing food motivation, anxiety-like behavior, or escape-mediated spatial memory. These HPCv meal-responsive neurons project to the lateral hypothalamic area (LHA) and single-nucleus RNA sequencing and in situ hybridization analyses indicate they are enriched in serotonin 2a receptors (5HT2aR). Either chemogenetic silencing of HPCv-to-LHA projections or intra-HPCv 5HT2aR antagonist yielded foraging-related spatial memory deficits, as well as alterations in caloric intake and the temporal sequence of spontaneous meal consumption. Collective results identify a population of HPCv neurons that dynamically respond to eating to encode meal-related memories.

Identifying genetic loci and phenomic associations of substance use traits: A multi-trait analysis of GWAS (MTAG) study

BACKGROUND AND AIMS

Genome-wide association studies (GWAS) of opioid use disorder (OUD) and cannabis use disorder (CUD) have lagged behind those of alcohol use disorder (AUD) and smoking, where many more loci have been identified. We sought to identify novel loci for substance use traits (SUTs) in both African- (AFR) and European- (EUR) ancestry individuals to enhance our understanding of the traits' genetic architecture.

DESIGN

We used multi-trait analysis of GWAS (MTAG) to analyze four SUTs in EUR subjects (OUD, CUD, AUD and smoking initiation [SMKinitiation]), and three SUTs in AFR subjects (OUD, AUD and smoking trajectory [SMKtrajectory]). We conducted gene-set and protein-protein interaction analyses and calculated polygenic risk scores (PRS) in two independent samples.

SETTING

This study was conducted in the United States.

PARTICIPANTS

A total of 5692 EUR and 4918 AFR individuals in the Yale-Penn sample and 29 054 EUR and 10 265 AFR individuals in the Penn Medicine BioBank sample.

FINDINGS

MTAG identified genome-wide significant (GWS) single nucleotide polymorphisms (SNPs) for all four traits in EUR: 41 SNPs in 36 loci for OUD; 74 SNPs in 60 loci for CUD; 63 SNPs in 52 loci for AUD; and 183 SNPs in 144 loci for SMKinitiation. MTAG also identified GWS SNPs in AFR: 2 SNPs in 2 loci for OUD; 3 SNPs in 3 loci for AUD; and 1 SNP in 1 locus for SMKtrajectory. In the Yale-Penn sample, the MTAG-derived PRS consistently yielded more significant associations with both the corresponding substance use disorder diagnosis and multiple related phenotypes than the GWAS-derived PRS.

CONCLUSIONS

Multi-trait analysis of genome-wide association studies boosted the number of loci found for substance use traits, identifying genes not previously linked to any substance, and increased the power of polygenic risk scores. Multi-trait analysis of genome-wide association studies can be used to identify novel associations for substance use, especially those for which the samples are smaller than those for historically legal substances.

GIP receptor agonism blocks chemotherapy-induced nausea and vomiting

OBJECTIVE

Nausea and vomiting remain life-threatening obstacles to successful treatment of chronic diseases, despite a cadre of available antiemetic medications. Our inability to effectively control chemotherapy-induced nausea and vomiting (CINV) highlights the need to anatomically, molecularly, and functionally characterize novel neural substrates that block CINV.

METHODS

Behavioral pharmacology assays of nausea and emesis in 3 different mammalian species were combined with histological and unbiased transcriptomic analyses to investigate the beneficial effects of glucose-dependent insulinotropic polypeptide receptor (GIPR) agonism on CINV.

RESULTS

Single-nuclei transcriptomics and histological approaches in rats revealed a topographical, molecularly distinct, GABA-ergic neuronal population in the dorsal vagal complex (DVC) that is modulated by chemotherapy but rescued by GIPR agonism. Activation of DVCGIPR neurons substantially decreased behaviors indicative of malaise in cisplatin-treated rats. Strikingly, GIPR agonism blocks cisplatin-induced emesis in both ferrets and shrews.

CONCLUSION

Our multispecies study defines a peptidergic system that represents a novel therapeutic target for the management of CINV, and potentially other drivers of nausea/emesis.

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.

GPR-160 Receptor Signaling in the Dorsal Vagal Complex of Male Rats Modulates Meal Microstructure and CART-Mediated Hypophagia

The g-protein coupled receptor GPR-160, recently identified as a putative receptor for the cocaine and amphetamine-regulated transcript (CART) peptide, shows abundant expression in the energy-balance control nuclei, including the dorsal vagal complex (DVC). However, its physiological role in the control of food intake has yet to be fully explored. Here, we performed a virally mediated, targeted knockdown (KD) of Gpr160 in the DVC of male rats to evaluate its physiological role in control of feeding. Our results indicate that DVC Gpr160 KD affects meal microstructure. Specifically, DVC Gpr160 KD animals consumed more frequent, but shorter meals during the dark phase and showed decreased caloric intake and duration of meals during the light phase. Cumulatively, however, these bidirectional effects on feeding resulted in no difference in body weight gain. We next tested the role of DVC GPR-160 in mediating the anorexigenic effects of exogenous CART. Our results show that DVC Gpr160 KD partially attenuates CART's anorexigenic effects. To further characterize Gpr160+ cells in the DVC, we utilized single-nucleus RNA sequencing data to uncover abundant GPR-160 expression in DVC microglia and only minimal expression in neurons. Altogether, our results suggest that DVC CART signaling may be mediated by Gpr160+ microglia, which in turn may be modulating DVC neuronal activity to control food intake.

Type I interferon response in astrocytes promotes brain metastasis by enhancing monocytic myeloid cell recruitment

Cancer metastasis to the brain is a significant clinical problem. Metastasis is the consequence of favorable interactions between invaded cancer cells and the microenvironment. Here, we demonstrate that cancer-activated astrocytes create a sustained low-level activated type I interferon (IFN) microenvironment in brain metastatic lesions. We further confirm that the IFN response in astrocytes facilitates brain metastasis. Mechanistically, IFN signaling in astrocytes activates C-C Motif Chemokine Ligand 2 (CCL2) production, which further increases the recruitment of monocytic myeloid cells. The correlation between CCL2 and monocytic myeloid cells is confirmed in clinical brain metastasis samples. Lastly, genetically or pharmacologically inhibiting C-C Motif Chemokine Receptor 2 (CCR2) reduces brain metastases. Our study clarifies a pro-metastatic effect of type I IFN in the brain even though IFN response has been considered to have anti-tumor effects. Moreover, this work expands our understandings on the interactions between cancer-activated astrocytes and immune cells in brain metastasis.

Unraveling Psychiatric Disorders through Neural Single-Cell Transcriptomics Approaches

The development of single-cell and single-nucleus transcriptome technologies is enabling the unraveling of the molecular and cellular heterogeneity of psychiatric disorders. The complexity of the brain and the relationships between different brain regions can be better understood through the classification of individual cell populations based on their molecular markers and transcriptomic features. Analysis of these unique cell types can explain their involvement in the pathology of psychiatric disorders. Recent studies in both human and animal models have emphasized the importance of transcriptome analysis of neuronal cells in psychiatric disorders but also revealed critical roles for non-neuronal cells, such as oligodendrocytes and microglia. In this review, we update current findings on the brain transcriptome and explore molecular studies addressing transcriptomic alterations identified in human and animal models in depression and stress, neurodegenerative disorders (Parkinson's and Alzheimer's disease), schizophrenia, opioid use disorder, and alcohol and psychostimulant abuse. We also comment on potential future directions in single-cell and single-nucleus studies.

An intensive longitudinal examination of topiramate treatment for alcohol use disorder: a secondary analysis of data from a randomized controlled trial

BACKGROUND AND AIMS

Previous findings have been equivocal as to whether a single-nucleotide polymorphism (rs2832407) in GRIK1, which encodes a glutamate receptor subunit, moderates the effects of topiramate treatment for drinking reduction. We leveraged intensive longitudinal data to provide greater precision and allow an examination of intermediate outcomes addressing this question. We used data from a randomized controlled trial (RCT) to test the hypotheses that topiramate treatment reduces daily heavy drinking, desire to drink and positive alcohol expectancies and that these effects are stronger in rs2832407*C-allele homozygotes.

DESIGN

Secondary data analysis of a randomized controlled trial.

SETTING

University of Pennsylvania Treatment Research Center in the United States.

PARTICIPANTS/CASES

Participants were 164 individuals (70.1% male, mean age = 51.42, 36.0% rs2832407*C-allele homozygotes) who sought to reduce or stop drinking.

INTERVENTION AND COMPARATOR

Participants were assigned to medication (topiramate or placebo), with stratification by genotype group (CC versus AA/AC) and treatment goal (reduce versus abstain).

MEASUREMENTS

During the 12-week treatment period, participants completed daily interactive voice response (IVR) surveys.

FINDINGS

On any given day during treatment, participants who received topiramate had lower odds of IVR-reported heavy drinking [odds ratio (OR) = 0.259, b (standard error, SE) = -1.351 (0.334), P < 0.001] and lower levels of desire to drink [b (SE) = -0.323 (0.122), P = 0.009] and positive alcohol expectancies [b (SE) = -0.347 (0.138), P = 0.013] than those who received placebo. Participants who received topiramate also reported greater reductions in positive alcohol expectancies during the first 2 weeks of treatment than those who received placebo [b (SE) = -0.028 (0.008), P = 0.001], but topiramate did not impact the daily rate of change in heavy drinking or desire to drink. Genotype did not moderate the effects of topiramate on any outcomes examined (P > 0.05).

CONCLUSIONS

Topiramate is an effective medication for individuals seeking to reduce heavy drinking. The effects are not moderated by the single-nucleotide polymorphism rs2832407.

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.

Cross-ancestry meta-analysis of opioid use disorder uncovers novel loci with predominant effects in brain regions associated with addiction

Despite an estimated heritability of ~50%, genome-wide association studies of opioid use disorder (OUD) have revealed few genome-wide significant loci. We conducted a cross-ancestry meta-analysis of OUD in the Million Veteran Program (N = 425,944). In addition to known exonic variants in OPRM1 and FURIN, we identified intronic variants in RABEPK, FBXW4, NCAM1 and KCNN1. A meta-analysis including other datasets identified a locus in TSNARE1. In total, we identified 14 loci for OUD, 12 of which are novel. Significant genetic correlations were identified for 127 traits, including psychiatric disorders and other substance use-related traits. The only significantly enriched cell-type group was CNS, with gene expression enrichment in brain regions previously associated with substance use disorders. These findings increase our understanding of the biological basis of OUD and provide further evidence that it is a brain disease, which may help to reduce stigma and inform efforts to address the opioid epidemic.

Identifying High-Risk Comorbidities Associated with Opioid Use Patterns Using Electronic Health Record Prescription Data

Introduction

Opioid use disorders (OUDs) constitute a major public health issue, and we urgently need alternative methods for characterizing risk for OUD. Electronic health records (EHRs) are useful tools for understanding complex medical phenotypes but have been underutilized for OUD because of challenges related to underdiagnosis, binary diagnostic frameworks, and minimally characterized reference groups. As a first step in addressing these challenges, a new paradigm is warranted that characterizes risk for opioid prescription misuse on a continuous scale of severity, i.e., as a continuum.

Methods

Across sites within the PsycheMERGE network, we extracted prescription opioid data and diagnoses that co-occur with OUD (including psychiatric and substance use disorders, pain-related diagnoses, HIV, and hepatitis C) for over 2.6 million patients across three health registries (Vanderbilt University Medical Center, Mass General Brigham, Geisinger) between 2005 and 2018. We defined three groups based on levels of opioid exposure: no prescriptions, minimal exposure, and chronic exposure and then compared the comorbidity profiles of these groups to the full registries and to those with OUD diagnostic codes.

Results

Our results confirm that EHR data reflects known higher prevalence of substance use disorders, psychiatric disorders, medical, and pain diagnoses in patients with OUD diagnoses and chronic opioid use. Comorbidity profiles that distinguish opioid exposure are strikingly consistent across large health systems, indicating the phenotypes described in this new quantitative framework are robust to health systems differences.

Conclusion

This work indicates that EHR prescription opioid data can serve as a platform to characterize complex risk markers for OUD using existing data.

Paternal nicotine taking elicits heritable sex-specific phenotypes that are mediated by hippocampal Satb2

Nicotine intake, whether through tobacco smoking or e-cigarettes, remains a global health concern. An emerging preclinical literature indicates that parental nicotine exposure produces behavioral, physiological, and molecular changes in subsequent generations. However, the heritable effects of voluntary parental nicotine taking are unknown. Here, we show increased acquisition of nicotine taking in male and female offspring of sires that self-administered nicotine. In contrast, self-administration of sucrose and cocaine were unaltered in male and female offspring suggesting that the intergenerational effects of paternal nicotine taking may be reinforcer specific. Further characterization revealed memory deficits and increased anxiety-like behaviors in drug-naive male, but not female, offspring of nicotine-experienced sires. Using an unbiased, genome-wide approach, we discovered that these phenotypes were associated with decreased expression of Satb2, a transcription factor known to play important roles in synaptic plasticity and memory formation, in the hippocampus of nicotine-sired male offspring. This effect was sex-specific as no changes in Satb2 expression were found in nicotine-sired female offspring. Finally, increasing Satb2 levels in the hippocampus prevented the escalation of nicotine intake and rescued the memory deficits associated with paternal nicotine taking in male offspring. Collectively, these findings indicate that paternal nicotine taking produces heritable sex-specific molecular changes that promote addiction-like phenotypes and memory impairments in male offspring.

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.

Post-treatment effects of topiramate on alcohol-related outcomes: A combined analysis of two placebo-controlled trials

Topiramate reduces drinking and alcohol-related problems and is increasingly being used to treat alcohol use disorder (AUD). In a randomized controlled trial (RCT) of topiramate, rs2832407, a single nucleotide polymorphism (SNP) in the GRIK1 gene moderated topiramate's effects (Study 1). However, a second RCT (Study 2) did not replicate the SNP's moderating effect during treatment. The current analysis combines data from these two studies to examine topiramate's effects on alcohol-related outcomes and on its pharmacogenetic moderation during a 6-month post-treatment period. This analysis includes 308 individuals with problematic alcohol use (67% male; mean age = 51.1; topiramate: 49%, placebo: 51%). It uses generalized linear mixed models to examine changes in self-reported alcohol consumption and alcohol-related problems and concentrations of the liver enzyme γ-glutamyltransferase. The report combines published 3- and 6-month follow-up data from Study 1 with similar, unpublished data from Study 2. Despite robust effects of topiramate on drinking during treatment, the overall multivariate medication effects on outcomes during 3- and 6-month follow-up were not significant (p = 0.08 and p = 0.26, respectively). The moderating effect of the SNP on primary treatment outcomes was also not significant during either follow-up period (p = 0.13 and p = 0.16, respectively). However, during the 3-month post-treatment period, drinks per day was significantly lower in the topiramate group than the placebo group in the rs2832407*CC-genotype group. The robust effects of topiramate on alcohol-related outcomes during treatment diminish substantially once the medication is discontinued. Research is needed both to determine the optimal treatment duration and to identify clinically useful pharmacogenetic moderators of topiramate for treating AUD.

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.

Single nuclei RNA sequencing of the rat AP and NTS following GDF15 treatment

Objective

Growth differentiation factor 15 (GDF15) is known to play a role in feeding, nausea, and body weight, with action through the GFRAL-RET receptor complex in the area postrema (AP) and nucleus tractus solitarius (NTS). To further elucidate the underlying cell type-specific molecular mechanisms downstream of GDF15 signaling, we used a single nuclei RNA sequencing (snRNAseq) approach to profile AP and NTS cellular subtype-specific transcriptomes after systemic GDF15 treatment.

Methods

AP and NTS micropunches were used for snRNAseq from Sprague Dawley rats 6 h following GDF15 or saline injection, and Seurat was used to identify cellular subtypes and cell type-specific alterations in gene expression that were due to the direct and secondary effects of systemic GDF15 treatment.

Results

Using the transcriptome profile of ∼35,000 individual AP/NTS nuclei, we identified 19 transcriptomically distinct cellular subtypes, including a single population Gfral and Ret positive excitatory neurons, representing the primary site of action for GDF15. A total of ∼600 cell type-specific differential expression events were identified in neurons and glia, including the identification of transcriptome alterations specific to the direct effects of GDF15 in the Gfral-Ret positive excitatory neurons and shared transcriptome alterations across neuronal and glial cell types. Downstream analyses identified shared and cell type-specific alterations in signaling pathways and upstream regulatory mechanisms of the observed transcriptome alterations.

Conclusions

These data provide a considerable advance in our understanding of AP and NTS cell type-specific molecular mechanisms associated with GDF15 signaling. The identified cellular subtype-specific regulatory mechanism and signaling pathways likely represent important targets for future pharmacotherapies.

•

GDF15 directly alters transcription in Gfral- and Ret-positive excitatory neurons.

•

GDF15 indirectly alters transcription in other neuronal and glial populations.

•

Cell type-specific expression changes identify regulatory and signaling mechanisms.

GIP Receptor Agonism Attenuates GLP-1 Receptor Agonist-Induced Nausea and Emesis in Preclinical Models

Glucagon-like peptide 1 receptor (GLP-1R) agonists decrease body weight and improve glycemic control in obesity and diabetes. Patient compliance and maximal efficacy of GLP-1 therapeutics are limited by adverse side effects, including nausea and emesis. In three different species (i.e., mice, rats, and musk shrews), we show that glucose-dependent insulinotropic polypeptide receptor (GIPR) signaling blocks emesis and attenuates illness behaviors elicited by GLP-1R activation, while maintaining reduced food intake, body weight loss, and improved glucose tolerance. The area postrema and nucleus tractus solitarius (AP/NTS) of the hindbrain are required for food intake and body weight suppression by GLP-1R ligands and processing of emetic stimuli. Using single-nuclei RNA sequencing, we identified the cellular phenotypes of AP/NTS cells expressing GIPR and GLP-1R on distinct populations of inhibitory and excitatory neurons, with the greatest expression of GIPR in γ-aminobutyric acid-ergic neurons. This work suggests that combinatorial pharmaceutical targeting of GLP-1R and GIPR will increase efficacy in treating obesity and diabetes by reducing nausea and vomiting.

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.

Analysis of genetic and clinical factors associated with buprenorphine response

BACKGROUND

Buprenorphine, approved for treating opioid use disorder (OUD), is not equally efficacious for all patients. Candidate gene studies have shown limited success in identifying genetic moderators of buprenorphine treatment response.

METHODS

We studied 1616 European-ancestry individuals enrolled in the Million Veteran Program, of whom 1609 had an ICD-9/10 code consistent with OUD, a 180-day buprenorphine treatment exposure, and genome-wide genotype data. We conducted a genome-wide association study (GWAS) of buprenorphine treatment response [defined as having no opioid-positive urine drug screens (UDS) following the first prescription]. We also examined correlates of buprenorphine treatment response in multivariable analyses.

RESULTS

Although no variants reached genome-wide significance, 6 loci were nominally significant (p < 1 × 10^-5), four of which were located near previously characterized genes: rs756770 (ADAMTSL2), rs11782370 (SLC25A37), rs7205113 (CRISPLD2), and rs13169373 (LINC01947). A higher maximum daily buprenorphine dosage (aOR = 0.98; 95 %CI: 0.97, 0.995), greater number of UDS (aOR = 0.97; 95 %CI: 0.96, 0.99), and history of hepatitis C (HCV) infection (aOR = 0.71; 95 %CI: 0.57, 0.88) were associated with a reduced odds of buprenorphine response. Older age (aOR: 1.01; 95 %CI: 1.000, 1.02) was associated with increased odds of buprenorphine response.

CONCLUSIONS

This study had limited statistical power to detect genetic variants associated with a complex human phenotype like buprenorphine treatment response. Meta-analysis of multiple data sets is needed to ensure adequate statistical power for a GWAS of buprenorphine treatment response. The most robust phenotypic predictor of buprenorphine treatment response was intravenous drug use, a proxy for which was HCV infection.

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.

Combined analysis of the moderating effect of a GRIK1 polymorphism on the effects of topiramate for treating alcohol use disorder

BACKGROUND

In an initial study, we reported that topiramate reduced heavy drinking among individuals who sought to reduce their drinking and that the effect was moderated by a single nucleotide polymorphism (SNP; rs2832407) in GRIK1, which encodes the kainate GluK1 receptor subunit (Kranzler et al., 2014). In a subsequent study that prospectively randomized patients to medication group based on their rs2832407 genotype, we replicated the main effect of topiramate but not the moderating effect of the SNP (Kranzler et al., 2021). Given the similar design of the two studies, here we combined the findings to provide greater statistical power to test the pharmacogenetic effect.

MATERIAL AND METHODS

This secondary analysis of two 12-week, randomized controlled trials of topiramate included a total of 292 European-ancestry individuals (67.1 % male; topiramate: 48.3 %, placebo: 51.7 %) with problematic alcohol use. Using MANOVA, we examined changes in self-reported alcohol consumption, problems resulting from alcohol use, and quality of life, and the biomarker γ-glutamyltransferase. To test the pharmacogenetic hypothesis, all patients were genotyped for rs2832407.

RESULTS

There was a significant overall effect of topiramate on the alcohol-related outcomes (partial η² = 0.134, p < 0.001), with follow-up analyses showing significant reductions in percent heavy drinking days (Cohen's d = 0.49), percent days abstinent (d = 0.23), drinks/day (d = 0.29) and alcohol-related problems (d = 0.45). Overall, the moderating effect of the SNP was non-significant (partial η² = 0.026, p = 0.37).

CONCLUSIONS

Although topiramate is an efficacious medication for reducing drinking and alcohol-related problems among patients with problematic alcohol use, rs2832407 does not appear to moderate its therapeutic effects. www.clinicaltrials.gov registrations: NCT00626925 and NCT02371889.

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.

Placental DNA methylation profiles in opioid-exposed pregnancies and associations with the neonatal opioid withdrawal syndrome

Opioid abuse during pregnancy can result in Neonatal Opioid Withdrawal Syndrome (NOWS). We investigated genome-wide methylation analyses of 96 placental tissue samples, including 32 prenatally opioid-exposed infants with NOWS who needed therapy (+Opioids/+NOWS), 32 prenatally opioid-exposed infants with NOWS who did not require treatment (+Opioids/-NOWS), and 32 prenatally unexposed controls (-Opioids/-NOWS, control). Statistics, bioinformatics, Artificial Intelligence (AI), including Deep Learning (DL), and Ingenuity Pathway Analyses (IPA) were performed. We identified 17 dysregulated pathways thought to be important in the pathophysiology of NOWS and reported accurate AI prediction of NOWS diagnoses. The DL had an AUC (95% CI) =0.98 (0.95-1.0) with a sensitivity and specificity of 100% for distinguishing NOWS from the +Opioids/-NOWS group and AUCs (95% CI) =1.00 (1.0-1.0) with a sensitivity and specificity of 100% for distinguishing NOWS versus control and + Opioids/-NOWS group versus controls. This study provides strong evidence of methylation dysregulation of placental tissue in NOWS development.

A Delta-Opioid Receptor Gene Polymorphism Moderates the Therapeutic Response to Extended-Release Buprenorphine in Opioid Use Disorder

BACKGROUND

Buprenorphine treatment is not equally effective in all patients with opioid use disorder (OUD). Two retrospective studies showed that, among African Americans (AAs), rs678849, a polymorphism in the delta-opioid receptor gene, moderated the therapeutic effect of sublingual buprenorphine.

METHODS

We examined rs678849 as a moderator of the response to an extended-release subcutaneous buprenorphine formulation (BUP-XR) in a 24-week OUD treatment study of 127 AAs and 327 European Americans (EAs). Participants were randomly assigned to receive: (1) BUP-XR as 2 monthly injections of 300 mg followed by either 300 mg monthly or 100 mg monthly for 4 months, or (2) monthly volume-matched placebo injections. Generalized estimating equations logistic regression analyses tested, per population group, the main and interaction effects of treatment (BUP-XR vs placebo) and genotype group (rs678849*CC vs CT/TT) on weekly urine drug screens (UDS).

RESULTS

Among AAs, the placebo group had higher rates of opioid-positive UDS than the BUP-XR group (log odds ratio = 1.67, 95% CI = 0.36, 2.98), but no genotype by treatment effect (P = .80). Among EAs, the placebo group also showed higher rates of opioid-positive UDS than the BUP-XR group (log odds ratio = 1.97, 95% CI = 1.14, 2.79) but a significant genotype by treatment interaction (χ 2(1) = 4.33, P = .04).

CONCLUSION

We found a moderating effect of rs678849 on the response to buprenorphine treatment of OUD in EAs, but not AAs. These findings require replication in well-powered, prospective studies of both AA and EA OUD patients treated with BUP-XR and stratified on rs678849 genotype.

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.

Assessment of Probable Opioid Use Disorder Using Electronic Health Record Documentation

IMPORTANCE

Electronic health records are a potentially valuable source of information for identifying patients with opioid use disorder (OUD).

OBJECTIVE

To evaluate whether proxy measures from electronic health record data can be used reliably to identify patients with probable OUD based on Diagnostic and Statistical Manual of Mental Disorders (Fifth Edition) (DSM-5) criteria.

DESIGN, SETTING, AND PARTICIPANTS

This retrospective cross-sectional study analyzed individuals within the Geisinger health system who were prescribed opioids between December 31, 2000, and May 31, 2017, using a mixed-methods approach. The cohort was identified from 16 253 patients enrolled in a contract-based, Geisinger-specific medication monitoring program (GMMP) for opioid use, including patients who maintained or violated contract terms, as well as a demographically matched control group of 16 253 patients who were prescribed opioids but not enrolled in the GMMP. Substance use diagnoses and psychiatric comorbidities were assessed using automated electronic health record summaries. A manual medical record review procedure using DSM-5 criteria for OUD was completed for a subset of patients. The analysis was conducted beginning from June 5, 2017, until May 29, 2020.

MAIN OUTCOMES AND MEASURES

The primary outcome was the prevalence of OUD as defined by proxy measures for DSM-5 criteria for OUD as well as the prevalence of comorbidities among patients prescribed opioids within an integrated health system.

RESULTS

Among the 16 253 patients enrolled in the GMMP (9309 women [57%]; mean [SD] age, 52 [14] years), OUD diagnoses as defined by diagnostic codes were present at a much lower rate than expected (291 [2%]), indicating the necessity for alternative diagnostic strategies. The DSM-5 criteria for OUD can be assessed using manual medical record review; a manual review of 200 patients in the GMMP and 200 control patients identifed a larger percentage of patients with probable moderate to severe OUD (GMMP, 145 of 200 [73%]; and control, 27 of 200 [14%]) compared with the prevalence of OUD assessed using diagnostic codes.

CONCLUSIONS AND RELEVANCE

These results suggest that patients with OUD may be identified using information available in the electronic health record, even when diagnostic codes do not reflect this diagnosis. Furthermore, the study demonstrates the utility of coding for DSM-5 criteria from medical records to generate a quantitative DSM-5 score that is associated with OUD severity.

Evolving Magnetically Levitated Plasma Proteins Detects Opioid Use Disorder as a Model Disease

There are several methods (e.g., enzyme-linked immunosorbent assay and liquid chromatography mass spectroscopy) that already use human plasma to detect a variety of possible diseases. However, this paper introduces the capabilities of magnetic levitation (Maglev) to detect disease (Opioid Use Disorder, used here as a model disease) by using levitation of human plasma proteins. The presented proof-of-concept findings revealed that the optical images of magnetically levitated plasma proteins carry important information about the health spectrum of plasma donors. In addition, the liquid chromatography mass spectroscopy analysis of the magnetically levitated plasma proteins demonstrated remarkable differences between the plasma of healthy individuals and patients with opioid use disorders. Overall, the presented method provides diagnostic value for disease detection using optical images of evolving magnetically levitated plasma proteins and/or proteomic information.

The importance of buprenorphine research in the opioid crisis

With the urgency to treat patients more effectively for opioid use disorder in the midst of the opioid epidemic, a key area for precision medicine is to improve individualized medication-assisted treatment for opioid use disorder. The expansion of medication-assisted treatment is a key to reducing illicit opioid use, preventing opioid overdose deaths, and reducing the comorbidities and societal impacts of opioid use disorder. The most common medication for opioid use disorder will soon be buprenorphine. Research to date shows the successful impact of buprenorphine treatment, including the pharmacogenomics of buprenorphine response and treatment efficacy. Buprenorphine is also a promising treatment for depression and anxiety, and neonatal opioid withdrawal syndrome (NOWS). However, the rates of success with medication-assisted treatment for opioid use disorder, particularly at the beginning of treatment, still show many individuals relapsing to illicit opioid use. With the scope of the opioid crisis, there is an urgent need for expansion of buprenorphine treatment research to provide critical information for improving outcomes of opioid use disorder. Implementing the best strategies for opioid use disorder treatment is of dire urgency and will save lives.

Genome-Wide Association Study Meta-Analysis of the Alcohol Use Disorders Identification Test (AUDIT) in Two Population-Based Cohorts

OBJECTIVE

Alcohol use disorders are common conditions that have enormous social and economic consequences. Genome-wide association analyses were performed to identify genetic variants associated with a proxy measure of alcohol consumption and alcohol misuse and to explore the shared genetic basis between these measures and other substance use, psychiatric, and behavioral traits.

METHOD

This study used quantitative measures from the Alcohol Use Disorders Identification Test (AUDIT) from two population-based cohorts of European ancestry (UK Biobank [N=121,604] and 23andMe [N=20,328]) and performed a genome-wide association study (GWAS) meta-analysis. Two additional GWAS analyses were performed, a GWAS for AUDIT scores on items 1-3, which focus on consumption (AUDIT-C), and for scores on items 4-10, which focus on the problematic consequences of drinking (AUDIT-P).

RESULTS

The GWAS meta-analysis of AUDIT total score identified 10 associated risk loci. Novel associations localized to genes including JCAD and SLC39A13; this study also replicated previously identified signals in the genes ADH1B, ADH1C, KLB, and GCKR. The dimensions of AUDIT showed positive genetic correlations with alcohol consumption (r_g=0.76-0.92) and DSM-IV alcohol dependence (r_g=0.33-0.63). AUDIT-P and AUDIT-C scores showed significantly different patterns of association across a number of traits, including psychiatric disorders. AUDIT-P score was significantly positively genetically correlated with schizophrenia (r_g=0.22), major depressive disorder (r_g=0.26), and attention deficit hyperactivity disorder (r_g=0.23), whereas AUDIT-C score was significantly negatively genetically correlated with major depressive disorder (r_g=-0.24) and ADHD (r_g=-0.10). This study also used the AUDIT data in the UK Biobank to identify thresholds for dichotomizing AUDIT total score that optimize genetic correlations with DSM-IV alcohol dependence. Coding individuals with AUDIT total scores ≤4 as control subjects and those with scores ≥12 as case subjects produced a significant high genetic correlation with DSM-IV alcohol dependence (r_g=0.82) while retaining most subjects.

CONCLUSIONS

AUDIT scores ascertained in population-based cohorts can be used to explore the genetic basis of both alcohol consumption and alcohol use disorders.

A review of opioid addiction genetics

Opioid use disorder (OUD) affects millions of people worldwide and the risk of developing the disorder has a significant genetic component according to twin and family studies. Identification of the genetic variants underlying this inherited risk has focused on two different methods: candidate gene studies and genome-wide association studies (GWAS). The most studied candidate genes have included the mu-opioid receptor (OPRM1), the delta-opioid receptor (OPRD1), the dopamine D2 receptor (DRD2), and brain-derived neurotrophic factor (BDNF). Variants in these genes have been associated with relatively small, but reproducible, effects on OUD risk. More recently, GWAS have identified potential associations with variants in KCNG2, KCNC1, CNIH3, APBB2, and RGMA. In total the genetic associations identified so far explain only a small portion of OUD risk. GWAS of OUD is still in the early stages when compared to studies of other psychiatric disorders, such as schizophrenia, which have found many relevant variants with small effect sizes only after large meta-analyses. Substantial increases in cohort sizes will likely be necessary in the OUD field to achieve similar results. In addition, it will be important for future studies of OUD to incorporate rare variants, epigenetics, and gene × environment interactions into models in order to better explain the observed heritability.

Pharmacogenetic analysis of opioid dependence treatment dose and dropout rate

BACKGROUND

Currently, no pharmacogenetic tests for selecting an opioid-dependence pharmacotherapy have been approved by the US Food and Drug Administration.

OBJECTIVES

Determine the effects of variants in 11 genes on dropout rate and dose in patients receiving methadone or buprenorphine/naloxone (ClinicalTrials.gov Identifier: NCT00315341).

METHODS

Variants in six pharmacokinetic genes (CYP1A2, CYP2B6, CYP2C19, CYP2C9, CYP2D6, CYP3A4) and five pharmacodynamic genes (HTR2A, OPRM1, ADRA2A, COMT, SLC6A4) were genotyped in samples from a 24-week, randomized, open-label trial of methadone and buprenorphine/naloxone for the treatment of opioid dependence (n = 764; 68.7% male). Genotypes were then used to determine the metabolism phenotype for each pharmacokinetic gene. Phenotypes or genotypes for each gene were analyzed for association with dropout rate and mean dose.

RESULTS

Genotype for 5-HTTLPR in the SLC6A4 gene was nominally associated with dropout rate when the methadone and buprenorphine/naloxone groups were combined. When the most significant variants associated with dropout rate were analyzed using pairwise analyses, SLC6A4 (5-HTTLPR) and COMT (Val158Met; rs4860) had nominally significant associations with dropout rate in methadone patients. None of the genes analyzed in the study was associated with mean dose of methadone or buprenorphine/naloxone.

CONCLUSIONS

This study suggests that functional polymorphisms related to synaptic dopamine or serotonin levels may predict dropout rates during methadone treatment. Patients with the S/S genotype at 5-HTTLPR in SLC6A4 or the Val/Val genotype at Val158Met in COMT may require additional treatment to improve their chances of completing addiction treatment. Replication in other methadone patient populations will be necessary to ensure the validity of these findings.

Reading LINEs within the cocaine addicted brain

INTRODUCTION

Long interspersed element (LINE)-1 (L1) is a type of retrotransposon capable of mobilizing into new genomic locations. Often studied in Mendelian diseases or cancer, L1s may also cause somatic mutation in the developing central nervous system. Recent reports showed L1 transcription was activated in brains of cocaine-treated mice, and L1 retrotransposition was increased in cocaine-treated neuronal cell cultures. We hypothesized that the predisposition to cocaine addiction may result from inherited L1s or somatic L1 mobilization in the brain.

METHODS

Postmortem medial prefrontal cortex (mPFC) tissue from 30 CA and 30 control individuals was studied. An Alexafluor488-labeled NeuN antibody and fluorescence activated nuclei sorting were used to separate neuronal from non-neuronal cell nuclei. L1s and their 3' flanking sequences were amplified from neuronal and non-neuronal genomic DNA (gDNA) using L1-seq. L1 DNA libraries from the neuronal gDNA were sequenced on an Illumina HiSeq2000. Sequences aligned to the hg19 human genome build were analyzed for L1 insertions using custom "L1-seq" bioinformatics programs.

RESULTS

Previously uncataloged L1 insertions, some validated by PCR, were detected in neurons from both CA and control brain samples. Steady-state L1 mRNA levels in CA and control mPFC were also assessed. Gene ontology and pathway analyses were used to assess relationships between genes putatively disrupted by novel L1s in CA and control individuals. L1 insertions in CA samples were enriched in gene ontologies and pathways previously associated with CA.

CONCLUSIONS

We conclude that neurons in the mPFC harbor L1 insertions that have the potential to influence predisposition to CA.

OPRD1 Genetic Variation and Human Disease

The OPRD1 gene encodes the delta-opioid receptor, which has multiple functions including regulating reward pathways. The gene contains more than 2,000 verified genetic variants but only 2 currently have evidence for specific functions: rs1042114 disrupts maturation of the receptor and rs569356 affects OPRD1 expression. These polymorphisms and others in the gene have been found to be associated with human diseases. The most reproducible data are associations between opioid addiction and three variants in intron 1 (rs2236861, rs2236857, and rs3766951), which have been described in a number of independent populations. Several publications also point toward an association between anorexia and a haplotype block containing rs569356 and rs533123. Unfortunately the mechanisms underlying these two effects are currently unknown. In contrast, rs1042114 has been linked to Alzheimer's disease through an increasingly well-defined mechanism by which the variant allele reduces production of the beta-amyloid plaques associated with the disease. Additional studies of OPRD1 variants are necessary to replicate current findings and to delineate the functional roles of relevant polymorphisms.

A polymorphism in the OPRM1 3'-untranslated region is associated with methadone efficacy in treating opioid dependence

The mu-opioid receptor (MOR) is the primary target of methadone and buprenorphine. The primary neuronal transcript of the OPRM1 gene, MOR-1, contains a ~13kb 3′ untranslated region with five common haplotypes in European-Americans. We analyzed the effects of these haplotypes on the percentage of opioid positive urine tests in European-Americans (n = 582) during a 24-week, randomized, open-label trial of methadone or buprenorphine/naloxone (Suboxone) for the treatment of opioid dependence. A single haplotype, tagged by rs10485058, was significantly associated with patient urinalysis data in the methadone treatment group. Methadone patients with the A/A genotype at rs10485058 were less likely to have opioid-positive urine drug screens than those in the combined A/G and G/G genotypes group (Relative Risk = 0.76, 95% confidence intervals = 0.73–0.80, p = 0.0064). Genotype at rs10485058 also predicted self-reported relapse rates in an independent population of Australian patients of European descent (n = 1215) who were receiving opioid substitution therapy (p = 0.003). In silico analysis predicted that miR-95-3p would interact with the G, but not the A allele of rs10485058. Luciferase assays indicated miR-95-3p decreased reporter activity of constructs containing the G, but not the A allele of rs10485058, suggesting a potential mechanism for the observed pharmacogenetic effect. These findings suggest that selection of a medication for opioid dependence based on rs10485058 genotype might improve outcomes in this ethnic group.

Abstract B36: Modifier genes and cancer susceptibility: Lessons from inbred mouse strains

The adenomatous polyposis coli (APC) tumor suppressor gene is mutated in Familial Adenomatous Polyposis (FAP), an inherited disorder that predisposes individuals to developing intestinal polyposis and eventually leads to cancer. In the Min mouse model, the genetic background of mice carrying a mutation in the murine homolog of the APC gene (ApcMin) is critical to the manifestation of tumor phenotypes, as inbred strains have been demonstrated to differ dramatically in their susceptibility to polyposis. Complex trait analyses from several laboratories, including our own, have identified modifier loci that alter intestinal tumor number and size in ApcMin/+ mice, but less than a handful of modifier genes have been identified to date. We review the results of several quantitative trait loci (QTL) studies that have exploited the genomic diversity among inbred mouse strains and report the identification of novel modifier loci that protect against polyposis. We compare the results with cancer susceptibility loci identified through GWAS and other studies in the human genome. The usefulness of using the mouse as a model organism for optimizing complex trait screens for cancer phenotypes and speeding the process of causative gene identification is discussed. Research supported in part by NCI grants to LDS.

Citation Format: Stephanie Nnadi, Xiang Wang, Richard C. Crist, Arthur M. Buchberg, Linda D. Siracusa. Modifier genes and cancer susceptibility: Lessons from inbred mouse strains. [abstract]. In: Proceedings of the AACR Special Conference: The Translational Impact of Model Organisms in Cancer; Nov 5-8, 2013; San Diego, CA. Philadelphia (PA): AACR; Mol Cancer Res 2014;12(11 Suppl):Abstract nr B36.

Pharmacogenetics of OPRM1

Pharmacogenetic research has the potential to explain the variation in treatment efficacy within patient populations. Understanding the interaction between genetic variation and medications may provide a method for matching patients to the most effective therapeutic options and improving overall patient outcomes. The OPRM1 gene has been a target of interest in a large number of pharmacogenetic studies due to its genetic and structural variation, as well as the role of opioid receptors in a variety of disorders. The mu-opioid receptor (MOR), encoded by OPRM1, naturally regulates the analgesic response to pain and also controls the rewarding effects of many drugs of abuse, including opioids, nicotine, and alcohol. Genetic variants in OPRM1, particularly the non-synonymous polymorphism A118G, have been repeatedly associated with the efficacy of treatments for pain and various types of dependence. This review focuses on the current understanding of the pharmacogenetic impact of OPRM1, primarily with regard to the treatment of pain and addiction.

Low frequency genetic variants in the μ-opioid receptor (OPRM1) affect risk for addiction to heroin and cocaine

The μ-opioid receptor (MOR) binds exogenous and endogenous opioids and is known to mediate the rewarding effects of drugs of abuse. Numerous genetic studies have sought to identify common genetic variation in the gene encoding MOR (OPRM1) that affects risk for drug addiction. The purpose of this study was to examine the contribution of rare coding variants in OPRM1 to the risk for addiction. Rare and low frequency variants were selected using the National Heart Lung and Blood Institute - Exome Sequencing Project (NHLBI-ESP) database, which has screened the exomes of over 6500 individuals. Two SNPs (rs62638690 and rs17174794) were selected for genotyping in 1377 European American individuals addicted to heroin and/or cocaine. Two different SNPs (rs1799971 and rs17174801) were genotyped in 1238 African American individuals addicted to heroin and/or cocaine. Using the minor allele frequencies from the NHLBI-ESP dataset as a comparison group, case-control association analyses were performed. Results revealed an association between rs62638690 and cocaine and heroin addiction in European Americans (p=0.02; 95% C.I. 0.47 [0.24-0.92]). This study suggests a potential role for rare OPRM1 variants in addiction disorders and highlights an area worthy of future study.

Case-control association analysis of polymorphisms in the δ-opioid receptor, OPRD1, with cocaine and opioid addicted populations

BACKGROUND

Addiction susceptibility and treatment responsiveness are greatly influenced by genetic factors. Sequence variation in genes involved in the mechanisms of drug action have the potential to influence addiction risk and treatment outcome. The opioid receptor system is involved in mediating the rewarding effects of cocaine and opioids. The μ-opioid receptor (MOR) has traditionally been considered the primary target for opioid addiction. The MOR, however, interacts with and is regulated by many known MOR interacting proteins (MORIPs), including the δ-opioid receptor (DOR).

METHODS

The present study evaluated the contribution of OPRD1, the gene encoding the DOR, to the risk of addiction to opioids and cocaine. The association of OPRD1 polymorphisms with both opioid addiction (OA) and cocaine addiction (CA) was analyzed in African American (OA n=336, CA n=503) and European American (OA n=1007, CA n=336) populations.

RESULTS

The primary finding of this study is an association of rs678849 with cocaine addiction in African Americans (allelic p=0.0086). For replication purposes, this SNP was analyzed in a larger independent population of cocaine addicted African Americans and controls and the association was confirmed (allelic p=4.53 × 10(-5); n=993). By performing a meta-analysis on the expanded populations, the statistical evidence for an association was substantially increased (allelic p=8.5 × 10(-7)) (p-values non-FDR corrected).

CONCLUSION

The present study suggests that polymorphisms in OPRD1 are relevant for cocaine addiction in the African American population and provides additional support for a broad role for OPRD1 variants in drug dependence.

A conserved tissue-specific homeodomain-less isoform of MEIS1 is downregulated in colorectal cancer

Colorectal cancer is one of the most common cancers in developed nations and is the result of both environmental and genetic factors. Many of the genetic lesions observed in colorectal cancer alter expression of homeobox genes, which encode homeodomain transcription factors. The MEIS1 homeobox gene is known to be involved in several hematological malignancies and solid tumors and recent evidence suggests that expression of the MEIS1 transcript is altered in colorectal cancer. Despite this potential connection, little is known about the role of the gene in the intestines. We probed murine gastrointestinal tissue samples with an N-terminal Meis1 antibody, revealing expression of two previously described isoforms, as well as two novel Meis1 products. A 32 kD Meis1 product was expressed in the nuclei of non-epithelial cells in the stomach and colon, while a 27 kD product was expressed in the cytoplasm of epithelial cells in the proximal colon. Our data suggest that the 27 kD and 32 kD Meis1 proteins are both forms of the Meis1d protein, a homeodomain-less isoform whose transcript was previously identified in cDNA screens. Both the MEIS1D transcript and protein were expressed in human colon mucosa. Expression of the MEIS1D protein was downregulated in 83% (10/12) of primary colorectal cancer samples compared to matched normal mucosa, indicating that MEIS1D is a biomarker of colorectal tumorigenesis. The decreased expression of MEIS1D in colon tumors also suggests that this conserved homeodomain-less isoform may act as a tumor suppressor in human colorectal cancer.

Identification of Mom12 and Mom13, two novel modifier loci of Apc (Min) -mediated intestinal tumorigenesis

Colorectal cancer is a heterogeneous disease resulting from a combination of genetic and environmental factors. The C57BL/6J (B6) Apc (Min/+) mouse develops polyps throughout the gastrointestinal tract and has been a valuable model for understanding the genetic basis of intestinal tumorigenesis. Apc (Min/+) mice have been used to study known oncogenes and tumor suppressor genes on a controlled genetic background. These studies often utilize congenic knockout alleles, which can carry an unknown amount of residual donor DNA. The Apc (Min) model has also been used to identify modifer loci, known as Modifier of Min (Mom) loci, which alter Apc (Min) -mediated intestinal tumorigenesis. B6 mice carrying a knockout allele generated in WW6 embryonic stem cells were crossed to B6 Apc (Min/+) mice to determine the effect on polyp multiplicity. The newly generated colony developed significantly more intestinal polyps than Apc (Min/+) controls. Polyp multiplicity did not correlate with inheritance of the knockout allele, suggesting the presence of one or more modifier loci segregating in the colony. Genotyping of simple sequence length polymorphism (SSLP) markers revealed residual 129X1/SvJ genomic DNA within the congenic region of the parental knockout line. An analysis of polyp multiplicity data and SSLP genotyping indicated the presence of two Mom loci in the colony: 1) Mom12, a dominant modifier linked to the congenic region on chromosome 6, and 2) Mom13, which is unlinked to the congenic region and whose effect is masked by Mom12. The identification of Mom12 and Mom13 demonstrates the potential problems resulting from residual heterozygosity present in congenic lines.

The armadillo repeat domain of Apc suppresses intestinal tumorigenesis

The adenomatous polyposis coli (APC) gene is known to act as a tumor suppressor gene in both sporadic and hereditary colorectal cancer by negatively regulating WNT signaling. Familial adenomatous polyposis (FAP) patients develop intestinal polyps due to the presence of a single germline mutation in APC. The severity of the FAP phenotype is a function of the position of the APC mutation, indicating a complex role for APC that extends beyond the canonical WNT pathway. APC encodes a large protein with multiple functional domains, including an armadillo repeat domain that has been linked to protein-protein interactions. To determine the effect of the armadillo repeat domain on intestinal tumorigenesis, we generated a congenic mouse line (Apc ( Δ242 )) carrying a gene trap cassette between exons 7 and 8 of the murine Apc gene. Apc ( Δ242/+) mice express a truncated Apc product lacking the armadillo repeat domain as part of a fusion protein with β-geo. Expression of the fusion product was confirmed by X-gal staining, ensuring that Apc ( Δ242 ) is not a null allele. In contrast, Apc ( Min/+) mice produce a truncated Apc product that contains an intact armadillo repeat domain. On the C57BL/6J background, Apc ( Δ242/+) mice develop more polyps than do Apc ( Min/+) mice along the entire length of the small intestine; however, polyps were significantly smaller in Apc ( Δ242/+) mice. In addition, polyp multiplicity in Apc ( Δ242/+) mice is affected by polymorphisms between inbred strains. These data suggest that the armadillo repeat domain of the Apc protein suppresses tumor initiation in the murine intestine while also promoting tumor growth.