Breaking the Chains: Advances in Substance Addiction Research through Single-Cell Sequencing, Epigenetics, and Epitranscriptomic

: Addiction is a complex brain disease influenced by genetic, environmental, and neurological factors. Psychostimulants, cocaine, and methamphetamine influence different cell types in different brain regions, with a focus on the neurons responsible for rewarding effects in the nucleus accumbens (NAc) and ventral tegmental area (VTA). Known markers for psychostimulant-induced neuronal plasticity in combination with droplet-based high-throughput single-cell sequencing divided the heterogeneity of cell populations in NAc and VTA into clusters, where all cells of the same type do not respond equally to exposure to psychostimulants. To explain psychostimulant-induced neuronal plasticity as changes in the amplitude and phase shifts of gene expression, we focused on epigenetic mechanisms of DNA and chromatin modifications, as well as DNA accessibility. We also comment on epitranscriptomics as a novel approach in the study of messenger RNA posttranslational modification, which regulates translation and potentially localized transcription in synapses in order to address the molecular chains that connect addiction from changes in gene expression to synaptic and, finally, neuronal plasticity.


Introduction
Behaviorally, addiction is defined as compulsive drug seeking and use despite negative consequences [1].The underlying molecular mechanisms primarily involve the brain reward system, associated with dopamine release [2][3][4].Acute drug exposure disrupts neurotransmitter signaling, localization, metabolism, and synthesis, resulting in neuroadaptive changes that contribute to the development of addiction [5,6].Prolonged drug exposure induces neuroplastic changes, altering the structure and function of neurons and influencing synaptic connections, thus impacting reward neuronal circuitry and resulting in long-term behavioral changes [7,8].
Drugs of abuse are a diverse group of compounds classified into various categories based on their pharmacological effects and molecular mechanisms of action.We have concentrated on the psychostimulants cocaine (COC) and methamphetamine (METH) because of the severity of the cognitive, behavioral, and physical disorders brought on by the damage to neuronal cells and the metabolic dysfunction in the central nervous system after its consumption [9].COC and METH are central nervous system stimulants that increase dopamine release in brain regions of the mesolimbic circuit, influencing reward and cognitive functions [10].Although COC and METH have similar behavioral and physiological effects, there are some major differences.COC is quickly removed and completely metabolized in the body, while METH has a longer duration of action due to slower metabolism and remains unchanged in the brain longer than COC, leading to prolonged stimulant effects [11].Both COC and METH are highly addictive stimulants that are widely recognized as one of the most abused drugs in the world, which is why this review focuses specifically on these two psychostimulants [12,13].Key brain regions involved in reward processing, influenced by COC and METH, include the nucleus accumbens (NAc) and ventral tegmental area (VTA).Long-term changes induced by psychostimulant exposure in these regions involve alterations in gene expression and cellular physiology [14,15].Psychostimulant-induced synaptic plasticity is linked to the known mechanism of phosphorylation of extracellular signal-regulated kinases (ERK1/2) and cyclic adenosine 3 ′ ,5 ′ -monophosphate (cAMP)-response element-binding protein (CREB), with cAMP-dependent protein kinase (PKA) also playing a role in drug-induced memory formation [16][17][18][19].However, the main challenge in addiction research remains the identification of cell types, gene/s, or different molecular markers capable of distinguishing natural from drug-induced alterations in the reward circuitry.
While some people only use drugs for experimentation and never develop an addiction, others develop an addiction after being exposed to psychoactive substances on a regular basis.Understanding the molecular basis for this phenomenon requires a complex interplay of genetic, environmental, and neurobiological factors.Individual variability and the dynamic nature of drug responses both influence the molecular mechanisms involved in the cascade of events that affect neurotransmitter release, receptor activation, and intracellular signaling pathways.Individual vulnerability and susceptibility to addiction were studied in a mouse model, and it was discovered that in wild-type populations, voluntary oral methamphetamine consumption can undergo bidirectional selective breeding, producing two strings with high and low preference for METH [20].Using a selective breeding approach, candidate genes can be identified in the absence of METH exposure.By using RNA-Seq to selectively breed low-and high-METH-preference mice, it was discovered that the trace amine-associated receptor 1 (Taar1) gene plays an important role [21].Indeed, it was confirmed that increasing TAAR1 function reduces drug self-administration and intake [22][23][24].In addition to genetic influence, non-genetic factors such as social isolation or enrichment have been shown to influence drug preference [25,26].Taken together, the complexity of the risk for addiction development as influenced by personal genetics and environment should be investigated further, because recent studies point to epigenetic changes caused by isolation that affect addiction development.
Previous studies on molecular changes associated with addiction have primarily focused on candidate genes identified by proteomic or genetic techniques using tissue extracts from drug-exposed individuals.With this approach, it is not possible to distinguish in which cells the genetic changes originate, whether they occur at the gene or protein level and whether the observed changes are time-dependent or stable.To explain the complexities of these influences, it is necessary to refine the accepted central dogma, foundational to our understanding of genetic information flow, by applying epigenetics and epitranscriptomics to gene regulation (Figure 1).Environmental factors influence gene expression by activating or deactivating genes via epigenetic modifications such as DNA methylation and histone acetylation [27].Moreover, recent studies show that environmental factors can lead to the modification of different RNAs, influencing RNA stability, splicing, translation, and other post-transcriptional processes that also have an impact on protein synthesis and/or function [28][29][30].
The primary focus of this review is on recent research using droplet-based singlecell sequencing methods and epigenetic and epitranscriptomic approaches in NAc and VTA brain regions after acute and chronic exposure to the psychostimulants COC and METH, regardless of the route of administration.The single-cell approach confirmed known and discovered new cell types involved in the response to psychostimulants, which were clustered based on changes in gene expression.In the field of epigenetics and epitranscriptomics, our attention has centered on the role of molecular players known as "readers", "writers", and "erasers", which contribute to the regulation of gene expression through modifications of DNA, histone proteins, and RNA.We have specifically addressed two crucial aspects of addiction development and treatment: first, we distinguished between the natural process of neuronal plasticity and psychostimulant-induced neuronal plasticity; second, we analyzed the common mechanisms and genes involved in obesity and addiction.Additionally, we discussed the systemic effects that are a consequence of addiction and psychostimulant use, including impaired oxidative regulation, immune response, and changes in the microbiome profile.These findings provide an alternative perspective on the long-term consequences of withdrawal and relapse, which are critical in addiction treatment.

Single-Cell Sequencing Techniques
Traditional classification divides cell types based on morphology rather than molecular features [31,32].Despite the fact that cells have almost identical genotypes, traditional RNA sequencing gives an average expression profile from a population of cells, ignoring cell-to-cell variability and transcriptome information from a subset of active genes [33][34][35].Single-cell RNA sequencing (scRNA-seq) is used for genetic profiling of cells that require prior separation of one cell type from a cell suspension using a fluorescence-activated cell sorting and flow cytometer [36].This approach requires a priori knowledge of a cell-typespecific marker, which is dependent on the availability of appropriate antibodies [37].Great progress in scRNA-seq was made by applying high-throughput droplet-based techniques, which do not require prior separation of the cells from the cell suspension, namely, In-Drop [38], Drop-seq [39], and 10× Chromium [40], supporting cost-effective capture and library production for thousands to millions of cells and allowing examination of gene expression heterogeneity among individual cells (Figure 2).Today, we have the possibility of simultaneous profiling of the transcriptome using 10× Chromium and epigenome using Assay for Transposase-Accessible Chromatin sequencing (scATAC-seq).The ATAC-seq technique is used for chromatin accessibility at the single-cell level by measuring the regions of DNA that are accessible to transposase enzymes, providing insights into the epigenetic landscape [41,42].Combined analysis of scRNA-seq and scATAC-seq data allows the simultaneous exploration of gene expression and chromatin accessibility within individual cells important in the study of addiction development [43,44].Finally, in the field of neurobiology, the Patch-seq approach stands out as a powerful adaptation of the patchclamp technique, allowing the simultaneous examination of individual neurons through a combination of electrophysiological, morphological, and transcriptomic analyses [45][46][47].These novel approaches provide a comprehensive understanding of neuronal properties by providing insights into both functional activities and molecular profiles' rising direction toward system biology's understanding of changes induced by drug taking, deprivation, and relapse.

Application of Single-Cell Sequencing in Addiction Research
By defining region-specific molecular signatures and neuronal circuits involved addictive behaviors [48][49][50], droplet-based high-throughput scRNA-seq enabled researc ers to compare gene expression profiles in different brain regions involved in addictio [51].At the molecular level, scRNA-seq of individual neurons identifies specific neuron subtypes associated with reward, reinforcement, and addiction processes [52], where transcriptome analysis identifies distinct transcriptional signatures associated with addi

Application of Single-Cell Sequencing in Addiction Research
By defining region-specific molecular signatures and neuronal circuits involved in addictive behaviors [48][49][50], droplet-based high-throughput scRNA-seq enabled researchers to compare gene expression profiles in different brain regions involved in addiction [51].At the molecular level, scRNA-seq of individual neurons identifies specific neuronal subtypes associated with reward, reinforcement, and addiction processes [52], whereas transcriptome analysis identifies distinct transcriptional signatures associated with addiction-related behaviors [53].Understanding how individual cells within a population exhibit different sensitivities or adaptations to drug use is crucial, while interactions between the nervous and immune systems must be integrated in this aspect.

Identifying Cell Populations by Molecular Clustering
scRNA-seq has been shown to be an effective tool for characterizing cellular diversity in brain regions associated with the reward system, such as the NAc and VTA [54,55].The NAc is a brain region primarily composed of two main types of cells, namely, Medium Spiny Neurons (MSNs) and interneurons.In NAc, clustering analysis revealed nine major cell populations, four neuronal and five non-neuronal [54].Using this approach, in the NAc, novel subpopulations of interneurons and MSNs were identified [52].scATAC-seq has been used to map cell-type-specific differences in chromatin accessibility in the NAc, providing insights into the epigenomic landscape of this brain region [55].Moreover, the characterization of γ-Aminobutyric Acid (GABA) MSNs and the discovery of notable variations in receptor expression patterns and MSN activation within NAc subterritories have highlighted the anatomical and functional heterogeneity of the NAc [56].The ventral tegmental area (VTA) is best known for containing dopaminergic neurons associated with reward and motivation.scRNA-seq identified in VTA selective markers for dopamine and combinatorial neurons revealed expression profiles for drugs of abuse receptors and population-specific enrichment of genes associated with brain disorders [57].This comprehensive molecular characterization highlights the heterogeneity of the NAc and VTA cell population, providing a valuable resource for future research into VTA and NAc gene expression and its implications for reward-related behaviors such as addiction.
By applying scRNAseq in the study of acute COC administration in rodents, distinct neuronal clusters within the NAc were discovered.These clusters exhibited known markers associated with two main types of MSNs: dopamine receptor D1-positive (Drd1-MSNs) and dopamine receptor D2-positive (Drd2-MSNs) [58].The same neuronal substrate that allows drugs of abuse to access Drd1 and Drd2 MSNs in NAc has been confirmed to augment and corrupt a shared pathway that normally serves physiological needs [59].Nevertheless, the investigation of genetically mediated variations in susceptibility to addiction-like behaviors was restricted by the use of inbred rodent strains.Furthermore, these studies primarily focused on acute drug treatments, ignoring the molecular changes associated with longterm addictive-like behaviors [60].
METH is another commonly abused psychostimulant that induces synaptic plasticity and pathological memory enhancement [61].Epigenetics plays an important role in regulating METH addiction [62].METH modulates dopamine (DA), norepinephrine (NE), serotonin, glutamate (Glu), and GABA neurotransmitters in the medial prefrontal cortex (mPFC), the VTA, and the NAc through histone acetylation, methylation, micro RNAs (miR-NAs), and ubiquitination.These epigenetic mechanisms do not regulate METH-induced addiction alone but rather collaborate with miRNA regulation of the Ubiquitin proteasome system [62].Cell clustering was performed on bulk RNA data using single-cell RNA data extracted from astrocytes for analysis of differently expressed genes in METH-exposed mice.The NF-κB signaling pathway, inflammation, and neurodegeneration were among the considerably enriched pathways found in the analysis.Immune infiltration analysis showed that the METH group had significantly lower neutrophil infiltration and significantly higher monocyte, T, and NK cell infiltration [63].Strong inflammatory responses occur early in METH withdrawal, and they decrease as withdrawal time increases.Cy-tokines, immune cells, complement, and immunoglobulin form a complex immune network that regulates immune responses after METH withdrawal, shifting the focus of METH addiction research from neurons to non-neuronal cell types [64].The novel approach in the study of molecular mechanisms is organoids.Brain organoids are generated through the differentiation of pluripotent stem cells within a three-dimensional culture environment, where neurodevelopmental signals are introduced to guide the process.Organoids closely mimic the cellular composition and anatomical structure of the brain [65].METH treatment of organoids resulted in transcriptional changes in astrocytes and neuronal progenitor cells, among other cell types.Additionally, METH induced novel networks of astrocytespecific gene expression that control cytokine responses and the inflammasome, resulting in modifications to the RNA and protein levels of neuroinflammatory and cytokine gene expression [66].After being exposed to METH, neutrospheres made from embryonic rat hippocampal tissue displayed abnormal cell differentiation in both neurons and astrocytes, as well as a decreased capacity for neutrosphere migration coupled with an increase in oxidative stress and apoptosis [67].Taken together, the scRNA-seq literature overview demonstrated that COC mainly influences neurons in VTA and NAc, while METH has a primary influence on the immune response.

Identifying Transcriptome Mechanism
By analyzing single-cell transcriptional responses in the prefrontal cortex (PFC) cells of mice undergoing COC self-administration, specific cell types that express genes associated with COC addiction were identified, namely, ∆FosB, Methyl CpG binding protein 2 (MeCP2) and brain-derived neurotrophic factor (BDNF) [68].Following COC administration, ∆FosB was reported to also be expressed in the VTA with dopamine neurons projecting to the NAc [69], MeCP2 was reported to be altered in the NAc [70], and BDNF expression played a role in the VTA-NAc pathway [71].These genes were selectively expressed in excitatory neurons, inhibitory neurons, and non-neuronal glial cells while showing varied expression levels over the stages of COC addiction [68].In addition to anatomic location, function and network connectivity have a significant impact on an excitatory neuron's transcriptional fate.COC, as rewarding stimuli, increases dopamine levels in NAc, which are responsible for COC's reinforcing effects.Dopamine dynamics is mediated by the activation of dopamine neurons in the VTA, which project to the PFC, hippocampus, and amygdala, separately from NAc cells [53].In the PFC of mice who self-administrated COC, some of the relevant genes had a specific expression pattern in different PFC cell clusters.Excitatory neurons were the predominant host for the expression of Drd1, Drd2, Gprotein-coupled receptors Gpr88 and Gpr6, and Regulator of G-protein signaling 9 (Rgs9), though Drd1, Gpr88, and Rgs9 were also detected in other cell types.Endothelial cells and microglia expressed the adenosine A2A receptor (Adora2a) the most, while inhibitory and excitatory neurons expressed it the least.According to computational analysis of publicly available scRNA-seq data, excitatory neurons in the PFC were the primary site of transcriptional reprogramming of relevant genes.This kind of transcriptional analysis of the addiction gene signature at cell-type-specific levels in the PFC may be useful in developing novel pharmacological strategies to combat addiction at the cellular level [53].Addiction has a significant impact on highly regulated gene networks in the Dorsum striatum and NAc as well as COC-related gene regulation.While these two brain regions exhibit distinct transcriptomes and addiction-related changes in gene expression, a subset of differentially expressed transcripts was identified that are regulated similarly by COC and METH [72].Further analysis supports the hypothesis that COC-responsive genes participate in the vulnerability to substance use disorders, as a convergence of up-regulated differentially expressed transcripts was found in individuals with confirmed addiction and genes associated with risk-taking behavior.It has been demonstrated that a significant amount of the striatal molecular pathology underlying COC addiction in humans can be replicated in mice that self-administer COC.This significant finding validates the use of mice models in the pathophysiological study of addiction.In particular, studies conducted on mice link aberrant gene expression patterns in humans with addiction to Drd1 MSNs.The study highlights that this connection was made using a human cohort that included a sizable number of Black individuals, who have not been adequately represented in previous transcriptional studies of addiction, even though there has been recent evidence that they are more susceptible to overdose deaths related to COC [72].
Reactive oxygen species and oxidative stress in striatal regions are known to be elevated by COC, and this effect is further compounded by elevated glutamate release and excessive dopamine levels.The potential involvement of both cell types in the regulation of conserved gene networks was revealed by integrating human transcriptomes with the Drd1 and Drd2 MSN transcriptome data from mice.Numerous transcription factors that are predicted to be upstream of the abnormalities in gene expression linked to addiction were identified through this analysis.The majority are shared between the NAc and dorsum striatum.Notably, several of these predicted upstream transcription factors have been implicated previously in COC addiction in rodent models; in particular, Activator protein 1 (AP-1) (FOS/JUN) family, EGR family, NFκB, E2F, and several nuclear hormone receptors cells [72].Simultaneously, a number of other transcriptional regulators that had not been previously linked were discovered to have mechanistic roles in COC addiction.Our knowledge of how COC and other addictive drugs cause transcriptomic plasticity in brain regions related to addiction has improved due to the use of animal models.Selfadministration procedures have, however, clearly demonstrated that abused drugs, such as COC, act as reinforcers in animal models.Nevertheless, it is still unclear whether the COCrelated molecular changes investigated in rodent models accurately reflect transcriptomic changes in human addiction.Preclinical animal models have been shown to be valid and valuable in providing important insight and mechanistic information across transcriptional networks, biological pathways, and even neuronal cell types in the striatum that lead to the harmful neuroadaptations of addiction, according to previous studies.To capture the unique molecular features of addiction, one must carefully consider treatment paradigms and experimental endpoints in these animal models.This is because chronic drug use will present differently at the transcriptional level depending on early versus late withdrawal times or relapse.Future evolution of single-cell and spatial sequencing technologies promises more direct insight into subtype-specific transcriptome alterations in human addiction cells [72].
In response to acute and chronic METH, mice hippocampi exhibit significant volumetric atrophy compared to controls.The genes involved in cytoskeleton organization and phagocytosis were downregulated in the acute METH-treated group compared to the control group.In the group receiving long-term METH treatment, genes linked to synaptic transmission, neuron differentiation regulation, and neurogenesis regulation were downregulated [73].Drd1 overstimulation after METH exposure induces metabolic changes and transcriptional pathways, switching gene expression and neuronal phenotype underlying addictive behavior [74].PKA phosphorylates voltage-dependent ion channels, GLUT receptors, transcription factors, and epigenetic enzymes involved in synaptic plasticity as naturally occurring in normal striatum.When the Drd1 is activated, PKA activates mitogen-activated protein kinases (MAPKs) and extracellular signal-regulated kinases 1/2 (ERK1/2) [75].Nuclear receptors, CREB, Elk-1, and H3 histones are all phosphorylated by ERK1/2 upon translocation to the nucleus, where they control the expression of certain genes [76].The DA-and cAMP-regulated phosphoprotein (DARPP-32) is a key substrate of DRD1/PKA signaling in the striatum [77].
METH consumption results in oxidative stress in the DA terminals due to excess of free DA undergoing oxidative metabolism and autoxidation, intracellularly and extracellularly.Along with hydrogen peroxide and reactive oxygen species (ROS) production, toxic DA metabolites like quinones and 3,4-Dihydroxyphenylacetaldehyde promote structural modifications of proteins.METH influences other cell organelles, endoplasmic reticulum, and mitochondria, leading to neurotoxic effects, while ROS accumulation leads to misfolded and insoluble proteins and severed organelles, decomposed by cell-clearing mechanisms of autophagy and ubiquitin proteasome [78].
ATACseq, a method used for mapping genome-wide chromatin accessibility, has been used for addiction research to explore the open chromatin regions associated with exposure to drugs of abuse [79,80] including METH [78].Changes in chromatin accessibility interact with DNA and histone modifications that are normally associated with transcription factor binding.Addictive drugs affect transcription factors like ∆FosB, EGRs, and MEF2, which regulate downstream target genes [81].
Patch-sequencing (Patch-seq), a variation of the patch-clamp technique, is a multimodal method that combines individual neuronal morphological, transcriptomic, and electrophysiological characterization.After electrophysiological recordings from individual neurons and morphologic reconstruction of the same cells, its cytoplasm is extracted and prepared for RNA-seq analysis.Approaches using Patch-seq have already made significant contributions to the identification of cell types and subtypes [45,47].This method was useful in the isolation of DA subtype cells that are influenced by METH consumption.DA neurons of the VTA-NAcLat circuit were identified and its gene expression was defined [82,83].

Epigenetic Mechanisms and Psychostimulant Addiction
Although single-cell sequencing technologies have transformed our understanding of cellular heterogeneity, single-cell data are challenging for downstream analysis [84].In addition, we have a process that alters gene activity without changing the DNA sequence, especially influenced by psychostimulants, which lead to neuronal plasticity [85].
Epigenetics refers to changes in the DNA molecule that do not alter the nucleotide sequence but control gene expression by inducing specific chemical modifications that regulate the activation or repression of genes and lead to different expression patterns [86].The predominant mechanisms of epigenetic modification in eukaryotes include DNA methylation, histone modifications, and non-coding RNA [87].In contrast to genetic modifications, epigenetic changes are reversible and respond to various environmental factors.DNA methylation is the most extensively studied epigenetic modification.It represents a stable gene silencing mechanism, which, in combination with modifications of histones and other chromatin-associated proteins, plays an important role in the regulation of gene expression and chromatin structure [88].In mammals, DNA methylation occurs primarily through covalent modification of cytosine residues in CpG dinucleotides-DNA sequences consisting of consecutive cytosine-guanine (CpG) nucleotide pairs.The precise DNA methylation patterns are generated and maintained by the cooperative activity of DNA methyltransferases (DNMTs).DNMT1 ensures the maintenance of existing methylation patterns, while DNMT3A and DNMT3B are responsible for the de novo establishment of new methylation patterns.Moreover, DNMTs are also responsible for the complex process of DNA demethylation together with the enzyme family of ten-eleven translocation (TET).Furthermore, the DNA is tightly wrapped around the histones, thus preventing the expression of specific genes.In particular, the N-terminal tails of histones can undergo post-translational covalent modifications at specific residues, influencing critical cellular processes such as transcription, replication, and DNA repair.Histone modifications primarily modulate the degree of chromatin condensation, thereby regulating gene expression based on the modified residues and the type of modifications [89,90].Distinct cell types exhibit specific patterns of histone modifications that contribute to cellular identity precisely by regulating gene expression.Enzymes that add or remove covalent modifications to histone proteins dynamically control these patterns.While histone deacetylases (HDACs) and histone demethylases (HDMs) remove acetyl and methyl groups, histone acetyltransferases (HATs) and histone methyltransferases (HMTs) add acetyl and methyl groups, respectively [90][91][92].Chromatin remodeling, transitioning from the condensed to the transcriptionally active state, controls gene expression: increased chromatin condensation impedes transcriptional and other factors from accessing the DNA molecule, suppressing, or disrupting transcription [90].
Mechanisms involved in chromatin remodeling include changes in histone composition within individual nucleosomes, various histone modifications, and the rearrangement of individual nucleosomes.All these mechanisms together influence chromatin compaction, thereby regulating transcriptional activity [87,89,93].Changes in the activity of enzymes that modify histones, known as writers, enzymes that remove modifications, known as erasers, and protein complexes that modulate gene expression, known as readers, have been linked to a variety of physiological changes and diseases [94,95].The number of known post-translational modifications (PTMs) of histones is increasing, supported by the development of technologies for their identification: antibody development, omics approaches in peptide analysis, and mass spectrometry [96,97].PTMs on individual histones are the focus of the current research, as they have been found to influence neuronal plasticity as the biological basis of learning and memory, but also the development of pathological forms of neuronal plasticity, such as in addiction [98].However, due to the complexity of the interactions between various epigenetic modifications and their multiplicity, developing an understanding of the mechanisms through which they change gene expression and related phenotypes is still challenging.In mammals, the administration of psychostimulants such as COC or METH causes a significant increase in dopamine signaling that is associated with a combinatorial cascade of histone PTMs.Such modifications in dopamine neurons in the reward-processing part of the brain result in behavioral consequences typical of the rewarding effects of psychostimulants.
The mechanisms of epigenetic regulation by RNA molecules are less explored compared to DNA methylation and histone modifications.Non-coding RNAs (ncRNAs) represent a diverse group of regulatory RNA molecules, some of which already have established functions, while others remain to be characterized [93].Based on their size, they are divided into two categories: long non-coding RNAs (lncRNAs), which comprise more than 200 nucleotides, and small ncRNAs, which comprise no more than 200 nucleotides [99].The ncRNAs relevant in the field of epigenetic regulation include lncRNAs, miRNAs (microRNA), siRNAs (small interfering RNA), and piRNAs (PIWI-interacting RNA) [93].

Cocaine-induced Epigenetic Modifications
Epigenetic modifications in the control of gene expression of specific brain regions in COC addiction are summarized in Table 1 and Figure 3.

DNA Methylation
Environmental stimuli in COC exposure can be translated to changes in gene expression and phenotypes that, through enzymatic modifications to DNA sequences, create long-lasting behavioral phenotypes.These methylated DNA regions are bound by methylbinding domain-containing proteins such as methyl CpG-binding protein 2 (MeCP2), which is important for recruiting co-repressors such as methyltransferases to the gene promoter [100].DNA methylation is critical for imprinting, X chromosome inactivation, and cell differentiation.Therefore, it is important that it can be altered at specific loci in germ cells by exposure to environmental factors such as toxins [101] and stress [102], and can therefore be inherited by offspring over multiple generations.Acute cocaine administration increases DNA methylation and DNMT3A and DNMT3B levels, while the binding of MeCP2 to specific gene promoters can decrease gene expression in the NAc [103].When cocaine is self-administered, MeCP2 is increased in both the striatum and PFC of rats during that response [104,105].Moreover, MeCP2 can be switched off in the striatum in association with an impaired COC-dependent increase in BDNF levels [104].It is demonstrated that following repeated COC exposure, including self-administration, in the NAc of mice and rats as well as in the NAc of COC-dependent humans, there is a decrease in protein-R-methyltransferase-6 (PRMT6) and its corresponding histone mark, asymmetric demethylation of R2 on histone H3 (H3R2me2a).First, PRMT6 down-regulation is observed specifically in NAc medium spiny neurons (MSNs) that express dopamine D2 receptors (D2-MSNs), whereas D1-MSNs exhibit the opposite regulation.

DNA Methylation
Environmental stimuli in COC exposure can be translated to changes in gene expression and phenotypes that, through enzymatic modifications to DNA sequences, create long-lasting behavioral phenotypes.These methylated DNA regions are bound by methyl-binding domain-containing proteins such as methyl CpG-binding protein 2 (MeCP2), which is important for recruiting co-repressors such as methyltransferases to the gene promoter [100].DNA methylation is critical for imprinting, X chromosome inactivation, and cell differentiation.Therefore, it is important that it can be altered at specific loci in germ cells by exposure to environmental factors such as toxins [101] and stress [102], and can therefore be inherited by offspring over multiple generations.Acute cocaine administration increases DNA methylation and DNMT3A and DNMT3B levels, while the binding of MeCP2 to specific gene promoters can decrease gene expression in the NAc [103].When cocaine is self-administered, MeCP2 is increased in both the striatum and PFC of rats during that response [104,105].Moreover, MeCP2 can be switched off in the striatum in association with an impaired COC-dependent increase in BDNF levels [104].It is demonstrated that following repeated COC exposure, including self-administration, in the NAc of mice and rats as well as in the NAc of COC-dependent humans, there is a decrease in protein-R-methyltransferase-6 (PRMT6) and its corresponding histone mark, asymmetric demethylation of R2 on histone H3 (H3R2me2a).First, PRMT6 down-regulation is observed specifically in NAc medium spiny neurons (MSNs) that express dopamine D2 receptors (D2-MSNs), whereas D1-MSNs exhibit the opposite regulation.

Histone Modifications Histones Methylation
Histones H3 and H4 can have their lysine or arginine residues methylated, with varying effects on transcription.While lysine methylation is connected to both transcriptional activation and repression depending on the methylation site, arginine methylation encourages transcriptional activation [106].This adaptability could be explained by the fact that, in contrast to acetylation, methylation has no effect on histone charge or histone-DNA interactions directly.H3 methylation increases the gene transcription, while lysine

Histone Modifications Histones Methylation
Histones H3 and H4 can have their lysine or arginine residues methylated, with varying effects on transcription.While lysine methylation is connected to both transcriptional activation and repression depending on the methylation site, arginine methylation encourages transcriptional activation [106].This adaptability could be explained by the fact that, in contrast to acetylation, methylation has no effect on histone charge or histone-DNA interactions directly.H3 methylation increases the gene transcription, while lysine (K) methylation on H3K9Me3 and H3K27Me3 at specific gene regions represses transcription [107].In addition to lysine (K) methylation, protein-R-methyltransferases (PRMT1 and 6) have the ability to methylate histone arginine (R), producing H3R2Me2 and H4R3Me2 in NAc following repeated cocaine exposure [108].

Histones Acetylation
COC administration leads to significant acetylation of histones H3 and H4 [109].H3K9, H3K14, H4K5, H4K8, H4K12, and H4K16 have emerged as the most extensively researched acetylation sites in the context of COC exposure [90].Histone acetylation in the NAc is associated with CREB-binding protein (CBP), a histone acetyltransferase [110].Repeated exposure to COC leads to a dual effect on histone acetylation, resulting in H3Ac/H4Acincreased and -decreased acetylation of gene promoters, attributed to the hypoacetylation of histones H3 and H4 [109].In the VTA during COC withdrawal, there is an observed enrichment of CBP in the promoter region of BDNF leading to an elevated level of histone acetylation, particularly H3K9/14Ac, at the Bdnf promoter [111].However, because the balance of acetylation and deacetylation is dynamic, increased HDAC5 deacetylase activity lowers COC preference in the Conditioned Place Preference (CPP) [112].A genetically modified version of HDAC5 that is expressed during COC administration, on the other hand, produces different results in various behavioral measures of COC effects [113][114][115].Furthermore, not only must the global level of acetylation and deacetylation be considered but also the location along the DNA where chromatin modulation occurs, i.e., the role of the PTM gene, which may act antagonistically on gene programs induced by psychostimulant administration [109].Acute COC exposure increases global histone H3 phospho-acetylation and histone H4 acetylation in the striatum [116,117].While long-term COC exposure increases the expression of HDAC5 target genes and decreases HDAC5 activity in the NAc, it also promotes histone acetylation [118].Ferguson et al. demonstrated that the overexpression of sirtuin (SIR1 and SIR2) in the NAc in response to chronic COC administration enhances the rewarding effects, identifying sirtuins as key mediators of cellular and molecular plasticity [119].Zell et al. revealed that m-opioid receptors may also regulate VTA neurotransmission by directly engaging glutamatergic neurons [120].Subclassing of the current dataset indicates that rather than directly modulating DA neurons, these effects might be mediated through GABAergic or glutamate/DA combinatorial populations.

Histones Ubiquitination
The ubiquitin-proteasome system (UPS), one of the epigenetic hallmarks, is a multifaceted network of ubiquitin ligases and proteasome structures that controls synaptic and epigenetic plasticity [121] and is also involved in memory processes and substance use disorders [122].E3 ubiquitin-protein ligase SMURF1 is a key mediator of neuroadaptations in the nucleus accumbens that follow COC exposure and mediates cue-induced COC seeking during withdrawal [123].All of these results point to the possibility that the SMURF1-SMAD1/5-RUNX2 pathway functions as a crucial transcriptional regulator that modifies plasticity after COC self-administration.The finding that RUNX2 and SMAD1/5 are upregulated in the NAc following COC self-administration allowed researchers to investigate RUNX2 binding to target genes that had previously been connected to COC plasticity.RUNX2 is found to bind at AP-1 sites and to the promoters of these genes following COC self-administration, indicating that it may be a master regulator of several pathways regulating COC-induced plasticity [124].There is a critical need for new research to identify novel potential therapeutic targets in order to develop effective treatments for COC use disorder.
In addition to histone acetylation, methylation, and ubiquitination, other histone modifications such as SUMOlation, propionylation, butyrylation, and the recently described dopa-mineylation and serotonylation play important roles in the regulation of psychostimulant-induced behavioral effects [80, 97,125].

ncRNAs
The most commonly investigated ncRNAs involved in epigenome modifications in COC addiction are miRNAs.They can regulate COC intake and potentially have an impact on developing compulsive consumption of the drug.The mechanisms that stand by miR-212 influences COC intake were also a point of exploration.The cAMP signaling cascade strongly activates the miR-212/132 gene cluster, with CREB upregulating both miRNAs.It was found that the pathway of miR-212 regulation dramatically boosts CREB signaling in cultured cells and also in the striatum of rats [104].miR-206 has been shown to negatively regulate BDNF expression, which is known to be important for the motivational effects of COC and other addictive substances [126].Additionally, it was demonstrated that the activation of miR-382 attenuates COC-induced increases in the transcription factor ∆FosB, which is known to be important in the rewarding effects of addictive drugs by negatively regulating the expression of the dopamine D1 receptor [127].

Methamphetamine-Induced Epigenetic Modifications
Epigenetic modifications in the control of gene expression of specific brain regions in METH addiction are summarized in Table 2 and Figure 4.

DNA Methylation
DNA methylation levels have been found to be changed in METH addicts [129], and moreover, this also occurs in their offspring [130].BDNF methylation was increased in the PFC of METH-addicted rats and patients but was decreased in the hippocampus of rats [131].The neurotoxic effects of METH exposure were partly caused by a decrease in BDNF expression [131,132].According to a report by Moszczynska et al., METH induction leads to impaired cognition and memory, accompanied by an increase in long-interspersed element-1 (LINE-1) activity in both the dentate gyrus of the hippocampus (Hip) and the dorsum striatum (DS) [133].METH improved spatial memory by downregulating the expressions of MeCP2 and DNMTs and reducing DNA methylation at the promoter of Synaptophysin (Syp) in the Hip.While METH increased the expression of MeCP2 and DNMTs, it also caused DNA hypermethylation at the Syp promoter in the PFC and decreased 3.2.1.DNA DNA methylation levels have been found to be changed in METH addicts [129], and moreover, this also occurs in their offspring [130].BDNF methylation was increased in the PFC of METH-addicted rats and patients but was decreased in the hippocampus of rats [131].The neurotoxic effects of METH exposure were partly caused by a decrease in BDNF expression [131,132].According to a report by Moszczynska et al., METH induction leads to impaired cognition and memory, accompanied by an increase in long-interspersed element-1 (LINE-1) activity in both the dentate gyrus of the hippocampus (Hip) and the dorsum striatum (DS) [133].METH improved spatial memory by downregulating the expressions of MeCP2 and DNMTs and reducing DNA methylation at the promoter of Synaptophysin (Syp) in the Hip.While METH increased the expression of MeCP2 and DNMTs, it also caused DNA hypermethylation at the Syp promoter in the PFC and decreased memory function.The rewarding effect of amphetamine was increased by the specific knockout of MeCP2 in the NAc [134].By reversing DNA methylation at Sy genes involved in METH-mediated changes in dendritic spines and synaptic transmission, oxytocin (OT) prevented METH-seeking behavior and relapse [135][136][137][138].

Histone Modifications Histones Methylation
According to recent research, epigenetic mechanisms that mediate drug-induced transcriptional and behavioral changes induced by METH consumption are mostly caused by histone modification [139].Indeed, it was shown that METH induced H3 methylation through increased tri-methylation of histone H3 at lysine 4 (H3K4me3) at the promoter site of chemokine receptor 2 associated with behavioral sensitization in mice [140].Iamjan et al. reported that rats given METH and humans with METH dependence, particularly METH-dependent psychosis, exhibit aberrant Brain-derived neurotrophic factor (BDNF) methylation [131].Salehzadeh et al. also confirmed this hypothesis and showed BDNF hypomethylation [132].METH consumption also influenced Long INterspersed Element (LINE-1) activity.Systemic METH administration boosts LINE-1 activity in two neurogenic niches of the adult rat brain in a promoter hypomethylation-independent manner.Neurological impairments can be triggered by altered levels of LINE-1 expression, indicating that LINE-1 induction may play a role in the development of cognitive impairments in METH users [133].

Histones Acetylation
Similar to COC, METH induces alterations in the acetylation levels of histones H3 and H4, influencing the expression of various enzymes in different brain regions [141].In NAc, METH exposure leads to differential acetylation changes on various histone lysine residues by regulating the protein levels of histone deacetylases [142].Acute METH exposure results in decreased H3 acetylation (H3K9Ac and H3K18Ac) and increased H4 acetylation (H4K5Ac, H4K8Ac, and H4K16Ac) [143].Chronic METH causes a reduction in histone H4 acetylation (H4K5Ac, H4K12Ac, and H4K16Ac) at glutamate (GLUT) receptor promoters, impacting the expression of alpha-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid (AMPA) and N-methyl-d-aspartate (NMDA) receptors, leading to oxidation, excitotoxicity, and neuroinflammation [144].However, differences in H3 and H4 acetylation patterns have been observed between single and repeated exposures, which may contribute to the escalating dose cravings observed in the drug-seeking behavior of individuals with METH addiction [18,145].Additionally, METH increased the Hip's expressions of H2BAc, H3K9Ac, and H4K12Ac as well as the DS' expression of H4K12Ac.In the medial prefrontal cortex (mPFC), repeated METH treatment increased H4AC enrichment at the promoters of D1, hypocretin (orexin) receptor 1 (HCRTR1), and NMDA 1, but decreased H3AC enrichment at the promoters of D2, HCRTR1, HCRTR2, histamine receptor H1 (HRH1), HRH3, and NMDA [135].In histone deacetylases (HDAC 1, 2, 3, 5, 6, 8, 9, and 11), a single treatment of METH reduces the mRNA level of HDAC3 in the NAc [146].On the other hand, rats' NAc and mPFC express more HDAC2 when given a large single dose of METH, while their HDAC1 expression is decreased.In mice, a single METH treatment lowers the amount of HDAC8 mRNA in the NAc [147].Similarly, the HDAC11 mRNA level decreased after acute METH exposure in the NAc in mice.METH induced activation of the D1 receptor (D1R) in the PFC, which affected HDAC1 and HDAC2 levels [135] and the miRNA 181a/d level in the ventral tegmentum area (VTA), increased the expressions of α-adrenergic receptors and the NMDA receptor subunit, and then regulated the function of dopamine receptors (DRs) [146,148].METH administration increased GluN2B, an NMDA receptor subunit expression and sequential activation of the ERK/CREB)/BDNF pathway.Histone acetylation has consequences at a given location and could be in and of during mediation of the targets in fully different neurons dopaminergic, glutamatergic, and GABAergic neurons (DA, Glu, and GABA).Therefore, enzymes controlling the transfer of acetyl groups are involved in METH addiction as well.ATF-2, a member of the ATF/CREB family, enhances CREB-dependent transcription and possesses intrinsic HAT activity on histone H4 [143].METH increases P300, one type of HAT that regulates Glu release induced by METH in human primary astrocytes [149].Another novel molecular mechanism that may affect METH-induced behavioral sensitivity is histone methylation.Methyltransferases (HMTs/KMTs) and demethylases (HDMs/KDMs)-such as KMT2A, an enzyme involved in H3K4me3-were linked to METH addiction.METH increased this enzyme, which is necessary for the development and maintenance of METH-associated memory.KDM5C is linked to METH because it demethylates H3K4.Using small interfering (siRNAs) in the NAc, some studies report decreased expressions of histone methyltransferase (HMT), histone lysine demethylase KDM5C, and mixed lineage leukemia 1 (Mll1).

Histones Ubiquitination
Neurotransmitter excitation and synaptic plasticity in brain disorders associated with drug addiction have been connected to the UPS, an enzymatic complex that controls proteolysis and turnover.Both pre-and post-synaptic neurons in the DA circuitry are impacted by UPS inhibition.Through endocytic internalization and degradation, the UPS reduces D1/D2-like DRs and Alphaamino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) [150].The UPS regulates the presynaptic release of Glu via both D1-like and D2-like receptors, which controls DA transmission [151], also influencing postsynaptic plasticity [152,153].As a result, the DA and Glu signaling pathways interact with UPS substrates [121].METH, despite being an extremely powerful DA releaser, impairs UPS activity, which is largely due to dopamine.As a result, pre-and post-synaptic neurons in the dopamine circuitry are extremely vulnerable to UPS inhibition [150].Parkin promotes the ubiquitination of substrate proteins, which aids in their degradation.Sharma et al. showed that studying the neuronal substrates underlying "resilience" or vulnerability to METH use disorder can be facilitated by using rats with excess or deficit PARKIN.Additionally, PARKIN may be a novel therapeutic target for the treatment of METH use disorder [154].After METH withdrawal, protein ubiquitination and E3 ubiquitin ligases are increased in the central amygdala.Research conducted by Cates et al. showed that ubiquitination of the central amygdala is linked to METH craving behavior [155].Synovial apoptosis inhibitor 1 (SYVN1) is an E3 ubiquitin ligase from the endoplasmic reticulum (ER).SYVN1 knockdown has been linked to METH-CPP by increasing GABAA1 in the DS [156].MiRNA-181a also regulates METH addiction via the ERAD pathway [136].

ncRNAs
Similarly, as in COC addiction, epigenetic modifications in METH addiction involving ncRNA predominantly include miRNAs.Using the KEGG pathway analysis, miRNAregulated genes were found to be involved in vesicular transport, METH addiction, the cyclic guanosine monophosphate cGMP-protein kinases G (PKG) signaling pathway, the dopaminergic synapse, and the GABAergic synapse [157].In genome-wide transcriptional profiling, the expression of multiple miRNAs is increased in the central amygdala alongside molecules related to METH addiction [155].METH increased the levels of miRNAs 237, 296, and 501 in the NAc.MiRNAs in the NAc regulate Wnt signaling and axon guidance genes [158].MiR-128 influenced METH-induced behavioral sensitization by altering synaptic plasticity-related molecules in the NAc [159].METH-induced locomotor sensitization is disrupted by Ago2-dependent miRNAs in the NAc.These Ago2/miR-3068-5p effects occur in conjunction with the glutamate receptor, GluN1/Grin1 [138].METH can increase the expression of miRNAs in the striatum, which harmed motor coordination and reduced striatal volume and dendritic length [160].

Epitranscriptomics and Psychostimulant Addiction
Epigenetic modifications to DNA and histones regulate gene transcription, whereas epitranscriptomic post-transcriptional RNA modifications influence gene expression [165].Although various types of RNA are involved in these changes, we have concentrated on messenger RNA (mRNA) modifications.The most common mRNA modifications are N1-methyladenosine (m 1 A), N6-methyladenosine (m 6 A), 5-methylcytosine (m 5 C), pseudouridine (Ψ), and others [166].m 6 A modification is one of the most abundant and reversible epitranscriptomic modifications, mediated by a set of proteins, which include methyltransferases ('writers'), demethylases ('erasers'), and m 6 A-binding proteins ('readers') [167].m 6 A methylation is associated with the control of mRNA metabolism, splicing, export, stability, translation, and degradation [168,169].Out of all tissues in the body, the brain has the highest abundance of m 6 A methylation, which is developmentally decreasing [168].In the brain, m 6 A methylation regulates neuronal transcripts and neuronal activity.Aside from its role in neuronal development [30], m 6 A modification is essential for the process of axon regeneration [170].Methyltransferases like-3 (METTL3) and like-14 (METTL14) [171][172][173], along with other proteins needed for m 6 A deposition, such as Wilms' tumor 1-associating protein (WTAP) [174] and RNA binding protein 15 (Rbm15) [175], form a stable protein complex that catalyzes m 6 A modification.Because they deposit RNA methylation modifications, these methyltransferases are collectively known as "writers".Fat mass obesity-associated protein (FTO) and alkB homologue 5 (ALKBH5) are two demethylases from the family of α-ketoglutarate dependent dioxygenases that can reverse m 6 A modification because it is dynamically regulated [176,177].Demethylases are known as "erasers" because they remove RNA methylation modifications.Posttranscriptional, site-specific adenosine-to-inosine base conversions, known as RNA editing, contribute to gene expression diversity and are catalyzed by Adenosine deaminases acting on RNA (ADARs) [178].Pseudouridine Synthase 7 (PUS7) is one of the major mRNA-modifying enzymes leading to pseudouridine (Ψ), a ubiquitous RNA modification [179].Both ADAR and PUS7 can lead to further mRNA modifications.m 6 A modifications are in direct connection to the so-called m 6 A "reader" proteins, which recognize the modified site.The proteins with YTH domains, which can specifically bind m 6 A through their YTH domain, are the most well-studied m 6 A readers.Fragile X mental retardation protein (FMRP) was also reported to be m 6 A reader and plays critical roles in synaptic plasticity and neuronal development.The identification of methylated nucleosides (m 6 A, m 5 C, m 1 A) is performed using immunoprecipitation and their variants using antibodies against methylated nucleosides or associated proteins methyltransferases, demethylases, and binding proteins.After fragmenting the RNA, fragments containing modified nucleosides are enriched before sequencing.The immunoprecipitate is analyzed using next-generation sequencing (NGS) to identify and map the modification [180].
The connection between epigenetic regulation and m 6 A RNA modification was associated with histone H3 trimethylation at Lys36 (H3K36me3), a marker for transcription elongation, which guides m 6 A deposition globally connected through METTL14 (DOI: 10.1038/s41586-019-1016-7).m 6 A modifications are mainly associated with neuronal plasticity in the brain, which is a consequence of learning and memory, and most of the literature is based on this line of research together with neurodegenerative disorders [176,181,182].Indeed, deficiency in m 6 A-dependent pathways significantly impairs neuronal function including dopamine signaling and dopamine-dependent learning.
Lowering neuronal m 6 A by overexpressing FTO or by adding m 6 A inhibitor led to the induction of N-methyl-d-aspartate (NMDA) receptor 1 expression, elevated oxidative stress, and Ca 2+ influx, resulting in dopaminergic neuron apoptosis [183].In addition, it was shown that the overexpression of FTO delays the dephosphorylation of CREB, increases the expression of the CREB, and targets neuropeptide receptor 1 (NPY1R) and BDNF known to regulate food intake and energy homeostasis [184].FTO affects dopamine (D2)-dependent responses to reward learning in meso-striato-prefrontal regions, suggesting a mechanism by which genetic predisposition alters reward processing not only in obesity but also in other disorders with altered D2R-dependent impulse control, such as addiction [185].Other than that, FTO was reported to be part of the regulation of BDNF processing [186].FTO demethylase, but not METTL3 and METTL14, was downregulated in the hippocampus following COC-induced CPP, leading to a higher level of m 6 A [187].Mice lacking the FTO gene exhibited lower body weight and decreased anxiety-and depression-like behaviors, mediated by changes in the gut microbiota [188].

Challenges of the Rewarding Molecular Pathway
Disrupting natural reward pathways, such as with food overeating leading to obesity and addictive drug usage leading to addiction, share many molecular regulatory sites in the central nervous system.Both influence neurotransmitters, leading to the dysregulation of the reward circuitry in the brain, and they have been studied as a complex interplay of genetic and environmental factors influencing their development.This review emphasizes the importance of the FTO gene, which, together with CREB and BDNF, regulates synaptic plasticity, but is also associated with cellular metabolism, oxidative stress in the brain, and the cycle of inflammation.Based on this, antioxidant and anti-inflammatory synthetic molecules have been implicated in the treatment of drug addiction.Future sc-RNA-seq studies as well as epigenetic and epitranscriptomic studies should focus on NAc D1R and D2R neurons, as these encode positive valence and reward responses to drugs of abuse and food and have a projection to the VTA.

Conclusions
In this review article, we covered single-cell sequencing research as well as epigenetic and epitranscriptomics methods applied to the brain regions VTA, NAc, and DS after both short-term and long-term exposure to the psychostimulants COC and METH.The single-cell method, which clustered cell types based on gene expression variations, was used to validate both known and unknown cell types implicated in the psychostimulant

Figure 1 .
Figure 1.Modifications to the Central Dogma of Molecular Biology: Gene Regulation Modulation by Environmental Factors.(A) Central Dogma: a theory stating that genetic information flows only in one direction, from DNA to RNA to protein, or RNA directly to protein.(B) Environmental Modulation of Central Dogma: Impact on DNA Unfolding, Transcription, and Translation Processes.Epigenome modification = DNA methylation or acetylation and post-translational modifications of histone tails that include phosphorylation, ubiquitination, acetylation, and methylation.Epitranscriptome modifications = RNA modifications primary messenger RNA (mRNA) modifications: N1-methyladenosine (m 1 A), N6-methyladenosine (m 6 A), 5-methylcytosine (m 5 C), pseudouridine (Ψ), and others.Created with BioRender.com(accessed on 11 December 2023).

Figure 2 .
Figure 2. Genomics, transcriptomics, and epigenomics technologies focused on the characterizati of individual cells from tissue samples.(A) Conventional massive population sequencing provid average expression signals for different cells, ignoring cell-to-cell variability.(B) By applying hig throughput droplet-based scRNA-seq to complex cell populations, it is possible to uncover differe cell types and interactions between cells, follow the development of distinct cell lineages, tra changes in gene expression, and identify regulatory relationships between genes.Created with B oRender.com (Accessed on 11 December 2023).

Figure 2 .
Figure 2. Genomics, transcriptomics, and epigenomics technologies focused on the characterization of individual cells from tissue samples.(A) Conventional massive population sequencing provides average expression signals for different cells, ignoring cell-to-cell variability.(B) By applying highthroughput droplet-based scRNA-seq to complex cell populations, it is possible to uncover different cell types and interactions between cells, follow the development of distinct cell lineages, track changes in gene expression, and identify regulatory relationships between genes.Created with BioRender.com(accessed on 11 December 2023).

Figure 3 .
Figure 3. Illustration of aberrant alterations in COC-induced epigenetic modifications.Green and red arrows indicate upregulation and downregulation, respectively.Created with BioRender.com(Accessed on 13 December 2023).

Figure 3 .
Figure 3. Illustration of aberrant alterations in COC-induced epigenetic modifications.Green and red arrows indicate upregulation and downregulation, respectively.Created with BioRender.com(accessed on 13 December 2023).

Figure 4 .
Figure 4. Illustration of aberrant alterations in METH-induced epigenetic modifications.Green and red arrows indicate upregulation and downregulation, respectively.Created with BioRender.com(Accessed on 13 December 2023).

Figure 4 .
Figure 4. Illustration of aberrant alterations in METH-induced epigenetic modifications.Green and red arrows indicate upregulation and downregulation, respectively.Created with BioRender.com(accessed on 13 December 2023).

Table 1 .
The role of epigenetic modifications in the control of gene expression of specific brain regions in cocaine addiction.