Processing math: 100%
Research article

Classification accuracy and functional difference prediction in different brain regions of drug abuser prefrontal lobe basing on machine-learning

  • Received: 25 April 2021 Accepted: 15 June 2021 Published: 23 June 2021
  • Taking different types of addictive drugs such as methamphetamine, heroin, and mixed drugs causes brain functional Changes. Based on the prefrontal functional near-infrared spectroscopy, this study was designed with an experimental paradigm that included the induction of resting and drug addiction cravings. Hemoglobin concentrations of 30 drug users (10 on methamphetamine, 10 on heroin, and 10 on mixed type) were collected. For these three types of individuals, the convolutional neural networks (CNN) was designed to classify eight Brodmann areas and the entire prefrontal area, and the average accuracy of the three classifications on each functional area was obtained. As a result, the classification accuracy was lower on the left side than on the right in the dorsolateral prefrontal cortex (DLPFC) of the drug users, while it was higher on the left than on the right in the ventrolateral prefrontal cortex (VLPFC), frontopolar prefrontal cortex (FPC), and orbitofrontal cortex (OFC). Then the differences in eight functional areas between the three types of individuals were statistically analyzed, and results showed significant differences in the right VLPFC and right OFC.

    Citation: Banghua Yang, Xuelin Gu, Shouwei Gao, Ding Xu. Classification accuracy and functional difference prediction in different brain regions of drug abuser prefrontal lobe basing on machine-learning[J]. Mathematical Biosciences and Engineering, 2021, 18(5): 5692-5706. doi: 10.3934/mbe.2021288

    Related Papers:

    [1] Haitao Jiang, Jiajia Guo, Hongwei Du, Jinzhang Xu, Bensheng Qiu . Transfer learning on T1-weighted images for brain age estimation. Mathematical Biosciences and Engineering, 2019, 16(5): 4382-4398. doi: 10.3934/mbe.2019218
    [2] Zhongwen Jia, Lihan Tang, Jidong Lv, Linhong Deng, Ling Zou . Depression-induced changes in directed functional brain networks: A source-space resting-state EEG study. Mathematical Biosciences and Engineering, 2024, 21(9): 7124-7138. doi: 10.3934/mbe.2024315
    [3] Sadia Anjum, Lal Hussain, Mushtaq Ali, Adeel Ahmed Abbasi, Tim Q. Duong . Automated multi-class brain tumor types detection by extracting RICA based features and employing machine learning techniques. Mathematical Biosciences and Engineering, 2021, 18(3): 2882-2908. doi: 10.3934/mbe.2021146
    [4] Yujie Kang, Wenjie Li, Jidong Lv, Ling Zou, Haifeng Shi, Wenjia Liu . Exploring brain dysfunction in IBD: A study of EEG-fMRI source imaging based on empirical mode diagram decomposition. Mathematical Biosciences and Engineering, 2025, 22(4): 962-987. doi: 10.3934/mbe.2025035
    [5] Zijian Wang, Yaqin Zhu, Haibo Shi, Yanting Zhang, Cairong Yan . A 3D multiscale view convolutional neural network with attention for mental disease diagnosis on MRI images. Mathematical Biosciences and Engineering, 2021, 18(5): 6978-6994. doi: 10.3934/mbe.2021347
    [6] Dongwei Liu, Ning Sheng, Tao He, Wei Wang, Jianxia Zhang, Jianxin Zhang . SGEResU-Net for brain tumor segmentation. Mathematical Biosciences and Engineering, 2022, 19(6): 5576-5590. doi: 10.3934/mbe.2022261
    [7] Gayathri Vivekanandhan, Mahtab Mehrabbeik, Karthikeyan Rajagopal, Sajad Jafari, Stephen G. Lomber, Yaser Merrikhi . Applying machine learning techniques to detect the deployment of spatial working memory from the spiking activity of MT neurons. Mathematical Biosciences and Engineering, 2023, 20(2): 3216-3236. doi: 10.3934/mbe.2023151
    [8] Yunyuan Gao, Zhen Cao, Jia Liu, Jianhai Zhang . A novel dynamic brain network in arousal for brain states and emotion analysis. Mathematical Biosciences and Engineering, 2021, 18(6): 7440-7463. doi: 10.3934/mbe.2021368
    [9] Yin Zhang, Xintong Wu, Jingwen Sun, Kecen Yue, Shuangshuang Lu, Bingjian Wang, Wenjia Liu, Haifeng Shi, Ling Zou . Exploring changes in brain function in IBD patients using SPCCA: a study of simultaneous EEG-fMRI. Mathematical Biosciences and Engineering, 2024, 21(2): 2646-2670. doi: 10.3934/mbe.2024117
    [10] Hassan Ali Khan, Wu Jue, Muhammad Mushtaq, Muhammad Umer Mushtaq . Brain tumor classification in MRI image using convolutional neural network. Mathematical Biosciences and Engineering, 2020, 17(5): 6203-6216. doi: 10.3934/mbe.2020328
  • Taking different types of addictive drugs such as methamphetamine, heroin, and mixed drugs causes brain functional Changes. Based on the prefrontal functional near-infrared spectroscopy, this study was designed with an experimental paradigm that included the induction of resting and drug addiction cravings. Hemoglobin concentrations of 30 drug users (10 on methamphetamine, 10 on heroin, and 10 on mixed type) were collected. For these three types of individuals, the convolutional neural networks (CNN) was designed to classify eight Brodmann areas and the entire prefrontal area, and the average accuracy of the three classifications on each functional area was obtained. As a result, the classification accuracy was lower on the left side than on the right in the dorsolateral prefrontal cortex (DLPFC) of the drug users, while it was higher on the left than on the right in the ventrolateral prefrontal cortex (VLPFC), frontopolar prefrontal cortex (FPC), and orbitofrontal cortex (OFC). Then the differences in eight functional areas between the three types of individuals were statistically analyzed, and results showed significant differences in the right VLPFC and right OFC.



    The brain function of people will be impaired after an extended period of drug abuse across a diverse variety. The classification accuracy of the eight encephalic regions (ERs) in the left and right hemispheres of the prefrontal lobe was analyzed on the basis of the fNIRS data acquired to identify the differences, if any. Statistical analysis was then performed on the eight ERs to find out which of them stand out in drug abusers. The study Innovations: 1) Drug abusers were classified; 2) The differences in classification accuracy were identified by classifying the drug abusers in conjunction with the ERs via CNN; 3) ERs with discernible distinctions in the drug abusers were identified through statistical analysis on the eight ERs in the left and right hemispheres of the prefrontal lobe.

    In recent years, near-infrared research is still a new field, and its advantages are gradually recognized by researchers [1]. NIRS in optical brain imaging is a spectral measurement method based on scalp detection, which measures the hemodynamic function of brain tissue using optical injection and detection points, and can record blood-oxygen level [2,3]. Compared with electroencephalography (EEG), a novel mode of brain functional imaging, it has the advantages of convenient wearing, high anti-interference ability and portability [4]. In human body, oxyhemoglobin (HbO2) and deoxyhemoglobin (Hbb) have specific absorption for near-infrared light at 700–900 nm. However, other biological tissues in the brain are relatively transparent in this wavelength range. Therefore, the changes in near-infrared light intensity at 760 and 850 nm emitted into the brain were measured, and the indirect changes in brain functional activity were based on the hemodynamic data transformed according to the Beer-Lambert law [5].

    Song et al, using fNIRS to extract oxy-Hb signals from schizophrenia patients and healthy people to construct a prefrontal brain network, and extract features from it to classify schizophrenia patients and healthy people, with a total accuracy rate of 85.5% [6]. Hu et al. propose a hybrid, multi-dimensional fusion structure of spatial and temporal segmentation model for the automated detection of thermography defects. In the model, creatively designed attention blocks the resources for local interaction along neighboring pixels, allowing for an adaptive recalibration of the feature maps. The performance of the two classifications is improved [7]. Kim et al. used NIRS technology to study the relationship between walking speed and blood oxygen activation in brain regions [8]. The proposed network of Koh et al. addresses the need in deep learning to match the data function of a time series with an appropriate network structure. This improves the classification accuracy of EEG and human activity signals [9]. Holtzer researched that changes in walking step length lead to activation of the prefrontal cortex, while Caliandro found that when the walking step length increases, the blood oxygen concentration of the lateral frontal cortex also shows an increasing trend [10,11]. Rodriguez-Rodriguez et al., in the process of monitoring type 1 diabetes mellitus (DM1), creatively used machine learning technology, developed six-feature selection techniques and four predictive algorithms, obtaining the best average performance and the best algorithm classification accuracy [12]. Li et al. calculated the normal people and schizophrenia patients in the 60 s task state [13]. Hennrich et al. collected human near-infrared signals and compared the classification performance of deep neural networks (DNN) for three mental tasks [14]. Abibullaev et al. used DNN to classify four mental tasks with an accuracy rate of 94% [15].

    fNIRS is used in the rehabilitation of stroke patients with motor imaging therapy, which can better monitor the changes in the treatment process [16]. Akiyama et al., using fNIRS to study the bilateral frontal cortex of depression patients [17]. fNIRS monitors the neonatal brain development stage, and its cerebral oxygen metabolism and hemodynamic characteristics are significantly different from those of adults [18].

    fNIRS is mainly aimed at evaluating and imaging the degree of brain damage in drug addicts, and detecting differences in the activation of the forehead brain area of people with different degree of addiction and Long-term methamphetamine users mid-near infrared detection showed abnormal activation of the prefrontal cortex, orbital frontal cortex, and the anterior cingulate gyrus [19,20]. Now fNIRS has been applied in many fields, such as the research of brain development, the research of perception and cognition, the research of motor control, the research of mental illness, the research of stroke and brain injury [21,22,23,24].

    According to the demographic scale, oral statement and medical examination, 30 drug addicts were selected. According to the experimental paradigm containing resting state and drug addiction induction, the near-infrared hemoglobin concentration data was collected. There are 30 subjects in total, classified into three classify. and Based on the CNN model, classification is performed.

    Study participant criteria: 1) Meet the diagnostic criteria for disorders caused by DSM-5 psychoactive substances; 2) Patients within six months of a withdrawal period; 3) Junior high school education and above; 4) Age 18–41 years; 5) Voluntary study participation and willingness to sign the informed consent form. Exclusion criteria: 1) Severe cognitive dysfunction, inability to cooperate and complete project-related assessment and testing; 2) Patients with severe physical diseases; 3) Patients with severe psychotic symptoms; 4) Current participation in other psychological interventions and treatments; 5) Current abuse of other substances that affect mental activities (except nicotine). The study was following with the declaration of Helsinki and was approved by the Ethics Committee of Shanghai University (Approval No. ECSHU2020-071).

    According to the demographic scale prepared in the early stage, data of 30 subjects (all men), such as drug use type, drug history, average drug dosage, and drug use frequency, were collected and recorded. Moreover, combining these data with personal health information collected from hospital records and experts' discussion and suggestions, three types of drugs were finally screened out, including methamphetamine, heroin, and mixed drugs, where ten subjects each were taking methamphetamine, heroin, and mixed drugs. Table 1 non-sensitive information of participants. METH: means Methamphetamine.

    Table 1.  Non-sensitive information of participants.
    N 30
    Sex (M/F) Male
    Age range (Year) 19–41
    Years of education 10 ± 2.72
    Years of drug abuse 8.2 ± 4.74
    Drug abuse per week: Most people 3–5 times a week.
    withdrawals 1.4 ± 0.98
    Reasons for taking drugs decompression needs; sex; emotion; curiosity; emotional frustration
    Subject No. 1 2 3 4 5 6 7 8 9 10
    Age: Year 41 28 25 35 35 39 29 35 33 32
    Type of drug METH METH METH METH METH METH METH METH METH METH
    Drug history 12 9 8 9 17 12 11 7 7 12
    Subject No. 11 12 13 14 15 16 17 18 19 20
    Age: Year 37 31 38 26 32 29 37 26 38 32
    Type of drug Heroin Heroin Heroin Heroin Heroin Heroin Heroin Heroin Heroin Heroin
    Drug history 14 6 13 8 4 5 12 13 17 17
    Subject No. 21 22 23 24 25 26 27 28 29 30
    Age: Year 37 31 31 24 40 34 35 40 39 36
    Type of drug Heroin, Cannabis, LSD, MDMA, Molly, Magu Pethidine hydrochloride, heroin, MDMA, ketamine, METH Cocaine, MDMA, Molly, ketamine, LSD, cannabis, Cannabis, Cocaine, LSD, trimethoprim, MDMA, Molly Heroin, Meth, MDMA METH, heroin, ketamine, MDMA METH, heroin, pethidine hydrochloride Heroin, METH, ketamine, MDMA, cannabis Heroin, METH, MDMA, ketamine, pethidine hydrochloride Heroin, METH, MDMA, ketamine, cannabis
    Drug history 16 11 7 6 19 13 11 20 20 19

     | Show Table
    DownLoad: CSV

    NIRSIT (OBELAB, Seoul, Korea) and a light source with dual wavelength vertical-cavity surface-emitting laser were used with the following technical spectrum: continuous wave, 780 nm and 850 nm; spatial resolution, 4 × 4 mm2; time resolution, 8.13 Hz; number of light points, 24; number of detection points, 32; and detection depth, 0.2–1.8 cm. The NIRS system used in the experiment allowed measurement of signals from four Source-Detector separations (15, 21.2, 30 and 33.5 mm) and hemodynamic reaction changes at different depths [25,26]. It is a functional NIRS device with 204 channels. Figure 1 NIRSIT wearing method in the experiment.

    Figure 1.  NIRSIT wearing method in the experiment.

    The NIRSIT channel and functional area are shown in Figure 3. The four advanced functional areas detected through the forehead near-infrared device are as follows: the dorsolateral prefrontal cortex (DLPFC), ventrolateral prefrontal cortex (VLPFC), frontopolar prefrontal cortex, and orbital frontal cortex (OFC). Regarding specific channel distribution, the right DLPFC has 1, 2, 3, 5, 6, 11, 17 and 18 channels. The left dorsolateral prefrontal lobe has 19, 20, 33, 34, 35, 38, 39 and 43 channels. The VLPFC of the left and right hemispheres have 4, 9, 10, 40, 44 and 45 channels. The left and right OFC have 14, 15, 16, 29, 30, 31, 32, 46, 47 and 48 channels. The frontopolar prefrontal cortex electrode arrangement has 7, 8, 12, 13, 21, 22, 23, 24, 25, 26, 27, 28, 36, 37, 41 and 42 channels.

    Figure 2.  NIRSTI channel and functional area division.
    Figure 3.  The example of drug pictures.

    When light passes through a uniform, non-scattering medium, only the absorption of photons by the medium is considered. According to Beer-Lambert law, the attenuation of light intensity is expressed as:

    OD=logII0=ε(λ)cdlge (1)

    Among them, I0is the incident light intensity, I is the output light intensity, ε(λ) is the extinction coefficient of the substance at the wavelength λ, which is determined by the absorbing medium and the wavelength of the light, c is the medium concentration, and the depth L of the penetration medium. Define the absorption coefficient as μa:

    μa(λ)=ε(λ)c (2)

    The total absorption coefficient in the medium can be expressed as the linear superposition of the absorption coefficient of each medium:

    μa(λ)=Niεi(λ)ci (3)

    The optical density can be expressed as:

    OD=logII0=Niεi(λ)ciL (4)

    The absorbance OD represents the degree of light attenuation after passing through the medium, which is related to the absorption coefficient μa of the medium, depth L of the penetration medium, and wavelength λ of the irradiated light. The available formula is as follows:

    OD=μa(λ)L (5)

    Among them, the absorption coefficient μa can be expressed as the product of the chromophore concentration c and its extinction coefficient ε. Considering the scattering problem, a differential path factor "DPF (λ)" is generally introduced, so it can be further expressed as follows:

    OD=cε(λ)LDPF(λ) (6)

    In the near-infrared band from 650 to 1000 nm, the main absorber of photons is hemoglobin. Although other chromophores also absorb light in this range, their contribution is assumed to be constant so that the change in the measured absorption signal over time can be attributed only to hemoglobin, and it can be expressed as follows:

    ΔOD(λ)=(Δ[Hbo]εHbo(λ)+Δ[Hbb]εHbb(λ))LDPF(λ) (7)

    In the formula, [Hbo] and [Hbb] represent the concentrations of oxyhemoglobin and deoxyhemoglobin, respectively. Changes in the concentrations of oxyhemoglobin and deoxyhemoglobin can be obtained by deformation:

    Δ[Hbo]=εHbb(λ1)ΔOD(λ2)DPF(λ2)εHbb(λ2)ΔOD(λ1)DPF(λ1)(εHbb(λ1)εHbo(λ2)εHbb(λ2)εHbo(λ1))L (8)
    Δ[Hbb]=εHbo(λ2)ΔOD(λ1)DPF(λ1)εHbo(λ1)ΔOD(λ2)DPF(λ2)(εHbb(λ1)εHbo(λ2)εHbb(λ2)εHbo(λ1))L (9)

    We used E-prime software package (Psychology Software Tools, Pittsburgh, PA) to write the experimental paradigm, with each picture numbered. The experimental paradigm consisted of three stages.

    Figures 3 and 4, are examples of drug abuse-related images and neutral images used in the experimental paradigm. Figure 5 The whole process of experimental paradigm. P means: drug picture; N means: Neutral Picture;

    Figure 4.  The example of neutral pictures.
    Figure 5.  The whole experimental structure.

    The first stage of the experiment, which lasted for 10 minutes, involved the participants closing their eyes for five minutes and then leaving their eyes open for five minutes.

    The second stage of the experiment, which lasted for six minutes, included the drug and neutral pictures. Each time block lasted for 10 seconds, for a total of 16 pictures, which were each displayed for 0.6 seconds. At the beginning, the first four pictures contained two drug images randomly displayed in a group, and the remaining 12 neutral images were displayed randomly. After a block had ended, a 4-second interval picture was shown, with a white background and a black cross.

    The third stage of the experiment, which lasted for 4.6 minutes, a completely neutral image was shown. Each time block lasted for 10 seconds and a total of 16 pictures were shown. The picture display time was 0.6 seconds, with a 4-second interval picture shown between each block.

    Methamphetamine easily causes intense excitement, which is difficult to eradicate after addiction [27,28].

    Heroin is a psychoanaesthetic drug. Once a person becomes addicted, their physiological reaction is intense, and they have a compulsion to seek medication [29,30].

    "Mixed drug abusers" refers to drug users who attempt to mix two or three drugs at a time. A major feature of mixed drug users is that they are the first to be exposed to new drugs. Mixed drugs have strong psychedelic properties, which causes the body to produce neurochemicals and form a dependence on them [31,32,33].

    Data on a total of 30 classified drug users were collected, each drug user containing 56 trials. Each trial corresponded to the channels of near-infrared data 0.625 seconds after the drug image appears. The whole CNN network has 20 layers, including 24 subjects in the training data, 3 in the validation data, and 3 in the testing data. Figure 6, CNN classification model of different types of drug abuse. The CNN model parameter are shown in Table 2.

    Figure 6.  CNN classification model of different types of drug abuse.
    Table 2.  CNN parameters.
    Layer (type) Output Shape Parameter #
    conv2d (Conv2D) (None, 5, 8, 32) 96
    activation (Activation) (None, 5, 8, 32) 0
    average_pooling2d (None, 5, 4, 32) 0
    conv2d_1 (Conv2D) (None, 5, 4, 16) 2064
    activation_1 (Activation) (None, 5, 4, 16) 0
    average_pooling2d_1 (None, 5, 3, 16) 0
    conv2d_2 (Conv2D) (None, 5, 3, 32) 1056
    activation_2 (Activation) (None, 5, 3, 32) 0
    average_pooling2d_2 (None, 5, 2, 32) 0
    conv2d_3 (Conv2D) (None, 5, 2, 32) 2080
    activation_3 (Activation) (None, 5, 2, 32) 0
    average_pooling2d_3 (None, 5, 1, 32) 0
    flatten (Flatten) (None, 160) 0
    dense1 (Dense) (None, 160) 25760
    dropout (Dropout) (None, 160) 0
    Dense2 (Dense) (None, 160) 25760
    dropout_1 (Dropout) (None, 160) 0
    dense3 (Dense) (None, 160) 25760
    dropout_2 (Dropout) (None, 160) 0
    dense4 (Dense) (None, 4) 644

     | Show Table
    DownLoad: CSV

    NIRS networks are designed with Butterworth filters. The expression of n-order Butterworth filter is:

    |H(f)|2=11+(ffc)2n=11+ϵ2(ffp)2n (10)

    where n is the order, fc is the cutoff frequency, and fp is the passband edge frequency.

    Convolutional layers are the core of convolutional neural networks. Its main role is to extract features from the input data. The calculation form is as follows:

    xlj=f(iMjxl1iklij+blj) (11)

    xlj is the jth feature of the layer l.klij is the jth feature of the layer l and the ith feature of the layer l1. blj is a bias parameter, f() is the activation function.

    The calculation form is as follows:

    xlj=f(βl1down(xl1j)+blj) (12)

    xlj is the jth feature of the layer l.βl1 is the Subsampling coefficient.blj is the bias parameter, down() is a sub-sampling function, f() is the activation function.

    The specific implementation process is as follows:

    Input: Values of x over a mini-batch:B={x1xm};
    Parameters to be learned:γ,β
    Output: {yi=BNγ,β(xi)}
    1. Calculate the mean of each mini-batch:
    μB=1mmi=1xi
    2. Calculate the variance of each mini-batch:
    σ2B1mmi=1(xiμB)2
    3. Normalize the data:
    ˆxi=xiμBσ2B+ε
    4. Scale and shift:
    yi=γ^xi+βBNγβ(xi)

     | Show Table
    DownLoad: CSV

    In this paper, the activation function uses a modified linear unit (ReLU), and the formula is as follows:

    f(x)=max(0,x) (13)

    The calculation formula is as follows:

    hw,b(x)=θ(wTx+b) (14)

    hw,b(x) is the output value of the neuron.x is the input feature vector of the neuron.w is the weight.b is the bias parameter θ() is the activation function.

    The softmax function formula is as follows:

    zi=Softmax(oi)=exp(oi)cexp(oc) (15)

    oi is the value of the output neuron corresponding to the Ith category

    Under the same CNN network structure and parameters, there are DLPFC, VLPFC, FPC and OFC prefrontal cortical regions in the left and right hemispheres of the brain. The average accuracy of 10 times obtained is shown in the Table 3.

    Table 3.  Average accuracy.
    Prefrontal cortex Channel Accuracy (%)
    Left-DLPFC 1, 2, 3, 17, 18, 5, 6, 11 51–61%
    Right-DLPFC 19, 20, 33, 34, 35, 38, 39, 43 56–72%
    Left-VLPFC 4, 9, 10 50–60%
    Right-VLPFC 40, 44, 45 40–50%
    Left-FPC 7, 8, 21, 22, 12, 13, 25, 26 55–70%
    Right-FPC 23, 24, 36, 37, 27, 28, 41, 42 50–60%
    Left-OFC 14, 15 16, 29, 30 50–62%
    Right-OFC 31, 32, 46, 47, 48 49–64%
    ALL-PFC 1, 2, 3, 17, 18, 5, 6, 11, 19, 20, 33, 34, 35, 38, 39, 43, 4, 9, 10, 40, 44, 45, 7, 8, 21, 22, 12, 13, 25, 26, 23, 24, 36, 37, 27, 28, 41, 42, 14, 15, 16, 29, 30, 31, 32, , 46, 47, 48 70–77%

     | Show Table
    DownLoad: CSV

    In this study, the accuracy rates of eight areas of the prefrontal cortex and of the entire brain were determined. Herein, the accuracy rates of all channels in terms of true-positive, true-negative, false-positive, and false-negative results are given as examples, with methamphetamine abusers labeled as 0, heroin abusers as 1, and mixed abusers as 2. Among the three categories of participants, the accuracy rate of identifying methamphetamine abuse is lower than that of heroin and mixed drugs. Further adjustments to the model can increase the uniformity of the accuracy rate. Figure 7, Confusion matrix of three types of drug abusers.

    Figure 7.  Confusion matrix of three types of drug abusers.

    According to the unified CNN framework, there are 8 functional areas in the prefrontal lobe of the brain, and the correct rate statistics are shown in Figure 8. According to the statistical results, the following conclusions are drawn that the accuracy rate obtained in Left-DLPFC is lower than the Right-DLPFC prefrontal cortex. However, In VLPFC, FPC, OFC functional areas, the accuracy of the left hemisphere of the brain is higher than that of the right hemisphere.

    Figure 8.  The accuracy of the forehead area.

    Heroin, methamphetamine, and mixed drug abusers, through one-way ANOVA analysis revealed three types of personnel differences. The analysis results of Table 4 are as follows: Right-VLPFC, Right-OFC, there are significant differences between the three types of personnel. The differences in other functional areas are not very obvious.

    Table 4.  statistical analysis of drug abusers in different functional areas.
    Heroin addicts Meth addicts Mixed drug addicts Left-DLPFC ANOVA summary F P-value R- square (P < 0.05)
    1.467 0.2327 0.01223 No
    ANOVA table SS DF MS F (DFn, DFd) P value
    Treatment (between columns) 0.001373 2 0.0006866 F (2,237) = 1.467 P = 0.2327
    Residual (within columns) 0.1109 237 0.0004680
    Total 0.1123 239
    Right-DLPFC ANOVA summary F P-value R- square (P < 0.05)
    2.794 0.0660 0.02267 No
    ANOVA table SS DF MS F (DFn, DFd) P value
    Treatment (between columns) 0.001266 2 0.0006332 F (2,237) = 2.749 P = 0.0660
    Residual (within columns) 0.05459 237 0.0002303 0.05459
    Total 0.05586 239
    Left-VLPFC ANOVA summary F P-value R- square (P < 0.05)
    0.2203 0.8032 0.01038 No
    ANOVA table SS DF MS F (DFn, DFd) P value
    Treatment (between columns) 0.001820 2 0.0009100 F (2, 42) = 0.2203 P = 0.8032
    Residual (within columns) 0.1735 42 0.004130
    Total 0.1753 44
    Right-VLPFC ANOVA summary F P-value R- square (P < 0.05)
    4.068 0.0243 0.1623 Yes
    ANOVA table SS DF MS F (DFn, DFd) P value
    Treatment (between columns) 0.000564 2 0.0002820 F (2, 42) = 4.068 P = 0.0243
    Residual (within columns) 0.002912 42 6.933e-005
    Total 0.003476 44
    Left-FPC ANOVA summary F P-value R- square (P < 0.05)
    1.492 0.2293 0.02486 No
    ANOVA table SS DF MS F (DFn, DFd) P value
    Treatment (between columns) 0.001555 2 0.0007773 F (2,117) = 1.492 P = 0.2293
    Residual (within columns) 0.06097 117 0.0005211
    Total 0.06252 119
    Right-FPC ANOVA summary F P-value R- square (P < 0.05)
    2.762 0.0657 0.02839 No
    ANOVA table SS DF MS F (DFn, DFd) P value
    Treatment (between columns) 0.000312 2 0.0001562 F (2,189) = 2.762 P = 0.0657
    Residual (within columns) 0.01069 189 5.655e-005
    Total 0.01100 191
    Left-OFC ANOVA summary F P-value R- square (P < 0.05)
    1.422 0.2446 0.01879 No
    ANOVA table SS DF MS F (DFn, DFd) P value
    Treatment (between columns) 0.000153 2 7.679e-005 F (2,147) = 1.422 P = 0.2446
    Residual (within columns) 0.007941 147 5.402e-005
    Total 0.008094 149
    Right-OFC ANOVA summary F P-value R- square (P < 0.05)
    6.930 0.0013 0.08780 Yes
    ANOVA table SS DF MS F (DFn, DFd) P value
    Treatment (between columns) 0.000685 2 0.0003425 F (2,144) = 6.930 P = 0.0013
    Residual (within columns) 0.007118 144 4.943e-005
    Total 0.007803 146

     | Show Table
    DownLoad: CSV

    The following results can be obtained:

    1) Long-term drug abuse of drug users has caused the function of the brain to become stronger or weaker.

    2) Long-term drug abuse by drug users on Right-VLPFC and Right-OFC can infer the functional damage of the cortical area.

    Due to the different pathogenic mechanisms of drugs, drug users who abuse different types of drugs suffer from different brain damages. In this paper, the blood oxygen levels of 30 drug users were obtained by high-density functional near-infrared spectroscopy. The eight Brodmann areas and the entire prefrontal area of the three types of individuals were identified and classified by the designed CNN to obtain the final results: the 3-class accuracy was lower on the left side than on the right in the DLPFC, while it was higher on the left than on the right in the VLPFC, FPC, and OFC. Statistical analysis on the eight functional areas showed significant differences in the right VLPFC and right OFC, and preliminarily predictions were for the brain damages caused by abusing different drugs.

    The authors thank the Shanghai Drug Rehabilitation Administration Bureau and Shanghai Qing Dong Drug Rehabilitation Center for their support, and all the people who contribute to the writing of the article. This project was supported by Key Research & Development Project of National Science and Technique Ministry of China (No. 2018YFC0807405, No. 2018YFC1312903), National Natural Science Foundation of China (No. 61976133). Defense Industrial Technology Development Program (JCKY2019413D002) and the Key Research and Development Project of Shandong (2019JZZY021010).

    The authors declare that they have no conflict of interest.



    [1] C. Bulgarelli, C. C. J. M. de Klerk, J. E. Richards, V. Southgate, A. Hamilton, A. Blasi, The developmental trajectory of fronto‐temporoparietal connectivity as a proxy of the default mode network: a longitudinal fNIRS investigation, Hum. Brain Mapp., 41 (2020), 2717-2740. doi: 10.1002/hbm.24974
    [2] A. S. Huhn, R. K. Brooner, M. M. Sweeney, S. W. Yip, K. E. Dunn, Increased neural activity in the right dorsolateral prefrontal cortex during a risky decision-making task is associated with cocaine use in methadone-maintained patients, Drug Alcohol. Depend., 205 (2019), 107650.
    [3] L. Huang, S. Y. Guo, Y. Wang, S. Wang, Q. B. Chu, L. Li, et al, Attention based residual network for medicinal fungi near infrared spectroscopy analysis, Math. Biosci. Eng., 16 (2019), 3003-3017.
    [4] L. Cai, T. Nitta, S. Yokota, T. Obata, E. Okada, Targeting brain regions of interest in functional near-infrared spectroscopy-Scalp-cortex correlation using subject-specific light propagation models, Hum. Brain Mapp., 42 (2021), 1969-1986. doi: 10.1002/hbm.25367
    [5] K. Saita, T. Morishita, K. Hyakutake, H. Fukuda, E. Shiota, Y. Sankai, et al, Combined therapy using botulinum toxin A and single joint hybrid assistive limb for upper-limb disability due to spastic hemiplegia, J. Neurol. Sci., 373 (2017), 182-187.
    [6] H. Song, L. Chen, R. Q. Gao, I. I. M. Bogdan, J. Yang, S. L. Wang, et al, Automatic schizophrenic discrimination on fNIRS by using complex brain network analysis and SVM, BMC Med. Inf. Decis., 17 (2017), 166.
    [7] B. Z, Hu, B. Gao, W. L. Woo, L. F. Ruan, J. K. Jin, Y. Yang, et al, A lightweight spatial and temporal multi-feature fusion network for defect detection, IEEE Trans. Image Process., 30 (2021), 472-486.
    [8] H. Y. Kim, S. P. Yang, G. L. Park, E. J. Kim, J. H. You, Best facilitated cortical activation during different stepping, treadmill, and robot-assisted walking training paradigms and speeds: A functional near-infrared spectroscopy neuroimaging study, NeuroRehabilitation, 38 (2016), 171-178.
    [9] B. H. D. Koh, C. L. P. Lim, H. Rahimi, W. L. Woo, B. Gao, Deep temporal convolution network for time series classification, Sensors, 21 (2021), 603.
    [10] R. Holtzer, R. Kraut, M. Izzetoglu, K. Ye, The effect of fear of falling on prefrontal cortex activation and efficiency during walking in older adults, GeroScience, 41 (2019), 89-100. doi: 10.1007/s11357-019-00056-4
    [11] P. Caliandro, M. Serrao, L. Padua, G. Silvestri, C. Iacovelli, C. Simbolotti, Prefrontal cortex as a compensatory network in ataxic gait: A correlation study between cortical activity and gait parameters, Restor. Neurol. Neurosci., 33 (2015), 177-187.
    [12] I. Rodriguez-Rodriguez, J. V. Rodriguez, W. L. Woo, B. Wei, D. J. Pardo-Quiles, A comparison of feature selection and forecasting machine learning algorithms for predicting glycaemia in type 1 diabetes mellitus, Comput. Inf. Sci., 11 (2021), 1742.
    [13] T. Liu, X. C. Liu, L. Yi, C. Z. Zhu, P. S. Markey, M. Pelowski, Assessing autism at its social and developmental roots: A review of autism spectrum disorder studies using functional near-infrared spectroscopy, Neuroimage, 185 (2019), 955-967. doi: 10.1016/j.neuroimage.2017.09.044
    [14] J. Hennrich, C. Herff, D. Heger, T. Schultz, Investigating deep learning for fNIRS based BCI, in Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), (2015), 2844-2847.
    [15] B. Abibullaev, J. An, J. I. Moon, Neural network classification of brain hemodynamic responses from four mental tasks, Int. J. Optomechatroni., 5 (2011), 340-359. doi: 10.1080/15599612.2011.633209
    [16] A. M. Atula, J. A. Mark, Y. E. Kim, H. Ayaz, Comparison of brain activation during motor imagery and motor movement using fNIRS, Comput. Intel. Neurosci., 2017 (2017), 5491296.
    [17] T. Akiyama, M. Koeda, Y. Okubo, M. Kimurab, Hypofunction of left dorsolateral prefrontal cortex in depression during verbal fluency task: A multichannel near-infrared spectroscopy study, J. Affect. Disord., 231 (2018), 83-90. doi: 10.1016/j.jad.2018.01.010
    [18] C. Issard, J. Gervain, Variability of the hemodynamic response in infants: Influence of experimental design and stimulus complexity, Dev. Cogn. Neurosci., 33 (2018), 182-193. doi: 10.1016/j.dcn.2018.01.009
    [19] A. C. May, R. L. Aupperle, J. L. Stewart, Dark times: The role of negative reinforcement in methamphetamine addiction, Front. Psychiatry, 11 (2020), 114.
    [20] K. Peng, M. A. Yucel, S. C. Steele, E. A. Bittner, C. M. Aasted, M. A. Hoeft, Morphine attenuates fNIRS signal associated with painful stimuli in the medial frontopolar cortex (medial BA 10), Front. Hum. Neurosci., 12 (2018), 394.
    [21] Y. L. Tseng, C. F. Lu, S. M. Wu, S. Shimada, T. Huang, G. Y. Lu, A functional near-infrared spectroscopy study of state anxiety and auditory working memory load, Front. Hum. Neurosci., 12 (2018), 313.
    [22] K. L. M. Koenraadt, E. G. J. Roelofsen, J. Duysens, N. L. W. Keijsers, Cortical control of normal gait and precision stepping: An fNIRS study, Neuroimage, 85 (2014), 415-422. doi: 10.1016/j.neuroimage.2013.04.070
    [23] A. C. Ehlis, S. Schneider, T. Dresler, A. J. Fallgatter, Application of functional near-infrared spectroscopy in psychiatry, Neuroimage, 85 (2014), 478-488. doi: 10.1016/j.neuroimage.2013.03.067
    [24] E. Jeong, H. Ryu, J. H. Shin, G. H. Kwon, G. Jo, J. Y. Lee, High oxygen exchange to music indicates auditory distractibility in acquired brain injury: An fNIRS study with a vector-based phase analysis, Sci. Rep., 8 (2018), 16767.
    [25] J. M. Kim, J. K. Choi, M. Choi, M. Ji, G. Hwang, S. B. Ko, Assessment of cerebral autoregulation using continuous-wave near-infrared spectroscopy during squat-stand maneuvers in subjects with symptoms of orthostatic intolerance, Sci. Rep., 8 (2018), 13257.
    [26] J. Shin, J. Kwon, J. Choi, C. H. Im, Performance enhancement of a brain-computer interface using high-density multi-distance NIRS, Sci. Rep., 7 (2017), 16545.
    [27] B. R. Ford, G. Bart, B. Grahan, R. D. Shearer, T. N. A. Winkelman, Associations between polysubstance use patterns and receipt of medications for opioid use disorder among adults in treatment for opioid use disorder, J. Addict. Med., 15 (2021), 159-162. doi: 10.1097/ADM.0000000000000726
    [28] J. C. Strickland, J. R. Havens, W. W. Stoops, A nationally representative analysis of "twin epidemics": Rising rates of methamphetamine use among persons who use opioids, Drug. Alcohol. Depen., 204 (2019), 107592.
    [29] A. Fugelstad, S. Bremberg, P. Hjelmstrom, I. Thiblin, Methadone-related deaths among youth and young adults in Sweden 2006-15, Addiction, 116 (2020), 319-327.
    [30] X. C. Duan, X. Z. Li, M. Martcheva, Dynamics of an age-structured heroin transmission model with vaccination and treatment, Math. Biosci. Eng., 16 (2019), 397-420. doi: 10.3934/mbe.2019019
    [31] R. van de Wetering, S. Schenk, Regional changes in increment FosB expression in rat brain following MDMA self-administration predict increased sensitivity to effects of locally infused MDMA, Addict. Biol., 25 (2019), e12814.
    [32] O. E. Hall, O. T. Hall, J. L. Eadie, J. Teater, J. Gay, M. Kim, Street-drug lethality index: A novel methodology for predicting unintentional drug overdose fatalities in population research, Drug. Alcohol. Depen., 221 (2021), 10836.
    [33] S. H. Ma, H. F. Huo, Global dynamics for a multi-group alcoholism model with public health education and alcoholism age, Math. Biosci. Eng., 16 (2019), 1683-1708. doi: 10.3934/mbe.2019080
  • This article has been cited by:

    1. Yun-Hsuan Chen, Jie Yang, Hemmings Wu, Kevin T. Beier, Mohamad Sawan, Challenges and future trends in wearable closed-loop neuromodulation to efficiently treat methamphetamine addiction, 2023, 14, 1664-0640, 10.3389/fpsyt.2023.1085036
    2. Yuxing Zhou, Xuelin Gu, Zhen Wang, Xiaoou Li, Identification of drug use degree by integrating multi-modal features with dual-input deep learning method, 2024, 1025-5842, 1, 10.1080/10255842.2024.2417206
  • Reader Comments
  • © 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(3070) PDF downloads(147) Cited by(2)

Figures and Tables

Figures(8)  /  Tables(4)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog