Skip to main content
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian
  • Log in
  • My Cart

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses

New Research In

Physical Sciences

Featured Portals

  • Physics
  • Chemistry
  • Sustainability Science

Articles by Topic

  • Applied Mathematics
  • Applied Physical Sciences
  • Astronomy
  • Computer Sciences
  • Earth, Atmospheric, and Planetary Sciences
  • Engineering
  • Environmental Sciences
  • Mathematics
  • Statistics

Social Sciences

Featured Portals

  • Anthropology
  • Sustainability Science

Articles by Topic

  • Economic Sciences
  • Environmental Sciences
  • Political Sciences
  • Psychological and Cognitive Sciences
  • Social Sciences

Biological Sciences

Featured Portals

  • Sustainability Science

Articles by Topic

  • Agricultural Sciences
  • Anthropology
  • Applied Biological Sciences
  • Biochemistry
  • Biophysics and Computational Biology
  • Cell Biology
  • Developmental Biology
  • Ecology
  • Environmental Sciences
  • Evolution
  • Genetics
  • Immunology and Inflammation
  • Medical Sciences
  • Microbiology
  • Neuroscience
  • Pharmacology
  • Physiology
  • Plant Biology
  • Population Biology
  • Psychological and Cognitive Sciences
  • Sustainability Science
  • Systems Biology
Research Article

Suppression of RGSz1 function optimizes the actions of opioid analgesics by mechanisms that involve the Wnt/β-catenin pathway

Sevasti Gaspari, Immanuel Purushothaman, Valeria Cogliani, Farhana Sakloth, Rachael L. Neve, David Howland, Robert H. Ring, Elliott M. Ross, Li Shen, and View ORCID ProfileVenetia Zachariou
PNAS February 27, 2018 115 (9) E2085-E2094; first published February 12, 2018; https://doi.org/10.1073/pnas.1707887115
Sevasti Gaspari
aFishberg Department of Neuroscience, Icahn School of Medicine at Mount Sinai, New York, NY 10029;
bFriedman Brain Institute, Icahn School of Medicine at Mount Sinai, New York, NY 10029;
cDepartment of Basic Sciences, University of Crete Medical School, 71003 Heraklion, Greece;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Immanuel Purushothaman
aFishberg Department of Neuroscience, Icahn School of Medicine at Mount Sinai, New York, NY 10029;
bFriedman Brain Institute, Icahn School of Medicine at Mount Sinai, New York, NY 10029;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Valeria Cogliani
aFishberg Department of Neuroscience, Icahn School of Medicine at Mount Sinai, New York, NY 10029;
bFriedman Brain Institute, Icahn School of Medicine at Mount Sinai, New York, NY 10029;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Farhana Sakloth
aFishberg Department of Neuroscience, Icahn School of Medicine at Mount Sinai, New York, NY 10029;
bFriedman Brain Institute, Icahn School of Medicine at Mount Sinai, New York, NY 10029;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Rachael L. Neve
dGene Delivery Technology Core, Massachusetts General Hospital, MA 01239;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
David Howland
eCure Huntington's Disease Initiative (CHDI) Foundation, Princeton, NJ 08540;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Robert H. Ring
fDepartment of Psychiatry, Icahn School of Medicine at Mount Sinai, New York, NY 10029;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Elliott M. Ross
gDepartment of Pharmacology, University of Texas Southwestern Medical Center, Dallas, TX 75390;
hGreen Center for Systems Biology, University of Texas Southwestern Medical Center, Dallas, TX 75390
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Li Shen
aFishberg Department of Neuroscience, Icahn School of Medicine at Mount Sinai, New York, NY 10029;
bFriedman Brain Institute, Icahn School of Medicine at Mount Sinai, New York, NY 10029;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Venetia Zachariou
aFishberg Department of Neuroscience, Icahn School of Medicine at Mount Sinai, New York, NY 10029;
bFriedman Brain Institute, Icahn School of Medicine at Mount Sinai, New York, NY 10029;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Venetia Zachariou
  • For correspondence: venetia.zachariou@mssm.edu
  1. Edited by Solomon H. Snyder, Johns Hopkins University School of Medicine, Baltimore, MD, and approved January 5, 2018 (received for review June 23, 2017)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Significance

Opioids are used to alleviate severe pain, but their long-term use leads to analgesic tolerance, dependence, and addiction. Here, we targeted specific intracellular pathways to dissociate the analgesic actions of opioids from addiction-related effects. Using genetically modified male and female mice in models of addiction and analgesia, we revealed a key role of an intracellular modulator of the mu opioid receptor, RGSz1, in opioid actions. We applied next-generation sequencing and biochemical assays to delineate the mechanism of RGSz1 action in the mouse periaqueductal gray. Findings from this work point to novel intracellular pathways that can be targeted to optimize the actions of opioids for the treatment of chronic pain.

Abstract

Regulator of G protein signaling z1 (RGSz1), a member of the RGS family of proteins, is present in several networks expressing mu opioid receptors (MOPRs). By using genetic mouse models for global or brain region-targeted manipulations of RGSz1 expression, we demonstrated that the suppression of RGSz1 function increases the analgesic efficacy of MOPR agonists in male and female mice and delays the development of morphine tolerance while decreasing the sensitivity to rewarding and locomotor activating effects. Using biochemical assays and next-generation RNA sequencing, we identified a key role of RGSz1 in the periaqueductal gray (PAG) in morphine tolerance. Chronic morphine administration promotes RGSz1 activity in the PAG, which in turn modulates transcription mediated by the Wnt/β-catenin signaling pathway to promote analgesic tolerance to morphine. Conversely, the suppression of RGSz1 function stabilizes Axin2–Gαz complexes near the membrane and promotes β-catenin activation, thereby delaying the development of analgesic tolerance. These data show that the regulation of RGS complexes, particularly those involving RGSz1-Gαz, represents a promising target for optimizing the analgesic actions of opioids without increasing the risk of dependence or addiction.

  • morphine
  • RNA sequencing
  • periaqueductal gray
  • G proteins
  • analgesia

Mu opioid receptor (MOPR) agonists, such as morphine, are prescribed for the treatment of severe acute or chronic pain conditions, but their long-term use is hindered by numerous side effects and by the development of dependence and addiction (1⇓⇓–4). Recent evidence highlights the rise of physical dependence and addiction in patients treated with opioid analgesics (5, 6). As research efforts aim to develop more efficient and safer morphine-like compounds (6, 7) and opioid alternatives, a more detailed understanding of the mechanisms mediating the analgesic effects of MOPR agonists is needed. Although information on the cellular and molecular adaptations associated with repeated exposure to opioids is emerging (8⇓⇓⇓⇓–13), the detailed brain region- and cell type-specific signal transduction events triggered by MOPR activation are not fully understood.

RGSz1, one product of the RGS20 gene (14), is a member of the Rz subfamily [composed of GAIP (RGS19), RGSz1, RGSz2 (RGS17)] of regulator of G protein signaling (RGS) proteins (15⇓–17) and is expressed in several brain regions (18⇓–20). RGSz1 accelerates the hydrolysis of GTP-bound Gαz more than 600-fold and may modulate the GTP hydrolysis of other Gαi subunits (18, 19). As with other RGS proteins, RGSz1 may serve as an effector antagonist for Gα subunits (21, 22). RGSz1 has been shown to modulate signal transduction from MOPRs and serotonin receptor 1A both in vitro (23, 24) and in vivo (25, 26) whereas knockout of the Gαz subunit accelerates the development of morphine tolerance (27, 28). Although a number of studies have documented brain region-specific effects of other RGS proteins (RGS9-2, RGS7, and RGS4) (29⇓⇓⇓–33) on the functional responses of MOPRs, the role of RGSz1 in opioid actions is not known. We hypothesized that the expression and activity of RGSz1 are altered by chronic opioid use and that the manipulation of RGSz1 expression will impact the development of tolerance to the drug and its analgesic and rewarding effects.

Here, we show that chronic morphine administration increases RGSz1 expression in the periaqueductal gray (PAG), a brain region that plays a key role in modulating endogenous analgesia and responses to MOPR agonists (34). Knockout of Rgsz1 (RGSz1KO) increased the analgesic efficacy of opioids and decreased the sensitivity to the rewarding and locomotor activating effects of morphine without impacting physical dependence. Global deletion or conditional knockdown of Rgsz1 in the ventrolateral (vl) PAG delayed the development of morphine tolerance in pain-free and chronic pain states. Coimmunoprecipitation (co-IP) assays indicated that, under states of morphine tolerance, there is an increased formation of complexes between Gαz and MOPR in the PAG. Next-generation RNA sequencing (RNA-Seq) revealed a robust down-regulation of Wnt/β-catenin signaling in the PAG of morphine-tolerant mice. Our Co-IP and Western blot analyses also suggested that RGSz1 competes with Axin2 (a member of the RGS protein family and a key component of the Wnt/β-catenin pathway that controls the nuclear translocation of β-catenin) for binding to Gαz. Using a morphine tolerance regimen, we showed that knockout of RGSz1 stabilizes Axin2/Gαz complexes, permitting the nuclear translocation of β-catenin and the transcription of target genes. Therefore, interfering with RGSz1 actions may be useful for promoting the analgesic effects of MOPR agonists while attenuating their addiction-related effects.

Results

Regulation of RGSz1 by Chronic Morphine.

We initially performed Western blot analyses to test the hypothesis that repeated morphine administration regulates RGSz1 levels in brain regions associated with analgesic responses. We applied a 4-d hot plate morphine tolerance paradigm and monitored analgesic responses for 4 consecutive days in groups of mice treated with morphine (15 mg/kg). RGSz1 protein levels in the PAG were up-regulated 2 h after the last morphine injection compared with that in saline-treated controls (Fig. 1A). No RGSz1 protein was detected in PAG samples from RGSz1KO mice. This morphine treatment regimen had no effect on RGSz1 expression in the nucleus accumbens (NAc), dorsal striatum (DS), or thalamus (THAL) (Fig. S1). Notably, morphine treatment up-regulated Gαz levels in the PAG of both RGSz1WT and RGSz1KO mice (Fig. 1B), but levels of other Gα subunits expressed in the PAG (Gαi3, Gαo, and Gαq) were not altered (Fig. 1 C–E). To further evaluate the association between RGSz1 and Gαz complexes in the PAG, we performed co-IP assays and found that morphine acutely (within 30 min) promoted the dissociation of Gαz from MOPRs (Fig. 1F). However, chronic morphine treatment induced the converse, promoting MOPR–Gαz complex formation (Fig. 1G).

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Morphine regulates RGSz1 and Gαz in the PAG. (A) RGSz1 protein levels were elevated in the PAG of RGSz1WT mice treated for 4 consecutive days with morphine [15 mg/kg; n = 16 per group; t test, t(30) = 2.144; *P < 0.05]. Chronic morphine also increased the levels of Gαz in the PAG (B), which occurred in both RGSz1WT and RGSz1KO mice [15 mg/kg; n = 11–12 per group; two-way ANOVA followed by Bonferroni’s post hoc test; F(1,42) = 11.94; *P < 0.05], but not the levels of Gαi3 [F(1,42) = 0.247] (C), Gαo [F(1,42) = 0.124] (D), or Gαq [F(1,43) = 0.366] (E) (dose: 15 mg/kg; n = 11–12 per group; two-way ANOVA followed by Bonferroni’s post hoc test; P > 0.05 for all three subunits). (F) Co-IP analysis revealed that Gαz dissociated from MOPRs 30 min after an acute morphine injection [15 mg/kg; n = 3–4 per group; two-way ANOVA followed by Bonferroni’s post hoc test, F(1,11) = 17.84; *P < 0.05]. (G) However, chronic treatment with morphine for 4 consecutive days promoted MOPR–Gαz complexes measured 30 min after the last injection [15 mg/kg; n = 4 per group; two-way ANOVA followed by Bonferroni’s post hoc test, F(1,12) = 32.94; *P < 0.05, ***P < 0.001]. M, morphine; S, saline; TL, total lysate.

RGSz1 Knockout Enhances the Analgesic Efficacy of MOPR Agonists and Delays Morphine Tolerance.

We next examined the role of RGSz1 in the analgesic actions of morphine and other MOPR agonists. In the hot plate test, male RGSz1KO mice show enhanced analgesia at lower doses of morphine compared with the responses of RGSz1WT controls (Fig. 2A). This phenotype was also observed in female RGSz1KO mice (Fig. 2A, Inset). The sensitivity to the analgesic effects of fentanyl and methadone was also increased in RGSz1KO mice compared with that of their RGSz1WT counterparts (Fig. 2B). We repeated the hot plate assay for 8 consecutive days to assess the impact of RGSz1 on the development of morphine analgesic tolerance. By day 3, RGSz1WT mice showed reduced responses to morphine due to the development of tolerance. Analgesic responses were significantly different between genotypes, indicating that the development of tolerance to morphine treatment (15 mg/kg, s.c.) is delayed in male (Fig. 2C) and female (Fig. 2C, Inset) RGSz1KO mice. We similarly assessed morphine tolerance under inflammatory pain conditions. Thermal hind paw hyperalgesia after intraplantar injection of complete Freund’s adjuvant (CFA) was assessed by the Hargreaves test (35). Similar to our observations of animals in pain-free states, morphine analgesia was reduced in RGSz1WT mice after 3 d of treatment whereas it was maintained in RGSz1KO mice throughout the testing period (Fig. 2D).

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Role of RGSz1 in morphine-induced analgesia. (A) Analgesic responses to morphine are enhanced in RGSz1KO mice compared with those in RGSz1WT controls [each dose was tested in a separate group of mice; n = 8–10 per group; two-way ANOVA followed by Bonferroni’s post hoc test, F(1,64) = 9.28; *P < 0.05, **P < 0.01]. A similar phenotype was observed in female mice [Inset: dose, 7.5 mg/kg; n = 10–13 per group; t test, t(21) = 2.885; **P < 0.001]. (B) RGSz1KO mice show increased analgesic responses compared with those in the RGSz1WT mice to fentanyl [0.125 mg/kg; n = 9–11 per genotype; t test; t(18) = 4.677; ***P < 0.001] and methadone [5 mg/kg; n = 11–12 per group; t test, t(21) = 2.519; *P < 0.05]. (C) The development of tolerance to 15 mg/kg morphine was delayed in male [n = 7 per group; two-way ANOVA repeated measures followed by Bonferroni’s post hoc test, F(1,12) = 20.99; ***P < 0.001] and female [Inset: n = 6–7 per group; two-way ANOVA repeated measures followed by Bonferroni’s post hoc test, F(1,11) = 22.89; ***P < 0.001] RGSz1KO mice. (D) Delayed development of morphine tolerance was also observed under inflammatory pain states. Withdrawal latencies of hind paws with CFA-induced inflammation were significantly higher in RGSz1KO mice after 3 and 4 d, demonstrating that morphine (3 mg/kg) analgesia was maintained [n = 4–5 per group; two-way ANOVA repeated measures followed by Bonferroni’s post hoc test, F(1,7) = 18.59; **P < 0.01, ***P < 0.001]. MPE, maximum possible effect.

Interestingly, whereas RGSz1 knockout enhanced the analgesic responses to morphine, it decreased the rewarding effects. Specifically, knockout of RGSz1 decreased the sensitivity to low morphine doses in the place preference paradigm (Fig. 3 A–C). When a higher dose of morphine was applied, both genotypes showed conditioned place preference (Fig. 3D). We also assessed locomotor activation, which is increased in rodents by repeated opioid administration (36). Knockout of RGSz1 lowered the locomotor activation in response to repeated morphine administration (Fig. 3E). Importantly, RGSz1KO mice showed normal responses to other reward-related behaviors, such as sucrose preference [sucrose consumption (percentage of total amount consumed): RGSz1WT, 60.4 ± 6.2; RGSz1KO, 69.7 ± 6.2].

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

RGSz1 modulates reward-like behaviors. (A–D) Knockout of RGSz1 affects the rewarding effect of low morphine doses. RGSz1KO mice did not show conditioned place preference with 1 mg/kg and spent less time in the drug-paired side when 3 mg/kg morphine was used. At higher doses (8 mg/kg), both genotypes showed conditioned place preferences to morphine [two-way repeated-measures ANOVA followed by Bonferroni’s post hoc test; 0.5 mg: n = 8 per group, F(1,14) = 1.12; 1 mg: n = 13 per group, F(1,24) = 0.33, **P < 0.01; 3 mg: n = 11–12 per group, F(1,21) = 2,94, **P < 0.01, ***P < 0.001; 8 mg: n = 16–19 per group, F(1,33) = 0.65, **P < 0.01, ***P < 0.001]. (E) Furthermore, RGSz1KO male mice showed lower locomotor activity in response to repeated morphine administrations than RGSz1WT control mice [10 mg/kg; n = 9–10 per group; two-way ANOVA repeated measures followed by Bonferroni’s post hoc test, F(1,17) = 11.16; *P < 0.05, **P < 0.01, ***P < 0.001]. CPP, conditioned place preference; ns, nonspecific.

RGSz1 in PAG Neurons Modulates Analgesic Tolerance.

We hypothesized that RGSz1 modulates the analgesic and rewarding effects of MOPR agonists via actions in distinct brain regions. To test this, we used viral approaches to conditionally knock down RGSz1 in specific brain regions of adult mice. On the basis of our initial Western blot analysis, we targeted the PAG. AAV2-CMV-CRE-GFP or an AAV2-CMV-GFP control virus was injected into the vlPAG of RGSz1fl/fl mice. Two weeks later, when maximal viral expression was achieved, mice were tested in the morphine tolerance paradigm. The development of morphine tolerance was significantly delayed by knockdown of RGSz1 in the vlPAG (Fig. 4A), corresponding to a 50% reduction of RGSz1 protein levels (Fig. 4B). As shown in Fig. 4 B and D and consistent with previous reports (37), the virus infects only neuronal cells and not astrocytes. Interestingly, the knockdown of RGSz1 in the NAc did not affect the development of analgesic tolerance (Fig. 4C). Notably, the knockdown of RGSz1 in the vlPAG of male mice did not affect analgesic responses to morphine in the hot plate assay {hot plate latencies [percentage of maximal possible effect (%MPE)] for 12 mg/kg morphine: AAV2-CMV-GFP, 50.3 ± 6.7; AAV2-CMV-CRE-GFP, 44.2 ± 8.4}. Furthermore, the actions of RGSz1 in the vlPAG did not affect the development of physical dependence. After treating mice for 5 d with escalating doses of morphine, withdrawal was precipitated with naloxone hydrochloride (1 mg/kg, s.c.). We immediately began monitoring several somatic withdrawal symptoms in AAV2-CMV-CRE–injected RGSz1fl/fl mice and their AAV2-CMV-GFP–injected controls, including jumps, wet-dog shakes (WDSs), tremor, ptosis, diarrhea, and weight loss for a period of 30 min (31). As shown in Fig. 4E, withdrawal symptoms were similar between the two groups of mice. Thus, RGSz1 actions in the vlPAG promote the development of analgesic tolerance without impacting the development of physical dependence.

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

Rgsz1 knockdown in vlPAG neurons suppresses morphine tolerance. (A) Conditional knockdown of Rgsz1 in neurons in the vlPAG of RGSz1fl/fl mice by stereotaxic infection of AAV2-CMV-CRE vectors maintained the analgesic effects of morphine in the hot plate assay over several days [15 mg/kg; n = 7 per group; two-way repeated-measures ANOVA followed by Bonferroni’s post hoc test, F(1,12) = 11.34; *P < 0.05, **P < 0.01, ***P < 0.001]. (B) Western blot analysis verifying RGSz1 down-regulation in the PAG of CRE virus-infected mice [n = 4–5 per group; t test, t(7) = 2.663; *P < 0.05] and verification of AAV2 viral vector distribution in the vlPAG. AAV2 vectors infect neurons (as indicated by NeuN staining) and not astrocytes (as indicated by GFAP staining). (C) The effect of Rgsz1 knockdown on the development of morphine tolerance was region-specific, and no effect was observed when neurons of the NAc were targeted [15 mg/kg; n = 6–7 per group; two-way repeated-measures ANOVA followed by Bonferroni’s post hoc test, F(1,11) = 0.08]. (D) Western blot analysis verifying RGSz1 down-regulation in the NAc of CRE virus-infected mice [n = 9–10 per group; t test, t(17) = 2.123; *P < 0.05] and verification of AAV2 viral vector distribution in the NAc. (E) Knockdown of RGSz1 in the vlPAG does not affect the intensity of naloxone-precipitated morphine withdrawal symptoms (n = 9–10 per group; t tests).

Morphine Distinctly Alters Gene Expression Patterns in the PAG of RGSz1KO and RGSz1WT Mice.

Repeated exposure to opioids promotes long-term adaptations in gene expression and neuronal plasticity (38, 39). We applied next-generation RNA-Seq to understand the impact of morphine tolerance and the influence of RGSz1 on gene expression in the PAG. PAG tissue samples were collected from mice on day 4 of the morphine hot plate assay, when analgesic tolerance emerges in RGSz1WT mice (Fig. 5A). After pooling tissue from two animals per sample, RNA from naive and morphine-treated RGSz1WT and RGSz1KO mice was analyzed. A heat map analysis comparing treated versus naive mice for each genotype showed that morphine tolerance altered gene expression patterns in RGSz1WT mice, but that a much larger number of genes were affected in the RGSz1KO group (Fig. 5B). The number of genes regulated by morphine in each genotype is shown in a Venn diagram in Fig. 5C. Our findings revealed that only one-quarter of the genes regulated by morphine were common between genotypes whereas three-quarters of the genes with altered expression were unique to RGSz1WT or RGSz1KO mice (Fig. 5C). A gene ontology analysis indicated that the majority of genes regulated by chronic morphine in the PAG of RGSz1WT and RGSz1KO mice are associated with transcription and cell proliferation/differentiation whereas genes associated with ion/transmembrane transport, cell adhesion, and behavior are also altered in RGSz1KO mice. We further investigated this observation using ingenuity pathway analysis (IPA), and, indeed, distinct pathways in the PAG were affected in each genotype, with Wnt/β-catenin signaling as the top pathway affected in RGSz1WT mice (P < 0.001) and serotonin receptor signaling targeted in RGSz1KO mice (P < 0.001) (Fig. 6A). Importantly, the pathway analysis of the preexisting gene expression adaptations in the PAG of RGSz1KO naive mice did not reveal alterations in any of the above-stated pathways (Fig. S3). Fig. 6B summarizes the changes in several components for each of these pathways after morphine treatment. We further validated the results by qPCR in separate cohorts of animals. Twenty selected genes per genotype were tested, and around 70% of the tested genes were validated. Differential gene expression patterns for a subset of the validated genes are presented in Fig. 6C.

Fig. 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 5.

RNA-Seq of PAG tissue from morphine-treated RGSz1WT and RGSz1KO mice. (A) Experimental time line showing that mice were treated for 4 consecutive days with 15 mg/kg morphine and tested in the hot plate assay to evaluate analgesic responses. On day 4 (when analgesic tolerance emerges in RGSz1WT mice), PAG tissue was collected 2 h after morphine administration. (B) Heat map analysis showing the effect of drug treatment on gene expression in RGSz1WT and RGSz1KO mice. The overall patterns of gene expression regulation were similar between the two genotypes, but a much larger group of genes were affected in the PAG of RGSz1KO mice. (C) Venn diagram depicting the actual numbers of genes affected in each group and their direction of regulation and associated Gene Ontology (GO) categories. The pie chart depicts results of a GO analysis of the 225 genes regulated by morphine in the PAG of both genotypes. Bars graphs show GO of genes uniquely regulated by morphine in the PAG of RGSz1WT (red) and RGSz1KO (light gray) mice.

Fig. 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 6.

Chronic morphine affects different intracellular pathways in the PAG of RGSz1KO and RGSz1WT mice. (A) IPA analysis of the genes regulated by chronic morphine revealed that distinct canonical pathways were affected in RGSz1WT and RGSz1KO mice (cutoff: P value < 0.003). (B) Selected genes in the canonical Wnt/β-catenin signaling or inflammatory response-related cascades and G protein-coupled/serotonin receptor signaling pathways and the direction of their regulation. (C) qPCR validation of a subset of genes in a separate cohort of animals by genotype [n = 7–8 per group; t tests, RGSz1WT: Ddc, t(13) = 2.135; Gch1, t(13) = 2.2; Pdyn, t(13) = 2.27; Bdkrb1, t(14) = 2.154; Wnt9a, t(14) = 3.4; Cxcl12, t(14) = 3.69; Ccr1, t(12) = 2.99; Il1r1, t(14) = 2.44; Cxcr4, t(14) = 6.08; RGSz1KO: Ddc, t(13) = 2.29; Tph2, t(13) = 2.4; Gch1, t(13) = 2.36; Pdyn, t(13) = 2.9; Wnt9a, t(14) = 0.87; Il1r1, t(12) = 5.49; NPY, t(14) = 2.86; Gabre, t(14) = 2.98; *P < 0.05, **P < 0.01, ***P < 0.001]. ns, not significant. MIF, macrophage migration inhibitory factor; PCP, planar cell polarity; TWEAK, TNF-related adaptor molecule.

Wnt/β-Catenin Signaling Is Essential for Analgesic Responses to Morphine.

Several key molecules in the Wnt/β-catenin signaling pathway were down-regulated by chronic morphine treatment in the PAG of RGSz1WT but not RGSz1KO mice (shown in blue in Fig. 7A). Furthermore, narrowing our lists of differentially regulated genes to known β-catenin targets (40) revealed that those genes in RGSz1WT mice were primarily down-regulated (Fig. 7B). We hypothesized that analgesic responses to morphine are maintained in RGSz1KO mice because β-catenin remains active, whereas signaling in this pathway is suppressed in RGSz1WT mice. Western blotting confirmed this hypothesis, as higher levels of active phospho-β-catenin (ser675) and inactive phospho-GSK3β (ser9) were observed in the PAG of RGSz1KO mice after 4 d of morphine treatment (Fig. 7 C and D). Furthermore, antagonism of β-catenin activity in the vlPAG of RGSz1KO mice by expressing a dominant-negative form of the protein (40) induced tolerance similar to that seen in RGSz1WT controls whereas overexpressing β-catenin itself in the vlPAG of RGSz1WT mice during morphine treatment was not sufficient to delay the development of tolerance (Fig. 7E; for viral validation, see Fig. S4). We therefore hypothesized that, in addition to β-catenin, other components of this pathway are necessary to obtain the prolonged analgesic response phenotype observed in RGSz1KO mice. It is also possible that the lack of an observed phenotype when β-catenin was overexpressed in the vlPAG may be related to compensatory changes in the expression of destruction complex components resulting from the high levels of transgene expression. In the absence of stimulation, β-catenin is bound to Axin2 in the cytoplasm, which targets it for proteasomal degradation (41). Wnt signaling releases β-catenin for translocation to the nucleus. Because Axin2 contains an RGS domain that has been shown to interact with Gα subunits (42⇓–44), we hypothesized that an interaction between Axin2 and Gαz indirectly affects Wnt/β-catenin signaling. This hypothesis was confirmed by our next set of biochemical assays. Although the total levels of Axin2 in the PAG were not affected by repeated morphine administrations (Fig. 8A), the interaction with Gαz in the synaptosomal fraction was significantly decreased in RGSz1WT mice, as indicated by coprecipitation of Gαz with Axin2 (Fig. 8B). However, this interaction was maintained in RGSz1KO mice that did not exhibit morphine tolerance (Fig. 8B). The dissociation of Axin2 from Gαz was accompanied by higher levels of Axin2 in the nuclear fractions of the samples, as shown in Fig. 8C. Furthermore, a correlation analysis revealed a significant negative association between the amount of Axin2 bound to Gαz in the synaptosomal fraction and the amount of Axin2 present in the nuclear fraction (Fig. 8D). The translocation of Axin2 to the nucleus was previously shown to block β-catenin–mediated transcription (45).

Fig. 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 7.

RGSz1 modulates the activity of the Wnt/β-catenin signaling pathway. (A) Schematic depiction of Wnt/β-catenin signaling pathway components (in blue) that are affected by chronic morphine in the PAG of RGSz1WT mice only. (B) Bar graph showing the subset of genes that are known β-catenin targets and the direction of their regulation after chronic morphine administration. Notably, the majority of genes in the RGSz1WT group were down-regulated. (C) Western blot analysis revealed regulation of active phospho-β-catenin (ser675) only in the PAG of RGSz1KO mice after 4 consecutive days of morphine administration [15 mg/kg; n = 14–15 per group; two-way ANOVA followed by Bonferroni’s post hoc test, F(1,53) = 7.37; *P < 0.05], with total β-catenin levels remaining unaffected [n = 11–12 per group; two-way ANOVA followed by Bonferroni’s post hoc test, F(1,40) = 0.18]. (D) Levels of inactive phospho-GSK3β (ser9) were also increased in the PAG of morphine-treated RGSz1KO mice [15 mg/kg; two-way ANOVA followed by Bonferroni’s post hoc test; for phospho-GSK3β: n = 14–15 per group; F(1,53) = 11.65, **P < 0.01; for total GSK3β: n = 11–12 per group; F(1,40) = 1.004]. (E) Antagonizing the activity of β-catenin in the vlPAG by overexpressing a dominant-negative form in RGSz1KO mice led to the development of analgesic tolerance to morphine [15 mg/kg; n = 7 for RGSz1KO-HSV-GFP and 12 for RGSz1KO-HSV-β-catDN; two-way repeated-measures ANOVA followed by Bonferroni’s post hoc test, F(1,51) = 9.9; *P < 0.05, ***P < 0.001]. However, overexpressing β-catenin in the vlPAG of RGSz1WT mice did not delay the development of analgesic tolerance [15 mg/kg; n = 9 for RGSz1WT-HSV-GFP and 10 for RGSz1WT-HSV-β-cat; two-way repeated-measures ANOVA followed by Bonferroni’s post hoc test, F(1,51) = 0.14]. M, morphine; S, saline.

Fig. 8.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 8.

RGSz1 knockout prevents dissociation of Axin2-Gαz complexes after chronic morphine. (A) Protein levels of Axin2 in the PAG were not affected by morphine treatment in RGSz1WT or RGSz1KO mice [15 mg/kg; n = 12 per group; two-way ANOVA followed by Bonferroni’s post hoc test, F(1,44) = 0.19]. (B) Co-IP assays on synaptosomal fractions revealed a reduced association of Axin2 with Gαz in RGSz1WT but not in RGSz1KO mice after repeated morphine administration [15 mg/kg; n = 3 per group; two-way ANOVA followed by Bonferroni’s post hoc test; F(1,8) = 7.6; **P < 0.01]. (C) Morphine tolerance in RGSz1WT also led to higher levels of Axin2 in the nucleus [15 mg/kg; n = 3 per group; two-way ANOVA followed by Bonferroni’s post hoc test; F(1,8) = 7; *P < 0.05]. (D) Correlation analysis revealed a negative association between the amount of Axin2 bound to Gαz in the synaptosomal fractions and the amount of Axin2 in the nuclear fractions of the same samples from RGSz1WT but not RGSz1KO mice. (E) Schematic depicting working hypothesis. In the PAG of tolerant mice, higher levels of RGSz1 and stabilized complexes between MOPR and Gαz lead to a reduced Axin2-Gαz interaction, thereby enabling the function of Axin2-containing destruction complexes, as well as the translocation of Axin2 to the nucleus and silencing β-catenin–mediated transcription. In the absence of RGSz1, Axin2 is maintained in the cytoplasm and forms complexes with Gαz whereas β-catenin may translocate to the nucleus to promote transcriptional activity. M, morphine; S, saline; TL, total lysate.

Discussion

Our findings provide evidence for a brain region-specific mechanism that modulates the functional responses of MOPRs. We show that constitutive deletion of Rgsz1 increases analgesic responses to morphine, methadone, and fentanyl and delays the development of analgesic tolerance in pain-free and chronic pain states. Importantly, the deletion of Rgsz1 decreases the rewarding effect of morphine. These features indicate that RGSz1 is a promising target for the development of adjunct medications to enhance the efficacy of typical opioid analgesics and lower their abuse potential. Our data also indicate that the modulatory roles of RGSz1 in morphine analgesia and tolerance are the same in both male and female mice, further supporting the value of this molecule as a potential new target for pain management.

By controlling G protein-coupled receptor (GPCR) activity, RGS proteins may dynamically affect the functions of ion channels, signal transduction cascades, epigenetic modifiers, and transcription factors (15⇓–17, 46). Several members of the RGS family have been shown to modulate the actions of MOPR agonists in vitro and in vivo (29⇓⇓⇓–33). In our previous studies, we found that RGS9-2 negatively modulates both morphine analgesia and reward-related behaviors (29, 31). Therefore, the inhibition of RGS9-2 not only improves the analgesic effects but would also likely increase the rewarding and locomotor effects of the drug and promote physical dependence. In addition, RGS9-2 has a ligand-dependent role on MOPR function, as it positively modulates the effects of fentanyl and oxycodone (30, 47). Another RGS protein, RGS7, similarly modulates both reward- and analgesia-related behaviors (32). RGS4 in the NAc has a less potent but significant effect on morphine reward, and, in the locus coeruleus, it modulates physical dependence; however, this molecule does not affect analgesia or tolerance to morphine (33). Although RGS4 and RGS7 are also expressed at moderate-to-high levels in the PAG (48), RGSz1 is unique in its differential modulation of opioid reward and analgesia.

Our study focused on understanding the mechanism by which RGSz1 modulates opioid analgesia and tolerance. Using conditional knockdown approaches, we demonstrate that the development of analgesic tolerance is fostered by RGSz1 in the vlPAG, a key area involved in the descending control of pain and morphine-induced analgesia (34). Importantly, the actions of RGSz1 in the vlPAG do not affect the development of physical dependence. A variety of studies document adaptations in this brain area in models of analgesic tolerance (49⇓⇓–52). However, the detailed mechanisms that contribute to MOPR signal transduction and desensitization in the PAG are not fully understood. Our results indicate that RGSz1 forms complexes with Gαz and MOPRs after chronic opioid use. These findings are in accordance with previous work by Hendry et al. (27) suggesting that compromised Gαz actions are associated with morphine tolerance. As mentioned above, several components of the MOPR desensitization machinery may affect the development of analgesic tolerance. This study highlighted the role of RGSz1, Gαz, and components of the β-catenin pathway in morphine tolerance.

Our group has used next-generation sequencing in knockout mouse models and with pharmacological treatment to monitor adaptations in gene expression under pain states and to identify the pathways involved (53, 54). We applied this methodology to determine the mechanism by which chronic morphine induces analgesic tolerance and the impact of RGSz1 in such adaptations. In agreement with previously reported data, our findings indicate that analgesic tolerance to morphine is linked to robust adaptations in genes associated with inflammatory cascades (55⇓⇓–58), involving 8 of the 15 most significantly affected pathways identified by IPA. Notably, most of these adaptations were observed only in RGSz1WT mice and therefore are likely to be linked to morphine tolerance. Furthermore, there was a robust up-regulation of genes associated with serotonin production and release in the PAG of RGSz1KO mice, highlighting the importance of the supraspinal serotonergic system in maintaining analgesic responses (59). We observed that several components of the Wnt/β-catenin signaling pathway were down-regulated in the RGSz1WT group at the time the mice become tolerant to morphine. Although there is no known direct link between Wnt/β-catenin signaling and opioid actions, earlier studies have shown that the inhibition of GSK3β, which promotes the degradation of β-catenin (41), delays the development of morphine tolerance (60, 61). Furthermore, there is evidence for Wnt/β-catenin signaling in the morphological effects of long-term exposure to MOPR agonists, such as reduced dendritic arborization, neurite outgrowth, and neurogenesis in several brain regions (62⇓⇓–65). Here, we confirm the involvement of this pathway by showing that morphine tolerance can develop in RGSz1KO mice when β-catenin actions in the vlPAG are directly antagonized. However, β-catenin overexpression is not sufficient to delay the development of tolerance in RGSz1WT mice, suggesting that additional factors control transcriptional activation as tolerance develops or that the robust overexpression of β-catenin affected the function of components of the destruction pathway. Our data suggest that states of morphine tolerance are associated with increased abundance of Axin2 in the nucleus, where it may repress β-catenin–mediated transcription (45).

Wnt/β-catenin signaling is among the most conserved intracellular pathways, playing pivotal roles in cell proliferation, migration, and homeostasis (41). Therefore, the pharmacological inhibition of Wnt components would broadly affect many tissues and cellular processes in the CNS and the periphery. The identification of molecules such as RGSz1 that modulate Wnt/β-catenin signal transduction in a brain region-specific manner may provide a novel avenue for developing adjunct medications to MOPR agonists for the management of chronic pain.

Additional research will be needed to continue to detail the actions of RGSz1 in the brain reward network and determine the key protein interactions and intercellular pathways by which RGSz1 promotes the rewarding effects of opioids. Finally, it will be valuable to investigate if RGSz1 modulates the analgesic effects of MOPR agonists in other brain regions and the intracellular pathways involved. Overall, our findings provide information on the signal transduction networks mediating the effects of MOPR agonists in male and female mice and identify RGSz1 as a target for optimizing the analgesic actions of MOPR agonists.

Materials and Methods

Animals.

For these studies, we used 2- to 3-mo-old male and female RGSz1 knockout (RGSz1KO) and WT (RGSz1WT) or RGSz1-floxed (RGSz1fl/fl) mice derived from homozygote breedings. The RGSz1KO and RGSz1fl/fl lines are on a DBA background, backcrossed to C57BL/6 mice for three and two generations, respectively. For the RGSz1 locus targeting strategy, see Fig. S2. The primers used for genotyping are listed in Table S1. Mice were housed with a 12-h dark/light cycle according to the Animal Care and Use Committee of the Icahn School of Medicine at Mount Sinai.

CFA Model.

Chronic inflammation was induced by an intraplantar injection of 25 μL of complete Freund’s adjuvant (diluted 1:1 in saline; Sigma-Aldrich) in the left hind paw (66). Morphine treatment (3 mg/kg) was initiated 24 h later and repeated for 4 consecutive days (day 1 to 4).

Hot Plate Analgesia Assays.

Analgesia was measured using a 52 °C hot plate apparatus (IITC Life Sciences) as previously described (29). Animals were habituated in the room for 1 h, and the baseline latencies to jump or lick the hind paw were measured. Morphine was administered s.c., and, 30 min later, mice were placed on the hot plate and the latencies to jump or lick the hind paw were monitored. A cutoff time of 40 s was used to avoid tissue damage and inflammation (31). For tolerance studies, this procedure was repeated for 4 to 8 consecutive days using a high dose of morphine (15 mg/kg, s.c.). Data are expressed as percentages of maximal possible effect [MPE = (latency − baseline)/(cutoff − baseline)].

Hargreaves Assay.

Morphine tolerance under inflammatory pain conditions was assessed by measuring thermal hyperalgesia in the Hargreaves test. For this assay, mice were placed in a Plexiglas box with a glass bottom, and the latency to withdraw the injured (CFA-injected) hind paw was measured after a high-intensity heat beam (30%) was applied to the midplantar area (IITC Life Sciences). Two measurements were obtained with a 10-min interval, and the average was defined as the thermal nociceptive threshold. A cutoff time of 20 s was used to avoid tissue damage (67).

Stereotaxic Virus Injection Surgery.

Conditional deletion of Rgsz1 was achieved by bilateral stereotaxic injections of the AAV2-CMV-CRE-EGFP (University of South Carolina Viral Core Facility) into the NAc or vlPAG of RGSz1fl/fl mice. Control animals received injections of AAV2-CMV-EGFP vectors. Stereotaxic coordinates for viral vector injections were as follows: NAc (with respect to bregma): anteroposterior (AP), +1.6 mm; anterolateral (AL), +1.5 mm; and dorsoventral (DV), −4.4 mm at 10° from the midline; vlPAG (with respect to lambda): AP, +0.6 mm; mediolateral (ML), +0.8 mm; and DV, −2.8 mm at 22° from the midline. HSV-β-catenin (HSV-β-cat), HSV-β-catenin dominant negative (HSV-β-catDN), and HSV-GFP vectors were generated by R.L.N. (40).

Conditioned Place Preference Test.

An unbiased place-conditioning procedure was performed using a two-chamber place-conditioning system (Med Associates, Inc.) (29, 31). Briefly, after monitoring baseline preferences for 20 min on day 1, animals were conditioned to the drug-paired or saline-paired side for 45 min on alternate days. After six conditioning sessions, animals were tested for 20 min, and the results are presented as the time spent in the drug-paired compartment at baseline compared with that on test day.

Locomotor Activation Paradigm.

Mice received a saline injection (s.c.) and were habituated to the locomotor apparatus (Med Associates, Inc.) for 30 min each day for 3 consecutive days. Following that, mice received injections of morphine (10 mg/kg, s.c.), and ambulatory activity was monitored for 30 min on 5 consecutive days (47).

Opiate Withdrawal Paradigm.

Mice were injected intraperitoneally (i.p.) with escalating morphine doses every 12 h for 4 consecutive days (day 1, 20 mg/kg; day 2, 40 mg/kg; day 3, 60 mg/kg; day 4, 80 mg/kg). On the 5th day, they were injected with 80 mg/kg morphine in the morning, and withdrawal was precipitated 3 h later using naloxone hydrochloride (1 mg/kg, s.c.; Sigma) (31). Withdrawal signs (jumps, wet-dog shakes, tremor, diarrhea, and weight loss) were monitored for 30 min after naloxone administration. Data are expressed as the percentages of the AAV2-GFP–injected control. For tremor and ptosis, we monitored the presence of signs at the beginning of each 5-min interval during the monitoring period (29).

Western Blotting, Subcellular Fractionations, and Co-IP Assays.

Two 14-gauge PAG punches (derived from a single animal) were used per sample for all assays. Western blotting was performed as previously described (53). For subcellular fractionation, samples were homogenized using a pestle in 50 μL of hypotonic cell lysis buffer (HCLB) plus protease and proteasome inhibitors, incubated on ice for 45 min, and centrifuged at 371 × g at 4 °C for 10 min, and supernatants were collected as the crude cytoplasmic fractions. The remaining nuclear pellets were washed twice with 100 μL of HCLB plus inhibitors, resuspended in 25 μL of nuclear lysis buffer with inhibitors, incubated for 2 h on ice, and centrifuged at 16,363 × g at 4 °C for 5 min, and the supernatants were collected as the crude nuclear fractions. To obtain the synaptosomal fractions, the crude cytoplasmic fractions were centrifuged twice to remove the remaining nuclei (371 × g at 4 °C for 10 min). The supernatants were collected and centrifuged at 9,361 × g at 4 °C for 10 min. Supernatants were collected as the cytoplasmic fractions, and the synaptosomal pellets were resuspended in 20 μL of HCLB.

Co-IPs for MOPR-Gαz were performed using whole lysate samples as previously reported (30). For co-IPs with Axin2, 1 μL of antibody was added to each synaptosomal fraction (co-IP with anti-GFP was used as the negative control), and samples were incubated overnight at 4 °C on a circular rotator. Each sample was then incubated for an additional 3 h with 10 μL of anti-rabbit Dynabeads (Thermo Fisher Scientific). Co-IP proteins were washed three times with 700 μL of blocking solution using a magnetic stand and then resuspended in 10 μL of radioimmunoprecipitation assay (RIPA) buffer. The following rabbit antibodies were used: anti-RGSz1 (19), anti-Gαz (no. 3904; Cell Signaling), anti-Gαi3 (Lilly Jiang, University of Texas Southwestern Medical Center, Dallas, TX) (68), anti-Gαo (no. 3975; Cell Signaling), anti-Gαq (Paul Sternweis, University of Texas Southwestern Medical Center, Dallas, TX) (69), anti-MOPR (no. 24216; Immunostar), anti–phospho-β-catenin (ser675, no. 4176; Cell Signaling), anti–β-catenin (no. 9562; Cell Signaling), anti–phospho-GSK3β (ser9 no. 9336; Cell Signaling), anti-GSK3β (no. 9315; Cell Signaling), anti-Axin2 (no. ab32197; Abcam), anti-GADPH (no. 5174; Cell Signaling), anti–β-actin (no. 4967; Cell Signaling), and anti-GFP (no. ab290; Abcam). Mouse anti–β-tubulin III (no. T8578; Sigma), GADPH, or β-actin antibodies were used as the loading control. Data are calculated as the optical density ratio of sample/loading control or IP/total lysate in the case of co-IPs and are expressed as the percentages of the control group. Bands were quantified using Image J software (53).

RNA-Seq and qPCR.

Four biological replicates per group were used for the RNA-Seq studies. PAG punches from two animals were pooled per biological replicate (four 14-gauge punches per sample). Total RNA was isolated with TRIzol, and the integrity was confirmed with an Agilent 2100 Bioanalyzer (53). mRNA-Seq libraries were prepared using a TruSeq RNA sample preparation kit v2 (Illumina). Sequencing was performed using the Illumina HiSeq 2000 apparatus. Read alignment, read counting, and differential analysis were performed using TopHat2 (70), HTSeq (71), and voom-limma R package (72), respectively. For all four gene lists, a cutoff of P value of <0.05 and log2(fold change) of <−0.5 or >0.5 were used to generate gene lists for further bioinformatics analysis. qPCR was performed using SYBR green and analyzed using the ΔΔCT method. Information on primers used for biological validations is provided in Table S1.

Bioinformatics Analysis.

Heat maps were generated using GENE-E, Venn diagrams were generated using VennPlex Version 1.0.0.2 (NIH), pathway analysis was conducted using IPA, and Gene Ontology (GO) analysis was conducted using the Database for Annotation, Visualization and Integrated Discovery (DAVID). For pathway analysis, a cutoff of P value of <0.003 was applied to the output pathways.

Statistical Analysis.

For the experiments monitoring behavior of the same group of mice over time (Figs. 2 C and D, 3, 4 A and C, and 7E and Fig. S2B), we used two-way repeated-measures ANOVAs followed by Bonferroni’s post hoc tests. For two-factor designs (Figs. 1 B–G, 2A, and 7 C and D and Fig. S2C), we used two-way ANOVAs followed by Bonferroni’s post hoc tests. For data containing a single independent variable (Figs. 1A, Upper, 2A, and 4B and Fig. S2 D and E), we used unpaired two-tailed t tests. For the data in graphs depicting multiple individual single-factor comparisons (Figs. 2B, 4E, and 6C), we used multiple t tests. Correlation analysis (Fig. 8D) was contacted using Spearman’s rho. Error bars are depicting ±SEM. F and t values for each dataset are provided in the figure legends.

Acknowledgments

This work was supported by National Institute on Drug Abuse Grant PPG-POIDAO8227 (to V.Z.), National Institute of Neurological Disorders and Stroke (NINDS) Grant NS086444 (to V.Z.), NINDS Grant NS093537 (to V.Z.), National Institute of General Medical Sciences Grant R01030355 (to E.M.R.), and National Institute of Deafness and Other Communication Disorders (NIDCD) Grant R21DC015744 (to E.M.R.).

Footnotes

  • ↵1To whom correspondence should be addressed. Email: venetia.zachariou{at}mssm.edu.
  • Author contributions: S.G., I.P., F.S., D.H., R.H.R., E.M.R., and V.Z. designed research; S.G., I.P., V.C., F.S., L.S., and V.Z. performed research; R.L.N., D.H., R.H.R., and E.M.R. contributed new reagents/analytic tools; S.G., I.P., V.C., F.S., D.H., R.H.R., L.S., and V.Z. analyzed data; and S.G., E.M.R., L.S., and V.Z. wrote the paper.

  • The authors declare no conflict of interest.

  • This article is a PNAS Direct Submission.

  • Data deposition: The data reported in this paper have been deposited in the Gene Expression Omnibus (GEO) database, https://www.ncbi.nlm.nih.gov/geo (accession no. GSE109165).

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1707887115/-/DCSupplemental.

Published under the PNAS license.

References

  1. ↵
    1. Dart RC, et al.
    (2015) Trends in opioid analgesic abuse and mortality in the United States. N Engl J Med 372:241–248.
    OpenUrlCrossRefPubMed
  2. ↵
    1. Ray WA,
    2. Chung CP,
    3. Murray KT,
    4. Hall K,
    5. Stein CM
    (2016) Prescription of long-acting opioids and mortality in patients with chronic noncancer pain. JAMA 315:2415–2423.
    OpenUrlCrossRefPubMed
  3. ↵
    1. Ling W,
    2. Mooney L,
    3. Hillhouse M
    (2011) Prescription opioid abuse, pain and addiction: Clinical issues and implications. Drug Alcohol Rev 30:300–305.
    OpenUrlCrossRefPubMed
  4. ↵
    1. Walwyn WM,
    2. Miotto KA,
    3. Evans CJ
    (2010) Opioid pharmaceuticals and addiction: The issues, and research directions seeking solutions. Drug Alcohol Depend 108:156–165.
    OpenUrlCrossRefPubMed
  5. ↵
    1. Walter C,
    2. Knothe C,
    3. Lötsch J
    (2016) Abuse-deterrent opioid formulations: Pharmacokinetic and pharmacodynamic considerations. Clin Pharmacokinet 55:751–767.
    OpenUrl
  6. ↵
    1. Skolnick P,
    2. Volkow ND
    (2016) Re-energizing the development of pain therapeutics in light of the opioid epidemic. Neuron 92:294–297.
    OpenUrl
  7. ↵
    1. Manglik A, et al.
    (2016) Structure-based discovery of opioid analgesics with reduced side effects. Nature 537:185–190.
    OpenUrlCrossRefPubMed
  8. ↵
    1. Cahill CM,
    2. Walwyn W,
    3. Taylor AMW,
    4. Pradhan AAA,
    5. Evans CJ
    (2016) Allostatic mechanisms of opioid tolerance beyond desensitization and downregulation. Trends Pharmacol Sci 37:963–976.
    OpenUrlCrossRef
  9. ↵
    1. Kreek MJ, et al.
    (2012) Opiate addiction and cocaine addiction: Underlying molecular neurobiology and genetics. J Clin Invest 122:3387–3393.
    OpenUrlCrossRefPubMed
  10. ↵
    1. Contet C,
    2. Kieffer BL,
    3. Befort K
    (2004) Mu opioid receptor: A gateway to drug addiction. Curr Opin Neurobiol 14:370–378.
    OpenUrlCrossRefPubMed
  11. ↵
    1. Bailey CP,
    2. Connor M
    (2005) Opioids: Cellular mechanisms of tolerance and physical dependence. Curr Opin Pharmacol 5:60–68.
    OpenUrlCrossRefPubMed
  12. ↵
    1. Trang T, et al.
    (2015) Pain and poppies: The good, the bad, and the ugly of opioid analgesics. J Neurosci 35:13879–13888.
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. Raehal KM,
    2. Bohn LM
    (2014) β-arrestins: Regulatory role and therapeutic potential in opioid and cannabinoid receptor-mediated analgesia. Handb Exp Pharmacol 219:427–443.
    OpenUrlCrossRefPubMed
  14. ↵
    1. Barker SA,
    2. Wang J,
    3. Sierra DA,
    4. Ross EM
    (2001) RGSZ1 and Ret RGS: Two of several splice variants from the gene RGS20. Genomics 78:223–229.
    OpenUrlCrossRefPubMed
  15. ↵
    1. Gerber KJ,
    2. Squires KE,
    3. Hepler JR
    (2016) Roles for regulator of G protein signaling proteins in synaptic signaling and plasticity. Mol Pharmacol 89:273–286.
    OpenUrlAbstract/FREE Full Text
  16. ↵
    1. Sjögren B,
    2. Blazer LL,
    3. Neubig RR
    (2010) Regulators of G protein signaling proteins as targets for drug discovery. Prog Mol Biol Transl Sci 91:81–119.
    OpenUrlCrossRefPubMed
  17. ↵
    1. Kimple AJ,
    2. Bosch DE,
    3. Giguère PM,
    4. Siderovski DP
    (2011) Regulators of G-protein signaling and their Gα substrates: Promises and challenges in their use as drug discovery targets. Pharmacol Rev 63:728–749.
    OpenUrlAbstract/FREE Full Text
  18. ↵
    1. Glick JL,
    2. Meigs TE,
    3. Miron A,
    4. Casey PJ
    (1998) RGSZ1, a Gz-selective regulator of G protein signaling whose action is sensitive to the phosphorylation state of Gzalpha. J Biol Chem 273:26008–26013.
    OpenUrlAbstract/FREE Full Text
  19. ↵
    1. Wang J, et al.
    (1998) RGSZ1, a Gz-selective RGS protein in brain. Structure, membrane association, regulation by Galphaz phosphorylation, and relationship to a Gz gtpase-activating protein subfamily. J Biol Chem 273:26014–26025.
    OpenUrlAbstract/FREE Full Text
  20. ↵
    1. Larminie C, et al.
    (2004) Selective expression of regulators of G-protein signaling (RGS) in the human central nervous system. Brain Res Mol Brain Res 122:24–34.
    OpenUrlCrossRefPubMed
  21. ↵
    1. Hepler JR,
    2. Berman DM,
    3. Gilman AG,
    4. Kozasa T
    (1997) RGS4 and GAIP are GTPase-activating proteins for Gq alpha and block activation of phospholipase C beta by gamma-thio-GTP-Gq alpha. Proc Natl Acad Sci USA 94:428–432.
    OpenUrlAbstract/FREE Full Text
  22. ↵
    1. Anger T,
    2. Zhang W,
    3. Mende U
    (2004) Differential contribution of GTPase activation and effector antagonism to the inhibitory effect of RGS proteins on Gq-mediated signaling in vivo. J Biol Chem 279:3906–3915.
    OpenUrlAbstract/FREE Full Text
  23. ↵
    1. Ajit SK, et al.
    (2007) RGSZ1 interacts with protein kinase C interacting protein PKCI-1 and modulates mu opioid receptor signaling. Cell Signal 19:723–730.
    OpenUrlCrossRefPubMed
  24. ↵
    1. Ghavami A, et al.
    (2004) Differential effects of regulator of G protein signaling (RGS) proteins on serotonin 5-HT1A, 5-HT2A, and dopamine D2 receptor-mediated signaling and adenylyl cyclase activity. Cell Signal 16:711–721.
    OpenUrlCrossRefPubMed
  25. ↵
    1. Creech RD,
    2. Li Q,
    3. Carrasco GA,
    4. Van de Kar LD,
    5. Muma NA
    (2012) Estradiol induces partial desensitization of serotonin 1A receptor signaling in the paraventricular nucleus of the hypothalamus and alters expression and interaction of RGSZ1 and Gαz. Neuropharmacology 62:2040–2049.
    OpenUrl
  26. ↵
    1. Garzón J,
    2. Rodríguez-Muñoz M,
    3. López-Fando A,
    4. García-España A,
    5. Sánchez-Blázquez P
    (2004) RGSZ1 and GAIP regulate mu- but not delta-opioid receptors in mouse CNS: Role in tachyphylaxis and acute tolerance. Neuropsychopharmacology 29:1091–1104.
    OpenUrlCrossRefPubMed
  27. ↵
    1. Hendry IA, et al.
    (2000) Hypertolerance to morphine in G(z alpha)-deficient mice. Brain Res 870:10–19.
    OpenUrlCrossRefPubMed
  28. ↵
    1. Yang J, et al.
    (2000) Loss of signaling through the G protein, Gz, results in abnormal platelet activation and altered responses to psychoactive drugs. Proc Natl Acad Sci USA 97:9984–9989.
    OpenUrlAbstract/FREE Full Text
  29. ↵
    1. Gaspari S, et al.
    (2014) Nucleus accumbens-specific interventions in RGS9-2 activity modulate responses to morphine. Neuropsychopharmacology 39:1968–1977.
    OpenUrlCrossRefPubMed
  30. ↵
    1. Psifogeorgou K, et al.
    (2011) A unique role of RGS9-2 in the striatum as a positive or negative regulator of opiate analgesia. J Neurosci 31:5617–5624, and erratum (2011) 31:7578.
    OpenUrlAbstract/FREE Full Text
  31. ↵
    1. Zachariou V, et al.
    (2003) Essential role for RGS9 in opiate action. Proc Natl Acad Sci USA 100:13656–13661.
    OpenUrlAbstract/FREE Full Text
  32. ↵
    1. Sutton LP, et al.
    (2016) Regulator of G-protein signaling 7 regulates reward behavior by controlling opioid signaling in the striatum. Biol Psychiatry 80:235–245.
    OpenUrlCrossRefPubMed
  33. ↵
    1. Han M-H, et al.
    (2010) Brain region specific actions of regulator of G protein signaling 4 oppose morphine reward and dependence but promote analgesia. Biol Psychiatry 67:761–769.
    OpenUrlCrossRefPubMed
  34. ↵
    1. Reichling DB,
    2. Kwiat GC,
    3. Basbaum AI
    (1988) Anatomy, physiology and pharmacology of the periaqueductal gray contribution to antinociceptive controls. Prog Brain Res 77:31–46.
    OpenUrlCrossRefPubMed
  35. ↵
    1. Hargreaves K,
    2. Dubner R,
    3. Brown F,
    4. Flores C,
    5. Joris J
    (1988) A new and sensitive method for measuring thermal nociception in cutaneous hyperalgesia. Pain 32:77–88.
    OpenUrlCrossRefPubMed
  36. ↵
    1. Robinson TE,
    2. Berridge KC
    (1993) The neural basis of drug craving: An incentive-sensitization theory of addiction. Brain Res Brain Res Rev 18:247–291.
    OpenUrlCrossRefPubMed
  37. ↵
    1. Taymans J-M, et al.
    (2007) Comparative analysis of adeno-associated viral vector serotypes 1, 2, 5, 7, and 8 in mouse brain. Hum Gene Ther 18:195–206.
    OpenUrlCrossRefPubMed
  38. ↵
    1. Robison AJ,
    2. Nestler EJ
    (2011) Transcriptional and epigenetic mechanisms of addiction. Nat Rev Neurosci 12:623–637.
    OpenUrlCrossRefPubMed
  39. ↵
    1. Russo SJ, et al.
    (2010) The addicted synapse: Mechanisms of synaptic and structural plasticity in nucleus accumbens. Trends Neurosci 33:267–276.
    OpenUrlCrossRefPubMed
  40. ↵
    1. Dias C, et al.
    (2014) β-catenin mediates stress resilience through Dicer1/microRNA regulation. Nature 516:51–55.
    OpenUrlCrossRefPubMed
  41. ↵
    1. Clevers H
    (2006) Wnt/beta-catenin signaling in development and disease. Cell 127:469–480.
    OpenUrlCrossRefPubMed
  42. ↵
    1. Egger-Adam D,
    2. Katanaev VL
    (2010) The trimeric G protein Go inflicts a double impact on axin in the Wnt/frizzled signaling pathway. Dev Dyn 239:168–183.
    OpenUrlCrossRefPubMed
  43. ↵
    1. Stemmle LN,
    2. Fields TA,
    3. Casey PJ
    (2006) The regulator of G protein signaling domain of axin selectively interacts with Galpha12 but not Galpha13. Mol Pharmacol 70:1461–1468.
    OpenUrlAbstract/FREE Full Text
  44. ↵
    1. Castellone MD,
    2. Teramoto H,
    3. Williams BO,
    4. Druey KM,
    5. Gutkind JS
    (2005) Prostaglandin E2 promotes colon cancer cell growth through a Gs-axin-beta-catenin signaling axis. Science 310:1504–1510.
    OpenUrlAbstract/FREE Full Text
  45. ↵
    1. Rennoll SA,
    2. Konsavage WM Jr,
    3. Yochum GS
    (2014) Nuclear AXIN2 represses MYC gene expression. Biochem Biophys Res Commun 443:217–222.
    OpenUrlCrossRefPubMed
  46. ↵
    1. Terzi D,
    2. Stergiou E,
    3. King SL,
    4. Zachariou V
    (2009) Regulators of G protein signaling in neuropsychiatric disorders. Prog Mol Biol Transl Sci 86:299–333.
    OpenUrlCrossRefPubMed
  47. ↵
    1. Gaspari S, et al.
    (2017) RGS9-2 modulates responses to oxycodone in pain-free and chronic pain states. Neuropsychopharmacology 42:1548–1556.
    OpenUrlCrossRefPubMed
  48. ↵
    1. Gold SJ,
    2. Ni YG,
    3. Dohlman HG,
    4. Nestler EJ
    (1997) Regulators of G-protein signaling (RGS) proteins: Region-specific expression of nine subtypes in rat brain. J Neurosci 17:8024–8037.
    OpenUrlAbstract/FREE Full Text
  49. ↵
    1. Macey TA, et al.
    (2009) Extracellular signal-regulated kinase 1/2 activation counteracts morphine tolerance in the periaqueductal gray of the rat. J Pharmacol Exp Ther 331:412–418.
    OpenUrlAbstract/FREE Full Text
  50. ↵
    1. Lane DA,
    2. Patel PA,
    3. Morgan MM
    (2005) Evidence for an intrinsic mechanism of antinociceptive tolerance within the ventrolateral periaqueductal gray of rats. Neuroscience 135:227–234.
    OpenUrlCrossRefPubMed
  51. ↵
    1. Bagley EE,
    2. Chieng BCH,
    3. Christie MJ,
    4. Connor M
    (2005) Opioid tolerance in periaqueductal gray neurons isolated from mice chronically treated with morphine. Br J Pharmacol 146:68–76.
    OpenUrlCrossRefPubMed
  52. ↵
    1. Bobeck EN,
    2. Chen Q,
    3. Morgan MM,
    4. Ingram SL
    (2014) Contribution of adenylyl cyclase modulation of pre- and postsynaptic GABA neurotransmission to morphine antinociception and tolerance. Neuropsychopharmacology 39:2142–2152.
    OpenUrlCrossRefPubMed
  53. ↵
    1. Mitsi V, et al.
    (2015) RGS9-2–Controlled adaptations in the striatum determine the onset of action and efficacy of antidepressants in neuropathic pain states. Proc Natl Acad Sci USA 112:E5088–E5097.
    OpenUrlAbstract/FREE Full Text
  54. ↵
    1. Descalzi G, et al.
    (2017) Neuropathic pain promotes adaptive changes in gene expression in brain networks involved in stress and depression. Sci Signal 10:eaaj1549.
    OpenUrlAbstract/FREE Full Text
  55. ↵
    1. Hutchinson MR, et al.
    (2011) Exploring the neuroimmunopharmacology of opioids: An integrative review of mechanisms of central immune signaling and their implications for opioid analgesia. Pharmacol Rev 63:772–810.
    OpenUrlAbstract/FREE Full Text
  56. ↵
    1. Johnston IN, et al.
    (2004) A role for proinflammatory cytokines and fractalkine in analgesia, tolerance, and subsequent pain facilitation induced by chronic intrathecal morphine. J Neurosci 24:7353–7365.
    OpenUrlAbstract/FREE Full Text
  57. ↵
    1. Watkins LR, et al.
    (2007) Norman Cousins Lecture. Glia as the “bad guys”: Implications for improving clinical pain control and the clinical utility of opioids. Brain Behav Immun 21:131–146.
    OpenUrlCrossRefPubMed
  58. ↵
    1. Mélik Parsadaniantz S,
    2. Rivat C,
    3. Rostène W,
    4. Réaux-Le Goazigo A
    (2015) Opioid and chemokine receptor crosstalk: A promising target for pain therapy? Nat Rev Neurosci 16:69–78.
    OpenUrlCrossRefPubMed
  59. ↵
    1. Tao R,
    2. Ma Z,
    3. Auerbach SB
    (1998) Alteration in regulation of serotonin release in rat dorsal raphe nucleus after prolonged exposure to morphine. J Pharmacol Exp Ther 286:481–488.
    OpenUrlAbstract/FREE Full Text
  60. ↵
    1. Parkitna JR, et al.
    (2006) Effects of glycogen synthase kinase 3beta and cyclin-dependent kinase 5 inhibitors on morphine-induced analgesia and tolerance in rats. J Pharmacol Exp Ther 319:832–839.
    OpenUrlAbstract/FREE Full Text
  61. ↵
    1. Liao W-W, et al.
    (2014) Coadministration of glycogen-synthase kinase 3 inhibitor with morphine attenuates chronic morphine-induced analgesic tolerance and withdrawal syndrome. J Chin Med Assoc 77:31–37.
    OpenUrl
  62. ↵
    1. Jin J, et al.
    (2010) Interaction of the mu-opioid receptor with GPR177 (Wntless) inhibits Wnt secretion: Potential implications for opioid dependence. BMC Neurosci 11:33.
    OpenUrlCrossRefPubMed
  63. ↵
    1. Reyes BA, et al.
    (2012) Opiate agonist-induced re-distribution of Wntless, a mu-opioid receptor interacting protein, in rat striatal neurons. Exp Neurol 233:205–213.
    OpenUrlPubMed
  64. ↵
    1. Li Y,
    2. Wang H,
    3. Niu L,
    4. Zhou Y
    (2007) Chronic morphine exposure alters the dendritic morphology of pyramidal neurons in visual cortex of rats. Neurosci Lett 418:227–231.
    OpenUrlCrossRefPubMed
  65. ↵
    1. Eisch AJ,
    2. Barrot M,
    3. Schad CA,
    4. Self DW,
    5. Nestler EJ
    (2000) Opiates inhibit neurogenesis in the adult rat hippocampus. Proc Natl Acad Sci USA 97:7579–7584.
    OpenUrlAbstract/FREE Full Text
  66. ↵
    1. Latremoliere A, et al.
    (2015) Reduction of neuropathic and inflammatory pain through inhibition of the tetrahydrobiopterin pathway. Neuron 86:1393–1406.
    OpenUrlCrossRefPubMed
  67. ↵
    1. Terzi D, et al.
    (2014) RGS9-2 modulates sensory and mood related symptoms of neuropathic pain. Neurobiol Learn Mem 115:43–48.
    OpenUrlCrossRefPubMed
  68. ↵
    1. Mumby SM,
    2. Gilman AG
    (1991) Synthetic peptide antisera with determined specificity for G protein alpha or beta subunits. Methods Enzymol 195:215–233.
    OpenUrlCrossRefPubMed
  69. ↵
    1. Gutowski S, et al.
    (1991) Antibodies to the alpha q subfamily of guanine nucleotide-binding regulatory protein alpha subunits attenuate activation of phosphatidylinositol 4,5-bisphosphate hydrolysis by hormones. J Biol Chem 266:20519–20524.
    OpenUrlAbstract/FREE Full Text
  70. ↵
    1. Kim D, et al.
    (2013) TopHat2: Accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol 14:R36.
    OpenUrlCrossRefPubMed
  71. ↵
    1. Anders S,
    2. Pyl PT,
    3. Huber W
    (2015) HTSeq–A Python framework to work with high-throughput sequencing data. Bioinformatics 31:166–169.
    OpenUrlCrossRefPubMed
  72. ↵
    1. Law CW,
    2. Chen Y,
    3. Shi W,
    4. Smyth GK
    (2014) voom: Precision weights unlock linear model analysis tools for RNA-seq read counts. Genome Biol 15:R29.
    OpenUrlCrossRefPubMed
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Suppression of RGSz1 function optimizes the actions of opioid analgesics by mechanisms that involve the Wnt/β-catenin pathway
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
RGSz1 modulates the actions of opioid analgesics
Sevasti Gaspari, Immanuel Purushothaman, Valeria Cogliani, Farhana Sakloth, Rachael L. Neve, David Howland, Robert H. Ring, Elliott M. Ross, Li Shen, Venetia Zachariou
Proceedings of the National Academy of Sciences Feb 2018, 115 (9) E2085-E2094; DOI: 10.1073/pnas.1707887115

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
RGSz1 modulates the actions of opioid analgesics
Sevasti Gaspari, Immanuel Purushothaman, Valeria Cogliani, Farhana Sakloth, Rachael L. Neve, David Howland, Robert H. Ring, Elliott M. Ross, Li Shen, Venetia Zachariou
Proceedings of the National Academy of Sciences Feb 2018, 115 (9) E2085-E2094; DOI: 10.1073/pnas.1707887115
Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley
Proceedings of the National Academy of Sciences: 115 (9)
Table of Contents

Submit

Sign up for Article Alerts

Article Classifications

  • Biological Sciences
  • Neuroscience

Jump to section

  • Article
    • Abstract
    • Results
    • Discussion
    • Materials and Methods
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Surgeons hands during surgery
Inner Workings: Advances in infectious disease treatment promise to expand the pool of donor organs
Despite myriad challenges, clinicians see room for progress.
Image credit: Shutterstock/David Tadevosian.
Setting sun over a sun-baked dirt landscape
Core Concept: Popular integrated assessment climate policy models have key caveats
Better explicating the strengths and shortcomings of these models will help refine projections and improve transparency in the years ahead.
Image credit: Witsawat.S.
Double helix
Journal Club: Noncoding DNA shown to underlie function, cause limb malformations
Using CRISPR, researchers showed that a region some used to label “junk DNA” has a major role in a rare genetic disorder.
Image credit: Nathan Devery.
Steamboat Geyser eruption.
Eruption of Steamboat Geyser
Mara Reed and Michael Manga explore why Yellowstone's Steamboat Geyser resumed erupting in 2018.
Listen
Past PodcastsSubscribe
Birds nestling on tree branches
Parent–offspring conflict in songbird fledging
Some songbird parents might improve their own fitness by manipulating their offspring into leaving the nest early, at the cost of fledgling survival, a study finds.
Image credit: Gil Eckrich (photographer).

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Special Feature Articles – Most Recent
  • List of Issues

PNAS Portals

  • Anthropology
  • Chemistry
  • Classics
  • Front Matter
  • Physics
  • Sustainability Science
  • Teaching Resources

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Librarians
  • Press
  • Site Map
  • PNAS Updates

Feedback    Privacy/Legal

Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490