ORIGINAL RESEARCH article

Front. Genet., 23 March 2021

Sec. Behavioral and Psychiatric Genetics

Volume 12 - 2021 | https://doi.org/10.3389/fgene.2021.653443

The Potential Regulatory Network of Glutamate Metabolic Pathway Disturbance in Chinese Han Withdrawal Methamphetamine Abusers

  • 1. Shanghai Mental Health Center, Shanghai Jiao Tong University School of Medicine, Shanghai, China

  • 2. Shanghai Key Laboratory of Psychotic Disorders, Shanghai, China

  • 3. CAS Center for Excellence in Brain Science and Intelligence Technology (CEBSIT), Chinese Academy of Sciences, Shanghai, China

Abstract

Objects:

To explore the long-term influence of methamphetamine abuse on metabolomics character, with gas chromatography-mass spectrometry (GS-MS) technology, and the potential regulatory network using the bioinformatics method.

Methods:

Forty withdrawal methamphetamine abusers (WMA) were recruited from Shanghai Gaojing Forced Isolation Detoxification Institute. Forty healthy controls (HC) were recruited from society. GS-MS technology was used to detect metabolic products in serum. A bioinformatics method was used to build a regulatory network. Q-PCR was used to detect the candidate gene expressions, and ELISA was used to detect the regulatory enzyme expressions.

Results:

Four pathways were significantly changed in the MA compared to the HC: (1) the arginine synthesis pathway, (2) alanine, aspartic acid and glutamate metabolic pathway, (3) cysteine and methionine metabolic pathway, and (4) the ascorbate and aldarate pathway (enrichment analysis p < 0.05, Impactor factor > 0.2). When focusing on the ‘Alanine, aspartate, and glutamate metabolism’ pathway, a regulatory network was established, and the expression of candidate regulatory genes and enzymes was verified. It was found that the expression of DLG2 (Discs large MAGUK scaffold protein 2), PLA2G4 (Phospholipase A2 group IVE), PDE4D (Phosphodiesterase 4D), PDE4B (Phosphodiesterase 4B), and EPHB2 (Ephrin type-B receptor 2) were significantly different between the two groups (p < 0.05), However, after adjusting for age and BMI, only DLG2, PLA2G4, and EPHB2 remained significant (p < 0.05). The expression of enzymes was not significantly different (p > 0.05).

Conclusion:

Methamphetamine abuse influences the metabolic process in the long term, and DLG2, PLA2G4, and EPHB2 may regulate the glutamate metabolism pathway.

Introduction

Methamphetamine (METH) is a widely abused addictive psychostimulant, which accounts for a considerable share of global disease burden. METH abuse is increasing rapidly (Yang X. et al., 2018), especially in East and South-East Asia, and Oceania. Long-term abuse of METH causes serious physical and mental damage. Compared to other psychostimulants, METH is better at penetrating the center nervous system (CNS) and has a longer duration of action (Yang X. et al., 2018), making METH more addictive while the damage to the CNS is more serious and persistent. METH addicts often suffer from permanent psychotic symptoms () and cognitive decline, even when deprived from drugs, bringing great difficulties for patients who try to return to society. These symptoms are mostly attributed to METH-induced long-term neurotoxicity and excitotoxicity (). The neurochemical mechanisms underlying neurotoxicity and excitotoxicity are complex and not well understood. Glutamate, as an important excitatory neurotransmitter is thought to activate various downstream signaling pathways, for example the apoptotic pathway, and eventually produce neurotoxicity (; Yang X. et al., 2018). A previous study found that methamphetamine may result in glutamine system disturbance (Tasic et al., 2017), and glutamate + glutamine in the right inferior frontal cortex is decreased in methamphetamine users (Zheng et al., 2016).

Metabolomics is a potentially powerful tool for understanding the global biological consequences of disease and drug administration. In recent years, with the improvement of metabolomics technology, metabolomics research in psychotic diseases has made great progress (; Zhang et al., 2016; ; ). However, in the field of substance addiction, metabolomics research is still in its initial stages. A few studies focus on methamphetamine addiction, using animals as study subjects. Energy metabolic disturbance is believed to be the biochemical base of toxicity and addictive property of MTHE (Zhang et al., 2016). Neurotransmitters in serum or plasma are also disrupted by METH, especially in the glutamate-glutamine-GABA system (). Nevertheless, metabolic character differences between animals and humans still exist. This study uses withdrawal methamphetamine abusers (WMA) as objects and explore the metabolomics characters and possible regulatory mechanisms that may be the basis of the long-term neurotoxicity of methamphetamine.

Materials and Methods

Ethics Statement

Informed consent of the participants was obtained after the nature of the procedures used in the study had been fully explained. The study was carried out in accordance with the last version of the Declaration of Helsinki, and the protocol was approved by the Ethics Committee of Shanghai Mental Health Center (SMHC IRB 2016-11R). All participants were identified using numbers, instead of names.

Research Participants

Forty male WMA were recruited from Shanghai Gaojing Forced Isolation Detoxification Institute in July 2017, and 40 male healthy controls were recruited from general society. The entire sample consisted of a Chinses Han population. The WMA accepted assessment from two senior psychiatric doctors and met the diagnosis criterion for methamphetamine abuse disorder according to the Diagnostic and Statistical Manual of mental disorders, fifth edition (DSM-5). Participants did not have the opportunity to use any illicit substance during the last 2 months. They completed a self-report form that included basic demographic information and accepted an interview of Addiction Severity Index (ASI). Exclusion criteria for both groups were: (1) history of craniocerebral trauma or epileptic seizure; (2) metal implants in the brain or an implanted pacemaker; (3) history of major physical diseases, organic mental disorders, and severe mental disorders. People in the control group did not have any history of illegal substance abuse. Venous blood was collected at 7:30–8:00 am before breakfast.

Gas Chromatography – Mass Spectroscopy (GC-MS)

After standing at room temperature for 30 min, venous blood was centrifuged at 4°C, 3000g, for 15 min. The upper layer of the serum was collected. Sample preparation, GC-MS experimental process, and GC-MS data processing were carried out using previously published methods (; ). Compound identification was performed by comparing the retention indices and mass spectral data with those previously generated from reference standards of known structures present in the JiaLib metabolite database using the proprietary software XploreMET. The current JiaLib comprises over 1,200 mammalian metabolites over a 15-year accumulation. The reference chemicals present in JiaLib were commercially purchased from Sigma-Aldrich (St. Louis, MO, United States), Santa Cruz (Dallas, TX, United States), Nu-Chek Prep (Elysian, MN, United States), and synthesized in the laboratory. XploreMET adopted a specific algorithm hypergeometric, which can enrich the information of metabolites, weight the contribution of metabolites, and focus on the key metabolic network. The P and impact value was used to measure the importance of metabolic pathway.

Bioinformatics Analysis

Based on the results of untargeted metabolomics, a glutamate metabolism pathway was noticed. Glutamate is a very important excitatory neurotransmitter, and a previous study has confirmed that glutamate + glutamine in the right inferior frontal cortex was decreased in methamphetamine users (Zheng et al., 2016). In serum or plasma, the glutamate-glutamine-GABA system was also disrupted by METH (). So, we chose the ‘alanine, aspartic acid and glutamate metabolic pathway’ to build the regulatory network.

The four differential expressed metabolites were treated as candidate metabolites. We collected 133 risk genes from previous genome-wide association studies (; ). The molecular network was constructed by connecting candidate metabolites, enzymes, intermediate proteins, and risk genes. First, we investigated the relative enzymes linked to these candidate metabolic compounds according to the KEGG database and the Human Metabolome Database (HMDB). To further explore correlation networks linked to the relative enzymes and risk genes, we searched the functional proteins network database (STRING v10.5) setting the requirement that the score must be more than 700 and that no more than one intermediate protein must be present between the relative enzyme and risk gene. Finally, we built a network including four enzymes: ASNS (Asparagine synthetase), AGT (Angiotensinogen), GAD1 (Glutamate decarboxylase 1), and GAD2 (Glutamate decarboxylase 2); 13 genes: DLG2 (Discs large MAGUK scaffold protein 2), NRG1 (Neuregulin 1), PLA2G4E (Phospholipase A2 group IVE), STX8 (Syntaxin 8), NEDD4L (NEDD4 like E3 ubiquitin protein ligase), PRKG1 (Protein kinase cGMP-dependent 1), PDE4D (Phosphodiesterase 4D), PDE4B (Phosphodiesterase 4B), SERPINA5 (Serpin family A member 5), PDE6C (Phosphodiesterase 6C), CD55 (CD55 molecule), EPHB2 (Ephrin type-B receptor 2), and RAPGEF5 (Rap guanine nucleotide exchange factor 5), and several intermediate proteins.

RNA Extraction and Gene Expression Analysis

RNA was extracted from the serum samples with Trizol reagent (Invitrogen, Carlsbad, CA, United States) according to the previous protocol (), and stored at −80°C for subsequent cDNA synthesis. Following RNA quantitation, 1ug of total RNA was used for cDNA synthesis (5 × All-In-One RT Master Mix, ABM, Canada), and the quantitative gene expression analysis was performed using TB Green Fast qPCR Mix (TaKaRa, Shiga, Japan), and Lightcycle® 96 real-time PCR machine (Roche Diagnostics GmbH, Germany). Primers are listed in Table 1. GAPDH was used as a housekeeper gene. The relative gene expression was calculated using the 2−△CT method (). The CT values of duplicates were averaged for the final analysis.

TABLE 1

genePrimers (F: Forward; R: Reverse)
GAPDHF AGGGCTGCTTTTAACTCTGGT R CCCCACTTGATTTTGGAGGGA
DLG2F CCCAATGGGATGGCAGACTT R AGTATTCCCACCTCCCTCCC
NRG1F GCAACTTTGTTTCCCGGGTC R TTGGGGGCAAAAGTCACACT
PLA2G4EF AATCAGTGCTCCCTTGAGCC R GCATTACCTGGATGGGGACC
STX8F CCCCTGGTTCTCCACATACG R TCTGTCTTCGGTCCCCTTCA
NEDD4LF AGGGTTAGCTTCCTGTTGGC R TGGTACGGGGTTGAGAATGC
PRKG1F ACTTGGTGCACCTCTCAACA R AGTCATTGCCAAGGTCCCAG
PDE4DF TCTTACAGCCCACGGGGATA R AGGCAGAATCAACCCATGCT
PDE4BF GTGTCGTTCACCGTGAGAGT R CGGGGTGAAGAGAGGAGGTA
CD55F CATCTTTCCTTCGGGTTGGC R CACCATCAACACCCCTGGTT
EPHB2F AGTCTGGGGAGGGACTCATC R CTGTGTGGCCATGGAAGCTA
RAPGEF5F GGGCTTTCACTGCTACCCAT R GATTGGCTGCGCCATTTAGG
SERPINA5F GCTATGGCCCATCTGTATGCT R TTCCCCAAGGCACCTGTATG
PDE6CF TTTTGGAGAGGCACCACCTG R AACTGTTTCAAACTGCCGCT

Primers used in Q-PCR.

Elisa

The serum levels of AGT, GAD1, GAD2, and ASNS were assessed using an ELISA kit (catalog number SEC904Hu, Cloud-Clone Corp, Shanghai, China for GAD1, and catalog numbers MBS726474, MBS2501551, MBS090144; Mybiosource, Southern California, San Diego, United States, for AGT, GAD2, and ASNS, respectively).

All assays were performed according to the manufacturer’s directions and were performed in duplicate. The average value of the duplicates was used as the final value.

Statistical Analysis

General data was presented as mean and standard error of the mean (Mean ± SD), and analyzed using SSPSS 22.0. Data was analyzed using the student t test, a was 0.05 and Cohen’s d was used to calculate the effect size of the t test. 2–△CT was used to calculate the relative expression of genes, and the expression ratio of the WMA group relative to the control group was calculated as previously described ().

Metabolomics data was managed with the XploreMET system (v3.0, Metabo-Profile, Shanghai, China), and R was used to analyze the data. The NIST 2017 database was used for metabolite certification and Metaboanalyst 3.0 was used for the metabolic pathway analysis. Partial least squares discriminant analysis (PSL-DA) was performed, and a t test was used after log transformation. Differential expressed metabolites were screened according to fold change (FC) and p value. FC > 1. 5 and p < 0. 05 showed that the metabolite was significantly differentially expressed. Based on the selected differential metabolites, a pathway analysis was performed using a hypergeometric test and KEGG pathways (Yang W. et al., 2018) and the p values were calculated as follows:

PtwHits is the hit number of differential metabolites in a certain metabolic pathway. TotalAll is the number of metabolites in all pathways of a certain species, such as Homo sapiens. PtwAll is the total number of metabolites in a certain pathway and TotalHits is the hit number of differential metabolites in all pathways. A significant p value means significantly changed metabolic pathway activities that need to be focused on. For multiple tests, a false discovery rate (FDR) was applied to correct the raw p values.

Results

Demography Characteristic

The demographic characteristic are shown in Table 2. Age, BMI, and education level between the two groups were not significantly different. However, WMA individuals have lower marriage rates and higher unemployment (p < 0.05).

TABLE 2

WMA (Mean ± SD)Control (Mean ± SD)t valuep value
Age (years)34.02 ± 6.8232.02 ± 8.571.130.260
BMI (kg/m2)24.44 ± 3.6423.54 ± 3.231.140.260
Education (years)5.95 ± 3.556.97 ± 3.141.140.256
Marital status16.5370.000
Married9 (23.1%)20 (52.6%)
Divorced/Separated15 (38.4%)1 (2.6%)
Unmarried15 (38.4%)17 (44.7%)
Employment8.2710.004
Unemployed18 (46.1%)6 (15.8%)
The age began to abuse (years)25.20 ± 7.52
Withdrawal time (month)3.20 ± 1.32
Duration of abuse (years)6.51 ± 3.26
Average dosage (g)0.64 ± 0.35
Frequency
Everyday24 (60.0%)
3–5 times/week12 (30.0%)
1 time/week2 (5.0%)
1–3 times/month2 (5.0%)

Demographic characteristics of subjects.

Metabolomics Profile and Pathways

Thirty-five differentially expressed metabolites were screened by t test (p < 0.05), in which 21 metabolites could be found by searching metabolic pathways in the KEGG information database (as shown in Table 3). Since age and BMI may influence the metabolic function, we adjusted age and BMI. After that, there were 19 metabolites differentially expressed in two groups, including two kinds of alcohols, seven kinds of amino acids, one kind of carbohydrate, one kind of fatty acid, three kinds of nucleotides, and five kinds of organic acids.

TABLE 3

Metabolites speciesMetabolitesMetabolic pathwaysFC valueBefore adjusting
After adjusting
Cohen’s d
p valuep value
AlcoholsGlycerinGalactose metabolism, glyceride metabolism0.60.0000.0000.898
AlcoholsInositolGalactose metabolism, inositol metabolism, inositol phosphate metabolism, phosphatidylinositol phosphate metabolism0.70.0010.0010.813
Amino acids1-MethylhistidineHistidine metabolism1.40.0010.001–0.802
Amino acidsL- Aspartic acidAmmonia recovery, arginine and proline metabolism, aspartic acid metabolism, B-alanine metabolism, malate-aspartic acid shuttle, transcription/translation, urea cycle1.50.0010.002–0.761
Amino acidsHomocysteineBetaine metabolism, catecholamine biosynthesis, glycine and serine metabolism, homocysteine degradation, methionine metabolism1.40.0040.004–0.692
Amino acidsL- MethionineBetaine metabolism, glycine and serine metabolism, methionine metabolism, spermidine and spermine biosynthesis, transcription/translation0.90.0050.0090.658
Amino acidsL- Glutamic acidAmino sugar metabolism, ammonia recovery, glutamate metabolism, phenyl acetate metabolism, purine metabolism, pyrimidine metabolism, transcription/translation, urea cycle1.10.0110.015–0.597
Amino acidsOrnithineMetabolism of arginine and proline, metabolism of glycine and serine, biosynthesis of spermine and spermine, urea cycle1.30.0130.020–0.582
Amino acidsL- CysteineCysteine Metabolism, Glutathione Metabolism, Glycine and Serine Metabolism, Methionine Metabolism, Pantothenate and CoA Biosynthesis, Tauurine and Tauurine Metabolism, Transcription/Translation0.90.0280.0160.517
Amino acidsL- IsoleucineDegradation, transcription/translation of valine, leucine and isoleucine1.10.0290.088–0.512
Amino acidsL- ProlineArginine and proline metabolism, transcription/translation1.20.0320.068–0.502
CarbohydrateSorbitolFructose and mannose degradation, galactose metabolism0.90.0250.0270.525
Fatty acidsLinoleic acidα-linolenic acid and linoleic acid metabolism0.80.0070.0080.643
NucleotidesHypoxanthinePurine metabolism1.40.0000.000–1.046
NucleotidesXanthinePurine metabolism1.20.0000.000–0.900
NucleotidesGuanosine hydratePurine metabolism0.70.0050.0060.658
Organic acidsGlyceric acidGlyceride metabolism, glycine and serine metabolism0.70.0000.0001.023
Organic acidsPyruvic acidAlanine metabolism, amino glucose metabolism, ammonia recovery, tricarboxylic acid cycle, cysteine metabolism, gluconeogenesis, glucose-alanine cycle, glycine and serine metabolism, glycolysis, pyruvic aldehyde degradation, pyruvate metabolism, acetyl mitochondrial transfer, urea cycle0.60.0000.0000.872
Organic acidsVanillylmandelic_acidTyrosine metabolism0.80.0020.0020.744
Organic acidsHydroxypropionic_acidPropionic acid metabolism0.90.0330.0270.497
Organic acidsFumaric acidArginine and proline metabolism, aspartic acid metabolism, citric acid cycle, mitochondrial electron transport chain, phenylalanine and tyrosine metabolism, tyrosine metabolism, urea cycle0.70.0370.0470.488

Differentially expressed metabolites in pathways.

According to the results of the Metabolic Pathway Enrichment Analysis (MPEA), four pathways were found to be significantly changed; (1) arginine synthesis pathway, (2) alanine, aspartic acid and glutamate metabolic pathway, (3) cysteine and methionine metabolic pathways, and (4) ascorbate and aldarate pathway (enrichment analysis p < 0.05,Impactor factor > 0.2) (Figure 1).

FIGURE 1

Glutamine Metabolites Regulatory Network

Previous research found that methamphetamine may result in glutamine system disturbance (Tasic et al., 2017), either in the brain and serum or in plasma (Zheng et al., 2016; ). This study also found a change in the glutamate pathway. We therefore focused on the ‘Alanine, aspartate, and glutamate metabolism’ pathway, and based on the bioinformatics method, the regulatory network was built. As shown in Figure 2, we found three metabolites (pyruvate, L-aspartate, and L-glutamine), which can be regulated by four enzymes, and which can be related to 13 candidate genes. Each candidate gene influences the regulatory enzyme through one intermediate protein (Figure 2).

FIGURE 2

Candidate Gene Expression

According to the regulatory network, we used Q-PCR to verify the 13 candidate genes expressed. 2–△CT was used to calculate the relative expression of genes () (Table 4). DLG2, PLA2G4, PDE4D, PDE4B, and EPHB2 were significantly differentially expressed between the two groups (p < 0.05). However, after adjusting for age and BMI, only DLG2, PLA2G4, and EPHB2 remained significant (p < 0.05). STX8, RAPGEF5, and PDE6C were not detected because of their low expression.

TABLE 4

WMAControlBefore adjusting
After adjusting
Cohen′sd
t valuep valuet valuep value
DLG20.250 ± 0.1040.190 ± 0.0832.8040.006–2.3480.0220.638
NRG10.220 ± 0.0710.210 ± 0.0661.2070.231–1.0020.3200.146
PLA2G40.050 ± 0.0280.020 ± 0.0244.2600.000–3.6720.0001.150
STX8
NEDD4L0.140 ± 0.0410.140 ± 0.053–0.5320.5960.5280.599–0.101
PDE4D0.140 ± 0.0380.120 ± 0.0502.0070.048–1.6650.1000.450
PDE4B0.170 ± 0.0390.150 ± 0.0402.1910.032–1.7220.0900.506
CD550.110 ± 0.0320.100 ± 0.0360.9120.365–0.6090.5450.294
EPHB20.280 ± 0.1370.190 ± 0.0653.7310.000–3.3240.0010.839
RAPGEF5
SERPINA50.220 ± 0.0760.190 ± 0.0521.9040.061–1.6390.1060.461
PDE6C

Expression of candidate genes.

For the DLG2, PLA2G4, PDE4D, PDE4B, and EPHB2 genes, the expression in the WMA group was 1.473, 4.226, 1.334, 1.153, and 1.595 times higher than that in control group, respectively (Figure 3).

FIGURE 3

Regulatory Enzyme Expression

ELISA was used to detect the regulatory enzyme expression. As Table 5 shows, no significant difference was found between the two groups (p > 0.05).

TABLE 5

WMAControlt valuep valueCohen’s d
ASNS2.28 ± 1.652.21 ± 1.440.180.860.045
AGT3.31 ± 1.233.52 ± 1.64–0.630.53–0.145
GAD10.49 ± 0.130.50 ± 0.12–0.340.74–0.080
GAD23.88 ± 0.923.71 ± 1.560.590.540.133

Expression of candidate enzymes.

Discussion

This study found that peripheral metabolites of WMA patients were mainly changed in amino acid metabolism. The most obvious pathways were; (1) the arginine synthesis pathway, (2) the alanine, aspartic acid and glutamate metabolic pathways, (3) the cysteine and methionine metabolic pathways, and (4) the purine metabolic pathway.

Previous studies on methamphetamine addicted animals found that abnormal energy metabolism was the biochemical basis of toxicity and addiction (; ). An increase in energy metabolism was found in the serum of amphetamine-administered rats. After 5 consecutive days of administration, the difference between the administration and control group decreased (). After high dose, high frequency, or long time administration, the energy reserve may be depleted (). The patients in this study were basically deprived from METH for about 2 months. This may explain why abnormal energy metabolism was not observed—the energy reserves may have been consumed.

In addition, amphetamines cause serum/plasma neurotransmitter disturbances mainly in the glutamate-glutamine-GABA system (). Glutamate-Glutamine disturbance was observed both in withdrawal animals and in the WMA patients (; ; Su et al., 2018), and this disturbance may be related to a continuous search for drugs and the high relapse risk (Zheng et al., 2014). Glutamate, aspartic acid, and methionine are elevated in rat neuroblastoma cells 48 h after methamphetamine treatment (). These studies showed that the effect of methamphetamine on the glutamate system could last for a long time. In this study, the “alanine, aspartic acid, and glutamate pathway” was found to be abnormal, and the expression of related metabolites was significantly changed, which is consistent with previous studies.

We chose the glutamate metabolic pathway and used bioinformatics methods to establish a regulatory network model, which predicts that GAD1, GAD2, ASNS, and AGT may directly regulate the glutamate pathway metabolites- pyruvate, L-aspartate, L-glutamine. The regulatory enzyme is regulated by DLG2, NRG1, PLA2G4E, STX8, NEDD4L, PRKG1, PDE4D, PDE4B, SERPINA5, PDE6C, CD55, EPHB2, and RAPGEF5 genes only through one intermediate protein. We verified the expression of the candidate genes and enzymes, and finally found five regulatory genes that were differentially expressed in the two groups, namely DLG2, PLA2G4, PDE4D, PDE4B, and EPHB2 genes. After adjusting for age and BMI, only DLG2, PLA2G4, and EPHB2 remained significant.

DLG2 plays a very important role in synaptic development and stability of plasticity (; Su et al., 2020). The function of the DLG2 gene is associated with complex cognitive and learning tasks (), and is also associated with developmental and intellectual disorders () such as Parkinson’s disease (). The DLG2 gene polymorphism was relevant for methamphetamine addiction (). This study also found that the expression of the DLG2 gene was up-regulated in WMA subjects. There was a long-term cognitive impairment in methamphetamine addicted populations (), as well as dysfunction of the glutamate system (). The DLG2 gene may play an important role in the development, plasticity, and stability of glutamatergic synapses and may further affect the cognitive function of WMA patients.

The protein coded by the PLA2G4E gene is involved in the regulation of membrane-mediated transport and affects phospholipase activity (Zhu et al., 2016). PLA2G4E is strongly expressed in neurons and regulates endogenous cannabinoids, by mobilizing intracellular calcium (Zhu et al., 2016). The endogenous cannabinoid system also regulates the function of the glutamate system. Cannabinoid 1 receptors affect glutamatergic and GABAergic regulation of anxiety and fear responses (Zheng et al., 2011; ). The PLA2G4E gene expression changed the most in this study, reaching 4.226 times that of the control group. That may influence cannabinoid, glutamate, and GABA system functions.

The PDE4D and PDE4B genes belong to the phosphodiesterase (PDE family) superfamily, and play an important role in the degradation of cAMP (). Since cAMP is an important second messenger, PDE4D/4B may play a role in signal transduction by regulating cAMP concentration. PDE4D/4B is widely expressed in the human brain. Previous studies have shown that PDE4D/4B is associated with neuropsychiatric disorders such as depression, schizophrenia (Wu et al., 2018), methamphetamine addiction (), alcohol-induced neuroinflammation (), and stroke (Zhong et al., 2016). Since signal transduction of neurotransmitters in the brain, such as dopamine, glutamate, and GABA, etc., requires the participation of second messenger cAMP, some researchers believe that the phosphodiesterase family (PDE) may be used as an important drug target for schizophrenia and affective disorder (Wu et al., 2018). This study found that the expression of PDE4D/4B in the peripheral blood of WMA patients was significantly higher than that of the control group. However, after adjusting for age and BMI, this difference was no longer significant. The expression ratios of the WMA group compared to the control group were also close to 1, showing that these two genes may not be important regulatory genes in the glutamine system of WMA.

The EPHB2 gene encodes Ephrin type-B receptor 2, which is a tyrosine kinase receptor involved in triggering and maturation. In mature brains, EPHB2 is highly abundant in large dendritic axes and in the frontal cortex and dendritic spines in the hippocampus, and is also involved in the regulation of glutamate receptor distribution and excitatory synapse formation (Su et al., 2017). EPHB2 also regulates synaptic plasticity and leads to anxiety behavior () and reversal of cognitive deficits (). EPHB2 reverses cognitive dysfunction in the Alzheimer’s disease model, possibly related with glutamate receptor-NMDA’s function (N-methyl-D-aspartic acid receptor), and EPHB2 is a key regulator of synaptic localization of NMDA receptors (). Knockout EPHB2 induces depression-like behavior and memory impairment in mice, which is mediated by NMDA receptor subtypes (Zhen et al., 2018). This study found a significant increase of the EPHB2 gene in WMA patients, which elicits that an abnormal expression of the EPHB2 gene may be a potential mechanism of glutamate system dysfunction and emotional or cognitive symptoms of methamphetamine addicts.

The expression of candidate regulatory enzymes was not differentially expressed between two groups. However, the activity of enzymes was not detected, which may influence the metabolites. Another complex regulatory mechanism could not be excluded.

Conclusion

MA use has a long-term and extensive impact on the human body, which can be reflected by the peripheral metabolites. The glutamate system metabolic disturbance may be regulated by the genes that are associated with multiple pathophysiological processes such as synaptic plasticity, signal transduction, cAMP degradation, phosphorylation, and synaptic localization, which may be related to clinical symptoms of withdrawal MA, such as long-term cognitive impairment, anxiety and fear responses, and so on.

There were some limitations in this study. First, only one time point was selected for the exploration of serum metabolomics in the addiction state of methamphetamine. If long-term follow-up can be performed during the acute phase and withdrawal period, the changes in metabolomics characteristics over time will be observed and the results will be more systematic. Second, the exploration of the mechanism is limited to biochemical and gene regulation mechanisms, and if it can be combined with the analysis of central neurotransmitter changes, it will be more in-depth. Third, the enzyme activity was not detected. Fourth, the sample size was limited, and all the enrolled participants were from the Shanghai area only. The conclusion can therefore only be extended to Shanghai Han nationality. A multi-center collaboration may therefore facilitate further discoveries.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors.

Ethics statement

Informed consent of the participants was obtained after the nature of the procedures used in the study had been fully explained. The study was carried out in accordance with the last version of the Declaration of Helsinki, and the protocol was approved by the Ethics Committee of Shanghai Mental Health Center (SMHC IRB 2016-11R). All participants were identified using numbers, instead of names.

Author contributions

SP: conceptualization, data curation, formal analysis, methodology, writing – original draft, review, and editing. HS: investigation, project administration, and supervision. TC: investigation and project administration. XL: investigation and project administration. JD: conceptualization, methodology, project administration, resources, and supervision. HJ: conceptualization, methodology, project administration, resources, and supervision. MZ: conceptualization, funding acquisition, methodology, resources, supervision, and writing – review and editing. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the National Natural Science Foundation (81771436), Shanghai Key Laboratory of Psychotic Disorders (13DZ2260500), Shanghai Municipal Science and Technology Major Project (2018SHZDZX05), Shanghai Clinical Research Center for Mental Health (19MC1911100), Shanghai Shenkang Hospital Development Center (SHDC2020CR3045B), and Shanghai Intelligent Engineering Technology Research Center for Addiction and rehabilitation (19DZ2255200).

Acknowledgments

We appreciate the support from Dr. Xingxing Jian, Shanghai Bioinformation Technology Research Center for the bioinformation analysis and Dr. Qingzhong Wang, Chinese Academy of Sciences for gene expression.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  • 1

    AshokA. H.MizunoY.VolkowN. D.HowesO. D. (2017). Association of stimulant use with dopaminergic alterations in users of cocaine, amphetamine, or methamphetamine a systematic review and meta-analysis.JAMA Psychiatry74511519. 10.1001/jamapsychiatry.2017.0135

  • 2

    BortolatoM.FrauR.PirasA. P.LuesuW.BiniV.DiazG.et al (2009). Methamphetamine induces long-term alterations in reactivity to environmental stimuli: correlation with dopaminergic and serotonergic toxicity.Neurotox. Res.15232245. 10.1007/s12640-009-9024-2

  • 3

    BuQ.LvL.YanG.DengP.WangY.ZhouJ.et al (2013). NMR-based metabonomic in hippocampus, nucleus accumbens and prefrontal cortex of methamphetamine-sensitized rats.Neurotoxicology361723. 10.1016/j.neuro.2013.02.007

  • 4

    CapestranoM.MariggioS.PerinettiG.EgorovaA. V.IacobacciS.SantoroM.et al (2014). Cytosolic phospholipase A2e drives recycling through the clathrin-independent endocytic route.J. Cell Sci.127977993. 10.1242/jcs.136598

  • 5

    ChenY. C.ChenC. K.WangL. J. (2015). Predictors of relapse and dropout during a 12-Week relapse prevention program for methamphetamine users.J. Psychoactive Drugs47317324. 10.1080/02791072.2015.1071447

  • 6

    ChristopherL. N.MichaelM.JosephB. W.KerstinA. F.WolfM. E. (2007). Behavioral sensitization to amphetamine is not accompanied by changes in glutamate receptor surface expression in the rat nucleus accumbens Christopher.J. Neurochem.86573579.

  • 7

    CourtneyK. E.RayL. A. (2014). Methamphetamine: an update on epidemiology, pharmacology, clinical phenomenology, and treatment literature.Drug Alcohol Depend.11121. 10.1016/j.drugalcdep.2014.08.003

  • 8

    CrockerC. E.BernierD. C.HanstockC. C.LakustaB.PurdonS. E.SeresP.et al (2014). Prefrontal glutamate levels differentiate early phase schizophrenia and methamphetamine addiction: a 1 H MRS study at 3 tesla.Schizophr. Res.157231237.

  • 9

    CservenkaA.RayL. A. (2017). Self-reported attentional and motor impulsivity are related to age at first methamphetamine use.Addict. Behav.65712. 10.1016/j.addbeh.2016.09.008

  • 10

    DarnaM.BeckmannJ. S.GipsonC. D.BardoM. T.DwoskinL. P. (2015). Effect of environmental enrichment on dopamine and serotonin transporters and glutamate neurotransmission in medial prefrontal and orbitofrontal cortex.Brain Res.1599115125. 10.1016/j.brainres.2014.12.034

  • 11

    GanH.ZhaoY.JiangH.ZhuY.ChenT.TanH.et al (2018). A research of methamphetamine induced psychosis in 1,430 individuals with methamphetamine use disorder: clinical features and possible risk factors.Front. Psychiatry9:551. 10.3389/fpsyt.2018.00551

  • 12

    HeikkinenN.KarkkainenO.LaukkanenE.KekkonenV.KaarreO.KivimakiP.et al (2018). Changes in the serum metabolite profile correlate with decreased brain gray matter volume in moderate-to-heavy drinking young adults.Alcohol758997. 10.1016/j.alcohol.2018.05.010

  • 13

    HuR.WeiP.JinL.ZhengT.ChenW.-Y.LiuX.-Y.et al (2017). Overexpression of EphB2 in hippocampus rescues impaired NMDA receptors trafficking and cognitive dysfunction in Alzheimer model.Cell Death Dis.8:e2717. 10.1038/cddis.2017.140

  • 14

    IkedaM.OkahisaY.AleksicB.WonM.KondoN.NaruseN.et al (2013). Evidence for shared genetic risk between methamphetamine-induced psychosis and schizophrenia.Neuropsychopharmacology3818641870. 10.1038/npp.2013.94

  • 15

    IrwinC.van ReenenM.MasonS.MienieL. J.WeversR. A.WesterhuisJ. A.et al (2018). The 1H-NMR-based metabolite profile of acute alcohol consumption: a metabolomics intervention study.PLoS One13:e0196850. 10.1371/journal.pone.0196850

  • 16

    KohnoM.OkitaK.MoralesA. M.RobertsonC. L.DeanA. C.GhahremaniD. G.et al (2016). Midbrain functional connectivity and ventral striatal dopamine D2-type receptors: link to impulsivity in methamphetamine users.Mol. Psychiatry2115541560. 10.1038/mp.2015.223

  • 17

    LiH.BuQ.ChenB.ShaoX.HuZ.DengP.et al (2014). Mechanisms of metabonomic for a gateway drug: nicotine priming enhances behavioral response to cocaine with modification in energy metabolism and neurotransmitter level.PLoS One9:e87040. 10.1371/journal.pone.0087040

  • 18

    MakerG. L.GreenT.MullaneyI.TrengoveR. D. (2018). Untargeted metabolomic analysis of rat neuroblastoma cells as a model system to study the biochemical effects of the acute administration of methamphetamine.Metabolites8:38. 10.3390/metabo8020038

  • 19

    McClayJ. L.AdkinsD. E.VunckS. A.BatmanA. M.VannR. E.ClarkS. L.et al (2013). Large-scale neurochemical metabolomics analysis identifies multiple compounds associated with methamphetamine exposure.Metabolomics9392402. 10.1007/s11306-012-0456-y

  • 20

    MellsopG.TapsellR.HolmesP. (2019). Mental health service users’ progression from illicit drug use to schizophrenia in New Zealand.Gen. Psychiatr.32:e100088. 10.1136/gpsych-2019-100088

  • 21

    MoallemN. R.CourtneyK. E.RayL. A. (2018). The relationship between impulsivity and methamphetamine use severity in a community sample.Drug Alcohol Depend.18717. 10.1016/j.drugalcdep.2018.01.034

  • 22

    MorrisC.O’GradaC.RyanM.RocheH. M.GibneyM. J.GibneyE. R.et al (2012). The relationship between BMI and metabolomic profiles: a focus on amino acids.Proc. Nutr. Soc.71634638. 10.1017/s0029665112000699

  • 23

    NithianantharajahJ.KomiyamaN. H.McKechanieA.JohnstoneM.BlackwoodD. H.St ClairD.et al (2013). Synaptic scaffold evolution generated components of vertebrate cognitive complexity.Nat. Neurosci.161624. 10.1038/nn.3276

  • 24

    ParsegianA.SeeR. E. (2014). Dysregulation of dopamine and glutamate release in the prefrontal cortex and nucleus accumbens following methamphetamine self-administration and during reinstatement in rats.Neuropsychopharmacology39811822. 10.1038/npp.2013.231

  • 25

    ProitsiP.KimM.WhileyL.SimmonsA.SattleckerM.VelayudhanL.et al (2017). Association of blood lipids with Alzheimer’s disease: a comprehensive lipidomics analysis.Alzheimers Dement.13140151.

  • 26

    ReggianiC.CoppensS.SekharaT.DimovI.PichonB.LufinN.et al (2017). Novel promoters and coding first exons in DLG2 linked to developmental disorders and intellectual disability.Genome Med.9:67.

  • 27

    RistM. J.RothA.FrommherzL.WeinertC. H.KrugerR.MerzB.et al (2017). Metabolite patterns predicting sex and age in participants of the Karlsruhe Metabolomics and Nutrition (KarMeN) study.PLoS One12:e0183228. 10.1371/journal.pone.0183228

  • 28

    Rodríguez CerdeiraC.Sánchez-BlancoE.Sánchez-BlancoB.González-CespónJ. L. (2017). Protein biomarkers of mood disorders.Int. J. Immunopathol. Pharmacol.30712.

  • 29

    RuehleS.ReyA. A.RemmersF.LutzB. (2012). The endocannabinoid system in anxiety, fear memory and habituation.J. Psychopharmacol.262339. 10.1177/0269881111408958

  • 30

    ShimaN.MiyawakiI.BandoK.HorieH.ZaitsuK.KatagiM.et al (2011). Influences of methamphetamine-induced acute intoxication on urinary and plasma metabolic profiles in the rat.Toxicology2872937. 10.1016/j.tox.2011.05.012

  • 31

    SuH.ChenT.JiangH.ZhongN.DuJ.XiaoK.et al (2020). Intermittent theta burst transcranial magnetic stimulation for methamphetamine addiction: a randomized clinical trial.Eur. Neuropsychopharmacol.31158161. 10.1016/j.euroneuro.2019.12.114

  • 32

    SuH.ZhongN.GanH.WangJ.HanH.ChenT.et al (2017). High frequency repetitive transcranial magnetic stimulation of the left dorsolateral prefrontal cortex for methamphetamine use disorders: a randomised clinical trial.Drug Alcohol Depend.1758491. 10.1016/j.drugalcdep.2017.01.037

  • 33

    SuM. F.LiuM. X.LiJ. Q.LappinJ. M.LiS. X.WuP.et al (2018). Epidemiological characteristics and risk factors of methamphetamine-associated psychotic symptoms.Front. psychiatry9:489. 10.3389/fpsyt.2018.00489

  • 34

    TasicL.PontesJ. G. M.CarvalhoM. S.CruzG.Dal MasC.SethiS.et al (2017). Metabolomics and lipidomics analyses by1H nuclear magnetic resonance of schizophrenia patient serum reveal potential peripheral biomarkers for diagnosis.Schizophr. Res.185182189. 10.1016/j.schres.2016.12.024

  • 35

    WuH. C.ChenC. M.ChenY. C.FungH. C.ChangK. H.WuY. R. (2018). DLG2, but not TMEM229B, GPNMB, and ITGA8 polymorphism, is associated with Parkinson’s disease in a Taiwanese population.Neurobiol. Aging64158.e1158.e6.

  • 36

    YangW.YangR.LuoJ.HeL.LiuJ.ZhangJ. (2018). Increased absolute glutamate concentrations and glutamate-to-creatine ratios in patients with methamphetamine use disorders.Front. psychiatry9:368. 10.3389/fpsyt.2018.00368

  • 37

    YangX.WangY.LiQ.ZhongY.ChenL.DuY.et al (2018). The main molecular mechanisms underlying methamphetamine- induced neurotoxicity and implications for pharmacological treatment.Front. Mol. Neurosci.11:186. 10.3389/fnmol.2018.00186

  • 38

    ZhangK.WangQ.JingX.ZhaoY.JiangH.DuJ.et al (2016). miR-181a is a negative regulator of GRIA2 in methamphetamine-use disorder.Sci. Rep.6:35691.

  • 39

    ZhenL.ShaoT.LuriaV.LiG.LiZ.XuY.et al (2018). EphB2 deficiency induces depression–like behaviors and memory impairment: involvement of NMDA 2B receptor dependent signaling.Front. Pharmacol.9:862. 10.3389/fphar.2018.00862

  • 40

    ZhengC.-Y.SeaboldG. K.HorakM.PetraliaR. S. (2011). MAGUKs, synaptic development, and synaptic plasticity.Neuroscientist17493512. 10.1177/1073858410386384

  • 41

    ZhengP.FangZ.XuX.-J.LiuM.-L.DuX.ZhangX.et al (2016). Metabolite signature for diagnosing major depressive disorder in peripheral blood mononuclear cells.J. Affect. Disord.1957581. 10.1016/j.jad.2016.02.008

  • 42

    ZhengT.LiuL.ShiJ.YuX.XiaoW.SunR.et al (2014). The metabolic impact of methamphetamine on the systemic metabolism of rats and potential markers of methamphetamine abuse.Mol. Biosyst.1019681977. 10.1039/c4mb00158c

  • 43

    ZhongN.JiangH.DuJ.ZhaoY.SunH.XuD.et al (2016). The cognitive impairments and psychological wellbeing of methamphetamine dependent patients compared with health controls.Prog. Neuropsychopharmacol. Biol. Psychiatry693137. 10.1016/j.pnpbp.2016.04.005

  • 44

    ZhuJ.ShangY.ZhangM. (2016). Mechanistic basis of MAGUK-organized complexes in synaptic development and signalling.Nat. Rev. Neurosci.17209223. 10.1038/nrn.2016.18

Summary

Keywords

methamphetamine abuse, metabolomics, bioinformatics, regulatory network, glutamate metabolic pathway

Citation

Peng S, Su H, Chen T, Li X, Du J, Jiang H and Zhao M (2021) The Potential Regulatory Network of Glutamate Metabolic Pathway Disturbance in Chinese Han Withdrawal Methamphetamine Abusers. Front. Genet. 12:653443. doi: 10.3389/fgene.2021.653443

Received

14 January 2021

Accepted

01 March 2021

Published

23 March 2021

Volume

12 - 2021

Edited by

Cunyou Zhao, Southern Medical University, China

Reviewed by

Jinsong Tang, Zhejiang University, China; Qiang Wang, Sichuan University, China

Updates

Copyright

*Correspondence: Min Zhao,

This article was submitted to Behavioral and Psychiatric Genetics, a section of the journal Frontiers in Genetics

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics