Skip to main content
Canna~Fangled Abstracts

Identification of key genes and therapeutic drugs for cocaine addiction using integrated bioinformatics analysis

By July 4, 2023July 20th, 2023No Comments


Abstract

Introduction

Cocaine is a highly addictive drug that is abused due to its excitatory effect on the central nervous system. It is critical to reveal the mechanisms of cocaine addiction and identify key genes that play an important role in addiction.

Methods

In this study, we proposed a centrality algorithm integration strategy to identify key genes in a protein–protein interaction (PPI) network constructed by deferential genes from cocaine addiction-related datasets. In order to investigate potential therapeutic drugs for cocaine addiction, a network of targeted relationships between nervous system drugs and key genes was established.

Results

Four key genes (JUN, FOS, EGR1, and IL6) were identified and well validated using CTD database correlation analysis, text mining, independent dataset analysis, and enrichment analysis methods, and they might serve as biomarkers of cocaine addiction. A total of seventeen drugs have been identified from the network of targeted relationships between nervous system drugs and key genes, of which five (disulfiram, cannabidiol, dextroamphetamine, diazepam, and melatonin) have been shown in the literature to play a role in the treatment of cocaine addiction.

Discussion

This study identified key genes and potential therapeutic drugs for cocaine addiction, which provided new ideas for the research of the mechanism of cocaine addiction.

Keywords: cocaine addiction, key genes, PPI network, biomarker, centrality algorithm

Introduction

Drug addiction is a chronic, recurrent disorder caused by the long-term effects of drugs on the brain (Leshner, ). Since 1985, cocaine, a highly addictive drug that has been abused due to its excitatory effects on the central nervous system, has become one of the world’s leading drugs, mostly in the Americas and Europe. According to the 2020 National Substance Use and Health Survey report released by the Substance Abuse and Mental Health Services Administration (SAMHSA), 1.9% of people 12 years of age or older in 2020 reported cocaine use in the past 12 months (NIDA, ). Cocaine abuse remains a major worldwide health problem (Richards and Le, ).

Numerous studies have shown that cocaine causes irreversible structural changes in organs such as the brain and heart (Riezzo et al., ; Dang et al., ). Research by Goertz et al. () found that cocaine increases dopaminergic neurons and motor activity through midbrain α1 adrenergic signaling. It is well known that the ventral tegmental area (VTA) is an area of the midbrain. In previous studies, the VTA was found to be associated with the addictive properties of many drugs, including cocaine (Cameron and Williams, ). Cocaine abuse results in significant adaptation of dopamine (DA) neurons in the VTA of the midbrain (Wolf et al., ; Stuber et al., ; Mameli and Lüscher, ). Therefore, studies based on the midbrain region could reveal the mechanisms of cocaine addiction.

A differential gene expression analysis is commonly used for the analysis of transcriptomic datasets to explore the underlying molecular mechanisms (Liu et al., ). The construction of the differential gene interaction network according to differential genes has become the main method for data analysis from the system level. Generally, centrality algorithms are mainly used to identify the role of specific nodes in a network and their impact on the network, and nodes with a high centrality ranking may affect other nodes and play an important role in the network. Using a variety of centrality algorithms to analyze the network, screening the most important key genes has become the main analysis method (Chaudhary et al., ; Ma et al., ; Bhattacharyya et al., ; Luan et al., ). In the study of Zhang et al. (), ten different centrality algorithms in cytoHubba were used to identify key genes in the protein–protein interaction (PPI) network, and it was finally verified that the key genes were potential biomarkers or therapeutic targets for opioid addiction. In Poisel et al. ()’s computational biology analysis of human postmortem brain tissues with cocaine addiction, a gene ontology (GO) enrichment analysis was carried out for addiction-related CpG sites. A PPI network analysis revealed several addiction-related genes as highly connected nodes, including CACNA1C, NR3C1, and JUN. Therefore, by identifying key genes in the network, the mechanisms of the addiction process were explored in depth at the system level to explain addiction.

To date, there are no FDA-approved drug treatments for cocaine addiction (Feng et al., ; Shang et al., ), so it is necessary to explore drugs to reduce the incidence and severity of cocaine abuse. In this study, we analyzed datasets related to cocaine addiction, constructed a cocaine addiction-related PPI network to identify potential biomarkers of cocaine addiction, and finally explored potential therapeutic drugs. This could provide new ideas for studying the mechanisms of cocaine addiction and potential cocaine addiction therapeutic drugs.

Materials and methods

The procedure of our study is shown in Figure 1, and the details are described in the following sections.

An external file that holds a picture, illustration, etc.
Object name is fnins-17-1201897-g0001.jpg

Workflow of our methodology. (A) Data. (B) Cocaine addiction-related PPI network. (C) Key gene screening. (D) Validation of key genes and identification of potential therapeutic drugs.

Data

Cocaine addiction-related data GSE54839 (Homo sapiens), GSE67281 (Homo sapiens), GSE186981 (Mus musculus), and GSE155313 (Mus musculus) were downloaded from the Gene Expression Omnibus (GEO) database (http://www.ncbi.nlm.nih.gov/geo). The dataset GSE54839 was chosen as our experimental set, which is based on the GPL6947 platform (Illumina HumanHT-12 V3.0 expression beadchip). This microarray-based study determined the profiles of midbrain gene expression in chronic cocaine abusers (n = 10) and well-matched drug-free control subjects (n = 10). Array-related procedures were performed in triplicate for each subject.

GSE67281 (Homo sapiens), GSE186981 (Mus musculus), and GSE155313 (Mus musculus) were chosen as our validation sets. GSE67281 is an expression profile in postmortem human midbrain specimens from chronic cocaine abusers (n = 11) and well-matched control subjects (n = 11). The GSE186981 is RNA-Seq data in hybrid mouse diversity panel (HMDP) mouse strains of nucleus accumbens (NAc) and prefrontal cortex brain regions (PFC). GSE155313 is the RNA-Seq data from the VTA region of mouse that underwent one of four commonly used paradigms: acute home cage injections of cocaine, chronic home cage injections of cocaine, cocaine-conditioning, or intravenous-self administration of cocaine.

Human PPI data were downloaded from the STRING database (https://string-db.org/) (Szklarczyk et al., ). With a combined score of >900 as the threshold, a total of 230,524 interactions between 11,763 genes were obtained.

Cocaine addiction-related PPI network

The dataset GSE54839 was differentially analyzed using the R package “limma” to obtain their differential genes. The p-value of <0.05 and |log2F>p20mm| > 0.2630344 (i.e., fold change ≥ 1.2 or fold change ≤ 0.8) were considered statistically significant.

To obtain the interaction relationships between DEGs, the downloaded PPI data were filtered using the DEGs obtained from the microarray data GSE54839. By using the gene interactions as edges and DEGs as nodes, a differential gene network was constructed. After removing scatters from the network, the core network was defined as a PPI network related to cocaine addiction.

Key gene identification

We proposed a centrality algorithm integration strategy to analyze genes in cocaine addiction-related PPI networks. The scores of the node in network under each centrality algorithm were calculated separately by applying a series of centrality measures, including degree, edge-percolated component (EPC), Laplacian centrality, maximum neighborhood component (MNC), Katz radiality, and semi-local centrality (SLC). The intersection of the top 10 genes of each centrality algorithm was considered the key gene.

In this study, G is the cocaine addiction-related PPI network we built, and V(G) is the collection of nodes in the network. For node x in G, N(x) is the set of direct neighbors of x in G. For collection A, |A| is used to represent the number of elements in the collection. The specific algorithms are as follows:

  1. Degree (Deg)

    ()=|()|.
  2. Edge percolated component (EPC)

    ()=1|()|=11000.

    Given a threshold of 0.5, 1,000 reduced networks were created by assigning each edge a random number between 0 and 1 and removing edges with associated random numbers less than the threshold. Let the Gk be the reduced network generated at the kth reduced process. If nodes x and y are connected in Gk, set  to 1; otherwise, =0 (Chin et al., ).

  3. Laplacian centrality (Qi et al., ):

    ()=()2+()+2()().
    (1)
  4. Maximum neighborhood component (MNC):

    ()=|((())|,
    (2)

    where m(N(x)) is a maximum connected component of the induced subgraph of G by N(x) (Lin et al., ).

  5. Katz centrality:

    ()==0=1|()|(),
    (3)

    where A is the adjacency matrix of the network G with eigenvalues λ, () is the number of paths from x to y with length k, α is a damping factor and 0<<1max. In all our experiments, we chose α = 0.1 (Wei et al., ).

  6. Radiality (rad):

    ()=()+1(,)|()|1,
    (4)

    where d is the diameter of the network G (Valente and Foreman, ).

  7. Semi-local centrality (SLC):

    ()=()()(),
    (5)

    where B(z) is the number of direct connections and two-step neighbors for node z (Chen et al., ).

Correlation analysis of key genes with cocaine addiction

The Human Protein Atlas (HPA; https://www.proteinatlas.org/) database creates a brain-centric knowledge resource on RNA and protein expression in three mammalian brains: human, pig, and mouse (Sjöstedt et al., ). The RNA expression of key genes in different brain regions in humans and mice was searched in the brain section of the HPA database.

The Comparative Toxicogenomics Database (CTD, https://ctdbase.org/) was used to obtain associations between key genes and cocaine addiction. In the CTD database, the inference score reflects the degree of similarity between the CTD chemical–gene–disease network and a similar scale-free random network (Davis et al., ). The higher the score, the higher the degree of association between the disease and the gene.

The role of key genes in the mechanism of cocaine addiction was identified by text mining in the PubMed database. The search keywords were “cocaine addiction” and the four key genes.

Three independent sets GSE67281 (human), GSE186981 (mouse), and GSE155313 (mouse) were used for pre-addiction and post-addiction differential expression analyses to verify key genes, and the threshold and differential analysis methods were consistent with those of the experimental set GSE54839.

The R package “homologene” was used to search for homologous genes between the human and the mouse. The “homologene” package is a package based on the NCBI HomoloGene (https://www.ncbi.nlm.nih.gov/homologene/) database. The HomoloGene database is a system that can automatically detect congeners in human and mouse genes (NCBI Resource Coordinators, ).

To investigate the possible molecular mechanisms of key genes for cocaine addiction, we used the Kyoto Encyclopedia of Genes and Genomes (KEGG; http://www.kegg.jp/ or http://www.genome.jp/kegg/) database for enrichment analysis. An adjusted p-value of <0.05 was considered to be statistically significant.

Potential therapeutic drug identification

In order to identify potential therapeutic drugs for cocaine addiction, a network of targeted relationships between nervous system drugs and key genes was constructed. Nervous system drug information was retrieved from the ATC classification system of the Drugbank (www.drugbank.ca) database (Wishart et al., ). The targeted effects of key genes with nervous system drugs were reflected in the CTD, where drugs that affect the expression level of genes were our screening criteria. According to the targeting relationship between key genes and nervous system drugs, a targeted relationship network between nervous system drugs and key genes was constructed. Cytoscape was used to visualize this network. Finally, the network was analyzed to screen for potential therapeutic drugs for cocaine addiction.

Results

Cocaine addiction-related PPI network

A total of 724 DEGs were identified from the GSE54839 dataset, including 409 up-regulated genes and 315 down-regulated genes (Supplementary Table 1). The volcano plot was plotted with the “ggplot2” package in R software to visualize the identified DEGs (Figure 2A).

An external file that holds a picture, illustration, etc.
Object name is fnins-17-1201897-g0002.jpg

DEGs and cocaine addiction-related PPI network: (A) A volcano plot of 724 DEGs. Red: upregulated genes; blue: downregulated genes; gray: unchanged genes. (B) Cocaine addiction-related PPI network. The higher the degree value, the larger the node.

The differential genes were mapped to the downloaded PPI data to obtain a differential gene network consisting of 236 nodes and 316 edges. Removing scatter points in the network, the core network had a total of 153 nodes and 263 edges, which was defined as a cocaine addiction-related PPI network. Cytoscape software was used to visualize the network (Figure 2B).

Key gene

The scores of each node in the cocaine addiction related PPI network were calculated separately using seven different centrality algorithms (Figure 3ASupplementary Table 2). Based on our proposed centrality algorithm integration strategy, the top ten genes scored by each algorithm were selected, and their intersections (FOS, IL6, EGR1, and JUN) were regarded as key genes (Figure 3B). The scores of the seven centrality algorithms for the four key genes are shown in Figure 3C.

An external file that holds a picture, illustration, etc.
Object name is fnins-17-1201897-g0003.jpg

Key gene identification: (A) Seven centrality algorithms calculate the distribution of scores for all nodes in the PPI network, and the red signal marks the top ten genes in the score. (B) The top 10 hub genes in the PPI network were identified by seven centrality algorithms and overlapped to obtain four key genes. (C) Seven centrality algorithm results for four key genes. (D) Co-expression analysis heat map of four key genes in samples from drug-free control subjects and chronic cocaine abusers. (E) Co-expression analysis heat map of four key genes in samples from drug-free control subjects and chronic cocaine abusers.

We performed correlation analyses for four key genes in non-drug control participants (Figure 3D) and chronic cocaine abusers (Figure 3E), respectively. The correlation heat map showed that all four key genes were positively correlated, and the correlation showed a significant increase in the cocaine group. It suggested that these four key genes might be more closely related to each other and had synergistic effects during addiction. In addition, the expression of the IL6 gene was the lowest of the four key genes.

Expression of key genes in brain regions

To investigate the expression of four key genes in the brain, we searched the HPA database for the expression of four key genes in different brain regions in the human and mouse (Figures 4A,B).B). The results showed that four key genes were expressed in all regions of the human brain, and IL6 was expressed in the human midbrain region lower than the other three genes, which is consistent with our findings. IL6 was not detected in the mouse midbrain, hypothalamus, pituitary gland, retina, pons, and medulla.

An external file that holds a picture, illustration, etc.
Object name is fnins-17-1201897-g0004.jpg

Brain regions heatmap of four key genes and histograms of expression of different brain regions obtained in the HPA database. (A) Expression of four key genes in human brain regions. (B) Expression of four key genes in mouse brain regions.

Correlation analysis of CTD databases

The correlation scores between each gene in the cocaine addiction-related PPI network and cocaine addiction were searched in CTD. The scores for the top 100 genes are shown in Figure 5. It showed a higher degree of association between the four key genes and cocaine addiction. CTD showed that FOS and EGR1 could be biomarkers of cocaine addiction or play a role in addiction, and EGR1 could be a gene for a therapeutic target in the treatment of cocaine addiction.

An external file that holds a picture, illustration, etc.
Object name is fnins-17-1201897-g0005.jpg

Top 100 correlation scores between all genes in the cocaine addiction-related PPI network and cocaine addiction. Red represents key genes and green represents other genes in the network.

Literature validation

The PubMed database showed that four key genes were all associated with cocaine addiction. Both FOS (Fos proto-oncogene) and JUN (Jun proto-oncogene) are members of the AP-1 transcription factor complex. Multiple studies have shown that cocaine affected the expression of FOS proteins (Todtenkopf et al., ; Imam et al., ; Larson et al., ; Lobo et al., ). Zhang et al. () prepared CPu extracts from D1 and D3 receptor mutant mice and wild-type control littermates at different time points after cocaine injection and found that ERK activation mediates acute cocaine-induced expression of c-fos (Fos). The study by Xu () found that c-fos might mediate cocaine-induced persistent changes by regulating the formation of AP-1 transcriptional complexes and gene expression. Previous studies have demonstrated that cocaine causes increased expression of the JUN protein (Malaplate-Armand et al., ; Paletzki et al., ). Cocaine affects the expression of the JUN protein (Imam et al., ).

There are some studies proving that EGR1 (early growth active protein 1) and c-fos expressions are reduced after cocaine induction (Helton et al., ; Ennulat et al., ). In experiments on mutant mice by Valjent et al. (), EGR1 was found to play a vital role in cocaine-related behavior. Humblot et al. () found that acute cocaine administration was effective in inducing c-FOS and EGR-1 direct early genes, and cocaine-induced EGR-1 and c-FOS expression was significantly reduced in brain regions of rats.

IL6 (interleukin-6) is a pro-inflammatory cytokine. The study by Halpern et al. () showed that men and women respond weakly to pro-inflammatory challenges to IL6 after intravenous cocaine. In experiments measuring changes in IL6 levels in crack cocaine-dependent adolescents after 21 days of withdrawal, it was found that IL6 was elevated in patients on admission compared to the control group (Pianca et al., ).

Independent set analysis

A differential expression analysis was performed on the human dataset GSE67281, which is the expression profile of human cocaine abusers in the midbrain region. A total of 200 DEGs were identified after annotation, including 110 upregulated genes and 90 downregulated genes. There were 20 intersecting genes in the datasets GSE54839 and GSE67281 (Figure 6), including the key genes JUN, FOS, and EGR1, all of which were downregulated in the addictive state (Table 1), while the IL6 gene was not annotated in GSE67281.

An external file that holds a picture, illustration, etc.
Object name is fnins-17-1201897-g0006.jpg

Venn diagram of differential genes in datasets GSE54839 and GSE67281.

Table 1

Differential expression of key genes in GSE54839 and GSE67281.

GSE54839 GSE67281
Gene Brain region logFC logFC
JUN Midbrain −0.842 −0.391
FOS Midbrain −0.978 −0.588
EGR1 Midbrain −1.090 −0.665
IL6 Midbrain −0.557

A differential expression analysis and a homology analysis were performed on the mouse dataset GSE155313 from the VTA brain region, and the Fos gene was identified as a downregulated differential gene under four different conditions. The degree of difference in the Fos gene was not the same between chronic and acute home cage injections of cocaine, and even greater in chronic home cage injections of cocaine condition (Table 2). It suggested that the Fos gene plays a crucial role in long-term addiction.

Table 2

Differential expression of key genes in GSE155313.

GSE155313
Gene logFC Paradigms Brain region
Fos −0.918 Cocaine/saline-conditioning VTA
Fos −0.750 Chronic home cage injections of cocaine/saline VTA
Fos −0.656 Acute home cage injections of cocaine/saline VTA
Fos −0.325 Chronic intravenous-self administration of cocaine/saline VTA

The JUN, FOS, and EGR1 genes were shown to be downregulated differential genes in the human validation set, and the FOS gene was also downregulated in the mouse validation set, which is the same as the experimental set. Analysis of the independent sets showed that the key genes we identified were well-confirmed.

Pathway verification

To reveal the roles of the key genes, we performed a KEGG enrichment analysis of all the genes in the cocaine addiction-related PPI network. A total of 75 KEGG pathways were enriched (Figure 7A). The four key genes were mainly enriched in the TNF signaling pathway, cocaine addiction, amphetamine addiction, IL-17 signaling pathway, MAPK signaling pathway, and Toll-like receptor signaling pathway. Previous studies have shown that these pathways were all linked to cocaine addiction (Northcutt et al., ; Lewitus et al., ; Brown et al., ; Ganguly et al., ; Montesinos et al., ; Bingor et al., ). The subnetwork associated with the key genes and their enrichment pathways (Figure 7B) showed that the four key genes were closely connected in the network, among which JUN, FOS, and IL6 were enriched into multiple pathways, and EGR1 was closely related to these pathways and played a very important synergy.

An external file that holds a picture, illustration, etc.
Object name is fnins-17-1201897-g0007.jpg

Functional enrichment analysis of key genes. (A) Chord diagram of four key genes vs. the top 20 KEGG pathways. (B) Subnetwork associated with key genes and addiction-related KEGG pathways.

The cocaine addiction pathway and the amphetamine addiction pathway are two enriched addiction-related pathways (Figures 8A,B).B). They have similar addiction mechanisms, and both have enhanced firing activity of dopamine neurons in the VTA of the midbrain, resulting in enhanced dopamine release from the NAc.

An external file that holds a picture, illustration, etc.
Object name is fnins-17-1201897-g0008.jpg

KEGG pathway map (Kanehisa and Goto, ; Kanehisa, ; Kanehisa et al., ). (A) Cocaine addiction pathway. (B) Amphetamine addiction pathway. Pink represents differential genes and red represents key genes.

With the stimulation of addictive drugs, the FOS gene induces and maintains an addictive state in the short term of addiction. In human who achieve long-term addiction after further drug use, the FOS gene causes long-term adaptive changes in the brain, and the JUN gene dimerizing with ΔFosB leads to an increased cocaine response. ΔFosB desensitizes c-fos mRNA induction after chronic amphetamine exposure (Renthal et al., ). Zhang et al. () have shown that FOS might mediate cocaine-induced persistent changes by regulating AP-1 transcriptional complexes and target gene expression. To sum up, both our key genes JUN and FOS played important roles in the addiction pathways.

Identification of potential therapeutic drugs for cocaine addiction

To find potential therapeutic drugs for cocaine addiction, a network of targeted relationships between nervous system drugs and key genes was constructed (Figure 9). Fourteen drugs affected four genes, eight drugs affected three genes, nineteen drugs affected two genes, and thirty-five drugs affected one gene. In particular, among drugs that affected four genes, disulfiram, cannabidiol, and dextroamphetamine have been used to mitigate the cocaine response. Many studies have shown that disulfiram might reduce cocaine use in patients with cocaine dependence (Petrakis et al., ; Gaval-Cruz and Weinshenker, ; De Mulder and Dom, ; Kosten et al., ). In the experiment conducted by Petrakis et al. (), disulfiram inhibited dopamine β-hydroxylase, resulting in dopamine overdose and decreased norepinephrine synthesis, possibly weakening cocaine cravings, leading to reduced cocaine use. Dextroamphetamine, a central nervous system stimulant, has been found to be a treatment for cocaine dependence (Grabowski et al., ; Shearer et al., ; Palis et al., ; Ndiaye et al., ). In experiments on rats by Chiodo and Roberts (), sustained dextroamphetamine treatment was found to weaken the potentiating effect of cocaine. Cannabidiol (CBD) is one of the main components of cannabis, and multiple studies have shown that CBD may act as a therapeutic drug for substance abuse (Katsidoni et al., ; Calpe-López et al., ; Anooshe et al., ; Karimi-Haghighi et al., ). Recent research showed that CBD can be effective in reducing the reward and reinforcement effects of addictive drugs (Galaj et al., ). Among the drugs that affected the three genes, diazepam and melatonin might be useful therapeutic agents for reducing cocaine abuse (Takahashi et al., ; Barbosa-Méndez et al., ; Sanchez et al., ).

An external file that holds a picture, illustration, etc.
Object name is fnins-17-1201897-g0009.jpg

A network of targeted relationships between nervous system drugs and key genes.

In twenty-two drugs affecting the expression of three or four key genes, five (disulfiram, dextroamphetamine, diazepam, cannabidiol, and melatonin) have been validated in the literature to reduce cocaine abuse and be used to treat cocaine addiction. They are distributed among drugs used in addictive diseases, psychoanaleptics, and antiepileptic drugs. We, therefore, speculated that seventeen drugs that affected the expression of three or four genes in these three classes might play the same role.

The Drugbank database was used to analyze the status of these seventeen drugs. All of them have been approved by the FDA for the treatment of other diseases, and their effects on cocaine addiction are still being studied. Twelve of these drugs, namely, disulfiram, nicotine, fluoxetine, donepezil, caffeine, amphetamine, cannabidiol, desipramine, valproic acid, dextroamphetamine, carbamazepine, and methylphenidate, are currently in clinical trials for the treatment of cocaine addiction.

The CTD was then used to analyze the relationship between seventeen drugs and cocaine addiction, and the results showed that disulfiram, dextroamphetamine, caffeine, fluoxetine, methylphenidate, desipramine, scopolamine, valproic acid, diazepam, haloperidol, donepezil, clozapine, carbamazepine, and cannabidiol were chemicals with known or potential therapeutic effects in cocaine addiction. Disulfiram, dextroamphetamine, nicotine, fluoxetine, caffeine, methylphenidate, desipramine, diazepam, scopolamine, amphetamine, haloperidol, and melatonin were chemical substances related to cocaine addiction or may play a role in the etiology of cocaine addiction. Although all seventeen drugs have been confirmed in the CTD to be associated with cocaine addiction, further research is needed to determine whether these drugs can be used to treat cocaine addiction.

Discussion

To study the mechanisms of chronic cocaine addiction, data on chronic cocaine abuse in the human midbrain region were used for analysis. Based on the differential expression analysis, a cocaine addiction-related PPI network was constructed, and seven different network centrality algorithms were used to calculate the scores of each gene in the network separately. Finally, four key genes were screened: FOS, IL6, JUN, and EGR1. Through CTD database correlation analysis, literature verification, independent dataset analysis, and enrichment analysis, we found that the four key genes were significantly associated with addiction, and they showed more significant changes under long-term addiction. The network of targeted relationships between nervous system drugs and key genes showed that seventeen drugs targeting three or four key genes were distributed among drugs used in addictive diseases, psychoanaleptics, and antiepileptic drugs, five of which have been shown to be associated with cocaine treatment in the literature. This suggested that key genes might serve as biomarkers for cocaine addiction and that potential therapeutic drugs for cocaine addiction could be found based on key genes.

In this study, the seven centrality algorithms we used were all calculated based on the attributes of the nodes themselves. To show the importance of identifying key genes, we also used these four centrality algorithms to analyze genes in the PPI network associated with cocaine addiction. The algorithm based on the shortest path (closeness, betweenness, EcCentricity, and stress) was not considered. The results of the four unused algorithms are shown in Tables 3,4.4. We obtained four key genes using seven centrality algorithms, of which three to four genes were also included in the results of these four unused algorithms.

Table 3

Top 10 results of three unused algorithms.

Closeness Betweenness Stress
JUN JUN JUN
IL6 IL6 IL6
FOS MAP2K2 MAP2K2
EGR1 VAMP2 VAMP2
CEBPB YWHAH YWHAH
CCL2 FGF2 EGR1
IL1B CEBPB IFITM2
MAP2K2 EGR1 BAG3
JUNB ASNS FGF2
FGF2 GOT1 CEBPB

Table 4

Top 20 results of EcCentricity algorithm.

Gene EcCentricity
GOT1 0.125
MAP2K2 0.125
YWHAH 0.125
IL6 0.125
EDN1 0.125
JUN 0.111
PAK1 0.111
TH 0.111
VAMP2 0.111
DDIT4 0.111
FOS 0.111
FGF2 0.111
DDC 0.111
TIMP1 0.111
CCL2 0.111
CCL20 0.111
CXCL10 0.111
IL1B 0.111
MDH1 0.111
CDK5 0.111

The HPA database showed that four key genes were expressed in multiple brain regions in humans and mice, so data from other brain regions in mice were used for analysis. The validation set GSE186981 was located in the NAc and PFC brain regions, and the difference analysis showed that the Fos and Egr1 genes were downregulated differential genes in both the NAc and PFC brain regions (Table 5). In the validation set of the two sets of mice, the Il6 and Jun genes were not differentially expressed genes. Human addiction to drug abuse is a long-term process, while animal model experiments are usually relatively short and may not fully mimic the process of long-term addiction in humans. In the human addiction pathway map, the FOS gene undergoes changes after acute drug administration, while the JUN gene undergoes changes after long-term addiction. The longest experimental period of the mouse validation set we used is only 7 days, which may not be enough to have formed long-term addiction, so there was no significant difference in the Jun gene. The expression of the IL6 gene was very low in both humans and mice, so we infer that it was too low to reach the difference.

Table 5

Differential expression of key genes in GSE186981.

GSE186981
Gene logFC Brain region
Fos −0.325 NAc
Fos −0.318 PFC
Egr1 −0.270 NAc
Egr1 −0.316 PFC

The genes ranked top in the network based on centrality algorithms were important since they were central in the cocaine addiction-related PPI network, so we analyzed the non-key genes in the top 10 genes of the seven centrality algorithms. The expression of these genes showed significant differences before and after addiction. To reveal their functions, we analyzed these genes using the KEGG database. The results showed that non-key genes were mainly enriched in the cocaine addiction pathway, amphetamine addiction pathway, MAPK signaling pathway, TNF signaling pathway, and synaptic vesicle circulation pathway, all of which were related to cocaine addiction. In the addiction pathway, long-term exposure to addictive drugs can induce a unique transcription factor, delta FosB, which can cause long-term adaptive changes in the brain. Research has confirmed that cocaine can induce the production of TNF (Kovalevich et al., ; Lewitus et al., ; Sil et al., ), thereby affecting non-key genes downstream of this pathway (Figure 10A). Upstream non-key genes (FGF2, VEGFA, and IL1B) affect the expression of downstream genes in this pathway (Figure 10B), while MAP2K2 regulates ERK through phosphorylation, thereby affecting cocaine addiction. Drug addiction is closely related to synapses, and non-key genes (SYT1, VAMP2, SNAP25, and STXBP1) play an important role in the synaptic vesicle circulation pathway (Figure 10C).

An external file that holds a picture, illustration, etc.
Object name is fnins-17-1201897-g0010.jpg

KEGG pathway map of non-key genes (Kanehisa and Goto, ; Kanehisa, ; Kanehisa et al., ). (A) MAPK signaling pathway. (B) TNF signaling pathway. (C) Synaptic vesicle cycle. Red represents key genes and yellow represents non-key genes in the top 10 of the seven centrality algorithms.

To date, there are too few datasets related to cocaine addiction in humans, and the sample size in our study is not very large. Therefore, the potential biomarkers and therapeutic targets of cocaine addiction identified in this study needed further experimental verification.

In summary, this study identified four key genes (FOS, IL6, EGR1, and JUN) that might be involved in cocaine addiction mechanisms and had potential roles as biomarkers and therapeutic targets for cocaine addiction. Our research provided new ideas for the study of the mechanism of cocaine addiction and was expected to help in the treatment of cocaine addiction.

Data availability statement

Publicly available datasets were analyzed in this study. This data can be found here: https://www.ncbi.nlm.nih.gov/geo/GSE54839GSE67281GSE155313, and GSE186981.

Author contributions

WL and LC: conceptualization, project administration, and supervision. XuW, SS, HC, BY, ZZ, XiW, and YW: data curation. XuW: formal analysis and investigation. WL: funding acquisition. XuW, WL, and LC: methodology. SS, HC, BY, ZZ, XiW, YW, JL, and YH: validation. XuW and WL: visualization and writing—original draft. LC: writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding Statement

This research was funded by the Natural Science Foundation of Heilongjiang Province (LH2021F043) and the National Natural Science Foundation of China (61702141).

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.

Publisher’s note

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2023.1201897/full#supplementary-material

Supplementary Table 1

Differential expression analysis results of GSE54839.

Supplementary Table 2

Scores of seven centrality algorithms.

References

  • Anooshe M., Nouri K., Karimi-Haghighi S., Mousavi Z., Haghparast A. (2021). Cannabidiol efficiently suppressed the acquisition and expression of methamphetamine-induced conditioned place preference in the ratBehav. Brain Res. 404, 113158. 10.1016/j.bbr.2021.113158 [PubMed] [CrossRef[]
  • Barbosa-Méndez S., Pérez-Sánchez G., Becerril-Villanueva E., Salazar-Juárez A. (2021). Melatonin decreases cocaine-induced locomotor sensitization and cocaine-conditioned place preference in ratsJ. Psychiatr. Res. 132, 97–110. 10.1016/j.jpsychires.2020.09.027 [PubMed] [CrossRef[]
  • Bhattacharyya N., Gupta S., Sharma S., Soni A., Bagabir S. A., Bhattacharyya M., et al.. (2022). CDK1 and HSP90AA1 appear as the novel regulatory genes in non-small cell lung cancer: a bioinformatics approachJ. Pers. Med. 12, 393. 10.3390/jpm12030393 [PMC free article] [PubMed] [CrossRef[]
  • Bingor A., Azriel M., Amiad L., Yaka R. (2021). Potentiated response of ERK/MAPK signaling is associated with prolonged withdrawal from cocaine behavioral sensitizationJ. Mol. Neurosci. 71, 2229–2236. 10.1007/s12031-021-01799-6 [PMC free article] [PubMed] [CrossRef[]
  • Brown K. T., Levis S. C., O’Neill C. E., Northcutt A. L., Fabisiak T. J., Watkins L. R., et al.. (2018). Innate immune signaling in the ventral tegmental area contributes to drug-primed reinstatement of cocaine seekingBrain Behav. Immun. 67, 130–138. 10.1016/j.bbi.2017.08.012 [PMC free article] [PubMed] [CrossRef[]
  • Calpe-López C., García-Pardo M. P., Aguilar M. A. (2019). Cannabidiol treatment might promote resilience to cocaine and methamphetamine use disorders: a review of possible mechanismsMolecules 24, 2583. 10.3390/molecules24142583 [PMC free article] [PubMed] [CrossRef[]
  • Cameron D. L., Williams J. T. (1994). Cocaine inhibits GABA release in the VTA through endogenous 5-HTJ. Neurosci. 14, 6763–6767. 10.1523/JNEUROSCI.14-11-06763.1994 [PMC free article] [PubMed] [CrossRef[]
  • Chaudhary R., Balhara M., Jangir D. K., Dangi M., Dangi M., Chhillar A. K. (2019). In silico protein interaction network analysis of virulence proteins associated with invasive aspergillosis for drug discoveryCurr. Top. Med. Chem. 19, 146–155. 10.2174/1568026619666181120150633 [PubMed] [CrossRef[]
  • Chen D. B., Lu L. Y., Shang M. S., Zhang Y. C., Zhou T. (2012). Identifying influential nodes in complex networksPhys. Stat. Mech. Appl. 391, 1777–1787. 10.1016/j.physa.2011.09.017 [CrossRef[]
  • Chin C. H., Chen S. H., Wu H. H., Ho C. W., Ko M. T., Lin C. Y. (2014). cytoHubba: identifying hub objects and sub-networks from complex interactomeBMC Syst. Biol. 8 (Suppl. 4), S11. 10.1186/1752-0509-8-S4-S11 [PMC free article] [PubMed] [CrossRef[]
  • Chiodo K. A., Roberts D. C. (2009). Decreased reinforcing effects of cocaine following 2 weeks of continuous D-amphetamine treatment in ratsPsychopharmacology 206, 447–456. 10.1007/s00213-009-1622-4 [PMC free article] [PubMed] [CrossRef[]
  • Dang J., Tao Q., Niu X., Zhang M., Gao X., Yang Z., et al.. (2022). Meta-analysis of structural and functional brain abnormalities in cocaine addictionFront. Psychiatry 13, 927075. 10.3389/fpsyt.2022.927075 [PMC free article] [PubMed] [CrossRef[]
  • Davis A. P., Wiegers T. C., Johnson R. J., Sciaky D., Wiegers J., Mattingly C. J. (2023). Comparative Toxicogenomics Database (CTD): update 2023Nucleic Acids Res. 51, D1257–d1262. 10.1093/nar/gkac833 [PMC free article] [PubMed] [CrossRef[]
  • De Mulder I., Dom G. (2012). Disulfiram as a treatment for cocaine dependencyTijdschr. Psychiatr. 54, 51–58. [PubMed[]
  • Ennulat D. J., Babb S., Cohen B. M. (1994). Persistent reduction of immediate early gene mRNA in rat forebrain following single or multiple doses of cocaineBrain Res. Mol. Brain Res. 26, 106–112. 10.1016/0169-328X(94)90080-9 [PubMed] [CrossRef[]
  • Feng H., Gao K., Chen D., Shen L., Robison A. J., Ellsworth E., et al.. (2022). Machine learning analysis of cocaine addiction informed by DAT, SERT, and NET-based interactome networksJ. Chem. Theory Comput. 18, 2703–2719. 10.1021/acs.jctc.2c00002 [PubMed] [CrossRef[]
  • Galaj E., Bi G. H., Yang H. J., Xi Z. X. (2020). Cannabidiol attenuates the rewarding effects of cocaine in rats by CB2, 5-HT(1A) and TRPV1 receptor mechanismsNeuropharmacology 167, 107740. 10.1016/j.neuropharm.2019.107740 [PMC free article] [PubMed] [CrossRef[]
  • Ganguly P., Honeycutt J. A., Rowe J. R., Demaestri C., Brenhouse H. C. (2019). Effects of early life stress on cocaine conditioning and AMPA receptor composition are sex-specific and driven by TNFBrain Behav. Immun. 78, 41–51. 10.1016/j.bbi.2019.01.006 [PMC free article] [PubMed] [CrossRef[]
  • Gaval-Cruz M., Weinshenker D. (2009). mechanisms of disulfiram-induced cocaine abstinence: antabuse and cocaine relapseMol. Interv. 9, 175–187. 10.1124/mi.9.4.6 [PMC free article] [PubMed] [CrossRef[]
  • Goertz R. B., Wanat M. J., Gomez J. A., Brown Z. J., Phillips P. E., Paladini C. A. (2015). Cocaine increases dopaminergic neuron and motor activity via midbrain α1 adrenergic signalingNeuropsychopharmacology 40, 1151–1162. 10.1038/npp.2014.296 [PMC free article] [PubMed] [CrossRef[]
  • Grabowski J., Rhoades H., Schmitz J., Stotts A., Daruzska L. A., Creson D., et al.. (2001). Dextroamphetamine for cocaine-dependence treatment: a double-blind randomized clinical trialJ. Clin. Psychopharmacol. 21, 522–526. 10.1097/00004714-200110000-00010 [PubMed] [CrossRef[]
  • Halpern J. H., Sholar M. B., Glowacki J., Mello N. K., Mendelson J. H., Siegel A. J. (2003). Diminished interleukin-6 response to proinflammatory challenge in men and women after intravenous cocaine administrationJ. Clin. Endocrinol. Metab. 88, 1188–1193. 10.1210/jc.2002-020804 [PubMed] [CrossRef[]
  • Helton T. E., Daunais J. B., McGinty J. F. (1993). Convulsant doses of cocaine alter immediate early gene and opioid peptide expression in rat limbic forebrainBrain Res. Mol. Brain Res. 20, 285–288. 10.1016/0169-328X(93)90054-S [PubMed] [CrossRef[]
  • Humblot N., Thiriet N., Gobaille S., Aunis D., Zwiller J. (1998). The serotonergic system modulates the cocaine-induced expression of the immediate early genes egr-1 and c-fos in rat brainAnn. N. Y. Acad. Sci. 844, 7–20. 10.1111/j.1749-6632.1998.tb08217.x [PubMed] [CrossRef[]
  • Imam S. Z., Duhart H. M., Skinner J. T., Ali S. F. (2005). Cocaine induces a differential dose-dependent alteration in the expression profile of immediate early genes, transcription factors, and caspases in PC12 cells: a possible mechanism of neurotoxic damage in cocaine addictionAnn. N. Y. Acad. Sci. 1053, 482–490. 10.1196/annals.1344.042 [PubMed] [CrossRef[]
  • Kanehisa M. (2019). Toward understanding the origin and evolution of cellular organismsProtein Sci28, 1947–1951. 10.1002/pro.3715 [PMC free article] [PubMed] [CrossRef[]
  • Kanehisa M., Furumichi M., Sato Y., Kawashima M., Ishiguro-Watanabe M. (2023). KEGG for taxonomy-based analysis of pathways and genomesNucleic Acids Res51, D587–d592. 10.1093/nar/gkac963 [PMC free article] [PubMed] [CrossRef[]
  • Kanehisa M., Goto S. (2000). KEGG: kyoto encyclopedia of genes and genomesNucleic Acids Res28, 27–30. 10.1093/nar/28.1.27 [PMC free article] [PubMed] [CrossRef[]
  • Karimi-Haghighi S., Razavi Y., Iezzi D., Scheyer A. F., Manzoni O., Haghparast A. (2022). Cannabidiol and substance use disorder: dream or realityNeuropharmacology 207, 108948. 10.1016/j.neuropharm.2022.108948 [PMC free article] [PubMed] [CrossRef[]
  • Katsidoni V., Anagnostou I., Panagis G. (2013). Cannabidiol inhibits the reward-facilitating effect of morphine: involvement of 5-HT1A receptors in the dorsal raphe nucleusAddict. Biol. 18, 286–296. 10.1111/j.1369-1600.2012.00483.x [PubMed] [CrossRef[]
  • Kosten T. R., Wu G., Huang W., Harding M. J., Hamon S. C., Lappalainen J., et al.. (2013). Pharmacogenetic randomized trial for cocaine abuse: disulfiram and dopamine β-hydroxylaseBiol. Psychiatry 73, 219–224. 10.1016/j.biopsych.2012.07.011 [PMC free article] [PubMed] [CrossRef[]
  • Kovalevich J., Yen W., Ozdemir A., Langford D. (2015). Cocaine induces nuclear export and degradation of neuronal retinoid X receptor-γ via a TNF-α/JNK- mediated mechanismJ. Neuroimmune Pharmacol. 10, 55–73. 10.1007/s11481-014-9573-x [PMC free article] [PubMed] [CrossRef[]
  • Larson E. B., Akkentli F., Edwards S., Graham D. L., Simmons D. L., Alibhai I. N., et al.. (2010). Striatal regulation of ΔFosB, FosB, and cFos during cocaine self-administration and withdrawalJ. Neurochem. 115, 112–122. 10.1111/j.1471-4159.2010.06907.x [PMC free article] [PubMed] [CrossRef[]
  • Leshner A. I. (1997). Addiction is a brain disease, and it mattersScience 278, 45–47. 10.1126/science.278.5335.45 [PubMed] [CrossRef[]
  • Lewitus G. M., Konefal S. C., Greenhalgh A. D., Pribiag H., Augereau K., Stellwagen D. (2016). Microglial TNF-α suppresses cocaine-induced plasticity and behavioral sensitizationNeuron 90, 483–491. 10.1016/j.neuron.2016.03.030 [PMC free article] [PubMed] [CrossRef[]
  • Lin C. Y., Chin C. H., Wu H. H., Chen S. H., Ho C. W., Ko M. T. (2008). Hubba: hub objects analyzer–a framework of interactome hubs identification for network biologyNucleic Acids Res. 36, W438–443. 10.1093/nar/gkn257 [PMC free article] [PubMed] [CrossRef[]
  • Liu Z., Ru L., Ma Z. (2021). Low expression of ADCY4 predicts worse survival of lung squamous cell carcinoma based on integrated analysis and immunohistochemical verificationFront. Oncol. 11, 637733. 10.3389/fonc.2021.637733 [PMC free article] [PubMed] [CrossRef[]
  • Lobo M. K., Covington H. E., 3rd, Chaudhury D., Friedman A. K., Sun H., Damez-Werno D., et al.. (2010). Cell type-specific loss of BDNF signaling mimics optogenetic control of cocaine rewardScience 330, 385–390. 10.1126/science.1188472 [PMC free article] [PubMed] [CrossRef[]
  • Luan H., Jian L., He Y., Zhang T., Zhou L. (2022). Exploration and validation of metastasis-associated genes for skin cutaneous melanomaSci. Rep. 12, 13002. 10.1038/s41598-022-17468-6 [PMC free article] [PubMed] [CrossRef[]
  • Ma Z., Xu J., Ru L., Zhu W. (2021). Identification of pivotal genes associated with the prognosis of gastric carcinoma through integrated analysisBiosci. Rep. 41, BSR20203676. 10.1042/BSR20203676 [PMC free article] [PubMed] [CrossRef[]
  • Malaplate-Armand C., Becuwe P., Ferrari L., Masson C., Dauça M., Visvikis S., et al.. (2005). Effect of acute and chronic psychostimulant drugs on redox status, AP-1 activation and pro-enkephalin mRNA in the human astrocyte-like U373 MG cellsNeuropharmacology 48, 673–684. 10.1016/j.neuropharm.2004.12.010 [PubMed] [CrossRef[]
  • Mameli M., Lüscher C. (2011). Synaptic plasticity and addiction: learning mechanisms gone awryNeuropharmacology 61, 1052–1059. 10.1016/j.neuropharm.2011.01.036 [PubMed] [CrossRef[]
  • Montesinos J., Castilla-Ortega E., Sánchez-Marín L., Montagud-Romero S., Araos P., Pedraz M., et al.. (2020). Cocaine-induced changes in CXCL1 and inflammatory signaling pathways in the hippocampus: association with IL1βNeuropharmacology 162, 107840. 10.1016/j.neuropharm.2019.107840 [PubMed] [CrossRef[]
  • NCBI Resource Coordinators (2014). Database resources of the National Center for Biotechnology InformationNucleic Acids Res. 42, D7–17. 10.1093/nar/gkt1146 [PMC free article] [PubMed] [CrossRef[]
  • Ndiaye N. A., Allain F., Samaha A. N. (2022). D-amphetamine maintenance therapy reduces cocaine use in female ratsPsychopharmacology 239, 3755–3770. 10.1007/s00213-022-06271-w [PubMed] [CrossRef[]
  • NIDA (2022). What is Cocaine? Available online at: https://nida.nih.gov/publications/research-reports/cocaine/what-cocaine (accessed February 1, 2023).
  • Northcutt A. L., Hutchinson M. R., Wang X., Baratta M. V., Hiranita T., Cochran T. A., et al.. (2015). DAT isn’t all that: cocaine reward and reinforcement require Toll-like receptor 4 signalingMol. Psychiatry 20, 1525–1537. 10.1038/mp.2014.177 [PMC free article] [PubMed] [CrossRef[]
  • Paletzki R. F., Myakishev M. V., Polesskaya O., Orosz A., Hyman S. E., Vinson C. (2008). Inhibiting activator protein-1 activity alters cocaine-induced gene expression and potentiates sensitizationNeuroscience 152, 1040–1053. 10.1016/j.neuroscience.2008.01.045 [PMC free article] [PubMed] [CrossRef[]
  • Palis H., MacDonald S., Jun J., Oviedo-Joekes E. (2021). Use of sustained release dextroamphetamine for the treatment of stimulant use disorder in the setting of injectable opioid agonist treatment in Canada: a case reportHarm Reduct. J. 18, 57. 10.1186/s12954-021-00500-9 [PMC free article] [PubMed] [CrossRef[]
  • Petrakis I. L., Carroll K. M., Nich C., Gordon L. T., McCance-Katz E. F., Frankforter T., et al.. (2000). Disulfiram treatment for cocaine dependence in methadone-maintained opioid addictsAddiction 95, 219–228. 10.1046/j.1360-0443.2000.9522198.x [PubMed] [CrossRef[]
  • Pianca T. G., Rosa R. L., Ceresér K. M. M., de Aguiar B. W., de Abrahão R. C., Lazzari P. M., et al.. (2017). Differences in biomarkers of crack-cocaine adolescent users before/after abstinenceDrug Alcohol Depend. 177, 207–213. 10.1016/j.drugalcdep.2017.03.043 [PubMed] [CrossRef[]
  • Poisel E., Zillich L., Streit F., Frank J., Friske M. M., Foo J. C., et al.. (2023). DNA methylation in cocaine use disorder-an epigenome-wide approach in the human prefrontal cortexFront. Psychiatry 14, 1075250. 10.3389/fpsyt.2023.1075250 [PMC free article] [PubMed] [CrossRef[]
  • Qi X. Q., Fuller E., Wu Q., Wu Y. Z., Zhang C. Q. (2012). Laplacian centrality: a new centrality measure for weighted networksInf. Sci. 194, 240–253. 10.1016/j.ins.2011.12.027 [CrossRef[]
  • Renthal W., Carle T. L., Maze I., Covington H. E., 3rd, Truong H. T., Alibhai I., et al.. (2008). Delta FosB mediates epigenetic desensitization of the c-fos gene after chronic amphetamine exposureJ. Neurosci. 28, 7344–7349. 10.1523/JNEUROSCI.1043-08.2008 [PMC free article] [PubMed] [CrossRef[]
  • Richards J. R., Le J. K. (2022). Cocaine Toxicity. FL: Treasure Island: StatPearls Publishing. [PubMed[]
  • Riezzo I., Fiore C., De Carlo D., Pascale N., Neri M., Turillazzi E., et al.. (2012). Side effects of cocaine abuse: multiorgan toxicity and pathological consequencesCurr. Med. Chem. 19, 5624–5646. 10.2174/092986712803988893 [PubMed] [CrossRef[]
  • Sanchez W. N., Pochapski J. A., Jessen L. F., Ellenberger M., Schwarting R. K., Robinson D. L., et al.. (2022). Diazepam attenuates the effects of cocaine on locomotion, 50-kHz ultrasonic vocalizations and phasic dopamine in the nucleus accumbens of ratsBr. J. Pharmacol. 179, 1565–1577. 10.1111/bph.15658 [PubMed] [CrossRef[]
  • Shang L., Jin Z., Wei H., Park S., Zhan C. G., Zheng F. (2023). Catalytic activities of a highly efficient cocaine hydrolase for hydrolysis of biologically active cocaine metabolites norcocaine and benzoylecgonineSci. Rep. 13, 640. 10.1038/s41598-022-27280-x [PMC free article] [PubMed] [CrossRef[]
  • Shearer J., Wodak A., van Beek I., Mattick R. P., Lewis J. (2003). Pilot randomized double blind placebo-controlled study of dexamphetamine for cocaine dependenceAddiction 98, 1137–1141. 10.1046/j.1360-0443.2003.00447.x [PubMed] [CrossRef[]
  • Sil S., Niu F., Tom E., Liao K., Periyasamy P., Buch S. (2019). Cocaine mediated neuroinflammation: role of dysregulated autophagy in pericytesMol. Neurobiol. 56, 3576–3590. 10.1007/s12035-018-1325-0 [PMC free article] [PubMed] [CrossRef[]
  • Sjöstedt E., Zhong W., Fagerberg L., Karlsson M., Mitsios N., Adori C., et al.. (2020). An atlas of the protein-coding genes in the human, pig, and mouse brainScience 367, eaay5947. 10.1126/science.aay5947 [PubMed] [CrossRef[]
  • Stuber G. D., Hopf F. W., Tye K. M., Chen B. T., Bonci A. (2010). Neuroplastic alterations in the limbic system following cocaine or alcohol exposureCurr. Top. Behav. Neurosci. 3, 3–27. 10.1007/7854_2009_23 [PubMed] [CrossRef[]
  • Szklarczyk D., Gable A. L., Nastou K. C., Lyon D., Kirsch R., Pyysalo S., et al.. (2021). The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement setsNucleic Acids Res. 49, D605–d612. 10.1093/nar/gkab835 [PMC free article] [PubMed] [CrossRef[]
  • Takahashi T. T., Vengeliene V., Spanagel R. (2017). Melatonin reduces motivation for cocaine self-administration and prevents relapse-like behavior in ratsPsychopharmacology 234, 1741–1748. 10.1007/s00213-017-4576-y [PubMed] [CrossRef[]
  • Todtenkopf M. S., Mihalakopoulos A., Stellar J. R. (2002). Withdrawal duration differentially affects c-fos expression in the medial prefrontal cortex and discrete subregions of the nucleus accumbens in cocaine-sensitized ratsNeuroscience 114, 1061–1069. 10.1016/S0306-4522(02)00272-5 [PubMed] [CrossRef[]
  • Valente T. W., Foreman R. K. (1998). Integration and radiality: measuring the extent of an individual’s connectedness and reachability in a networkSoc. Networks 20, 89–105. 10.1016/S0378-8733(97)00007-5 [CrossRef[]
  • Valjent E., Aubier B., Corbill,é A. G., Brami-Cherrier K., Caboche J., Topilko P., et al.. (2006). Plasticity-associated gene Krox24/Zif268 is required for long-lasting behavioral effects of cocaineJ. Neurosci. 26, 4956–4960. 10.1523/JNEUROSCI.4601-05.2006 [PMC free article] [PubMed] [CrossRef[]
  • Wei P. J., Wu F. X., Xia J., Su Y., Wang J., Zheng C. H. (2020). Prioritizing cancer genes based on an improved random walk methodFront. Genet. 11, 377. 10.3389/fgene.2020.00377 [PMC free article] [PubMed] [CrossRef[]
  • Wishart D. S., Feunang Y. D., Guo A. C., Lo E. J., Marcu A., Grant J. R., et al.. (2018). DrugBank 5.0: a major update to the DrugBank database for 2018Nucleic Acids Res. 46, D1074–d1082. 10.1093/nar/gkx1037 [PMC free article] [PubMed] [CrossRef[]
  • Wolf M. E., Sun X., Mangiavacchi S., Chao S. Z. (2004). Psychomotor stimulants and neuronal plasticityNeuropharmacology 47 (Suppl. 1), 61–79. 10.1016/j.neuropharm.2004.07.006 [PubMed] [CrossRef[]
  • Xu M. (2008). c-Fos is an intracellular regulator of cocaine-induced long-term changesAnn. N. Y. Acad. Sci. 1139, 1–9. 10.1196/annals.1432.049 [PMC free article] [PubMed] [CrossRef[]
  • Zhang J., Zhang L., Jiao H., Zhang Q., Zhang D., Lou D., et al.. (2006). c-Fos facilitates the acquisition and extinction of cocaine-induced persistent changesJ. Neurosci. 26, 13287–13296. 10.1523/JNEUROSCI.3795-06.2006 [PMC free article] [PubMed] [CrossRef[]
  • Zhang L., Lou D., Jiao H., Zhang D., Wang X., Xia Y., et al.. (2004). Cocaine-induced intracellular signaling and gene expression are oppositely regulated by the dopamine D1 and D3 receptorsJ. Neurosci. 24, 3344–3354. 10.1523/JNEUROSCI.0060-04.2004 [PMC free article] [PubMed] [CrossRef[]
  • Zhang X., Yu H., Bai R., Ma C. (2020). Identification and characterization of biomarkers and their role in opioid addiction by integrated bioinformatics analysisFront. Neurosci. 14, 608349. 10.3389/fnins.2020.608349 [PMC free article] [PubMed] [CrossRef[]

Articles from Frontiers in Neuroscience are provided here courtesy of Frontiers Media SA

Leave a Reply