Skip to main content

ORIGINAL RESEARCH article

Front. Psychiatry, 04 February 2021
Sec. Neuroimaging
This article is part of the Research Topic Neurobiological Biomarkers for Developing Novel Treatments of Substance and Non-Substance Addiction View all 45 articles

Characterizing the Structural Pattern of Heavy Smokers Using Multivoxel Pattern Analysis

\nYufeng Ye,&#x;Yufeng Ye1,2Jian Zhang&#x;Jian Zhang3Bingsheng Huang,,Bingsheng Huang1,4,5Xun CaiXun Cai5Panying Wang,,Panying Wang6,7,8Ping Zeng,Ping Zeng5,6Songxiong Wu,Songxiong Wu5,6Jinting MaJinting Ma5Han HuangHan Huang5Heng LiuHeng Liu9Guo Dan
Guo Dan10*Guangyao Wu,,
Guangyao Wu6,7,8*
  • 1Department of Radiology, Panyu Central Hospital, Guangzhou, China
  • 2Medical Imaging Institute of Panyu, Guangzhou, China
  • 3Health Science Center, Shenzhen University, Shenzhen, China
  • 4Shenzhen University Clinical Research Center for Neurological Diseases, Shenzhen, China
  • 5Medical AI Lab, School of Biomedical Engineering, Health Science Center, Shenzhen University, Shenzhen, China
  • 6Department of Radiology, Shenzhen University General Hospital, Shenzhen, China
  • 7Shenzhen University International Cancer Center, Shenzhen, China
  • 8Zhongnan Hospital of Wuhan University, Wuhan, China
  • 9Medical Imaging Center of Guizhou Province, Department of Radiology, The Affiliated Hospital of Zunyi Medical University, Zunyi, China
  • 10School of Biomedical Engineering, Health Science Center, Shenzhen University, Shenzhen, China

Background: Smoking addiction is a major public health issue which causes a series of chronic diseases and mortalities worldwide. We aimed to explore the most discriminative gray matter regions between heavy smokers and healthy controls with a data-driven multivoxel pattern analysis technique, and to explore the methodological differences between multivoxel pattern analysis and voxel-based morphometry.

Methods: Traditional voxel-based morphometry has continuously contributed to finding smoking addiction-related regions on structural magnetic resonance imaging. However, voxel-based morphometry has its inherent limitations. In this study, a multivoxel pattern analysis using a searchlight algorithm and support vector machine was applied on structural magnetic resonance imaging to identify the spatial pattern of gray matter volume in heavy smokers.

Results: Our proposed method yielded a voxel-wise accuracy of at least 81% for classifying heavy smokers from healthy controls. The identified regions were primarily located at the temporal cortex and prefrontal cortex, occipital cortex, thalamus (bilateral), insula (left), anterior and median cingulate gyri, and precuneus (left).

Conclusions: Our results suggested that several regions, which were seldomly reported in voxel-based morphometry analysis, might be latently correlated with smoking addiction. Such findings might provide insights for understanding the mechanism of chronic smoking and the creation of effective cessation treatment. Multivoxel pattern analysis can be efficient in locating brain discriminative regions which were neglected by voxel-based morphometry.

Introduction

Tobacco smoking in the form of cigarettes continues to be the leading cause of preventable illness and mortality in the world (1). China, the largest producer of tobacco, is estimated to contain 311 million individuals who are current smokers, with 295 million men and 16 million women, respectively (2). Chronic smoking is known to correlate with a series of diseases including stroke, lung cancer, hepatocellular carcinoma, and vascular dysfunctions (37). Related neuroimaging studies also suggested that the numerous toxic chemicals contained in a cigarette, especially nicotine, could promote potential brain afflictions in chronic cigarette smokers (8, 9).

Apart from these serious public health problems, 78% of smokers who expressed willingness to quit smoking reported a relapse situation in China, and the percentage in America is currently 80% (10, 11), indicating the ineffectiveness of existing cessation treatments. Poor treatment outcome for smoking cessation may result from a lack of awareness of the mechanism behind smoking addiction and the available biomarkers that characterize heavy smokers (12, 13).

Structural magnetic resonance imaging (sMRI), which visualizes the central neural system with high-resolution in a non-invasive way, provides indispensable spatial information in the procedure of identifying such biomarkers (i.e., discriminative brain regions). For the past decade, univariate approaches such as voxel-based morphometry (VBM) have continuously been applied on sMRI studies (14, 15). In these studies, morphological abnormalities are consistently discovered in heavy smokers in some brain regions, including the prefrontal cortex, anterior cingulate cortex, thalamus, and the insula (9, 1418).

Despite the accumulative results reported, VBM, as a traditional univariate approach, has its inherent limitations in identifying spatial patterns that exist in a certain population (e.g., a group of heavy smokers). Generally, VBM has the multiple comparisons problem (19), which causes a loss of sensitivity and overlooks the dependency of the focal set of voxels in localizing informative regions relevant to specific brain abnormalities (20). Multi-voxel pattern analysis (MVPA), which utilizes a multi-variate technique and is driven by machine learning algorithms, provides a sensitive and different approach in identifying group-wise differences. MVPA takes multiple voxels into account and considers patterns across a group of voxels that may respond weakly but consistently differently between conditions (21, 22). Therefore, MVPA can be sensitive in distinguishing different experimental conditions. MVPA has become favorable in neuroimaging research for its ability to detect subtle anatomical differences (2325).

The field of utilizing MVPA to discover structural abnormalities in heavy smokers remains less active. Notably, in an sMRI study based on support vector machine (SVM), the authors (26) used the mean gray matter volumes (GMVs) of 1,024 self-defined brain regions as input features to the SVM to identify the most discriminative regions by finding the most informative features in the SVM. Using average brain region GMVs as features may be an over simplistic way to characterize the anatomical structure of the neural system. The subtle difference within the same self-defined region may be overlooked. In this study, in order to discover smaller regions that reflect the experimental conditions of interest, we used a searchlight algorithm by moving a searchlight window through the volume of the brain, to sample the gray matter values of voxels as the input features to our linear classifier (SVM). A searchlight algorithm is able to preserve subtle differences between a group of voxels while a linear SVM can detect such subtle differences by efficiently defining a boundary that maximally separates two classes (e.g.,: heavy smokers and healthy non-smokers) of samples.

With regard to the peer studies on smoking addiction and the limitations of the current univariate approach, in this study, we aimed to: (1) explore the most discriminative gray matter regions between heavy smokers and otherwise healthy controls (HCs) with a data-driven MVPA technique and sMRI data; and (2) explore the methodological difference between MVPA and VBM. We hypothesized that MVPA would be a more sensitive method in locating group-wise brain differences than VBM.

Methods

Participants

The study adhered to the Declaration of Helsinki and was approved by the Medical Ethics Review Board of Zhongnan Hospital of Wuhan University. After a complete description of the study, written informed consent was obtained from the subjects.

Sixty-eight right-handed subjects (39 heavy smokers and 29 healthy non-smoker control subjects) were recruited via advertisement flyers and referrals for MRI studies. All subjects were screened for psychiatric and non-psychiatric medical disorders using the Mini International Neuropsychiatric Interview (27). All participants underwent an interview session followed by an inclusion procedure. The inclusion criteria for both groups were: the absence of any history of medical (e.g., cardiac disease) or neurological (e.g., stroke) disorders, intellectual disability, drug abuse or dependence (other than nicotine dependence for heavy smokers), or psychiatric disease; none of the subjects reported the daily consumption of alcohol or experienced social consequences secondary to alcohol intake. Heavy smokers were defined as those who met the DSM-IV criteria for nicotine dependence, had smoked at least 20 cigarettes daily in the past 5 years, and had no period of smoking abstinence longer than 3 months. The severity of the heavy smokers' nicotine addiction was measured using the Fagerström Test for Nicotine Dependence (28). The healthy control subjects in this study were defined as those who had a history of smoking no more than five cigarettes in a lifetime.

Data from two heavy smokers and one non-smoker were excluded for not meeting the above criteria. As a result, 37 heavy smokers (male: 30; female: 7) and 28 healthy non-smoking control subjects (male: 21; female: 7) were included in this study. More detailed information on the subjects in each group is presented in Table 1. Demographic data were compared between the two groups by using two-sample t-tests or a Chi-square test in the Statistical Package for Social Science, version 19 (SPSS Inc., USA). The threshold level of significance was set as p < 0.05.

TABLE 1
www.frontiersin.org

Table 1. Demographic and clinical characteristics of the heavy smokers and HCs.

MR Acquisition

All participants underwent a high-resolution 3dn-MTC T1-weighted structural scan using a Siemens Trio 3.0-Tesla MR scanner (Erlangen, Germany) with a standard birdcage head coil. sMRI images were obtained using a MPRAGE pulse sequence with the following parameters: repetition time = 25 ms; echo time = 4.51 ms; flip angle = 25; acquisition matrix = 256 × 256; and slice thickness = 1 mm with a 1 mm gap.

Image Preprocessing

The raw DICOM images of sMRI were converted to the NIFTI format using MRIcron (University of South Carolina, Columbia, SC, USA, http://www.mricro.com). The following preprocessing steps were then performed to obtain gray matter (GM) maps using VBM8 toolbox (http://dbm.neuro.unijena.de/vbm) in SPM8 (Version 6313, Wellcome Department of Imaging Neuroscience, London, UK, http://www.fil.ion.ucl.ac.uk/spm) on MATLAB R2013a. Firstly, the sMRI images were registered to the Montreal Neurological Institute (MNI) stereotactic space and resampled to a 1.5 mm isotropic voxel spacing. Secondly, the co-registered images were segmented into three types of tissues, namely GM, white matter, and cerebrospinal fluid. Thirdly, a study-specific template was created using the high-dimensional Diffeomorphic Anatomical Registration Through Exponentiated Lie Algebra (DARTEL) algorithm and with the predefined templates in the VBM8 toolbox. Next, in order to preserve the total volume of each brain tissue, the segmented images were modulated using non-linear deformation which can compensate for the effect of spatial normalization. This step multiplied the spatially normalized gray matter by its relative volume before and after spatial normalization (29). At last, the images were smoothed with an 8-mm full-width-half-maximum Gaussian kernel.

MVPA

An MVPA technique combining a searchlight algorithm and a linear SVM was used to classify prominent regions that distinguished heavy smokers from HCs. Generally, by moving a searchlight region through the brain volume, one can continuously map the information content regarding the experimental conditions of interest in the brain (30). The procedure of our MVPA method is as follows. The smoothed GM maps computed from the data preprocessing step were used as inputs in the MVPA. Firstly, at each voxel (Vi) of the GM images in the normalized space, a three-dimensional cubical region size of 3 × 3 × 3 centering at Vi was identified. For each subject, the gray matter volume values of all 27 voxels (at a specific voxel position) were extracted and converted into a high-dimensional vector and used to construct the feature matrix. To train and test the SVM, a leave-one-out (LOO) cross-validation (CV) strategy was adopted, which excludes one subject as a testing set each time and trains the classifier using the remaining subjects. As a result, two feature matrices MF1*S and MF2*S were retrieved as the training set and testing set (F1 = 64, F2 = 1, S = 27), respectively. F1 and F2 indicate the number of subjects in the two sets respectively, and S indicates the number of voxels of the feature matrix. Next, the training set was fed into the SVM implemented in the LIBSVM toolbox (http://www.csie.ntu.edu.tw/~cjlin/Libsvm). In each training set, a nested 5-fold CV was applied to determine the optimized parameter C (regularization) and g (gamma for radius basis function) for testing. That is, each time one-fifth of the training set was selected as a testing sample and the classifier was built upon the remaining data. Parameters that produced the highest accuracy across these 5-folds was identified as the optimized C and g. The identified C and g were then used in the corresponding testing procedure. Finally, the accuracy of the trained classifier was assigned to the chosen voxel Vi. After repeating this procedure on every voxel, a three-dimensional accuracy map denoting the classification ability between heavy smokers and HCs for all voxels was obtained.

To evaluate the statistical significance of the experimental results, we converted the accuracy map to a p-value map under an assumption of binomial distribution. Detailed information about the conversion procedure can be found in another publication (25). A connected domain algorithm was conducted on the p-value map to produce clusters with significant predictive power, whose threshold was set at p < 0.001 with more than 50 adjacent voxels (24).

Post hoc Analysis of MVPA Analysis Results

To further explore the difference of GMVs between heavy smokers and HCs, a post hoc analysis was carried out, with which a further VBM analysis was conducted within the brain clusters detected by MVPA. The analysis was corrected for multiple comparisons using family wise error (FWE) at the cluster level (p < 0.05). To do this, firstly the GMVs of the brain regions detected by the MVPA were extracted using the MarsBar toolbox (http://www.mrc-cbu.cam.ac.uk/Imaging/marsbar.html). Then, with age and gender as co-variates, voxel-wise two independent samples t-tests were further performed in these regions to determine the significant differences between heavy smokers and HCs (p < 0.05, FWE-corrected, two-tailed).

Whole Brain VBM Analysis

In order to compare the methodological differences between MVPA and VBM, a whole brain VBM analysis was conducted. With age and gender as co-variates, the whole brain GM maps of heavy smokers and HCs were fed to VBM. The same threshold values of post hoc analysis were applied in order to conduct a comparison between whole brain VBM analysis and MVPA. Discriminative regions discovered by the whole brain VBM were defined as clusters with a cluster size > 50 voxels and a p-value < 0.001.

ROC Analysis

To evaluate the ability of MVPA or VBM in distinguishing the smokers from the healthy controls, we performed ROC analysis based on the brain regions with significant differences. Firstly, we used the average gray matter values of brain regions as features to establish an SVM model using 5-fold cross validation. Then the ROC analysis was performed with the predicted probability of each participant in the procedure of cross validation. The area under the ROC curve (AUC) denotes the ability of classifying heavy smokers from healthy controls.

Results

Demographic Characteristics

Heavy smokers did not differ significantly with HCs in terms of age (p = 0.96), years of education (p = 0.36), or gender ratio (p = 0.77). For heavy smokers, the average FTND score was 8.89 ± 0.68 (range, 8–10; median, 9) and the average number cigarettes per day was 35.13 ± 10.70, which indicated a relatively high dependence on cigarettes. Detailed demographic information for both groups can be found in Table 1.

MVPA and post hoc Analysis

The discriminative regions recognized by our MVPA technique without covariates regression before classifier training are shown in Figure 1. Our proposed technique yielded a voxel-wise accuracy of at least 81% for classifying heavy smokers from HCs. Several cortical and subcortical regions demonstrated a strong classification ability with GM differences between heavy smokers and HCs. These regions were primarily located at the temporal cortex and prefrontal cortex, occipital cortex, thalamus (bilateral), insula (left), anterior and median cingulate gyri (ACG and MCG), and precuneus (left).

FIGURE 1
www.frontiersin.org

Figure 1. Brain regions with high classification accuracy identified by MVPA. The color bar indicates the classification accuracy values of the whole brain GM voxels. The image is displayed in the neurologic convention, with the left side corresponding to the left brain.

Two-sample t-tests (post hoc analysis) of GMVs revealed that 16 out of the 18 brain regions identified by MVPA demonstrated a significant GMV decrease in heavy smokers as compared to healthy controls (p < 0.05, FWE-corrected, two-tailed). Notably, the inferior temporal gyrus (left) and cerebellum were found with a significantly increased GMV in heavy smokers than in HCs. The peak accuracy values of these clusters and the corresponding t-values which were derived from the post hoc analysis are reported in Table 2. The peak t-value was defined as the maximum t-value found in the extracted cluster. A positive t-value indicated that the GM volume in HS was significantly smaller than HC in a specific region while a negative value indicated the opposite. The brain regions recognized by VBM are shown in Figure 2.

TABLE 2
www.frontiersin.org

Table 2. Brain regions with high predictive accuracy identified by MVPA and the corresponding post-hoc (VBM) results, which was corrected for multiple comparisons using family wise error (FWE) at the cluster level (p < 0.05).

FIGURE 2
www.frontiersin.org

Figure 2. Brain regions identified by whole brain VBM analysis. The color bar indicates the peak t-values of the maximum voxel found in the extracted cluster. A positive t-value indicates that the GM volume in HS is significantly smaller than HC in a specific region while a negative value indicates the opposite. The image is displayed in the neurologic convention, with the left side corresponding to the left brain.

The ROC analysis results are shown in Figure 3. The ROC results indicate that MVPA outperformed VBM in differentiating the heavy smokers from the HCs (AUC, 0.81 vs. 0.78; sensitivity, 0.82 vs. 0.79; specificity, 0.77 vs. 0.77; positive predictive value (PPV), 0.72 vs. 0.71; negative predictive value (NPV), 0.85 vs. 0.82).

FIGURE 3
www.frontiersin.org

Figure 3. The ROC analysis results of MVPA (left) and VBM (right). The results (AUC, 0.81 vs. 0.78; sensitivity, 0.82 vs. 0.79; specificity, 0.77 vs. 0.77; PPV, 0.72 vs. 0.71; NPV, 0.85 vs. 0.82) indicate that MVPA outperformed VBM in differentiating the smokers from healthy controls by using the GMVs of the detected brain regions as the features.

VBM Analysis

As summarized in Table 3, whole brain VBM analysis revealed that six brain regions demonstrated a significant decrease or increase of GMV in heavy smokers than in HC. These regions were located in the occipital, precuneus lingual, and cerebrum. The peak t-value was defined as the maximum t-value found in the extracted cluster. A positive t-value indicated that the GM volume in HS was significantly smaller than HC in a specific region while a negative value indicated the opposite.

TABLE 3
www.frontiersin.org

Table 3. Brain regions detected by VBM with a significant decrease or increase GMV in heavy smokers than HC.

Discussion

Traditional univariate studies prior to this research have constantly discovered that heavy smokers shared similar GM alterations in the cingulum, thalamus, cerebellum, prefrontal gyrus, and precuneus (9, 1518, 31). In this study, these regions were identified by MVPA as discriminative regions between heavy smokers and HCs. In a following post hoc analysis on these particular regions, our analysis results also demonstrated significant GM decrease in several areas including the temporal cortex, prefrontal cortex, thalamus, ACG and MCG, and precuneus in heavy smokers, indicating a strong correlation between GM alteration in regional GM deficiency and chronic smoking. In addition, our MVPA result also revealed that regions (insula, cerebellum) which were seldomly reported in traditional VBM analysis, may participate in the chronic smoking mechanism. The resulting contrast between MVPA and VBM may indicate the methodological differences between these two methods.

In our study, the bilateral thalamus were identified by MVPA with high accuracy. Particularly, in the post hoc analysis, the right thalamus was found to have significant GMV loss in heavy smokers. The observation of GMV deficiency in heavy smokers may primely result in the cognitive impairment reported in several studies (8, 17), as the thalamus relays information between the cerebral and different subcortical cortexes (32). In addition to cognitive impairment, the thalamus has been identified as a brain region with a relatively high density of nicotine acetycholine receptors (α4β2* nAChRs) (23, 33, 34). Frequent nicotine binding activities could be the leading reason for the GMV decrease in the thalamus. Such a finding also provides insight for formulating more effective withdrawal treatment. In a related Positron Emission Computed Tomography (PET) study (33), the author suggested that maintaining α4β2* nAChRs in the desensitization state may be a possible way to alleviate withdrawal symptoms.

The cingulum, bidirectionally connected with the medial temporal lobes, was another region widely reported to be associated with smoking addiction in univariate studies. In this study, the bilateral cingulum could discriminate heavy smokers from HCs with at least an accuracy of 75%. In the corresponding post hoc analysis, both the ACG and MCG demonstrated significant GM loss in heavy smokers. As a region connected to sites repeatedly implicated in cognitive control, the ACG is involved in executive function behaviors, such as inhibitory control and conflict resolution (35), which may play an important part in the process of quitting smoking.

The prefrontal area was another discriminative region identified by MVPA with high accuracy. Especially, in the post hoc analysis, the inferior and middle areas were found with significant GMV loss in heavy smokers as compared to HCs. As the prefrontal cortex is associated with concentration, emotion, and other higher brain functions, our findings with the prefrontal cortex may partially explain the decreased attention and impairment of working memory reported in heavy smokers (3639).

Among other regions detected by MVPA, the insula presented relatively high predictive accuracy and stood out as a discriminative region. Acting as a critical neural substrate for addiction including nicotine dependence (40), the insula plays a potential role in the decision-making task that is associated with relapse to drug use (41) and cue-induced drug urges (14, 42, 43). However, the insula was not reported in our whole brain VBM analysis nor in other univariate sMRI studies. In contrast, our finding of the insula was in accordance with another multi-variate sMRI study deploying SVM (26). The fact that the insula along with the ACG are part of the salience network (SN) (44), a network that mediates one's subjective feelings, is thought to be one of the solutions to developing effective treatment plans. A previous fMRI study discovered that damage to the insula could disrupt addiction to cigarette smoking which in turn may help improve smoking cessation outcomes (40). Besides, evidence of the insula's role in mediating addiction urges may partially justify why the insula stands out as a discriminative region.

Notably, nicotine binding activity in the cerebellum, a region that was identified by our MVPA method with high accuracy (82%) has also been reported in several biochemical studies (45, 46). Heavy smokers showed a greater density of nicotine binding in the cerebellum as compared to non-smoking individuals (47), and brain blood flow in the cerebellum was increased by cigarette smoking (45, 46). This evidence lends support to the fact that decreased cerebellar GMV is associated with smoking addiction. However, like the insula and other regions, significantly lower GMV was not found in the cerebellum in heavy smokers more than in healthy controls in the whole brain VBM analysis. Moreover, as indicated by sMRI studies deploying ML techniques (26, 48), the insula and cerebellum were found to participate in the mechanism of chronic smoking. Such differences may be attributed to the fundamental difference between MVPA and traditional univariate analysis.

To be specific, MVPA fixes what was considered a disadvantage in traditional univariate methods. Conventional univariate analysis methods try to separately find voxels that show a statistically significant response to the experimental/physical conditions by deriving the average of the chosen voxels in all subjects, which assumes that the covariance across neighboring voxels is not informative about the experimental conditions under examination (49). Such assumption leads to the following inadequacies. Firstly, low response voxels that still carry important information could therefore be dismissed. Secondly, a consistent spatial pattern is neglected by simply averaging all voxels. MVPA uses a different solution to boost accuracy. Instead of uniformly averaging all voxels in VBM, MVPA assigns a weighted average in different conditions. This operation allows MVPA to discover spatial patterns neglected by VBM. Further, MVPA tries to optimize these weights by involving data-driven machine learning techniques. The evidence that the regions showing correlation to chronic smoking also implied that MVPA could be more sensitive than conventional univariate methods. In this view, MVPA can be better than VBM in discovering specific patterns in the chronic smoking population.

The study still presents the following limitations. Firstly, the sample size in the current study is relatively small. In a future study, more data should be collected in order to construct a more robust classifier and to verify our findings. Secondly, a more comprehensive dataset combining different modalities, such as fMRI and diffusion tensor imaging, would be helpful in exploring the mechanism of discriminative ability of such regions in heavy smokers.

In this study, MVPA was applied to structural MRI to identify the brain regions in discriminating heavy smokers from healthy individuals. The anatomical deficiency in heavy smokers was mainly discovered in the insula, ACG, MCG, prefrontal cortex, precuneus, and cerebellum, which are highly involved in the addiction of chronic smoking. These findings in accordance with previous results (45, 47, 48, 50) might provide insights for understanding the mechanism of chronic smoking and effective cessation treatment. Such insights also indicate the potential of using MVPA in future neuroimaging research. Moreover, the comparison between VBM and MVPA revealed that MVPA can be efficient in locating brain discriminative regions which were neglected by VBM.

Data Availability Statement

The data analyzed in this study is subject to the following licenses/restrictions: The data for this study are not publicly available because Zhongnan Hospital of Wuhan University, the center from which the data were collected, does not agree to make the data publicly accessible. Requests to access these datasets should be directed to huangb@szu.edu.cn.

Ethics Statement

The studies involving human participants were reviewed and approved by Medical Ethics Review Board of Zhongnan Hospital of Wuhan University. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author Contributions

YY, JZ, GD, and GW: conception and design. GD and GW: administrative support. GW: provision of study materials or patients. YY, JZ, GD, GW, and SW: collection and assembly of data. YY, JZ, GW, and SW: data analysis and interpretation. All authors: manuscript writing. All authors final approval of manuscript.

Funding

This study was jointly funded by The National Natural Science Foundation of China (Nos. 61001057, 61973220, 61973220), National Key Program of China (Nos. (2016YFC1304702), The Shenzhen Science and Technology Project (JCYJ20200109114014533, JCYJ20190808175413552, JCYJ20180507182040213, and JCYJ20170818141853626), The Guangzhou Science and Technology Project (201804010111), Shenzhen University Top Ranking Project (860/000002100108), Shenzhen-Hong Kong Institute of Brain Science-Shenzhen Fundamental Research Institutions (2019SHIBS0003), and Shenzhen University Presidential Fund (No. 85706-0000040544).

Conflict of Interest

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

Acknowledgments

The authors want to thank the Zhongnan Hospital of Wuhan University and the participants who took part in the study and data collection.

References

1. Lin F, Wu G, Zhu L, Lei H. Heavy smokers show abnormal microstructural integrity in the anterior corpus callosum: a diffusion tensor imaging study with tract-based spatial statistics. Drug Alcohol Depend. (2013) 129:82–7. doi: 10.1016/j.drugalcdep.2012.09.013

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Li S, Meng L, Chiolero A, Ma C, Xi B. Trends in smoking prevalence and attributable mortality in China, 1991–2011. Preventive Med. (2016) 93:82–7. doi: 10.1016/j.ypmed.2016.09.027

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Kannel W, Higgins M. Smoking and hypertension as predictors of cardiovascular risk in population studies. J Hypertens Suppl. (1990) 8:S3–8.

PubMed Abstract | Google Scholar

4. Howard MO, Kivlahan D, Walker RD. Cloninger's tridimensional theory of personality and psychopathology: applications to substance use disorders. J Stud Alcohol. (1997) 58:48–66. doi: 10.15288/jsa.1997.58.48

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Longstreth W, Bernick C, Manolio TA, Bryan N, Jungreis CA, Price TR. Lacunar infarcts defined by magnetic resonance imaging of 3660 elderly people: the Cardiovascular Health Study. Arch Neurol. (1998) 55:1217–25. doi: 10.1001/archneur.55.9.1217

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Islami F, Torre LA, Jemal A. Global trends of lung cancer mortality and smoking prevalence. Transl Lung Cancer Res. (2015) 4:327. doi: 10.3978/j.issn.2218-6751.2015.08.04

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Abdel-Rahman O, Helbling D, Schöb O, Eltobgy M, Mohamed H, Schmidt J, et al. Cigarette smoking as a risk factor for the development of and mortality from hepatocellular carcinoma: An updated systematic review of 81 epidemiological studies. J Evid Based Med. (2017) 10:245–54. doi: 10.1111/jebm.12270

PubMed Abstract | CrossRef Full Text

8. Franklin TR, Wetherill RR, Jagannathan K, Johnson B, Mumma J, Hager N, et al. The effects of chronic cigarette smoking on gray matter volume: influence of sex. PloS ONE. (2014) 9:e104102. doi: 10.1371/journal.pone.0104102

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Fritz H-C, Wittfeld K, Schmidt CO, Domin M, Grabe HJ, Hegenscheid K, et al. Current smoking and reduced gray matter volume—a voxel-based morphometry study. Neuropsychopharmacology. (2014) 39:2594. doi: 10.1038/npp.2014.112

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Yang T, Mao A, Feng X, Jiang S, Wu D, Bottorff JL, et al. Smoking cessation in an urban population in China. Am J Health Behav. (2014) 38:933–41. doi: 10.5993/AJHB.38.6.15

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Wetherill RR, Rao H, Hager N, Wang J, Franklin TR, Fan Y. Classifying and characterizing nicotine use disorder with high accuracy using machine learning and resting-state fMRI. Addict Biol. (2018) 24:811–21. doi: 10.1111/adb.12644

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Sutherland MT, Mchugh MJ, Pariyadath V, Stein EA. Resting state functional connectivity in addiction: lessons learned and a road ahead. Neuroimage. (2012) 62:2281–95. doi: 10.1016/j.neuroimage.2012.01.117

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Fedota JR, Stein EA. Resting-state functional connectivity and nicotine addiction: prospects for biomarker development. Ann N Y Acad Sci. (2015) 1349:64–82. doi: 10.1111/nyas.12882

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Brody AL, Mandelkern MA, London ED, Childress AR, Lee GS, Bota RG, et al. Brain metabolic changes during cigarette craving. Arch Gen Psychiatry. (2002) 59:1162–72. doi: 10.1001/archpsyc.59.12.1162

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Liao Y, Tang J, Liu T, Chen X, Hao W. Differences between smokers and non-smokers in regional gray matter volumes: a voxel-based morphometry study. Addict Biol. (2012) 17:977–80. doi: 10.1111/j.1369-1600.2010.00250.x

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Brody AL, Mandelkern MA, Jarvik ME, Lee GS, Smith EC, Huang JC, et al. Differences between smokers and nonsmokers in regional gray matter volumes and densities. Biol Psychiatry. (2004) 55:77–84. doi: 10.1016/S0006-3223(03)00610-3

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Gallinat J, Meisenzahl E, Jacobsen LK, Kalus P, Bierbrauer J, Kienast T, et al. Smoking and structural brain deficits: a volumetric MR investigation. Eur J Neurosci. (2006) 24:1744–50. doi: 10.1111/j.1460-9568.2006.05050.x

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Zhang X, Salmeron BJ, Ross TJ, Geng X, Yang Y, Stein EA. Factors underlying prefrontal and insula structural alterations in smokers. Neuroimage. (2011) 54:42–8. doi: 10.1016/j.neuroimage.2010.08.008

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Adluru N, Hanlon BM, Lutz A, Lainhart JE, Alexander AL, Davidson RJ. Penalized likelihood phenotyping: unifying Voxelwise analyses and multi-voxel pattern analyses in neuroimaging. Neuroinformatics. (2013) 11:227–47. doi: 10.1007/s12021-012-9175-9

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Haynes JD, Rees G. Decoding mental states from brain activity in humans. Nat Rev Neurosci. (2006) 7:523–534. doi: 10.1038/nrn1931

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Haynes JD, Rees G. Predicting the stream of consciousness from activity in human visual cortex. Curr Biol. (2005) 15:1301–7. doi: 10.1016/j.cub.2005.06.026

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Norman KA, Polyn SM, Detre GJ, Haxby JV. Beyond mind-reading: multi-voxel pattern analysis of fMRI data. Trends Cogn Sci. (2006) 10:424–30. doi: 10.1016/j.tics.2006.07.005

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Paterson D, Nordberg A. Neuronal nicotinic receptors in the human brain. Prog Neurobiol. (2000) 61:75–111. doi: 10.1016/S0301-0082(99)00045-3

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Uddin LQ, Menon V, Young CB, Ryali S, Chen T, Khouzam A, et al. Multivariate searchlight classification of structural magnetic resonance imaging in children and adolescents with autism. Biol Psychiatry. (2011) 70:833–41. doi: 10.1016/j.biopsych.2011.07.014

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Zhang J, Cao W, Wang M, Wang N, Yao S, Huang B. Multivoxel pattern analysis of structural MRI in children and adolescents with conduct disorder. Brain Imaging Behav. (2018) 13:1–8. doi: 10.1007/s11682-018-9953-6

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Ding X, Yang Y, Stein EA, Ross TJ. Multivariate classification of smokers and nonsmokers using SVM-RFE on structural MRI images. Hum Brain Mapp. (2015) 36:4869–79. doi: 10.1002/hbm.22956

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Sheehan DV, Lecrubier Y, Sheehan KH, Amorim P, Janavs J, Weiller E, et al. The Mini-International Neuropsychiatric Interview (M.I.N.I.): the development and validation of a structured diagnostic psychiatric interview for DSM-IV and ICD-10. J Clin Psychiatry. (1998) 59(Suppl 20):22–33; quiz 34–57.

PubMed Abstract | Google Scholar

28. Heatherton TF, Kozlowski LT, Frecker RC, Fagerstrom KO. The Fagerstrom Test for Nicotine Dependence: a revision of the Fagerstrom Tolerance Questionnaire. Br J Addict. (1991) 86:1119–27. doi: 10.1111/j.1360-0443.1991.tb01879.x

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Mechelli A, Price CJ, Friston KJ, Ashburner J. Voxel-based morphometry of the human brain: methods and applications. Curr Med Imaging. (2005) 1:105–113. doi: 10.2174/1573405054038726

CrossRef Full Text | Google Scholar

30. Oosterhof NN, Wiestler T, Downing PE, Diedrichsen J. A comparison of volume-based and surface-based multi-voxel pattern analysis. Neuroimage. (2011) 56:593–600. doi: 10.1016/j.neuroimage.2010.04.270

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Almeida OP, Garrido GJ, Lautenschlager NT, Hulse GK, Jamrozik K, Flicker L. Smoking is associated with reduced cortical regional gray matter density in brain regions associated with incipient Alzheimer disease. Am J Geriatr Psychiatry. (2008) 16:92–8. doi: 10.1097/JGP.0b013e318157cad2

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Gazzaniga MS, Ivry RB, Mangun GR. Cognitive Neuroscience: The Biology of the Mind. New York, NY: W.W. Norton (2009).

Google Scholar

33. Brody AL, Mandelkern MA, London ED, Olmstead RE, Farahi J, Scheibal D, et al. Cigarette smoking saturates brain α4β2 nicotinic acetylcholine receptors. Arch Gen Psychiatry. (2006) 63:907–14. doi: 10.1001/archpsyc.63.8.907

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Picard F, Sadaghiani S, Leroy C, Courvoisier DS, Maroy R, Bottlaender M. High density of nicotinic receptors in the cingulo-insular network. Neuroimage. (2013) 79:42–51. doi: 10.1016/j.neuroimage.2013.04.074

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Volkow ND, Fowler JS, Wang GJ, Swanson JM. Dopamine in drug abuse and addiction: results from imaging studies and treatment implications. Mol Psychiatry. (2004) 9:557–69. doi: 10.1038/sj.mp.4001507

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Snyder FR, Davis FC, Henningfield JE. The tobacco withdrawal syndrome: performance decrements assessed on a computerized test battery. Drug Alcohol Depend. (1989) 23:259–66. doi: 10.1016/0376-8716(89)90090-2

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Ernst M, Heishman SJ, Spurgeon L, London ED. Smoking history and nicotine effects on cognitive performance. Neuropsychopharmacology. (2001) 25:313–9. doi: 10.1016/S0893-133X(01)00257-3

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Slawecki CJ, Ehlers CL. Lasting effects of adolescent nicotine exposure on the electroencephalogram, event related potentials, and locomotor activity in the rat. Brain Res Dev Brain Res. (2002) 138:15–25. doi: 10.1016/S0165-3806(02)00455-8

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Jacobsen LK, Krystal JH, Mencl WE, Westerveld M, Frost SJ, Pugh KR. Effects of smoking and smoking abstinence on cognition in adolescent tobacco smokers. Biol Psychiatry. (2005) 57:56–66. doi: 10.1016/j.biopsych.2004.10.022

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Naqvi NH, Rudrauf D, Damasio H, Bechara A. Damage to the insula disrupts addiction to cigarette smoking. Science. (2007) 315:531–4. doi: 10.1126/science.1135926

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Paulus MP, Tapert SF, Schuckit MA. Neural activation patterns of methamphetamine-dependent subjects during decision making predict relapse. Arch Gen Psychiatry. (2005) 62:761–8. doi: 10.1001/archpsyc.62.7.761

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Wang GJ, Volkow ND, Fowler JS, Cervany P, Hitzemann RJ, Pappas NR, et al. Regional brain metabolic activation during craving elicited by recall of previous drug experiences. Life Sci. (1999) 64:775–84. doi: 10.1016/S0024-3205(98)00619-5

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Bonson KR, Grant SJ, Contoreggi CS, Links JM, Metcalfe J, Weyl HL, et al. Neural systems and cue-induced cocaine craving. Neuropsychopharmacology. (2002) 26:376–86. doi: 10.1016/S0893-133X(01)00371-2

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Seeley WW, Menon V, Schatzberg AF, Keller J, Glover GH, Kenna H, et al. Dissociable intrinsic connectivity networks for salience processing and executive control. J Neurosci. (2007) 27:2349–56. doi: 10.1523/JNEUROSCI.5587-06.2007

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Domino EF, Ni L, Xu Y, Koeppe RA, Guthrie S, Zubieta JK. Regional cerebral blood flow and plasma nicotine after smoking tobacco cigarettes. Prog Neuropsychopharmacol Biol Psychiatry. (2004) 28:319–27. doi: 10.1016/j.pnpbp.2003.10.011

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Zubieta JK, Heitzeg MM, Xu Y, Koeppe RA, Ni L, Guthrie S, et al. Regional cerebral blood flow responses to smoking in tobacco smokers after overnight abstinence. Am J Psychiatry. (2005) 162:567–77. doi: 10.1176/appi.ajp.162.3.567

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Court JA, Lloyd S, Thomas N, Piggott MA, Marshall EF, Morris CM, et al. Dopamine and nicotinic receptor binding and the levels of dopamine and homovanillic acid in human brain related to tobacco use. Neuroscience. (1998) 87:63–78. doi: 10.1016/S0306-4522(98)00088-8

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Yu R, Zhao L, Lu L. Regional grey and white matter changes in heavy male smokers. PLoS ONE. (2011) 6:e27440. doi: 10.1371/journal.pone.0027440

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Mahmoudi A, Takerkart S, Regragui F, Boussaoud D, Brovelli A. Multivoxel pattern analysis for fmri data: a review. Comput Math Methods Med. (2012) 2012:14. doi: 10.1155/2012/961257

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Girshick R. Fast r-cnn. In: Proceedings of the IEEE international conference on computer vision. Santiago (2015).

Google Scholar

Keywords: smoking addiction, multivoxel pattern analysis, voxel-based morphometry, machine learning, structural magnetic resonance imaging

Citation: Ye Y, Zhang J, Huang B, Cai X, Wang P, Zeng P, Wu S, Ma J, Huang H, Liu H, Dan G and Wu G (2021) Characterizing the Structural Pattern of Heavy Smokers Using Multivoxel Pattern Analysis. Front. Psychiatry 11:607003. doi: 10.3389/fpsyt.2020.607003

Received: 16 September 2020; Accepted: 16 December 2020;
Published: 04 February 2021.

Edited by:

Yanhui Liao, Sir Run Run Shaw Hospital, China

Reviewed by:

Kiyotaka Nemoto, University of Tsukuba, Japan
Dahua Yu, Inner Mongolia University of Science and Technology, China

Copyright © 2021 Ye, Zhang, Huang, Cai, Wang, Zeng, Wu, Ma, Huang, Liu, Dan and Wu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Guo Dan, danguo@szu.edu.cn; Guangyao Wu, wuguangy2002@163.com

These authors have contributed equally to this work

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