Next Article in Journal
A Metabolomics-Based Screening Proposal for Colorectal Cancer
Next Article in Special Issue
In Vitro Metabolic Study of Four Synthetic Cathinones: 4-MPD, 2-NMC, 4F-PHP and bk-EPDP
Previous Article in Journal
A New Look at Novel Cardiovascular Risk Biomarkers: The Role of Atherogenic Lipoproteins and Innovative Antidiabetic Therapies
Previous Article in Special Issue
The Application of Metabolomics in Forensic Science with Focus on Forensic Toxicology and Time-of-Death Estimation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Postmortem Metabolomics Reveal Acylcarnitines as Potential Biomarkers for Fatal Oxycodone-Related Intoxication

1
Department of Forensic Genetics and Forensic Toxicology, National Board of Forensic Medicine, 587 58 Linköping, Sweden
2
Department of Biomedical and Clinical Sciences, Division of Clinical Chemistry and Pharmacology, Linköping University, 581 83 Linköping, Sweden
*
Author to whom correspondence should be addressed.
Submission received: 14 December 2021 / Revised: 18 January 2022 / Accepted: 20 January 2022 / Published: 25 January 2022
(This article belongs to the Special Issue Metabolomics in Forensic Sciences)

Abstract

:
Postmortem metabolomics has recently been suggested as a potential tool for discovering new biological markers able to assist in death investigations. Interpretation of oxycodone concentrations in postmortem cases is complicated, as oxycodone tolerance leads to overlapping concentrations for oxycodone intoxications versus non-intoxications. The primary aim of this study was to use postmortem metabolomics to identify potential endogenous biomarkers that discriminate between oxycodone-related intoxications and non-intoxications. Ultra-high performance liquid chromatography-quadrupole time-of-flight mass spectrometry data from 934 postmortem femoral blood samples, including oxycodone intoxications and controls positive and negative for oxycodone, were used in this study. Data were processed and evaluated with XCMS and SIMCA. A clear trend in group separation was observed between intoxications and controls, with a model sensitivity and specificity of 80% and 76%. Approximately halved levels of short-, medium-, and long-chain acylcarnitines were observed for oxycodone intoxications in comparison with controls (p < 0.001). These biochemical changes seem to relate to the toxicological effects of oxycodone and potentially acylcarnitines constituting a biologically relevant biomarker for opioid poisonings. More studies are needed in order to elucidate the potential of acylcarnitines as biomarker for oxycodone toxicity and their relation to CNS-depressant effects.

Graphical Abstract

1. Introduction

Metabolomics strives at quantifying as many low weight molecules as possible within a well determined biological sample, i.e., the metabolome [1]. Metabolomics has proven suitable both for biomarker discovery and for generating new biological hypotheses, and since its introduction has been applied within several of different scientific fields [2]. So far, the number of metabolomics applications used for death investigations are limited, and most efforts have focused on using metabolomics for predicting time since death [3]. However, we recently demonstrated on autopsy cases with pneumonia that postmortem metabolomics could potentially aid in determining the cause of death [4]. Postmortem metabolomics offers a new and unique possibility within the field of forensic toxicology as the postmortem metabolome probably reflects the events leading up to death, the events after death, and potentially also the cause of death.
Overconsumption of opioids can potentially produce life-threatening respiratory depression through their action on μ-opioid receptors. Studies indicates that opioids decrease the human response to hypoxia, which might lead to irregular breathing and, for some cases, complete cessation of rhythmic respiratory function [5]. Other CNS-depressants such as benzodiazepines and ethanol can lead to harmful synergistic effects, further increasing the risk of respiratory depression [6,7]. Oxycodone is a semi-synthetic opioid that is commonly prescribed due to its analgesic properties. Between 2001 and 2015, there was a significant increase in prescription opioid overdose deaths [8,9]. Compared to other Nordic countries, Sweden has experienced a consistent increase in the prevalence of prescribed oxycodone [10]. Expanding the use of oxycodone has led to a higher count of unintentional overdose deaths, and simultaneous use of alcohol and other CNS-depressants may further add to this risk. However, interpretation of oxycodone concentrations in postmortem cases is complicated, as oxycodone tolerance leads to overlapping concentrations for oxycodone intoxications and non-intoxications [11,12].
Hence, the aim of this study was to investigate if postmortem metabolomics could be used to identify metabolic differences between intoxication cases positive for oxycodone in comparison to two control groups, positive and negative for oxycodone. Furthermore, this study aimed to link potential identified biomarkers to hypoxia and the toxicological effect of oxycodone. We found a significant decrease in femoral blood for 25 acylcarnitines in intoxication cases in comparison to control cases. The homeostatic imbalance of acylcarnitines could potentially be linked to respiratory depression and the effect of oxycodone intoxication.

2. Results

Postmortem metabolomics provides a biochemical overview of the postmortem metabolome. This overview might provide vital insight into the agonal period and the cause of death. In this study, the postmortem metabolome of oxycodone positive intoxication cases were compared with two control groups, positive and negative for oxycodone. The primary aims were to find potential biomarkers for oxycodone intoxications and to investigate if these biomarkers could be linked to respiratory function and hypoxia.

2.1. Study Cohort

A demographic overview of the 934 autopsy cases included in this study are presented in Table 1. A skewed demography was observed between oxycodone intoxications and controls (Table 1). In general, the oxycodone-related intoxication cases included more females, were younger, and had a higher body mass in comparison to the two control groups. In addition, other CNS-depressants, such as benzodiazepines, were more common for oxycodone-related intoxications in comparison to the two control groups. In Supplementary Table S1, all analytical findings of other opioids and common prescribed CNS-depressants are presented.

2.2. Multivariate Models

Out of 934 autopsy cases, 625 were used in a training set for model building. A total number of 1397 chromatographic peaks with specific accurate masses and retention times were included in the PCA model, as seen in Supplementary Figure S1. The PCA model showed no clear outliers or unwanted trends related to gender, age, or BMI. Furthermore, no trend related to acquisition date could be observed, indicating no or limited batch-to-batch variation. In addition to this, the internal standards in the blank whole blood samples analyzed over the whole study period (n = 4065) showed an absolute area variation of <26% (CV), a maximum mass accuracy deviation of 5.7 ppm, and a standard deviation for the retention times less than 2.5 s, altogether indicating adequate data quality.
All chromatographic peaks in the PCA model except chromatographic peaks belonging to oxycodone, noroxycodone, oxycodol, and noroxycodol (including different forms of adducts, isotopes and insource fragments) were then used in an OPLS-DA model in order to identify potential markers able to differentiate oxycodone-related intoxication from oxycodone nonrelated intoxication. In the final OPLS-DA model, without any oxycodone-related peaks, an apparent group separation was observed between oxycodone intoxications and positive control cases, as seen in Figure 1a. The OPLS-DA model described the data well, and the cross-validation score was adequate, indicating that the model was reproducible (R2 = 0.41, Q2 = 0.21).
The remaining 309 autopsy cases were used in a validation set enabling a full validation of the multivariate model, thereby minimizing the risk of overfitting models. The group separation for the validation set was comparable to the training set, as seen in Figure 1b. For the validation set, the model’s sensitivity and specificity were 80% and 74%, while the sensitivity and specificity were 80% and 80% for the training set. In depth analysis of false positives (positive controls classified as intoxications by the multivariate model) showed that several autopsy cases had high levels of oxycodone (or other opioids), and/or ethanol or benzodiazepines even though the primary cause of death was not intoxication (e.g., hanging or drowning). In depth analysis of false negatives (oxycodone-related intoxications classified as positive controls by the multivariate model) showed that several autopsy cases had low levels of oxycodone but high levels of other non-CNS-depressant drugs such as acetaminophen, as seen in Supplementary Table S2.

2.3. Metabolites of Importance

Due to the skewed demographics, several supervised OPLS-models were built in order to ensure interpretation of the results was not affected by the systematic differences. Within each separate group, no or limited trends were observed for the supervised OPLS-models based on gender and BMI (Q2 < 0.028, data not shown). However, several chromatographic peaks correlated with age, and an OPLS model based on age as a y-variable for the negative control is presented in the supplementary material, Supplementary Figure S2. In addition, to further evaluate the influence of age, a new age matched OPLS-DA model was developed. In this model, two thirds of the samples with an age <56 years in the intoxication and two thirds of the samples with an age >56 years in the positive controls were excluded. In both groups, a mean of 56 years was acquired and an age mean comparison with student’s t-test gave a p-value of 0.92. The age matched OPLS-DA score plot showed similar results to the model using the entire training set but with a reduced predictive power. The corresponding volcano plot of variables of importance over p(corr) showed that acylcarnitines were the most prominent metabolites. The OPLS-DA score plot and volcano plot are presented in the supplementary material, as seen in Figure S3. Even so, in order not to neglect the skewed demographics, only features important for discriminating intoxications from both control groups were considered, and a criteria of p(corr) > 0.2, or p(corr) < −0.2 was used. This is illustrated in a shared and unique structure plot (SUS-plot), as seen in Figure 2. Only chromatographic peaks significantly different between intoxication and both control groups were considered, which are highlighted as red squares in Figure 2.
Furthermore, to ensure that postmortem interval (PMI) differences were not neglected, an age and PMI compensated model was developed, as seen in Figure S4. An increased group separation was observed between intoxication cases versus positive controls but with a reduced model robustness. In addition, a postmortem degradation modified OPLS-model was developed. In this model, all autopsy cases where the forensic pathologist had documented a beginning decomposition, (e.g., some/medium/extensive decomposition or maggots), were excluded. A similar trend in group separation was observed, and no apparent differences in those chromatographic peaks that were important for group discrimination were observed, as seen in Supplementary Figure S5.
All chromatographic peaks identified and important for group separation belonged to the biochemical class acylcarnitines, as seen in Table 2. All acylcarnitines were significantly decreased for intoxication cases compared to control cases. For example, propionylcarnitine, heptanoylcarnitine, and hexadecadienoylcarnitine that belong to short-, medium-, and long-chain acylcarnitines, respectively, all had a significant difference in fold change of approximately ½ for intoxication in comparison to positive controls, as seen in Figure 3. One unidentified chromatographic peak with m/z 307.122 was significantly higher in the intoxication group in comparison with the controls.

3. Discussion

In simple terms, postmortem metabolomics gives a vast overview of the biochemical composition in a given metabolome after death. That overview has primarily been used for postmortem interval (PMI) determination [3]. However, postmortem metabolomics might also provide vital information about the agonal period as well as the biological processes leading up to death [4]. Interpretation of oxycodone concentrations in postmortem cases is complicated due to tolerance and overlapping concentrations for fatal and non-fatal levels. Therefore, this study aimed at identifying biomarkers able to discriminate between intoxication and non-intoxications and possibly link these changes to respiratory depression and hypoxia.

3.1. Study Strengths and Limitations

Due to the sheer number of autopsy cases, this study was able to use a two-sets study design, with a training set that was cross-validated and an untouched validation set enabling an external validation. Using a two-sets study design to fully validate the multivariate model gave a strong foundation for this postmortem metabolomics investigation.
Postmortem samples are bound to show high inter-individual variations, as death itself and time since death most likely result in extensive biochemical changes that might be unrelated to the cause of death [14]. The PMI was only known for a small portion of the included autopsy cases, making it difficult to interpret PMI’s impact on the multivariate models. However, as suggested by Chigine et al., a PMI control group was created, as seen in Figure S4 [14]. In this age and PMI-compensated model, an improved group separation was observed. Unfortunately, the model robustness decreased, probably owing to the reduced number of samples. In addition, the multivariate model only including samples without decomposition showed similar results to the previous models. The discriminative power of the OPLS models with intoxications versus controls were highly intriguing, and similar results were observed both for the training set and the validation set, confirming that no model overfitting occurred. In addition, as neither acquisition date nor PMI or postmortem degradation seemed to skew or conceal any results, the study data were considered robust and reliable.
It is important to mention that all supervised OPLS models assume that the group classification is perfect. For this study, perfect diagnostic accuracy by the forensic pathologist is assumed. We have not found any studies that estimate the error rate of forensic pathological intoxication diagnosis with regard to the correct classification of involved substances or with regard to other potential differential diagnoses. However, in a clinical setting, it has been estimated that approximately 10–15% of diagnoses are incorrect [15]. While these percentages might not be comparable to a forensic pathological setting, it can also not be assumed that the forensic pathologist is correct by default. Therefore, it is very difficult to evaluate the model’s sensitivity (number of false negative) and the specificity (number of false positives). In addition, as the model is based on biochemical changes, some misclassifications by the model could be due to the same biochemical processes being present but to a lesser extent. In order to evaluate this aspect, a small set of false negatives and false positives were inspected in depth. Several of the positive controls, falsely classified as intoxications, had intoxication as a secondary contributing cause of death. For these cases, a high concentration of oxycodone or other opioids or the presence of other CNS depressants were observed, even though the primary cause of death was something else (e.g., hanging or drowning). Looking at the false negatives, some were intoxications with non-opioid substances (for example acetaminophen). In this aspect, it could thereby be argued that these are not true false positives or false negatives as the metabolic pattern was explainable.
It is also important to mention that the model includes all intoxication where oxycodone is present, even though the concentration might be rather low and the contribution to the primary cause of death could be questionable. The diagnostic accuracy might have been improved by manually evaluating the intoxication cases and only include cases where the contribution from oxycodone is certain. This approach, however, would have been very tedious, and the final multivariate model would only be able to handle a very small and homogeneous set of samples. In addition, a stringent selection might miss vital results, as it is well known that a combination of oxycodone together with other CNS-depressants further increases the risk of respiratory depression [6,7]. However, it is important in future studies to include more CNS-depressants in order to investigate their general effect on acylcarnitine homeostasis. In any case, some error is to be expected as neither a model nor a forensic pathologist can be expected to have complete diagnostic accuracy. With the current sensitivity and specificity, the model is able to highlight cases that merit a closer review in order to confirm or exclude an intoxication diagnosis.

3.2. Acylcarnitines Potential as Biomarker for Oxycodone-Related Intoxications

Acylcarnitines have been linked to a number of age-related diseases where most commonly elevated levels of acylcarnitines are observed. For example, elevated levels of acylcarnitines have been suggested as a marker for cardiovascular disease [16,17] and diabetes [18,19]. However, Jarrell et al. recently demonstrated that long-chain and very long-chain acylcarnitines increase with age [20]. These results are troubling as there is an age-related systematic difference between the oxycodone-related intoxications and positive controls in this study. In order to ensure the validity of the results in this postmortem metabolomics investigation, several multivariate models were developed in order to investigate a possible age bias. In the age-matched OPLS-DA model, an apparent group separation was still observed, and acylcarnitines showed high model impact and correlation, as seen in Supplementary Figure S3. In addition, the age-based OPLS model showed that acylcarnitines had a poor model impact and a low model correlation in comparison to the OPLS-DA models of intoxications versus positive controls, as seen in Figure S4. To bypass the age-related differences, only results differentiating oxycodone intoxications from both control groups are presented. This minimized the confounding effects of age, as there were no age-related differences between intoxications and negative controls.
Acylcarnitines are considered important metabolic intermediates in mitochondrial metabolism [21], and the decreased levels of acylcarnitines observed in our study could therefore be an effect of oxygen depletion due to the presumed respiratory depression for intoxication cases. The are several relevant in vivo studies that report that hypoxia and ischemia affect the homeostasis of acylcarnitines, thereby supporting our results [22,23]. Bjorndalen et al. uses an animal model with a TTP diet (known inhibitor of β-oxidation and mitochondrial function) with the primary aim to evaluate mitochondrial function and respiratory activity. Interestingly, Bjorndalen et al. observes a hepatic lipid accumulation together with a reduction of short chain acylcarnitines in plasma for the rats on a TTP diet [22]. Similarly, Hal et al. observed drastic decreases in levels of acylcarnitines by measuring total carnitine during neonatal hypoxia [23]. Bruder and Raff also observe similar results in neonatal rats (Postnatal Day 2) during acute hypoxia; however, increased levels of plasma acylcarnitines were observed in older rats (Postnatal Day 60) [24]. Whitmer et al. found decreased concentrations of acylcarnitines together with an accumulation of long-chain acylcarnitines in ischemic and hypoxic rat hearts [25]. Friolet et al. demonstrated increased concentrations of short-chain acylcarnitines in plasma after exhaustive exercise under both hypoxic and normoxic conditions [26]. The observed decreases of acylcarnitines in our study might reflect a prolonged hypoxic state and an increase in oxidative stress, resulting in a depletion of acylcarnitine reserves.
Furthermore, several in vitro studies regarding ischemia in cerebral rat cells show that acylcarnitines (primarily acetylcarnitine) have a neuroprotective effect. Pretreatment with acylcarnitines before induced ischemia showed increased cell survival, increased cytochrome c oxidase activity, and reduced oxidative stress compared to controls, among other protective parameters. Even though the mechanism is not fully understood, there are suggested theories that include maintenance of mitochondrial proteins, protection against oxidative stress, and inhibition of apoptosis [27,28,29,30,31]. As high levels of acylcarnitines are protective in ischemia, perhaps the low level observed in this study with a terminal outcome reflects a sustained ischemia in which the acylcarnitines are consumed.
Even though controversy exists regarding which acylcarnitines are increased, decreased, or unaltered in whole blood or plasma during hypoxia and ischemia, the potential of postmortem metabolomics for biomarker discovery is unquestionable. Further studies are needed to fully elucidate the potential of acylcarnitines as biomarkers for hypoxia and oxycodone toxicity and acylcarnitines’ relation to age and other CNS-depressants.

4. Materials and Methods

All autopsy cases admitted to the Swedish National Board of Forensic Medicine between late June 2017 until end of October 2020 with femoral blood and with a toxicological screening using high-resolution mass-spectrometry were considered for this study (n = 17,008). All intoxications with ICD9-codes 965, 967, 969, 970, 977, 980, and 995 and positive for oxycodone were included in the intoxication group (oxycodone intoxication, n = 375). All other autopsy cases positive for oxycodone were grouped into the oxycodone-positive control group (Positive Controls, n = 364). A third, oxycodone-negative control group with the inclusion criteria ICD9-codes 994, 958, 933 869, 861, 852, and 804; BMI 18.5–30.0; age 20–75; and the exclusion criteria of intoxications (including insulin, cyanide, and carbon monoxide), diseases potentially affecting the metabolic fingerprint (such as cancer, cardiovascular diseases, sepsis and pneumonia), and signs of alcohol damage (Negative Controls, n = 195). Each parameter in the demographic overview was compared between the three groups with the statistical methods χ2 and one-way ANOVA.
All autopsy cases were analyzed with a standardized procedure described elsewhere [32]. In short, each femoral blood sample was prepared by protein precipitation (MeCN:EtOH, 90:10), including an addition of three internal standards (amphetamine-d8, diazepam-d5 and mianserin-d3). All samples were injected on a UHPLC-ESI-QToF system (Agilent 6540 QTOF with a Jet Stream interface and an Agilent 1290 Infinity LC instrument). Separation was performed on a C18 column (Waters Acquity, HSS T3 column; 150 mm × 2.1 mm, 1.8 μm) using gradient elution. MS data were collected in positive mode, and the total acquisition time for each sample was 12 min. Each analytical run included a blank drug-free bovine whole blood sample (purchased from a local slaughterhouse) also containing the three internal standards, analyzed in the beginning and at the end of each run.
The raw LC/MS data from the selected autopsy cases were exported to mzData-files using Masshunter. The exported files were loaded into R (2.14) and the XCMS package for peak detection and retention time alignment [33]. In XCMS, the centWave algorithm was used for peak detection using the following parameters: Δm/z of 15 ppm, minimum peak width of 3 s, maximum peak width of 20 s, and signal-to-noise threshold of 3 with noise variable set to 1500. Retention time correction was performed using the obiwarp function, and for the grouping, an m/z width of 0.05, base width of 3, and minimum fraction of 0.8 were used.
All chromatographic peaks before 90 s and after 660 s were excluded, and the remaining features were normalized in Excel using the probabilistic quotient normalization. All autopsy cases (n = 934) were ordered after class and acquisition date, and every third case was assigned to a validation set (n = 309) used for external validation, while the remaining were used in a training set (n = 625) for model building. All variables were scaled with unit variance, log transformed, and subjected to multivariate analysis using SIMCA 15.0.2 (Umetrics, Umeå, Sweden). Principal component analysis (PCA) was used to give an overview of the data, enabling identification of outliers and observation of trends. PLS models for each variable (i.e., age, gender, weight, length, and BMI) were created to investigate systematic differences in the metabolic profiles. Orthogonal partial least-squares discriminant analysis (OPLS-DA) was used to identify features directly related to oxycodone, i.e., isotopes, different adducts, in-source fragments, and oxycodone metabolites. All oxycodone-related features were excluded, and two OPLS-DA models of intoxication (n = 249) vs. positive controls (n = 244) and intoxication (n = 249) vs. negative controls (n = 132) were developed. The OPLS-DA model was used to identify features important for group separation with a potential link to oxycodone toxicity. Variable importance for the projection plots (VIP) and shared unique structure plots (SUS) were used to identify potential biomarkers important for discrimination of oxycodone intoxication. Metabolites with a p(corr) > 0.2 or <−0.2 were putatively identified and annotated by matching molecular weight (±5 ppm) with the Human Metabolome Database, ref [34] as well as the METLIN database [35]. For visualization, each metabolites’ normalized mean intensity was also statistically evaluated with univariate analysis using student’s t-test. All presented p-values were Bonferroni-corrected (Microsoft Office Excel 2013). Several OPLS and OPLS-DA models were developed to check and correct for systematic differences related to age, gender, BMI, PMI, grade of decomposition, and acquisition date. For model specifics, the reader is referred to the Figure legends in the supplemental material.
Experimental reproducibility was assessed by determination of the coefficients of variation (CV) for the isotope-labeled standards in blank whole blood samples and by visual examination of the score plot from the PCA. All multivariate models has been validated using k-fold cross-validation, and the OPLS-DA model of oxycodone intoxications versus positive controls was also evaluated using external validation. For the cross-validation, samples were divided into seven cross-validation groups (k = 7). Each group was, in turn, kept out of model development, and the kept out part was then predicted by the model. After this, the predictions of the kept out part were compared with the actual values. These steps were repeated until all parts (samples) had been kept out. For all included components, SIMCA computes an overall Q2 = 1 − PRESS/SS, where the prediction error sum of squares (PRESS) is the squared differences between observed and predicted values for the Y-data kept out of the model fitting and SS is the sum of squares of Y. For the external validation, left-out autopsy cases were predicted using the corresponding OPLS-DA model. The model predictability was evaluated by determining the model’s sensitivity and specificity. False positives and negatives in the validation set with a predicted score value above 5 or below −5 were investigated in depth in order to elucidate the usability of the multivariate model.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/metabo12020109/s1, Figure S1: PCA score plot, Figure S2: Age-based OPLS model, Figure S3: Age-matched OPLS-DA model, Figure S4: Age and PMI-matched OPLS-DA model, Figure S5. Modified OPLS-DA model with regard to decomposition, Table S1: Distribution of commonly prescribed CNS-depressants, Table S2: Cause of death and analytical findings for false positives and false negatives (t > 5 or t < −5).

Author Contributions

Conceptualization, A.E., C.S., G.J., R.K. and H.G.; formal analysis, A.E., C.S. and H.G.; investigation, A.E., C.S., G.J., R.K. and H.G.; methodology, A.E., C.S., G.J., R.K. and H.G.; software, A.E.; visualization, A.E.; writing—original draft, A.E. and C.S.; writing—review and editing, A.E., C.S., G.J., R.K. and H.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

This study was approved by the Swedish Ethical Review Authority (Dnr 2019-04530).

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fiehn, O. Metabolomics—The link between genotypes and phenotypes. Plant Mol. Biol. 2002, 48, 155–171. [Google Scholar] [CrossRef] [PubMed]
  2. Newgard, C.B. Metabolomics and Metabolic Diseases: Where Do We Stand? Cell Metab. 2017, 25, 43–56. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Dawidowska, J.; Zanowska, M.K.; Markuszewski, M.J.; Kaliszan, M. The Application of Metabolomics in Forensic Science with Focus on Forensic Toxicology and Time-of-Death Estimation. Metabolites 2021, 11, 801. [Google Scholar] [CrossRef] [PubMed]
  4. Elmsjö, A.; Vikingsson, S.; Söderberg, C.; Kugelberg, F.C.; Green, H. Post-Mortem Metabolomics: A Novel Approach in Clinical Biomarker Discovery and a Potential Tool in Death Investigations. Chem. Res. Toxicol. 2021, 34, 1496–1502. [Google Scholar] [CrossRef]
  5. Pattinson, K.T.S. Opioids and the control of respiration. Br. J. Anaesth. 2008, 100, 747–758. [Google Scholar] [CrossRef] [Green Version]
  6. Van Der Schrier, R.; Roozekrans, M.; Olofsen, E.; Aarts, L.; Van Velzen, M.; De Jong, M.; Dahan, A.; Niesters, M. Influence of Ethanol on Oxycodone-induced Respiratory DepressionA Dose-escalating Study in Young and Elderly Individuals. Anesthesiology 2017, 126, 534–542. [Google Scholar] [CrossRef]
  7. Pirnay, S.O.; Mégarbane, B.; Borron, S.W.; Risède, P.; Monier, C.; Ricordel, I.; Baud, F.J. Effects of Various Combinations of Benzodiazepines with Buprenorphine on Arterial Blood Gases in Rats. Basic Clin. Pharmacol. Toxicol. 2008, 103, 228–239. [Google Scholar] [CrossRef]
  8. Rookey, B.D. Drugged Driving in the Opioid Era: Spatial Patterns of Oxycodone Use in Fatal Traffic Crashes, 2001–2016. Subst. Use Misuse 2018, 53, 1288–1298. [Google Scholar] [CrossRef]
  9. Simonsen, K.W.; Kriikku, P.; Thelander, G.; Edvardsen, H.M.E.; Thordardottir, S.; Andersen, C.U.; Jönsson, A.K.; Frost, J.; Christoffersen, D.J.; Delaveris, G.J.M.; et al. Fatal poisoning in drug addicts in the Nordic countries in 2017. Forensic Sci. Int. 2020, 313, 110343. [Google Scholar] [CrossRef]
  10. Muller, A.E.; Clausen, T.; Sjøgren, P.; Odsbu, I.; Skurtveit, S. Prescribed opioid analgesic use developments in three Nordic countries, 2006–2017. Scand. J. Pain 2019, 19, 345–353. [Google Scholar] [CrossRef]
  11. Nedahl, M.; Johansen, S.S.; Linnet, K. Postmortem Brain–Blood Ratios of Codeine, Fentanyl, Oxycodone and Tramadol. J. Anal. Toxicol. 2021, 45, 53–59. [Google Scholar] [CrossRef] [PubMed]
  12. Drobnik, W.; Liebisch, G.; Audebert, F.-X.; Fröhlich, D.; Glück, T.; Vogel, P.; Rothe, G.; Schmitz, G. Plasma ceramide and lysophosphatidylcholine inversely correlate with mortality in sepsis patients. J. Lipid Res. 2003, 44, 754–761. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Sumner, L.W.; Samuel, T.; Noble, R.; Gmbh, S.D.; Barrett, D.; Beale, M.H.; Hardy, N. Proposed minimum reporting standards for chemical analysis Chemical Analysis Working Group (CAWG) Metabolomics Standards Initiative (MSI). Metabolomics 2007, 3, 211–221. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Chighine, A.; Locci, E.; Nioi, M.; D’Aloja, E. Looking for Post-Mortem Metabolomic Standardization: Waiting for Godot—The Importance of Post-Mortem Interval in Forensic Metabolomics. Chem. Res. Toxicol. 2021, 34, 1946–1947. [Google Scholar] [CrossRef]
  15. Graber, M.L. The incidence of diagnostic error in medicine. BMJ Qual. Saf. 2013, 22, ii21–ii27. [Google Scholar] [CrossRef] [Green Version]
  16. Kalim, S.; Clish, C.B.; Wenger, J.; Elmariah, S.; Yeh, R.W.; Deferio, J.J.; Pierce, K.; Deik, A.; Gerszten, R.E.; Thadhani, R.; et al. A Plasma Long-Chain Acylcarnitine Predicts Cardiovascular Mortality in Incident Dialysis Patients. J. Am. Heart Assoc. 2013, 2, e000542. [Google Scholar] [CrossRef] [Green Version]
  17. Koh, A.S.; Gao, F.; Liu, J.; Fridianto, K.T.; Ching, J.; Tan, R.S.; Wong, J.-I.; Chua, S.J.; Leng, S.; Zhong, L.; et al. Metabolomic profile of arterial stiffness in aged adults. Diabetes Vasc. Dis. Res. 2018, 15, 74–80. [Google Scholar] [CrossRef] [Green Version]
  18. Bouchouirab, F.-Z.; Fortin, M.; Noll, C.; Dubé, J.; Carpentier, A.C. Plasma Palmitoyl-Carnitine (AC16:0) Is a Marker of Increased Postprandial Nonesterified Incomplete Fatty Acid Oxidation Rate in Adults with Type 2 Diabetes. Can. J. Diabetes 2018, 42, 382–388. [Google Scholar] [CrossRef]
  19. Adams, S.H.; Hoppel, C.L.; Lok, K.H.; Zhao, L.; Wong, S.W.; Minkler, P.E.; Hwang, D.H.; Newman, J.W.; Garvey, W.T. Plasma Acylcarnitine Profiles Suggest Incomplete Long-Chain Fatty Acid β-Oxidation and Altered Tricarboxylic Acid Cycle Activity in Type 2 Diabetic African-American Women. J. Nutr. 2009, 139, 1073–1081. [Google Scholar] [CrossRef]
  20. Jarrell, Z.R.; Smith, M.R.; Hu, X.; Orr, M.; Liu, K.H.; Quyyumi, A.A.; Jones, D.P.; Go, Y.-M. Plasma acylcarnitine levels increase with healthy aging. Aging 2020, 12, 13555–13570. [Google Scholar] [CrossRef]
  21. Dambrova, M.; Zuurbier, C.J.; Borutaite, V.; Liepinsh, E.; Makrecka-Kuka, M. Energy substrate metabolism and mitochondrial oxidative stress in cardiac ischemia/reperfusion injury. Free Radic. Biol. Med. 2021, 165, 24–37. [Google Scholar] [CrossRef]
  22. Bjørndal, B.; Alterås, E.K.; Lindquist, C.; Svardal, A.; Skorve, J.; Berge, R.K. Associations between fatty acid oxidation, hepatic mitochondrial function, and plasma acylcarnitine levels in mice. Nutr. Metab. 2018, 15, 1–13. [Google Scholar] [CrossRef] [PubMed]
  23. Çam, H.; Yildirim, B.; Aydin, A.; Say, A. Carnitine Levels in Neonatal Hypoxia. J. Trop. Pediatr. 2005, 51, 106–108. [Google Scholar] [CrossRef] [PubMed]
  24. Bruder, E.D.; Raff, H. Cardiac and plasma lipid profiles in response to acute hypoxia in neonatal and young adult rats. Lipids Health Dis. 2010, 9, 3–6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Whitmer, J.T.; Idell-Wenger, J.A.; Rovetto, M.J.; Neely, J.R. Control of fatty acid metabolism in ischemic and hypoxic hearts. J. Biol. Chem. 1978, 253, 4305–4309. [Google Scholar] [CrossRef]
  26. Friolet, R.; Hoppeler, H.; Krähenbühl, S. Relationship between the coenzyme A and the carnitine pools in human skeletal muscle at rest and after exhaustive exercise under normoxic and acutely hypoxic conditions. J. Clin. Investig. 1994, 94, 1490–1495. [Google Scholar] [CrossRef]
  27. Ferreira, G.C.; McKenna, M.C. l-Carnitine and Acetyl-l-carnitine Roles and Neuroprotection in Developing Brain. Neurochem. Res. 2017, 42, 1661–1675. [Google Scholar] [CrossRef]
  28. Bagetta, V.; Barone, I.; Ghiglieri, V.; Di Filippo, M.; Sgobio, C.; Bernardi, G.; Calabresi, P.; Picconi, B. Acetyl-l-Carnitine selectively prevents post-ischemic LTP via a possible action on mitochondrial energy metabolism. Neuropharmacology 2008, 55, 223–229. [Google Scholar] [CrossRef]
  29. Zanelli, S.A.; Solenski, N.J.; Rosenthal, R.E.; Fiskum, G. Mechanisms of Ischemic Neuroprotection by Acetyl-L-carnitine. Ann. N. Y. Acad. Sci. 2005, 1053, 153. [Google Scholar] [CrossRef] [Green Version]
  30. Barhwal, K.; Hota, S.K.; Prasad, D.; Singh, S.B.; Ilavazhagan, G. Hypoxia-induced deactivation of NGF-mediated ERK1/2 signaling in hippocampal cells: Neuroprotection by acetyl-L-carnitine. J. Neurosci. Res. 2008, 86, 2705–2721. [Google Scholar] [CrossRef]
  31. Al-Majed, A.A.; Sayed-Ahmed, M.M.; Al-Omar, F.A.; Al-Yahya, A.A.; Aleisa, A.M.; Al-Shabanah, O.A. Carnitine esters prevent oxidative stress damage and energy depletion following transient forebrain ischaemia in the rat hippocampus. Clin. Exp. Pharmacol. Physiol. 2006, 33, 725–733. [Google Scholar] [CrossRef] [PubMed]
  32. Roman, M.; Ström, L.; Tell, H.; Josefsson, M. Liquid chromatography/time-of-flight mass spectrometry analysis of postmortem blood samples for targeted toxicological screening. Anal. Bioanal. Chem. 2013, 405, 4107–4125. [Google Scholar] [CrossRef] [PubMed]
  33. Smith, C.A.; Want, E.J.; O’Maille, G.; Abagyan, R.; Siuzdak, G. XCMS: Pocessing Mass Spectrometry Data for Metabolite Profiling Using Nonlinear Peak Alignment, Matching, an Identification. Anal. Chem. 2006, 78, 779–787. [Google Scholar] [CrossRef] [PubMed]
  34. Wishart, D.S.; Jewison, T.; Guo, A.C.; Wilson, M.; Knox, C.; Liu, Y.; Djoumbou, Y.; Mandal, R.; Aziat, F.; Dong, E.; et al. HMDB 3.0—The Human Metabolome Database in 2013. Nucleic Acids Res. 2013, 41, D801–D807. [Google Scholar] [CrossRef]
  35. Smith, C.A.; O’Maille, G.; Want, E.J.; Qin, C.; Trauger, S.A.; Brandon, T.R.; Custodio, D.E.; Abagyan, R.; Siuzdak, G. METLIN: A metabolite mass spectral database. Ther. Drug Monit. 2005, 27, 747–751. [Google Scholar] [CrossRef]
Figure 1. OPLS-DA score plots. An apparent group separation is observed between oxycodone intoxications (▲) and the positive controls (■) for both the training set (a) and the validation set (b) (R2 = 0.41, Q2 = 0.21).
Figure 1. OPLS-DA score plots. An apparent group separation is observed between oxycodone intoxications (▲) and the positive controls (■) for both the training set (a) and the validation set (b) (R2 = 0.41, Q2 = 0.21).
Metabolites 12 00109 g001
Figure 2. Shared and unique structure plot (SUS-plot). Shared metabolites between the two OPLS-DA models of intoxication vs. positive controls and intoxication vs. negative controls with a p(corr) > 2 or p(corr) < −2 are highlighted as red squares.
Figure 2. Shared and unique structure plot (SUS-plot). Shared metabolites between the two OPLS-DA models of intoxication vs. positive controls and intoxication vs. negative controls with a p(corr) > 2 or p(corr) < −2 are highlighted as red squares.
Metabolites 12 00109 g002
Figure 3. Boxplots of normalized area intensities with first, second, and third quartile. The three identified metabolites illustrates that short-, medium-, and long-chain acylcarnitines are affected in a similar manner. * p-value < 0.001 and n.s. non-significant.
Figure 3. Boxplots of normalized area intensities with first, second, and third quartile. The three identified metabolites illustrates that short-, medium-, and long-chain acylcarnitines are affected in a similar manner. * p-value < 0.001 and n.s. non-significant.
Metabolites 12 00109 g003
Table 1. Demographic overview of study cohort.
Table 1. Demographic overview of study cohort.
Oxycodone IntoxicationsPositive ControlsNegative ControlsStatistics
n375364195
Females/males160/215123/24144/151p < 0.001 1
Age (yrs)48 (35–60)65 (53–74)47 (30–59)p < 0.001 2
Body weight (kg)86 (72–100)76 (62–92)74(65–83)p < 0.001 2
Body Height (cm)173 (165–181)173 (165–179)176 (170–182)p < 0.001 2
Body Mass Index (kg/m2)29 (25–33)26 (22–30)24 (22–26)p < 0.001 2
Data are presented as median with quartile range 25–75% in parentheses. 1 p-value calculated with Chi22)-test, 2 p-values calculated between the three groups with one-way ANOVA.
Table 2. Metabolite information for identified metabolites.
Table 2. Metabolite information for identified metabolites.
Metabolites 1Chain Length 2Identifier 3Mean m/z 4Exact m/z 5Δ 6% 7p-Value 8
AcetylcarnitineC2M204T119204.123204.12300.00.704.6 × 10−4
PropionylcarnitineC3M218T126218.139218.1387−1.40.673.5 × 10−12
Butyrylcarnitine C4M232T167_1232.155232.1543−3.00.642.8 × 10−7
(Iso)valerylcarnitineC5M246T226246.170246.17000.00.511.6 × 10−19
HexanoylcarnitineC6M260T289260.186260.1856−1.50.698.7 × 10−8
HeptanoylcarnitineC7M274T349274.201274.20131.10.532.6 × 10−15
OctanoylcarnitineC8M288T405288.217288.2169−0.30.557.4 × 10−12
NonanoylcarnitineC9M302T456_2302.232302.23262.00.693.6 × 10−9
DecenoylcarnitineC10:1M314T459314.232314.23261.90.598.8 × 10−3
DecanoylcarnitineC10M316T504316.248316.24820.60.467.1 × 10−8
HydroxyhexadecadiencarnitineC16:2-OHM412T568412.304412.30574.10.642.7 × 10−4
TetradecadiencarnitineC14:2M368T572368.279368.27951.40.557.1 × 10−4
DodecanoylcarnitineC12M344T573344.279344.27951.50.596.4 × 10−6
HydroxyhexadecenoylcarnitineC16:1-OHM414T588414.321414.32141.00.644.0 × 10−14
TetradecenoylcarnitineC14:1M370T591370.295370.29520.50.591.3 × 10−7
HexadecadienoylcarnitineC16:2M396T601396.311396.3108−0.50.555.8 × 10−7
Hydroxyhexadecanoyl carnitineC16-OHM416T606416.337416.33710.20.793.7 × 10−8
TetradecenoylcarnitineC14M372T608372.311372.3108−0.50.642.9 × 10−8
LinolenylcarnitineC18:3M422T611422.327422.3265−1.20.691.9 × 10−6
HydroxyoctadecanoylcarnitineC18-OHM442T617442.353442.3527−0.70.652.3 × 10−10
LinoleylcarnitineC18:2M424T628424.343424.3421−2.10.781.9 × 10−4
UnidentifiedNaM307T132307.122NANA1.593.6 × 10−10
1 Putative identities according to MSI [13], 2 number of carbons, saturations, and alcohols on the ester-group on the acylcarnitines, 3 specific identifier for each chromatographic peak, including mass and retention time, 4 mean measured accurate mass for [H]+ over all samples, 5 theoretical monoisotopic mass 6 ppm difference between measured and theoretical mass, 7 % of median of the normalized peak intensities of the intoxication group divided by the positive control group including both training and validation set, and 8 Bonferroni-corrected p-values for the log transformed normalized intensities in the intoxication group vs. positive controls.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Elmsjö, A.; Söderberg, C.; Jakobsson, G.; Green, H.; Kronstrand, R. Postmortem Metabolomics Reveal Acylcarnitines as Potential Biomarkers for Fatal Oxycodone-Related Intoxication. Metabolites 2022, 12, 109. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo12020109

AMA Style

Elmsjö A, Söderberg C, Jakobsson G, Green H, Kronstrand R. Postmortem Metabolomics Reveal Acylcarnitines as Potential Biomarkers for Fatal Oxycodone-Related Intoxication. Metabolites. 2022; 12(2):109. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo12020109

Chicago/Turabian Style

Elmsjö, Albert, Carl Söderberg, Gerd Jakobsson, Henrik Green, and Robert Kronstrand. 2022. "Postmortem Metabolomics Reveal Acylcarnitines as Potential Biomarkers for Fatal Oxycodone-Related Intoxication" Metabolites 12, no. 2: 109. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo12020109

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop