国际医学放射学杂志ＩｎｔＪＭｅｄＲａｄｉｏｌ２０21May鸦44穴3雪：254-258；276
论著
早期城市化对健康青年人脑灰质体积和
执行功能的影响
刘嘉玮1郭力宁1李魏薇1马林1满玮琦1，2徐佳圆1刘风1于春水1
【摘要】目的探讨早期城市环境暴露对健康青年人脑灰质体积（GMV）及执行功能的影响，以加深对城市化与人脑及认知发育之间关系的理解。方法在天津地区前瞻性收集1021名18~30岁的汉族健康青年人，平均（23.76±2.34）岁；平均受教育（16.64±1.91）年。采集受试者早期（0~15岁）的城市化分数、全脑T
WI数据及Go/No-go
1执行功能数据，记录Go条件的正确率（GA）、No-go条件的正确率（NGA）及Go条件的反应时间（GRT）。采用SPM12软件进行影像数据分析，找出与早期城市化分数显著相关的GMV脑区，利用Spearman等级相关分析早期城市化分数、显著团块GMV、执行功能数据之间的相关性。构建早期城市
化分数、显著团块GMV、执行功能的中介模型，在控制中介变量GMV的条件下，分析早期城市化分数对执行功能的直接效应，以及早期城市化分数通过GMV对执行功能的中介效应。结果早期城市化分数与背侧前扣带皮质和内侧额上回组成团块的GMV呈负相关（Ｐ<0.05，ＦWE 校正）。早期城市化分数与背侧前扣带回和内侧额上回形成团块的GMV、GRT呈负相关（r=-0.17，r=-0.18，P<0.001），与GA呈正相关（r=0.18，P<0.001）。背侧前扣带回和内侧额上回形成团块的GMV与GA、NGA呈负相关（r=-0.10，P=
0.001，r=-0.06，P=0.05），与GRT呈正相关（r=0.08，P=0.01）。背侧前扣带回和内侧额上回形成团块的GMV在早期城
市化分数与GRT之间起着部分中介的作用，早期城市化分数对GRT的直接效应值为-0.17，显著团块的GMV在早期城市化分数与GRT间的中介效应值为-0.01，中介效应占总效应的6.74%。结论早期城市环境暴露对青年人前额叶灰质发育及执行功能具有促进作用。
【关键词】磁共振成像；城市化；灰质体积；执行功能；中介效应分析
中图分类号：R322.8；R455.2文献标志码：A
The effects of early urbanization on brain grey matter volume and executive function of healthy youn
g adults LIU Jiawei1,GUO Lining1,LI Weiwei1,MA Lin1,MAN Weiqi1,2,XU Jiayuan1,LIU Feng1,YU Chunshui1.1Department of Radiology and Tianjin Key Laboratory of Functional Imaging,Tianjin Medical University General Hospital,Tianjin300052, China;2Academy of Medical Engineering and Translational Medicine,Tianjin University.Corresponding author:YU Chunshui,E-mail:******************
【Abstract】Objective To explore the effects of early exposure to urban environment on grey matter volume(GMV) and executive function of healthy young adults,and to improve our understanding of the association between urbanization and developments of brain and cognitive.Methods A total of1021healthy young adults aged18through30years old were collected in Tianjin in this prospective study,with an average age of23.76±2.34years,the average years of education was
16.64±1.91.The early urbanization scores(the sum of1-15years of age),whole brain T1WI data,and executive function scores
derived from Go/no-GO test were collected,and the Go trial accuracy(GA),No-go trial accuracy(NGA),and Go trial reaction time(GRT)were recorded.Voxel-wi multiple regression was conducted to identify the brain regions where GMV were correlated with early urbanization scores ba
d on SPM12,and Spearman correlation was performed to test the relationship between the early urbanization scores and the GMV of significant brain regions or executive function.Mediation model was ud to test the mediation effect of the GMV of significant brain regions on the associations between the early urbanization scores and executive function.The direct effect of the early urbanization scores on GRT under control of the GMV of cluster, and the mediating effect of the GMV of cluster between early urbanization and GRT were analyzed.Results Early 作者单位：1天津医科大学总医院医学影像科，天津市功能影像重点实验室，天津300052；2天津大学医学工程与转化医学研究院通信作者：于春水，E-mail：******************
基金项目：国家重点研发计划（2018YFC1314300）；国家自然科学基金重点项目（82030053）；天津市教委科研计划项目（2019KJ195）DOI:10.19300/j.2021.L18803
254
国际医学放射学杂志ＩｎｔｅｒｎａｔｉｏｎａｌＪｏｕｒｎａｌｏｆＭｅｄｉｃａｌＲａｄｉｏｌｏｇｙ２０21May鸦44穴3雪
urbanization scores were negatively correlated with GMV in a cluster consisting of dorsal anterior cingulate cortex(dACC) and medial superior frontal gyrus(mSFG)(P<0.05,FWE correction).Early urb
anization scores were negatively correlated with GMV in a cluster consisting of dACC and mSFG,and GRT(r=-0.17,r=-0.18,P<0.001),and positively correlated with GA (r=0.18,P<0.001).The GMV in a cluster consisting of dACC and mSFG were negatively correlated with GA,NGA(r=-0.10, P=0.001and r=-0.06,P=0.05,respectively),and positively correlated with GRT(r=0.08,P=0.01).GMV in a cluster consisting of dACC and mSFG showed a partial mediation effect between the early urbanization scores and GRT.The direct effect value of the early urbanization scores on GRT was-0.17,the mediating effect value of the GMV of cluster between the early urbanization and GRT was-0.01,and the mediating effect accounted for6.74%of the total effect.Conclusion The results suggest that early exposure to urban environment may facilitate prefrontal grey matter development and executive function in young adults.
【Keywords】Magnetic resonance imaging;Urbanization;Grey matter volume;Executive function;Mediation effect analysis
ＩｎｔＪＭｅｄＲａｄｉｏｌ，２０21，44（3）：254-258；276
环境无处不在且影响着人类的精神健康，精神类疾病中50%的风险因素可归因于环境[1-2]。国家统计局数据表明，中国的城镇人口比例从1978年的17.92%提高到2019年的60.6%[3-4]，不断增长的比例
强调了研究城市化程度对精神健康影响的重要性。近期研究[5]发现，城市化与大脑处理社会压力相关区域的活动存在关联，这种大脑结构或功能的变化可能参与调节城市化对精神健康或认知功能的影响，而目前较少有研究关注三者之间的关联。在生命的早期阶段，即出生到青春期（15岁之前），人脑处于快速发育时期，认知功能也在不断增强，人脑及认知功能更容易受到城市化程度的影响[5]。本研究旨在探索早期阶段的城市环境暴露对健康青年人脑灰质体积及执行功能的影响。
1资料与方法
1.1一般资料前瞻性在天津市招募健康青年人。所有受试者在进行实验前均详细了解实验内容并签署知情同意书。纳入标准：①18～30周岁，身体健康；②汉族人；③右利手。排除标准：①吸烟史，过度饮酒史，药物、毒品滥用或依赖史；②妊娠期或月经期女性；③MRI禁忌证；④有神经精神疾病家族史；
⑤实验前一晚睡眠不足7h。最终纳入1021名受试者，男467名，女554名；年龄18~30岁，平均（23.76±2.34）岁；受教育程度9~22年，平均（16.64±1.91）年。本项研究得到了天津医科大学总医院伦理委员会批准。
凉拌藕怎么做
1.2早期城市化分数采用城市化分数[5-6]评价受试者早期阶段（0~15岁）的城市暴露水平。首先，基于问卷调查获得受试者15岁以前的居住地址，根据国家统计局2008年8月1日起施行的《关于统计上划
分城乡的规定》（国函〔2008〕60号）把居住地址分为3个等级：1代表农村，2代表城镇，3代表城市。然后，将受试者居住地的等级乘以受试者在该地区居住的时间（年），如有住址迁移则将结果累加计算获得前15年的累积值，即为该受试者的早期城市化分数。早期城市化分数为15～45分；15分代表一直生活在农村，45分代表一直生活在城市；分数越高表示城市环境暴露时间越长。
1.3MRI数据采集及处理
1.3.1数据采集采用GE Discovery MR7503.0T MR扫描设备，头部8通道相控阵线圈。对全部受试者先行常规T
2
WI检查以确定无颅内器质性病变，
随后行高分辨矢状面三维T
1
WI全脑数据采集。扫描参数：TR8.14ms，TE3.17ms，TI450ms；翻转角12°；视野256mm×256mm；矩阵256×256；层厚1mm，层间距0；层数188。
1.3.2数据预处理采用SPM12软件中的插件VBM8进行数据预处理。具体步骤：①根据先验模板采用“联合分割”算法将全脑结构像分割为灰质、白质和脑脊液，即可生成灰质概率图，它代表每个体素隶属于灰质的概率；②采用VBM8内置的DARTEL高阶配准法将每名受试者的灰质概率图配准到MNI标准空间以保证每个体素在受试者之间相对应，并将每个体素重采样到1.5mm×1.5mm×1.5mm；③使用标准化过程中产生的非线性变形场参数对灰质概率图进行调制得到GMV图，每个体素内的值代表绝对GMV；④将GMV图进行平滑处理（高斯核为8mm×8mm×8mm）以提高图像的信噪比和数据正态性，并且减少配准不齐的影响。
1.4执行功能数据采集采用Go/No-go任务评估
255
国际医学放射学杂志ＩｎｔＪＭｅｄＲａｄｉｏｌ２０21May 鸦44穴3雪
被试执行抑制能力，使用E-Prime 2.0软件呈现并收集结果。在实验中，屏幕上会呈现字母X 或者Y ，如果当前呈现的字母和前一次出现的字母不同，被试需要快速按键做出反应，这种条件为Go 条件；如果当前呈现字母和前一次呈现字母相同，被试此时不能按键，如果按键则记为错误，这种条件为No-go 条件，占整个测试的10%[7]。任务包括练习和正式测试，练习通过后开始正式测试。正式测试分为2组，每组210个试次，要求在字母显示后尽可能快速和准确地做出反应。记录Go 条件的正确率
（Go trial
accuracy,GA ）、
No -go 条件的正确率（No -go trial accuracy,NGA ）、Go 条件的反应时间（Go trial reaction time,GRT ）。1.5
统计学方法
采用SPM12软件进行影像数据
分析。采用基于体素的多元线性回归方法分析早期城市化分数与体素水平GMV 之间的关联性，采用
整体错误率（family wi error ，FWE ）方法对结果进
行多重比较校正（P <0.05），并提取通过校正的脑区团块的平均GMV 值。采用Spearman 等级相关分析早期城市化分数值、显著团块GMV 、执行功能之间的相关性。采用SPSS24.0软件PROCESS 程序中的Model4构建早期城市化分数值、
显著团块GMV 、执行功能之间的中介模型，并分析总效应、直接效应、中介效应。其中，总效应由
直接效应和中介效应共同构成，直接效应是指在控制中介变量GMV 时，早期城市化分数对执行功能的效应；中介效应是早期城市化分数对GMV 效应和GMV 对执行功能效应的乘积。用Bootstrap 重采样法抽样5000次来检验
中介效应的显著性，效应值以标准误差、
95%置信区间的下限和上限表示，如果置信区间不包含零，则
认为中介效应显著。上述分析均将年龄、性别、教育年限以及全脑灰质体积作为控制变量予以控制。
2结果
2.1关联性分析
缺点英语2.1.1早期城市化分数与体素水平GMV 之间的回归分析早期城市化分数与背侧前扣带皮质（dorsal anterior cingulate cortex,dACC ）和内侧额上回（medial superior frontal cortex,mSFG ）组成团块的GMV 呈负相关（MNI 峰值坐标：0，35，39；体素2124，峰值T =
-6.22）
（图1），即早期城市化分数值越高，该团块的GMV 越小。
2.1.2各变量之间的相关性早期城市化分数与dACC 及mSFG 组成团块的平均GMV 呈负相关（r =-0.17，P <0.001）。早期城市化分数与GA 呈正相关（r =0.18，P <0.001），与GRT 呈负相关（r =-0.18，
P <0.001）。dACC 和mSFG 组成团块的平均GMV 与NGA 、GA 呈负相关（GMV 与NGA ：r=-0.06，P =0.05；GMV 与GA ：r=-0.10，P =0.001），与GRT 呈正相关（r =0.08，P=0.01）
。2.2GMV 在早期城市化分数与执行功能关联中的
中介效应早期城市化分数能负向预测显著团块的GMV 值（βa =-0.14，P <0.001）和GRT （βc ’=-0.17，P <0.001）。显著团块的GMV 能正向预测GRT （βb =
0.10，P <0.05），详见图2。早期城市化分数可经中介路径和直接路径预测GRT ，中介路径（早期城市化
分数→显著团块GMV →GRT ）和直接路径（早期城市化分数→GRT ）的Bootstrap 效应值的95%置信区间（CI ）均不包含0，表明中介效应和直接效应均达
到显著水平，详见表1。因此，dACC 和mSFG 团块的平均GMV 值在早期城市化分数与GRT 的关联中
图1GMV 与早期城市化分数值显著相关的脑区（P <0.05，FWE 校正）。色带中T 值的绝对值越大，回
归系数越显著，
相关性越强。
256
才士
国际医学放射学杂志ＩｎｔｅｒｎａｔｉｏｎａｌＪｏｕｒｎａｌｏｆＭｅｄｉｃａｌＲａｄｉｏｌｏｇｙ２０21May 鸦44穴3雪
表1
早期城市化分数对执行功能影响的中介效应路径及效应值分析
执行功能效应路径
效应值标准化效应值
g姓标准误95%CI
下限上限GA
早期城市化分数→显著团块GMV →GA 0.00010.01890.00000.00000.0001早期城市化分数→GA 0.00040.15070.00010.00020.0007总效应
0.0005
0.1696
0.0000
0.0000
0.0010
NGA 早期城市化分数→显著团块GMV →NGA 0.00010.03180.0001-0.00010.0003早期城市化分数→NGA 0.00080.06220.0005-0.00010.0018总效应
0.0009
0.0940
0.0000
0.0000
0.0020
GRT 早期城市化分数→显著团块GMV →GRT -0.0790-0.01260.0411-0.1646-0.0030早期城市化分数→GRT -0.9880-0.17430.1814-1.3485-0.6344总效应
-1.0670
-0.1869
0.1770
-1.4110
-0.7180
起到部分中介作用，且部分中介效应占总效应的比
例为6.74%。3讨论
近年来，中国经历了快速的城市化过程。由于不断增长的城市化水平带来的环境变化与精神健康息息相关，故明确早期城市环境暴露与大脑结构及认知的关联有助于更好地理解城市化水平对精神健康的影响。本研究使用城市化分数来评估受试者早期所在城市环境暴露水平，通过关联分析建立了早期城市环境暴露与受试者背侧前扣带皮质和内侧额上回GMV 以及执行功能之间的关联，并通过中介效应模型发现，该脑区的GMV 在早期城市化与执行功能关联中具有显著的中介效应。
既往研究[8-9]显示，城市化水平较高地区居民比较低地区的有更好的认知功能表现，尤其是在信息处理速度、执行功能和口语学习等方面。本研究结果提示早期城市环境暴露越高，健康青年人的执行功能越强，与既往研究结果相符。这可能是多因素共同作用的结果。首先，城市化水平的提升会带来
更好的生活条件，如完善的医疗、教育、营养条件等[10-11]，从而促进认知的发展。其次，城市生活会增加周围环境的复杂度，研究[12]表明丰富的外界环境对认知功能的发展有促进作用。综上，城市化发展在一定程度上促进个体的认知发展，尤其是执行功能。
本研究还发现，早期城市环境暴露水平较高的个体dACC 及mSFG 的GMV 较小。dACC 和mSFG
具有很强的解剖与功能连接[13-14]，提示这2个脑区参与了相似的功能。dACC 是执行控制功能的经典脑区，参与冲突监控、错误检测和奖励评估等多种执行功能[14-17]。mSFG 通过与dACC 形成丰富的解剖联系，也参与执行控制功能[18-19]。本研究同时发现，dACC 和mSFG 的GMV 越小，执行功能越强。一个脑区GMV 大小与认知功能的关系与年龄因素有关。例如，老年人前额叶脑区的GMV 越小提示认知功能越差[20]。但是，青年人并非如此，前额叶皮质是成熟时间最晚的脑区，完全成熟可能要持续到25岁左右[21]。前额叶成熟之前主要经历了经验依赖性突触修剪过程，该过程将通过去除冗余的突触，使得前额叶皮质变薄、体积变小，而认知处理能力增强[22-24]。本研究纳入的受试者年龄为18~30岁，正处于前额叶成熟的后期阶段，因此dACC 和mSFG 的GMV 与执行功能呈负相关。
中介效应分析发现，个体dACC 和mSFG 的GMV 减小介导了早期城市环境暴露水平与执行功能之间的正性效应。这一发现提示，早期城市环境暴露水平越高，受试者的经济条件、营养水平、教育及医疗条件等社会环境因素越好，
这有助于青春期
图2显著团块平均GMV 在早期城市化分数对GRT 影响的中介效应模型。*P <0.05,**P <0.001。
小学德育原则
257
国际医学放射学杂志ＩｎｔＪＭｅｄＲａｄｉｏｌ２０21May 鸦44穴3雪
（下转第276页）
前额叶的突触修剪过程，进而导致前额叶GMV 减小，而更加高效的突触连接使得执行功能增强。这一解释也与以往的研究[25-26]发现相一致，如青春期大量饮酒作为个体特殊环境暴露可以干扰突触修剪过程导致前额叶GMV 增大，进而造成执行功能减低。以上均提示环境暴露可能通过影响青春期突触修剪过程进而影响GMV 及执行功能。
本研究还存在一些局限性。首先，早期城市化分数指标只是对受试者15岁前生活成长地区居住时间的反映，而未进一步细化家庭居住地址和学校所在地的不同，可能对结果造成影响。其次，本研究只探讨了城市环境暴露对人脑及行为的影响，而未能考虑遗传因素及遗传环境交互作用的影响。这些因素的效应需要在今后研究中进行探讨。再次，脑成像指标中只选择了GMV ，而且GMV 只能解释一部分执行功能的变异。今后研究中需要纳入更多的脑成像指标，以全面认识早期城市化对人脑结构、功能
与连接的影响，以及这些指标与执行功能之间的关系。最后，本研究只考虑了青年人15岁前城市化程度的累积效应，并没有进行敏感期分析，无法明确哪个年龄阶段的城市环境暴露对人脑及认知功能的影响最大。尽管存在这些局限性，本研究基于较大样本发现了早期城市化程度对青年人脑灰质发育以及执行功能的正性效应，增进了对城市环境暴露与人脑健康及行为关系的理解。
参考文献：
[1]
于春水.人脑影像表型的基因组与暴露组广度关联研究[J/OL].国际医学放射学杂志，2021,44:249-253.YU Chunshui.Genome -and exposome -wide association studies of human brain imaging phenotypes [J/OL].Int J Med Radiol,2021,44:249-253.DOI:10.19300/j.2021.S18881.
[2]van Os J,Kenis G,Rutten BPF.The environment and schizophrenia [J].Nature,2010,468:203-212.DOI:10.1038/nature09563.
[3]
中华人民共和国统计局.人口总量平稳增长人口素质显著提升———新中国成立70周年经济社会发展成就系列报告之二十[EB/OL].(2019-08-22)[2021-03-22].v/tjsj/zxfb/201908/t20190
822_1692898.html.
[4]中华人民共和国统计局.中国统计年鉴-2020[M].北京：中国统计出版社,2020.National Bureau of Statistics of China.China statistical yearbook 2020[M].Beijing:China Statistics Press,2020.
最强记忆力[5]Lederbogen F,Kirsch P,Haddad L,et al.City living and urban upbringing affect neural social stress processing in humans [J].Nature,2011,474:498-501.DOI:10.1038/nature10190.
[6]Haddad L,Sch 覿fer A,Streit F,et al.Brain structure correlates of urban upbringing,an environmental risk factor for schizophrenia[J].Schizophr Bull,2015,41:115-122.DOI:10.1093/schbul/sbu072.
[7]Kaufman JN,Ross TJ,Stein EA,et al.Cingulate hypoactivity in cocaine urs during a GO-NOGO task as revealed by event-related functional magnetic resonance imaging [J].J Neurosci,2003,23:7839-7843.DOI:10.1523/JNEUROSCI.23-21-07839.2003.
[8]Cassarino M,O ’
Sullivan V,Kenny RA,et al.Environment and cognitive aging:a cross -ctional study of place of residence and cognitive performance in the Irish longitudinal study on aging [J].Neuropsychology,20
维生素a作用16,30:543-557.DOI:10.1037/neu0000253.
[9]Gupta S,Vaida F,Riggs K,et al.Neuropsychological performance in mainland China:the effect of urban/rural residence and lf -reported daily academic skill u [J].J Int Neuropsychol Soc,2011,17:163-173.DOI:10.1017/S1355617710001384.
[10]Phillips MR,Zhang J,Shi Q,et al.Prevalence,treatment,and
associated disability of mental disorders in four provinces in China during 2001-05:an epidemiological survey [J].Lancet,2009,373:2041-2053.DOI:10.1016/S0140-6736(09)60660-7.
[11]曹雷，郑景山，曹玲生，等.中国2012年国家免疫规划疫苗常规免
疫报告接种率分析[J].中国疫苗和免疫,2013,19:389-396.Cao L,Zhang JS,Cao LS,Analysis on the coverage of national immunization program vaccines reported in China [J].Chin J
Vaccines Immunizaion,2013,19:389-396.DOI:1006-916X （2013）05-0389-09.
[12]Smail MA,Smith BL,Nawreen N,et al.Differential impact of stress
and environmental enrichment on corticolimbic circuits [J].Pharmacol Biochem Behav,2020,197:172993.DOI:10.1016/j.pbb.2020.172993.
[13]Li W,Qin W,Liu H,et al.Subregions of the human superior frontal
gyrus and their connections [J].Neuroimage,2013,78:46-58.DOI:10.uroimage.2013.04.011.
[14]Heilbronner SR,Hayden BY.Dorsal anterior cingulate cortex:a
bottom -up view [J].Annu Rev Neurosci,2016,39:149-170.DOI:10.1146/annurev-neuro-070815-013952.
白玉虎皮兰[15]Botvinick MM,Braver TS,Barch DM,et al.Conflict monitoring and
cognitive control[J].Psychol Rev,2001,108:624-652.DOI:10.1037/0033-295x.108.3.624.[16]Botvinick
MM.
Conflict
monitoring
and
decision
making:
reconciling two perspectives on anterior cingulate function [J].Cogn Affect Behav Neurosci,2007,7:356-366.DOI:10.3758/cabn.7.4.356.[17]Holroyd CB,Coles M.The neural basis of human error processing:
reinforcement learning,dopamine,and the error -related negativity [J].Psychol Rev,2002,109:679-709.DOI:10.1037/0033-295X.109.4.679.
[18]Domenech P,Koechlin E.Executive control and decision-making in
the prefrontal cortex[J].Curr Opin Behav Sci,2015,1:101-106.DOI:10.beha.2014.10.007.
[19]Mace RA,Waters AB,Sawyer KS,et al.Components of executive
function model regional prefrontal volumes [J].Neuropsychology,2019,33:1007-1019.DOI:10.1037/neu0000563.
[20]Ruscheweyh R,Deppe M,Lohmann H,et al.Executive performance
258
