Jon Eric Sprague
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Research Interests
Jon Eric Sprague has a total of 119 co-authors. Here are the top co-authors by number of shared publications:
Jeremy R. Canfield
Department of Pharmacology, University of Michigan Medical School, Ann Arbor, MI, 48109, USA
9 shared publicationsRay Larsen
The Department of Biological Sciences, Bowling Green State University, Bowling Green, OH, 43403, USA
6 shared publicationsSudhan Pachhain
Department of Biological Sciences, Bowling Green State University, OH, USA
6 shared publicationsAmal Aburahma
The Ohio Attorney General's Center for the Future of Forensic Science, USA
5 shared publicationsSayantan Roy Choudhury
The Department of Biological Sciences, Bowling Green State University, Bowling Green, OH, USA
5 shared publications- Caroline E. Freiermuth
Center for Addiction Research, University of Cincinnati College of Medicine, Cincinnati, OH 45267, USA
4 shared publications
Latest Publications
- Open Access
Genome-Wide Association Studies and Deep-Learning Functional Annotation of Opioid Use Disorder across Three Ancestries in the All of Us Research Program
Background: Opioid use disorder (OUD) is heritable, yet most genome-wide association studies (GWAS) have focused on European populations, leaving the genetic architecture of OUD in non- European populations underexplored. Methods: We conducted GWAS of OUD across three ancestries using electronic health records and genomic data from 52,357 All of Us Research Program participants (8,912 cases; 43,445 matched opioid-exposed controls; 48.5% female). Participants were stratified into European (EUR), African (AFR), and Admixed American (AMR) ancestry groups for logistic regression GWAS, with independent replication in the Million Veteran Program. We then applied the deep- learning model AlphaGenome to predict the tissue-specific transcriptomic and splicing consequences of top risk variants across 13 reward-pathway brain regions. Results: We identified and replicated a novel DDX6 risk locus, alongside established OPRM1 and FURIN signals. AlphaGenome predicted the DDX6 regulatory allele downregulates the stress-resistance gene FOXR1 in the nucleus accumbens, while the protective OPRM1 variant (rs1799971) upregulates OPRM1 expression across reward networks. Other signals of interest included IL6R and SHISA9 (EUR); GHR (AFR); and ASTN2 (AMR). Conclusions: This study identifies DDX6 as a novel OUD risk locus, replicates associations with OPRM1 and FURIN , and highlights biologically plausible ancestry-specific signals in AFR and AMR populations. We also replicated top variants in an independent population. Finally, integrating GWAS with deep-learning annotations provides specific, localized biological hypotheses to guide future experimental validation and targeted therapeutics.
Open AccessGenome-Wide Association Studies and Deep-Learning Functional Annotation of Opioid Use Disorder across Three Ancestries in the All of Us Research Program
Background: Opioid use disorder (OUD) is heritable, yet most genome-wide association studies (GWAS) have focused on European populations, leaving the genetic architecture of OUD in non- European populations underexplored. Methods: We conducted GWAS of OUD across three ancestries using electronic health records and genomic data from 52,357 All of Us Research Program participants (8,912 cases; 43,445 matched opioid-exposed controls; 48.5% female). Participants were stratified into European (EUR), African (AFR), and Admixed American (AMR) ancestry groups for logistic regression GWAS, with independent replication in the Million Veteran Program. We then applied the deep- learning model AlphaGenome to predict the tissue-specific transcriptomic and splicing consequences of top risk variants across 13 reward-pathway brain regions. Results: We identified and replicated a novel DDX6 risk locus, alongside established OPRM1 and FURIN signals. AlphaGenome predicted the DDX6 regulatory allele downregulates the stress-resistance gene FOXR1 in the nucleus accumbens, while the protective OPRM1 variant (rs1799971) upregulates OPRM1 expression across reward networks. Other signals of interest included IL6R and SHISA9 (EUR); GHR (AFR); and ASTN2 (AMR). Conclusions: This study identifies DDX6 as a novel OUD risk locus, replicates associations with OPRM1 and FURIN , and highlights biologically plausible ancestry-specific signals in AFR and AMR populations. We also replicated top variants in an independent population. Finally, integrating GWAS with deep-learning annotations provides specific, localized biological hypotheses to guide future experimental validation and targeted therapeutics.
- Open Access
- 5 May 2026
- Published by Elsevier in Pharmacology Biochemistry and Behavior
- Vol. 265, 174201
The abuse-related effects of illicitly manufactured fentanyls (IMFs): A mini review
- 5 May 2026
- Published by Elsevier in Pharmacology Biochemistry and Behavior
- Vol. 265, 174201
Open AccessThe abuse-related effects of illicitly manufactured fentanyls (IMFs): A mini review
- Open Access
- 13 August 2025
- Published by Taylor & Francis in Pharmacogenomics
- Vol. 26 (10-12) , 391-425
Scoping review of associations between cytochrome P450 3A4/5 single nucleotide polymorphisms and risk factors for fentanyl overdose
Fentanyl overdose is a public health crisis in the United States, as fentanyl was implicated in nearly 70% of drug overdose deaths in 2023. To provide insight into genetic factors that may influence risk of fentanyl overdose, we conducted a scoping review of associations between cytochrome P450 3A4 (CYP3A4) and 3A5 (CYP3A5) genetic variants and relevant phenotypes. We searched databases through August 2025 for peer-reviewed studies of human subjects or postmortem samples, and integrated evidence from 64 genetic association studies that analyzed single nucleotide polymorphisms in CYP3A4 or CYP3A5. We considered a diverse range of phenotypes relevant to fentanyl overdose, including opioid overdose, fentanyl pharmacokinetics and pharmacodynamics, opioid use (disorder), and pharmacotherapy response. Evidence from 64 studies suggested that the no-function CYP3A5 * 3 (rs776746) allele contributes to increased fentanyl overdose risk, with strongest support coming from studies of fentanyl pharmacokinetics in clinical settings. There was less robust evidence for the role of CYP3A4 variants (e.g. rs2242480). Future research should prioritize prospective genotyping of at-risk populations, development of models that integrate pharmacogenetics with psychiatric genetics, and large-scale harmonization of relevant datasets.
- 13 August 2025
- Published by Taylor & Francis in Pharmacogenomics
- Vol. 26 (10-12) , 391-425
Open AccessScoping review of associations between cytochrome P450 3A4/5 single nucleotide polymorphisms and risk factors for fentanyl overdose
Fentanyl overdose is a public health crisis in the United States, as fentanyl was implicated in nearly 70% of drug overdose deaths in 2023. To provide insight into genetic factors that may influence risk of fentanyl overdose, we conducted a scoping review of associations between cytochrome P450 3A4 (CYP3A4) and 3A5 (CYP3A5) genetic variants and relevant phenotypes. We searched databases through August 2025 for peer-reviewed studies of human subjects or postmortem samples, and integrated evidence from 64 genetic association studies that analyzed single nucleotide polymorphisms in CYP3A4 or CYP3A5. We considered a diverse range of phenotypes relevant to fentanyl overdose, including opioid overdose, fentanyl pharmacokinetics and pharmacodynamics, opioid use (disorder), and pharmacotherapy response. Evidence from 64 studies suggested that the no-function CYP3A5 * 3 (rs776746) allele contributes to increased fentanyl overdose risk, with strongest support coming from studies of fentanyl pharmacokinetics in clinical settings. There was less robust evidence for the role of CYP3A4 variants (e.g. rs2242480). Future research should prioritize prospective genotyping of at-risk populations, development of models that integrate pharmacogenetics with psychiatric genetics, and large-scale harmonization of relevant datasets.
- Open Access
- 20 July 2025
- Published by Elsevier in Data in Brief
- Vol. 62, 111901
Data on ethanol drinking and microbiome in alcohol preferring and non-preferring rats after a fecal microbiota transplant
- 20 July 2025
- Published by Elsevier in Data in Brief
- Vol. 62, 111901
Open AccessData on ethanol drinking and microbiome in alcohol preferring and non-preferring rats after a fecal microbiota transplant
- Open Access
- 14 May 2025
- Published by Elsevier in Drug and Alcohol Dependence
- Vol. 272, 112710
Enhanced pharmacokinetic and pharmacodynamic effects of combining p-fluorofentanyl and fentanyl
- 14 May 2025
- Published by Elsevier in Drug and Alcohol Dependence
- Vol. 272, 112710
Open AccessEnhanced pharmacokinetic and pharmacodynamic effects of combining p-fluorofentanyl and fentanyl
- Open Access
- 17 October 2024
- Published by Springer Nature in Archives of Toxicology
- Vol. 99 (1) , 287-297
In vivo pharmacokinetic, pharmacodynamic and brain concentration comparison of fentanyl and para-fluorofentanyl in rats
In 2022, para-fluorofentanyl (pFF) rose to the 6th most reported drug and the most reported fentanyl analog in the United States according to the Drug Enforcement Administration (DEA). pFF differs from fentanyl by the addition of a single fluorine group. To date, pFF has not been extensively evaluated in vivo and is frequently seen in combination with fentanyl. In the present study, the pharmacodynamic (PD) and pharmacokinetic (PK) properties and brain region-specific concentrations of pFF were evaluated in male Sprague–Dawley rats and compared to fentanyl. A 300 μg/kg subcutaneous dose of fentanyl or pFF was administered to assess PD and PK parameters as well as brain region concentrations. PD parameters were evaluated via a tail flick test to evaluate analgesia and core body temperature to measure hypothermia, a surrogate marker of overall opioid toxicity. Fentanyl and pFF were found to be equally active at the tested dose in terms of tail flick response with both compounds producing an analgesic response that lasted up to 240 min post-drug treatment. pFF induced a significantly greater hypothermic effect compared to fentanyl with a maximum temperature decrease of −5.6 ℃. Plasma PK parameters (T1/2, AUC, etc.) did not differ between fentanyl and pFF. However, pFF concentrations in the medulla, hippocampus, frontal cortex and striatum were more than two times the fentanyl concentrations. The increase in brain concentrations and greater hypothermic effect suggests that pFF is potentially more dangerous than fentanyl.
- 17 October 2024
- Published by Springer Nature in Archives of Toxicology
- Vol. 99 (1) , 287-297
Open AccessIn vivo pharmacokinetic, pharmacodynamic and brain concentration comparison of fentanyl and para-fluorofentanyl in rats
In 2022, para-fluorofentanyl (pFF) rose to the 6th most reported drug and the most reported fentanyl analog in the United States according to the Drug Enforcement Administration (DEA). pFF differs from fentanyl by the addition of a single fluorine group. To date, pFF has not been extensively evaluated in vivo and is frequently seen in combination with fentanyl. In the present study, the pharmacodynamic (PD) and pharmacokinetic (PK) properties and brain region-specific concentrations of pFF were evaluated in male Sprague–Dawley rats and compared to fentanyl. A 300 μg/kg subcutaneous dose of fentanyl or pFF was administered to assess PD and PK parameters as well as brain region concentrations. PD parameters were evaluated via a tail flick test to evaluate analgesia and core body temperature to measure hypothermia, a surrogate marker of overall opioid toxicity. Fentanyl and pFF were found to be equally active at the tested dose in terms of tail flick response with both compounds producing an analgesic response that lasted up to 240 min post-drug treatment. pFF induced a significantly greater hypothermic effect compared to fentanyl with a maximum temperature decrease of −5.6 ℃. Plasma PK parameters (T1/2, AUC, etc.) did not differ between fentanyl and pFF. However, pFF concentrations in the medulla, hippocampus, frontal cortex and striatum were more than two times the fentanyl concentrations. The increase in brain concentrations and greater hypothermic effect suggests that pFF is potentially more dangerous than fentanyl.
- 17 September 2024
- Published by Elsevier in Journal of Substance Use and Addiction Treatment
- Vol. 167, 209506
State program enables the identification of factors associated with naloxone awareness, self-efficacy, and use for overdose reversal: A cross-sectional, observational study in an urban emergency department population
- 17 September 2024
- Published by Elsevier in Journal of Substance Use and Addiction Treatment
- Vol. 167, 209506
State program enables the identification of factors associated with naloxone awareness, self-efficacy, and use for overdose reversal: A cross-sectional, observational study in an urban emergency department population
- Open Access
- 12 September 2024
- Published by Elsevier in Chemico-Biological Interactions
- Vol. 403, 111242
Association of genetic variants in CYP3A5, DRD2 and NK1R with opioid overdose
- 12 September 2024
- Published by Elsevier in Chemico-Biological Interactions
- Vol. 403, 111242
Open AccessAssociation of genetic variants in CYP3A5, DRD2 and NK1R with opioid overdose
- Open Access
- 1 August 2024
- Published by Springer Nature in The Pharmacogenomics Journal
- Vol. 24 (4) , 1-6
Opioid use disorder risk alleles in self-reported assigned African American/Afro-Caribbean and European biogeographical genetic ancestry groups and in males and females
The influence of genetic variants related to opioid use disorder (OUD) was evaluated using multiple logistic regression analysis in self-reported assigned African American/Afro-Caribbean and European biogeographical ancestry groups (BGAGs) and by sex. From a sample size of 1301 adult patients (>18 years of age) seen in emergency departments of three medical centers in Ohio, six variants were found to be associated with OUD. Two of the variants, rs2740574 (CYP3A4) and rs324029 (DRD3), were included in the analysis having met criteria of at least five subjects for each BGAG, variant carrier status, and OUD status combinations. Variant carriers in the African/Afro-Caribbean BGAG had slightly lower predicted probabilities of OUD. Variant carriers in the European BGAG had slightly higher predicted probabilities of OUD. Relative to sex, all the six variants met evaluation criteria (five subjects for all sex, variant, and OUD status combinations). No statistically significant interactions were found between a given variant, BGAGs and sex. Findings suggest variant testing relative to OUD risk can be applied across BGAGs and sex, however, studies in larger populations are needed.
- 1 August 2024
- Published by Springer Nature in The Pharmacogenomics Journal
- Vol. 24 (4) , 1-6
Open AccessOpioid use disorder risk alleles in self-reported assigned African American/Afro-Caribbean and European biogeographical genetic ancestry groups and in males and females
The influence of genetic variants related to opioid use disorder (OUD) was evaluated using multiple logistic regression analysis in self-reported assigned African American/Afro-Caribbean and European biogeographical ancestry groups (BGAGs) and by sex. From a sample size of 1301 adult patients (>18 years of age) seen in emergency departments of three medical centers in Ohio, six variants were found to be associated with OUD. Two of the variants, rs2740574 (CYP3A4) and rs324029 (DRD3), were included in the analysis having met criteria of at least five subjects for each BGAG, variant carrier status, and OUD status combinations. Variant carriers in the African/Afro-Caribbean BGAG had slightly lower predicted probabilities of OUD. Variant carriers in the European BGAG had slightly higher predicted probabilities of OUD. Relative to sex, all the six variants met evaluation criteria (five subjects for all sex, variant, and OUD status combinations). No statistically significant interactions were found between a given variant, BGAGs and sex. Findings suggest variant testing relative to OUD risk can be applied across BGAGs and sex, however, studies in larger populations are needed.
- Open Access
- 24 June 2024
- Published by Springer Nature in Scientific Reports
- Vol. 14 (1) , 1-11
Bile acids and the gut microbiome are involved in the hyperthermia mediated by 3,4-methylenedioxymethamphetamine (MDMA)
Hyperthermia induced by phenethylamines, such as 3,4–methylenedioxymethamphetamine (MDMA), can lead to life-threatening complications and death. Activation of the sympathetic nervous system and subsequent release of norepinephrine and activation of uncoupling proteins have been demonstrated to be the key mediators of phenethylamine-induced hyperthermia (PIH). Recently, the gut microbiome was shown to also play a contributing role in PIH. Here, the hypothesis that bile acids (BAs) produced by the gut microbiome are essential to PIH was tested. Changes in the serum concentrations of unconjugated primary BAs cholic acid (CA) and chenodeoxycholic acid (CDCA) and secondary BA deoxycholic acid (DCA) were measured following MDMA (20 mg/kg, sc) treatment in antibiotic treated and control rats. MDMA-induced a significant hyperthermic response and reduced the serum concentrations of three BAs 60 min post-treatment. Pretreatment with antibiotics (vancomycin, bacitracin and neomycin) in the drinking water for five days resulted in the depletion of BAs and a hypothermic response to MDMA. Gut bacterial communities in the antibiotic-treated group were distinct from the MDMA or saline treatment groups, with decreased microbiome diversity and alteration in taxa. Metagenomic functions inferred using the bioinformatic tool PICRUSt2 on 16S rRNA gene sequences indicated that bacterial genes associated to BA metabolism are less abundant in the antibiotic-MDMA treated group. Overall, these findings suggest that gut bacterial produced BAs might play an important role in MDMA-induced hyperthermia.
- 24 June 2024
- Published by Springer Nature in Scientific Reports
- Vol. 14 (1) , 1-11
Open AccessBile acids and the gut microbiome are involved in the hyperthermia mediated by 3,4-methylenedioxymethamphetamine (MDMA)
Hyperthermia induced by phenethylamines, such as 3,4–methylenedioxymethamphetamine (MDMA), can lead to life-threatening complications and death. Activation of the sympathetic nervous system and subsequent release of norepinephrine and activation of uncoupling proteins have been demonstrated to be the key mediators of phenethylamine-induced hyperthermia (PIH). Recently, the gut microbiome was shown to also play a contributing role in PIH. Here, the hypothesis that bile acids (BAs) produced by the gut microbiome are essential to PIH was tested. Changes in the serum concentrations of unconjugated primary BAs cholic acid (CA) and chenodeoxycholic acid (CDCA) and secondary BA deoxycholic acid (DCA) were measured following MDMA (20 mg/kg, sc) treatment in antibiotic treated and control rats. MDMA-induced a significant hyperthermic response and reduced the serum concentrations of three BAs 60 min post-treatment. Pretreatment with antibiotics (vancomycin, bacitracin and neomycin) in the drinking water for five days resulted in the depletion of BAs and a hypothermic response to MDMA. Gut bacterial communities in the antibiotic-treated group were distinct from the MDMA or saline treatment groups, with decreased microbiome diversity and alteration in taxa. Metagenomic functions inferred using the bioinformatic tool PICRUSt2 on 16S rRNA gene sequences indicated that bacterial genes associated to BA metabolism are less abundant in the antibiotic-MDMA treated group. Overall, these findings suggest that gut bacterial produced BAs might play an important role in MDMA-induced hyperthermia.