BDNF基因、同伴关系与青少年早期抑郁：基于动态发展视角

/BDNF基因、同伴关系与青少年早期抑郁：基于动态发展视角【摘要】青少年抑郁是遗传基因与环境动态交互作用的结果，但是现有研究忽视了抑郁遗传效应的发展动态性。本研究通过对1086名青少年(平均年龄12.32，50％女生)进行3年的追踪，分别从遗传效应的年龄差异以及遗传效应影响抑郁发展轨迹的角度，考察BDNF基因与同伴关系对青少年抑郁的动态影响。结果显示：(1)在3个时间点上，BDNF基因与同伴拒绝交互影响青少年抑郁，但其作用模式存在年龄差异：12岁时，MetMet基因型携带者对环境敏感性高于ValMet基因型携带者；13岁时，MetMet和ValVal基因型携带者对环境的敏感性均高于ValMet基因型携带者；14岁时，ValVal基因型携带者对环境的敏感性高于ValMet基因型携带者。(2)青少年早期抑郁呈线性增长趋势，但是抑郁初始水平与增长速度无关。(3)BDNF基因与同伴拒绝交互预测青少年抑郁的初始水平，相比ValMet基因型，携带MetMet基因型的青少年在经历同伴拒绝后抑郁初始水平更高。(4)BDNF基因显著预测青少年抑郁增长速度，相比ValMet基因型携带者，携带MetMet和ValVal基因型的青少年抑郁增长速度更快。【关键词】青少年抑郁；BDNF基因；同伴拒绝；年龄差异；潜变量增长模型

1问题提出 近20年来，抑郁的分子遗传研究取得了系列重要突破。既有研究探索了众多抑郁的风险基因及其与环境的交互作用机制，揭示了基因通过内表型影响抑郁的潜在过程(如Caoetal.，2021；Caspietal.，2003)。然而，多数分子遗传研究仍然采用静态视角，对抑郁遗传效应的动态发展问题关注不足。特别是在该领域一系列关键问题上，譬如“抑郁的遗传效应是否随年龄增长发生变化?”“风险基因如何在抑郁发展中发挥作用”等，已有研究尚不能给出确定的答案。本研究采用动态发展视角，考察脑源性神经营养因子(brain-derivedneurotrophicfactor，BDNF)基因与同伴关系交互作用的年龄差异及其对青少年抑郁发展轨迹的影响。 1.1BDNF基因、同伴关系与青少年抑郁 根据抑郁的神经营养假说(neurotrophinhypothesis)，中枢BDNF功能缺陷会导致海马和前额叶皮质神经元受损，引发抑郁障碍(Martinowichetal.，2007)。此外，抑郁的下丘脑-垂体-皮质醇假说(hypothalamic-pituitary-cortisolhypothesis)指出，应激反应失调是引发抑郁的重要机制(Belmaker&Agam，2008)。BDNF基因通过调节BDNF的表达，不仅在突触可塑性和神经元连接性中发挥重要作用(Kowian'skietal.，2018)，而且影响个体对压力的反应性(Notaras&vandenBuuse，2020)，因而是抑郁的“易感”基因。BDNF基因位于11号染色体p13区，其第11号外显子区的Va166Met(rs6265)单核苷酸多态性与BDNF的细胞内转运和活性依赖性分泌密切相关。该多态性上C到A突变导致缬氨酸-蛋氨酸(Val-Met)替换，使得Met(A)等位基因与BDNF分泌减少与活性降低相关(Gonzalez-Castroetal.，2015)。基于上述神经生理基础，诸多研究发现了BDNF基因与抑郁的关联(Hiltetal.，2007)，并且在抑郁相关内表型——杏仁核和海马体积、功能研究中获得了支持(Caseyetal.，2009)。 在探索抑郁易感基因的同时，研究者亦关注影响抑郁的环境风险，尤其是青少年期的同伴关系。研究显示，同伴关系(如同伴接纳和同伴拒绝)与抑郁的发生发展密切相关，其在青少年期的重要性超过了其他环境风险(Rudolphetal.，2000)。同伴接纳是抑郁的保护性因素，高水平的同伴接纳不仅是青少年自尊和归属感的重要来源，并且与较低的皮质醇浓度有关，因而能够降低抑郁风险(Baumeister&Leary，1995；Blackhartetal.，2007)。同伴拒绝则会导致青少年形成消极的自我图式，导致应激功能失调，增加抑郁风险(Blackhartetal.，2007；Plattetal.，2013)。尽管同伴接纳和同伴拒绝具有中等程度的负相关，但其对抑郁的影响和作用模式存在差异(Guyeretal.，2012；Hsuetal.，2015)，因此本研究分别考察同伴接纳和同伴拒绝在抑郁发展中的作用。 更重要的是，BDNF基因与环境交互影响抑郁。譬如，Gottfredson等(2015)对青少年的研究发现，在经历同伴侵害后，携带ValVal基因型的青少年比携带Met等位基因的青少年具有更高的抑郁水平。此外，Chen等(2013)研究显示，BDNF基因与压力性生活事件、母亲教养行为交互影响青少年抑郁，相比Met等位基因携带者，Val等位基因携带者对积极或消极环境因素更加敏感，表现出更少或更多的抑郁症状。有关儿童(Meyeretal.，2018)和成人(Aguileraetal.，2009)的研究虽然也发现了BDNF基因与环境的交互作用，却显示Met等位基因携带者对环境的敏感性更高。譬如，Meyer等(2018)发现BDNF基因与消极教养对错误相关负波(error-relatednegativity，ERN，一种抑郁的内表型)具有交互作用，相比ValVal基因型携带者，携带Met等位基因的幼儿在经历高消极教养后具有更高的ERN。尽管既有研究在抑郁的易感基因上(Valvs.Met)存在分歧，但这些不一致的研究结果提示BDNF基因及其与环境的交互作用可能存在年龄差异。 1.2BDNF基因及其与环境因素的交互作用：年龄差异的证据 遗传与环境交互作用的发展动态性问题是发展心理学和行为遗传学的重要理论问题。根据Lerner(1978)的天性与教养的动态交互观，遗传基因与环境因素的交互效应会随年龄或发展阶段而动态变化。最近，Boyce等(2020，2021)的“基因-环境-时间”交互理论框架(theinteractivegene-environment-timeframework)，进一步指出基因(DNA序列、表观遗传过程等)、环境(积极、消极经历等)和时间(年龄、发展阶段、关键期等)三个因素交互影响个体的心理健康适应结果。 近年来，在抑郁、反思(rumination)、线索学习(cuedlearning)等多个领域的研究均发现BDNF基因的直接效应可能是随年龄变化的函数。譬如，Hilt等(2007)发现BDNFVal66Met多态性的ValVal基因型与青少年期起始的抑郁和更高的反思有关，而ValMet基因型与成年期起始的抑郁有关。Casey等(2009)发现在童年期，相比ValVal基因型携带者，Met等位基因携带者对条件刺激具有更长的反应时(线索学习的测量指标)，但在青少年时期，两种基因型间的反应时差异消失。 据我们所知，目前尚未有研究直接考察BDNF基因×环境对抑郁的年龄差异，但是有关大脑皮质体积的研究为该年龄差异提供了支持(Caseyetal.，2009)。研究显示，在无压力环境下(控制组)，Met等位基因和ValVal基因型携带者间的皮质体积差异随年龄增长逐渐减小，在BDNF含量达到顶峰时，其差异基本消失；但是在压力环境下，BDNF在海马等关键脑区的含量降低，并且使得Met等位基因携带者具有更小的皮质体积，其与ValVal基因型携带者间的差距随年龄逐渐增大。类似的，Dalton等(2014)关于BDNF基因、5-HTTLPR基因累加遗传得分与家庭环境质量对抑郁的研究发现，在15岁前，累加基因得分与家庭环境的交互作用模式符合不同易感性模式；在15岁后，累加基因得分与环境的交互模式符合素质-压力模式。这些研究不仅为BDNF基因及其与环境交互作用的年龄差异提供了证据，更重要的是，这些研究提示青少年期可能是BDNF基因效应产生动态变化的关键时期。 关于脑内BDNF浓度的发展研究和表观遗传研究为理解BDNF基因效应在青少年期发生动态变化提供了启示。一方面，在青少年期，脑内BDNF浓度发生重要变化，这可能是引发BDNF基因效应呈现年龄差异的重要原因。动物研究显示，脑内BDNF含量呈倒U型发展模式，在青少年期达到峰值(Ivanova&Beyer，2001)。关于人类背外侧前额叶的BDNFmRNA含量变化的研究也显示从童年期到青少年期，BDNF表达持续增加(Websteretal.，2002)。Casey等(2009)进一步提示BDNF基因与脑内BDNF浓度共同影响BDNF功能。基于此，在前青少年期，Met等位基因能够加剧较低的BDNF浓度导致的功能失调，呈现风险性；但是在青少年期，脑内BDNF含量的升高反而使得Met等位基因的风险性消失。 另一方面，青少年的环境经历也发生重要变化(Steinberg，2022)，不同时期环境经历的差异也会导致BDNF遗传效应的动态变化(Kageyamaetal.，2022；Surietal.，2013)。动物研究显示，不同年龄段的高果糖摄入会导致BDNF基因甲基化差异，仅在童年期和青少年期的高果糖摄入会导致小鼠BDNF基因启动子区甲基化增加，但是在成年早期和晚期小鼠中并未观察到DNA甲基化的变化(Kageyamaetal.，2022)。此外，早期消极经历对小鼠BDNF基因甲基化的短期和长期影响亦存在差异，Suri等(2013)的研究显示早期母子分离经历会导致成年早期小鼠的BDNFIV启动子区甲基化降低进而增加BDNF含量，但是会增加成年中期小鼠的BDNF甲基化。基于动物表观遗传研究结果，不同时期的环境经历或者早期环境经历对不同发展阶段BDNF基因甲基化水平存在差异性影响，可能在不影响基因型的情况下，使得BDNF基因在不同年龄阶段表现出不同的功能，进而呈现出年龄差异的特点。 值得指出的是，青少年早期可能是观察BDNF基因功能变化的关键期。一方面，有关BDNF神经生理功能的研究显示，相比其他发展阶段，BDNF基因调节恐惧相关神经环路的关键期是青少年早期(Dinchevaetal.，2016)，并且从青少年早期到青少年中期，BDNF含量及其TrkB受体在眶额叶和内侧前额叶的表达水平发生明显的变化(Shapiroetal.，2017)。恐惧神经环路、眶额叶和内侧前额叶等均是与抑郁密切相关的关键脑区(Hare&Duman，2020；Ho&King，2021；Subramaniametal.，2018)，因而青少年早期及其向青少年中期的过渡阶段是观察BDNF基因效应变化的关键时期。另一方面，一项4～15岁的追踪研究发现，BDNF基因及其与母亲产前焦虑的交互作用模式对内化问题的影响在4～13岁时相对稳定，但是在青少年中期(15岁)BDNF基因的作用模式与4～13岁时存在差异(O'Donnelletal.，2014)。这也提示，BDNF基因在童年期与青少年早期的效应差异可能没有青少年早期与青少年中期的差异显著。基于此，本研究假设青少年早期是观测BDNF基因功能变化的窗口期。 1.3BDNF基因与环境交互模式的发展动态性：素质-压力、优势敏感与不同易感模型 BDNF基因与环境的交互作用模式可能被三种理论模型所解释：“素质-压力”(diathesis-stress)、“优势敏感”(vantagesensitivity)与“不同易感性”(differentialsusceptibility)模型。“素质-压力”模型指出携带易感基因的个体对消极环境更加敏感，因而在经历不利环境时才会产生心理或行为问题(Monroe&Simons，1991)；“优势敏感”模型则认为携带易感基因的个体对积极环境更加敏感，在支持性环境中则会发展得更好(Pluess&Belsky，2013)；“不同易感性”模型则是前两种模型的整合，认为所谓的易感基因事实上是一种遗传可塑性，更具可塑性的个体对积极和消极环境均更加敏感，因而相比不携带可塑性基因的个体，具有可塑性的个体在积极环境中发展“更好”，但在消极环境中发展“更坏”(Belskyetal.，2007)。 然而，这三种理论模型并非相互对立。近期，研究者指出个体的环境敏感性会受早期环境经历影响而产生变化，携带环境敏感性基因的个体在生命早期对积极和消极环境均具有敏感性(即不同易感性模式)，但当这些个体在偏积极环境中成长时，其对积极环境的敏感性则更显著(即优势敏感模式)；当其在偏消极环境中成长时，则对消极环境的敏感性更强(即素质-压力模式；Pluess，2015)。由此，基因与环境的交互作用模式可能也会随个体发展而动态变化，其变化模式取决于个体的环境经历。如前所述，Dalton等(2014)研究显示在15岁前后，基因-环境交互作用从不同易感性模式向素质-压力模式变化。因此，本研究在考察BDNF基因与环境交互作用年龄差异的基础上，进一步检验其交互作用模式随年龄的动态变化。 综上，青少年期的个体面临一系列生理、环境经历的变化，这不仅增加了个体的抑郁易感性，并且可能使得BDNF基因及其与环境交互作用模式呈现随年龄的动态变化。本研究目的之一是追踪童年期向青少年期过渡的个体，考察BDNF基因与同伴接纳、同伴拒绝对青少年抑郁影响的年龄差异。基于脑内BDNF含量的倒U型发展模式，本研究假设在青少年早期，Met等位基因对抑郁的易感性可能逐渐减弱，而Val等位基因的易感性逐渐显现，且其交互作用模式会随年龄产生动态变化。 1.4BDNF基因与环境因素对青少年抑郁发展轨迹的影响 青少年早期不仅是遗传效应动态变化的时期，也是青少年抑郁发展的关键期。纵向研究显示，抑郁发端于童年期，但是直到青少年早期(12岁左右)，抑郁发病率才开始迅速上升，并在青少年中期(15岁左右)达到顶峰(Avenevolietal.，2015；Hankinetal.，2015)。由此可见，遗传效应动态变化与抑郁发展变化的关键期相重叠，这提示青少年期遗传效应动态变化可能与抑郁的发展轨迹密切相关。 已有研究显示，遗传基因及其与环境的交互作用既影响青少年抑郁的初始水平，也可以解释抑郁、孤独感等情绪障碍随年龄的发展变化趋势(如Nelemansetal.，2021；vanRoekeletal.，2011)。譬如，Nelemans等(2021)发现重性抑郁风险基因累加得分与父母批评交互影响青少年抑郁的初始水平，也预测抑郁增长速度。具体表现为，携带更多风险基因并且经历更多父母批评的青少年不仅具有更高的抑郁症状，并且其抑郁增长速度更快。此外，Petersen等(2012)发现5-HTTLPR低转录效能基因型携带者比高转录效能基因型携带者在16～17岁时具有更高的抑郁加速发展模式，而与其他年龄段的抑郁发展速度和加速度无关。据我们所知，目前仅有一项研究考察BDNF基因与青少年孤独感发展轨迹的关系(Verhagenetal.，2014)。追踪研究发现，孤独感不仅能够预测个体的抑郁风险，并且与抑郁具有相似的易感因素(Limetal.，2020)，由此可推测BDNF基因及其与环境的交互作用亦可能解释抑郁发展轨迹。 发展级联模型(developmentalcascades)能够为理解BDNF基因与环境交互影响抑郁发展轨迹提供启示。发展级联模型指出，多种发展系统间存在跨时间的交互作用和相互作用，随着年龄的增长，不同水平(如生理水平和行为水平)、同水平不同领域(如认知领域和社会性领域)以及不同系统(如家庭系统和同伴系统)或世代(如父母和子女)的心理社会因素间发生的多重相互作用并产生累积效应，引起跨水平、跨领域甚至跨系统的扩散效应(Masten&Cicchetti，2010)。根据发展级联模型，在生命早期发展阶段，携带BDNF易感基因型的个体会具有更高的抑郁初始水平(如Hiltetal.，2007)；其高水平的抑郁症状将进一步唤起其他影响情绪发展的消极经历，如同伴侵害、同伴拒绝等(Beesonetal.，2020)；消极环境经历不仅会直接增加个体的抑郁症状，并且会继续增加个体对消极环境的生物敏感性(Pluess，2015)，个体在后续发展过程中则更易受到不利环境的消极影响，增加抑郁的风险。由此，随着时间的推移，具有遗传易感性的个体在发展过程中抑郁的多种风险因素不断累积，其与不携带遗传易感性的个体在消极经历、消极环境敏感性及抑郁水平间的差异也越来越大，进而表现出更高的抑郁增长速度。 基于此，BDNF基因及其与环境的交互作用不仅会影响青少年抑郁的初始水平，并且影响抑郁的发展速度。本研究目的之二是考察BDNF基因与同伴关系对青少年抑郁发展轨迹的影响。基于发展级联模型，本研究假设，携带BDNF风险基因的个体不仅具有更高的抑郁初始水平，并且具有更高的抑郁增长速度，且该差异在个体经历高风险环境时进一步扩大。 总体来说，本研究基于动态发展视角，考察BDNF基因与同伴接纳、同伴拒绝对青少年抑郁的动态影响，主要探讨以下问题：(1)BDNF基因与同伴关系对青少年抑郁影响的年龄差异模式；(2)BDNF基因及其与同伴关系的交互作用对青少年抑郁初始水平和发展速度的影响。 2研究方法 2.1被试 本研究为国内一项儿童青少年发展大型追踪项目的子研究。由于本研究关注青少年早期这一遗传效应动态变化的关键阶段，基于Steinberg(2022)对青少年期各阶段的划分，选取该追踪项目中的青少年早期及其向青少年中期过渡的3个时间点(12～14岁)。被试为在山东省济南市14所小学整群抽取的六年级学生1090名，在每年春季学期进行一次问卷测查，连续追踪3年，其中青少年抑郁和同伴关系在3个时间点均进行测评，在第三年采集被试唾液样本进行基因分型。其中，4人DNA含量过低无法分型而缺失基因数据。在1086名具有基因数据的被试中，被试平均年龄12.32±0.50岁，男生和女生各半，每年因转学、请假等流失被试分别为12人、8人和34人，Little'sMCAR检验不显著(=58.91，df=49，p=0.16)，且流失被试和具有全部数据的被试在性别(=3.61，df=1，p=0.06)、年龄(t(1063)=0.27，p=0.78)、民族(=0.20，df=1，p=0.65)、基因型(=2.85，df=2，p=0.24)等人口学变量上亦不存在显著差异，表明被试为随机缺失。 2.2研究工具 2.2.1抑郁 采用儿童抑郁量表(Children'sDepressionInventory，CDI；Kovacs，1992)测量被试抑郁症状。该量表共包含27个项目，要求被试报告近两周内的抑郁症状，分别采用0、1、2记分，均分越高则表明青少年抑郁症状越多。该量表在正常青少年群体中应用广泛(Zhangetal.，2015)，具有良好的测量学特性。本研究中3个时间点抑郁的Cronbach'sα系数均为0.88。 2.2.2同伴关系 采用同伴提名法测评青少年的同伴拒绝和同伴接纳。要求被试写出班里最不喜欢的三位同学(消极提名)和最喜欢的三位同学(积极提名)。对每个被试的被提名次数进行班级内标准化，以便于不同班级青少年提名分数的比较。消极提名的标准分是同伴拒绝的指标，积极提名的标准分是同伴接纳的指标。 2.2.3DNA提取与分型 采用Sequenom(SanDiego，CA，USA)芯片基质辅助激光解吸/电离飞行时间(MALDI-TOF)质谱平台从唾液样本中提取DNA，并对BDNF基因rs6265多态性位点进行基因分型。PCR引物为：forward5’-TCAAGAGGCTTGACATCATTGG-3’；reverse5’-GCCGAAC3TTTCTGGTCCTCAT-3’。PCR反应条件：94℃15min；94℃20s，56℃30s，72℃1min，共45个循环；最终72℃=3min。随后经历单碱基延伸反应，基因分型采用MassARRAYTyper3.4软件系统进行分析。本研究所使用的检测平台和技术具有较高的可靠性(基因分型有效率＞97％)。 2.3研究程序 研究经过山东师范大学伦理委员会审核批准。首先，项目组将问卷测查、唾液样品采集、DNA提取与分型流程等信息告知施测学校负责人、监护人和青少年，获得三方的知情同意后进行了数据收集。其次，以班级为单位，由经过严格培训的教师和研究生担任主试，采用现场问卷调查的方式测评青少年抑郁症状、同伴拒绝和同伴接纳，施测结束后问卷当场收回。第三，在第三次问卷测评结束后一周内以班级为单位采集被试唾液样本，每人采集2ml以上。 2.4数据处理与分析 本研究中同伴关系(同伴接纳和同伴拒绝)与青少年抑郁在T1～T3三个时间点上进行追踪测评，在T3时，进行被试唾液采集并对BDNF基因进行分型。基于纵向数据，本研究统计分析如下： 第一，对T1～T3的同伴接纳、同伴拒绝和青少年抑郁进行描述统计、相关分析，并对BDNF基因进行Hardy-Weinberg平衡检验。后续分析中对BDNF基因型进行虚拟编码，以ValMet为参照组，BDNF1=MetMetvs.ValMet.BDNF2=ValValvs.ValMet。 第二，采用SPSS23.0进行分层线性回归分析。分别以T1～T3的抑郁为因变量，以性别为控制变量，以BDNF基因型(以ValMet基因型为参照组形成的两个虚拟变量)、各时间点的同伴关系、BDNF基因型与各时间点同伴关系的交互项为预测变量，在3个时间点考察BDNF基因与同伴关系对青少年抑郁的同时性影响并比较3个时间点上结果的差异，若任一时间点的基因和环境的交互项显著，则通过简单斜率检验进一步分析其交互作用模式。为进一步检验本研究发现的年龄差异效应，参照Odgers等(2008)研究，合并3个时间点的数据，将时间点作为预测变量(T1=0，T2=1，T3=2)，考察基因、环境与时间点的三者交互作用。 第三，采用再参数化回归(Widamanetal.，2012)检验BDNF基因与同伴关系的交互作用符合何种理论模型(“不同易感性”“素质-压力”和“优势敏感”模型)。再参数化回归公式如下，D1代表ValMet基因型(1=ValMet，0=MetMet和ValVal)，D2代表MetMet基因型(1=MetMet，0=ValMet和ValVal)，D3代表ValVal基因型(1=ValVal，0=MetMet和ValMet)，Peer代表同伴关系，Gender代表性别，C代表交叉点，BO代表截距，B1、B2、B3分别代表三种基因型中同伴关系的回归系数，B4代表性别的回归系数。 第四，采用Mplus6.0，构建无条件潜变量增长模型(LatentGrowthCurveModeling，LGCM)分析抑郁的一般发展趋势以及抑郁起始水平和发展速度是否存在显著的个体差异。 第五，在无条件模型(Model1)基础上，逐次增加性别(Model2)、BDNF基因型(Model3)、T1同伴关系(Model4)、BDNF基因与T1同伴关系的交互项(Model5)为预测变量，考察其对抑郁发展截距和斜率的预测作用。上述潜变量增长模型采用全息最大似然估计，并分别构建同伴拒绝和同伴接纳模型，以区分两者的差异。 第六，进行敏感性分析，检验结果的可靠性。其一，采用同伴地位()为环境指标重复上述检验。其二，采用平行潜变量增长模型，考察同伴关系的变化与青少年抑郁变化间的关系。 3结果 3.1描述统计 各变量均值、标准差和相关系数见下页表1。青少年抑郁呈上升趋势，青少年抑郁与同伴接纳在不同时间点均呈显著负相关，青少年抑郁与同伴拒绝在不同时间点均呈显著正相关，抑郁、同伴接纳和同伴拒绝均具有中等程度的稳定性。重复测量方差分析显示，青少年抑郁随年龄呈增长趋势(F(2，1036)=68.80，p＜0.001)，同伴拒绝呈下降趋势(F(2，1081)=4.43，p=0.01)，但同伴接纳保持稳定(F(2，1081)=1.04，p=0.35)。独立样本t检验显示，在3个时间点上，女生均具有更低的同伴拒绝水平(ts(1083)≥4.52，ps＜0.001)，但在抑郁和同伴接纳水平上不存在性别差异(ts(1072)≤1.74，ps≥0.08)。BDNF基因型分布为MetMet(252人，23.2％)，ValMet(533人，49.1％)，ValVal(301人，27.7％)；基因型分布符合Hardy-Weinberg平衡(=0.29，df=2，p=0.59)。方差分析显示，不同基因型在各时间点同伴接纳、拒绝水平上均无显著差异(Fs(2，1081)≤2.40，ps≥0.09)。 3.2BDNF基因与同伴关系对青少年抑郁的交互作用及其年龄差异 以性别、BDNF基因、同伴关系以及基因与同伴关系交互项为预测变量，分别对T1～T3的抑郁进行同时性预测(见表2和表3)，以检验基因与环境交互作用的年龄差异。结果显示，在3个时间点上，BDNF基因型的主效应均不显著，同伴拒绝均能显著正向预测青少年抑郁。BDNF基因型与同伴拒绝的交互作用对青少年抑郁的同时性影响均显著，但是其交互作用模式存在差异。T1(12岁)时，相比ValMet基因型，在MetMet基因型携带者中，同伴拒绝对青少年抑郁的预测力更强；T2(13岁)时，相比ValMet基因型，在MetMet和ValVal基因型携带者中，同伴拒绝对青少年抑郁的预测力更强；T3(14岁)时，相比ValMet基因型，在ValVal基因型携带者中，同伴拒绝对青少年抑郁的预测力更强。在同伴接纳模型中，在3个时间点上，同伴接纳均能显著负向预测青少年抑郁，但是BDNF基因型的主效应及其与同伴接纳的交互作用均不显著。 进一步的简单斜率检验显示(见图1)，T1时，在MetMet(b=0.06，t=3.96，p＜0.001)和ValVal基因型(b=0.03，t=2.48，p=0.01)携带者中，同伴拒绝显著正向预测青少年抑郁，但在ValMet基因型中，同伴拒绝对青少年抑郁的预测作用不显著(b=0.01，t=1.16，p=0.24)。T2时，在MetMet(b=0.09，t=5.45，p＜0.001)、ValVal(b=0.06，t=4.52，p＜0.001)和ValMet基因型(b=0.02，t=2.06，p=0.04)携带者中，同伴拒绝均显著正向预测青少年抑郁。T3时，在MetMet(b=0.04，t=2.34，p=0.02)和ValVal基因型(b=0.07，t=4.08，p＜0.001)携带者中，同伴拒绝显著正向预测青少年抑郁，但在ValMet基因型中，同伴拒绝对青少年抑郁的预测作用不显著(b=0.004，t=0.41，p=0.68)。 图1BDNF基因与同伴拒绝对青少年抑郁的交互作用图(T1～T3) 同时性回归分析显示，BDNF基因与同伴拒绝交互作用模式在不同的年龄存在差异，并且从简单斜率检验可知12～14岁，ValVal基因型和MetMet基因型对环境的敏感性存在较大的变化，为进一步比较MetMet和ValVal对环境的敏感性差异随时间的变化。本研究以MetMet为参照组重编码基因，参照Odgers等(2008)研究，将时间点作为预测变量，考察基因、环境与时间点的三者交互作用。结果显示，BDNF基因型(ValValvs.)、同伴拒绝与时间点的三者交互作用显著(b=0.03，t=2.01，p=0.045)，这提示ValVal和MetMet基因型对环境的敏感性差异随年龄而变化。但是BDNF基因型(MetMetvs.)、同伴拒绝与时间点的三者交互(b=-0.01，t=-0.81，p=0.42)，BDNF基因型(ValValvs.)、同伴拒绝与时间点的三者交互(b=0.02，t=1.58，p=0.11)未达到显著水平。 3.3基因与环境交互模式检验：素质-压力、优势敏感与不同易感模型 基于回归分析结果，采用再参数化回归检验BDNF基因与同伴拒绝的交互作用符合何种理论模型。T1同伴拒绝与BDNF基因对抑郁的再参数化分析(表4)显示，交叉点C及其95％CI均在同伴拒绝全距内([-0.74，6.74])，限定ValMet回归系数为0后，Modela与Modelb在解释率上不存在显著差异，且Modela的AIC和BIC值最小，强不同易感模型为最优模型。T2同伴拒绝与BDNF基因对抑郁的再参数化分析(见表5)显示，交叉点C及其95％CI均在同伴拒绝全距内([-1.01，6.01])，限定ValMet基因型回归系数为0后，Modela与Modelb在解释率上存在显著差异，且Modelb的AIC值最小，弱不同易感模型为最优模型。T3同伴拒绝与BDNF基因对抑郁的再参数化分析(见表6)显示，交叉点C及其95％CI超出同伴拒绝全距内([-0.80，6.04])，该基因-环境交互不符合不同易感模型，优势敏感模型解释率与不同易感模型不存在显著差异，且Modele的AIC和BIC更小，强优势敏感模型为最优模型。 3.4抑郁的发展趋势：BDNF基因与同伴关系的交互作用 首先，对青少年抑郁构建无条件潜变量增长模型，如表7(68页)所示，无条件模型拟合良好。从12至14岁，青少年抑郁呈线性增长趋势，截距(=0.03，p＜0.001)和斜率(=0.01，p＜0.001)的变异均呈现出明显的个体间差异。截距和斜率之间的相关不显著(r=-0.16，p=0.053)，表明抑郁增长速度与初始水平无关。 在无条件模型基础上，以截距和斜率为结果变量，逐次增加性别(Model2)、基因型(Model3)、T1同伴关系(Model4)，BDNF基因与T1同伴关系的交互项(Model5)等预测变量，建立条件模型。结果显示，BDNF1(MetMetvs.ValMet)和BDNF2(ValValvs.ValMet)对斜率的预测作用显著，相比ValMet基因型，Met和Val纯合子携带者的抑郁增长速度更快(见图2)。 图2不同BDNF基因型青少年抑郁的发展变化 如表7和图3所示，同伴拒绝模型(Model5a)拟合良好，BDNF1与T1拒绝的交互项对截距的预测显著。如表7和图4所示，同伴接纳模型(Model5b)拟合良好，但该模型未发现BDNF基因与接纳交互作用对截距或斜率具有显著预测作用。此外，多群组比较显示，条件模型无性别差异(△=22.03，△df=15，p=0.12)。 图3BDNF基因与同伴拒绝对青少年早期抑郁发展模式的影响(Model5a) 图4BDNF基因与同伴接纳对青少年早期抑郁发展模式的影响(Model5b) 3.5敏感性分析 第一，采用同伴地位(Z同伴接纳-Z同伴拒绝)为环境指标重复上述检验。同时性回归分析(见原刊网络版附录表A1)显示，在3个时间点上，BDNF基因型与同伴地位的交互作用对青少年抑郁的同时性影响均显著。12～13岁时，相比ValMet基因型，在MetMet基因型携带者中，同伴拒绝对青少年抑郁的预测力更强；14岁时，相比ValMet基因型，在ValVal基因型携带者中，同伴拒绝对青少年抑郁的预测力更强。潜变量增长模型(见表7)显示，同伴地位模型(Model5c)拟合良好，仅BDNF1与T1同伴地位的交互项对截距的预测显著。在T1～T2时，BDNF基因与同伴地位的交互作用更符合不同易感性模型，而在T3时，优势敏感模型能够很好地解释基因×环境交互作用模式(见网络版附录表A2～A4)。 第二，采用平行潜变量增长模型，考察同伴关系的变化与青少年抑郁变化间的关系。结果显示(见网络版附录表A5)，同伴拒绝、同伴接纳和同伴地位的初始水平均显著预测抑郁的初始水平，而非抑郁的增长速度，同伴关系的增长速度与青少年抑郁的增长速度无关。多群组比较分析显示，不同基因型携带者仅在同伴拒绝初始水平到青少年抑郁初始水平这一路径上存在显著差异(△=16.01，△df=2，p＜0.001)，与上述分析结果相一致。 4讨论 本研究基于动态发展视角，考察了BDNF基因与同伴关系对青少年抑郁交互作用的年龄差异，以及该基因及其与环境交互作用对青少年抑郁一般发展轨迹的影响。结果显示，12～14岁，BDNF基因与同伴拒绝的交互作用均显著同时性预测青少年抑郁，但是其交互作用模式存在年龄差异。具体表现为，12岁时，MetMet基因型携带者对环境敏感性高于ValMet基因型携带者；13岁时，MetMet和ValVal基因型携带者对环境的敏感性均高于ValMet基因型携带者；14岁时，ValVal基因型携带者对环境的敏感性高于ValMet基因型携带者。此外，BDNF基因与同伴拒绝交互预测青少年抑郁的初始水平，相比ValMet基因型携带者，携带MetMet基因型的青少年在经历同伴拒绝后抑郁初始水平更高；BDNF基因显著预测青少年抑郁的增长速度，相比ValMet基因型，MetMet和ValVal基因型携带者的抑郁增长速度更快。 各时间点的同时性回归分析均发现BDNF基因与同伴拒绝交互预测青少年抑郁，携带易感基因的青少年在经历同伴拒绝时报告了更多抑郁症状。在12～14岁，个体对环境的易感基因型从MetMet基因型逐渐转变为ValVal基因型。该研究结果能够为以往研究在BDNF基因与环境交互结果上的分歧提供启示。如前所述，以成年人或儿童为被试的研究显示Met等位基因携带者具有更高的抑郁水平并且对环境更加敏感(Aguileraetal.，2009；Meyeretal.，2018)，而以青少年为被试的研究却显示Val等位基因携带者表现出对环境更高的敏感性(Chenetal.，2013)。针对该分歧，已有研究经常提及但鲜少直接检验BDNF基因效应的年龄差异(Caseyetal.，2009)，本研究则为这一发展动态性模式提供了支持。此外，除BDNF基因外，已有关于5-HTYLPR(Hankinetal.，2015)、DRD2基因(Zhangetal.，2015)的研究亦发现了类似的年龄差异。譬如，Zhang等(2015)的研究考察了DRD2基因与母亲教养行为影响抑郁的年龄差异，在个体11和12岁时，发现携带A1等位基因的个体对消极教养行为更加敏感，但在13岁时，该基因与教养的交互作用不显著。这提示基因及其与环境交互作用的发展动态性并不罕见，尤其是在发展过渡的关键阶段(如本研究的青少年早期)更容易观测到基因及其与环境交互作用的变化。但是，本研究基于较短的年龄跨度且在没有外部验证样本的情况下，所观测到BDNF基因与同伴拒绝交互作用随时间的变化仍可能是一种假阳性结果，未来仍需要大量研究重复验证该结果的可靠性。 如前所述，BDNF基因及其与环境交互作用随年龄的动态变化既可能与BDNF含量的发展变化有关，也可能与环境因素的动态变化有关。第一，来自动物和人类研究的发现均显示，从童年期到青少年期，BDNF浓度呈增长趋势(如Ivanova&Beyer，2001；Websteretal.，2002)。因而，不同年龄阶段基底BDNF浓度使得BDNF基因功能存在动态变化。在童年晚期或前青少年期，基底BDNF浓度较低，高活性Val等位基因能够缓冲BDNF低浓度的风险性，而低活性Met等位基因则表现出风险性；而随着年龄增长基底BDNF含量上升，使得Met与Val基因型在BDNF含量上的差异被掩盖，因而低活性Met等位基因的风险性逐渐减弱并消失。更重要的是，除了Met等位基因的风险性随年龄增长而消失，本研究发现Val等位基因风险性逐渐显现出来。这提示BDNF功能与抑郁间的关系可能呈现非线性趋势。换言之，过高或过低的BDNF功能均可能增加抑郁的风险性。因而，不能仅仅依据BDNF基因活性判断等位基因是“好”或“坏”，还需要考虑BDNF基因在BDNF功能与抑郁U型函数上所处的位置。第二，在童年期和青少年期，压力易感的神经生理机制差异亦可能影响BDNF基因与环境交互作用的年龄差异。研究显示，海马区发育成熟时间早于杏仁核(Andersen&Teicher，2008)，因而童年期海马对压力更加敏感，青少年期压力易感脑区是杏仁核(Lupienetal.，2009)。暴露于压力环境会导致海马区BDNF含量降低，而杏仁核BDNF含量上升(Changetal.，2021)。更重要的是，BDNF在海马和杏仁核的功能截然相反，海马区较低的BDNF含量与较高的抑郁水平有关(Duman&Monteggia，2006)，而较高的BDNF含量会导致较大的杏仁核体积进而增加抑郁风险(Thomasetal.，2001)。基于此，可推测在童年期，低活性Met等位基因与压力共同导致海马体更低的BDNF含量，导致抑郁风险增加；而在青少年期，高活性Val等位基因与压力共同导致杏仁核更高的BDNF含量，增加抑郁风险。第三，本研究被试处于“小升初”的生态过渡期，青少年面临着从熟悉同伴群体过渡到陌生同伴群体，这也是个体同伴地位重组的关键时期。本研究也发现随着年龄增长，青少年同伴拒绝水平呈现下降趋势，如前所述环境变化能够改变BDNF甲基化水平(Kageyamaetal.，2022)，由此可推测环境变化引起的易感基因DNA甲基化会使该基因失效从而掩盖风险等位基因的消极影响，呈现出BDNF基因与同伴拒绝交互作用的动态变化。 基因-环境交互作用的发展动态性不仅体现在易感基因型的变化上，也体现在其基因与环境交互作用模式的变化上。再参数化分析显示，在12～13岁时，BDNF基因与同伴拒绝的交互作用模式更符合不同易感模型，即携带易感基因型(ValVal和MetMet基因型)的个体对积极和消极环境均更加敏感，而表现出更低或更高的抑郁水平。在14岁时，相比不同易感模型，优势敏感模型能够更好地解释BDNF基因与同伴拒绝的交互模式。如前所述，个体的环境敏感性会受早期环境经历影响而产生动态变化，携带环境敏感性基因的个体在生命早期对积极和消极环境均具有敏感性(即不同易感性模式)，但当这些个体在积极或消极环境中成长时，则其环境敏感性则会向优势敏感或素质压力模式发展(Pluess，2015)。基于该理论假设，本研究发现青少年同伴拒绝水平随年龄逐渐降低，而同伴接纳水平保持稳定，即个体的积极同伴经历更为显著，因而在相对积极环境中成长起来的青少年，其基因-环境交互作用模式呈现从不同易感性向优势敏感性的转变。值得指出的是，BDNF基因与同伴拒绝而非同伴接纳存在交互作用，该基因-环境交互作用模式的转变可能受到同伴拒绝有限取值范围的影响，因而本研究补充分析了同伴地位(同伴社会偏好)的结果，亦发现14岁时优势敏感模型能够很好地解释基因-环境交互模式。这一结果也提示早期的环境经历可能会影响个体对未来环境的反应，但是由于本研究仅在消极环境中发现了个体的环境敏感性从不同易感模式向优势敏感模型的变化，应该对本研究结果保持谨慎，未来研究应该进一步关注个体环境敏感性的发展变化，尤其是早期环境经历对个体环境敏感性的塑造作用。 无论是易感基因型还是基因-环境交互模式随年龄的变化均为“天性与教养的动态交互观”(Lerner，1987)和“基因-环境-时间”交互理论框架(Boyceetal.，2020，2021)提供了支持。尽管在个体发展过程中DNA序列不会发生变化，但是遗传基因的效应并非是静态的，生理成熟因素、环境经历等均会通过复杂的相互作用机制影响基因功能的表达，并且提示未来研究应该深入考察基因-环境的动态交互的基础理论问题。首先，本研究仅采用年龄作为时间(发展)因素，但是发展时程(timing，如青春期发展阶段)、神经环路发展关键期(不同脑区的发展关键期存在差异)、世代变迁(如代际变迁、进化)等均是揭示基因-环境的动态变化的时间因素(Boyceetal.，2020)。不同类型的时间因素可能揭示了不同的动态发展机制，譬如，环境因素(如营养和社会情绪因素)会延长或缩短特定发展阶段(如青春期)所需的时间(Ellisetal.，2011)，使得发育阶段和年龄无关，因而产生不同的基因-环境-时间交互模式。与此一致，Hankin等(2015)有关青少年抑郁的研究发现基因-环境-年龄存在交互作用，但青春期发展阶段对基因-环境交互无调节作用。由此，未来研究应该进一步扩展对时间(发展)因素的探讨以更全面描述“基因-环境-时间”交互理论框架。其次，探索基因与环境交互作用在时间轴上的动态展开机制是未来研究的重点问题。根据Cole(2009)的递归发展(recursivedevelopment)模型，表观遗传过程可能是解释遗传-环境交互效应发展动态性的关键机制。该模型指出某一时间点的环境因素可以通过中枢神经系统影响个体的行为表型，并通过表观遗传机制影响RNA转录进而塑造出下一个时间点新的神经反应功能，新塑造的神经反应功能对环境因素的反应则会与上一时间点存在差异，进而导致新的行为表型和新的RNA转录过程，以此循环递归产生动态变化。事实上，这一模型与发展级联模型的观点相似，均提示不同水平的影响因素跨时间的相互作用是解释个体动态发展的关键机制，但该模型更明确地指出基因表达机制在其中的作用。由此，未来研究应该在发展级联模型或递归发展模型框架下，采用表观遗传研究范式，在基因表达水平上，揭示遗传-环境交互作用动态变化的潜在机制。 本研究并未发现BDNF基因与同伴接纳间存在显著的交互作用，这可能与社会接纳和拒绝在脑内产生的不同影响有关。研究显示，社会接纳会导致伏隔核的内源性阿片释放增加，而社会拒绝则会导致调节压力、情绪和动机脑区的内源性阿片释放降低(Hsuetal.，2015)。同时，动物研究显示注射δ-阿片受体激动剂(其效应类似于阿片含量增加)的大鼠额叶皮质、海马和基底外侧杏仁核中BDNFmRNA表达增加(Torregrossaetal.，2004)。由此，可推测同伴接纳会增加内源性阿片含量，通过增加BDNFmRNA表达，提高脑内BDNF浓度而产生“天花板效应”，使得Val和Met等位基因间的差异消失，因而无法观察到基因-环境交互作用。此外，这也可能与青少年对消极环境的敏感性高于积极环境有关。研究发现，在9～17岁儿童青少年中，面对同伴拒绝时个体瞳孔扩张反应(反映情绪调节脑区活动程度的指标)比面对同伴接纳时更大，并且该差异随年龄增长而增加(Silketal.，2012)。 与以往研究结果(Avenevolietal.，2015)一致，无条件模型显示青少年早期抑郁呈线性增长趋势。该时期抑郁的迅速增加可能源于青少年早期个体面临一系列生理、认知、社会角色过渡，不仅使青少年具有更高的压力易感性，并且伴随着压力性经历增多，加剧了青少年罹患抑郁的风险(Hyde&Mezulis，2020)。此外，本研究还发现抑郁的初始水平与其发展速度并无显著相关。这提示，无论个体进入青春期时的初始抑郁水平高还是低，均可能具有更高或更低的抑郁增长速度。因此，初始抑郁水平较低并不意味着个体具有低抑郁风险，高抑郁初始水平也可能因保护因素而降低其抑郁的增长速度。更重要的是，本研究显示青少年抑郁的初始水平和增长速度存在明显的个体差异，BDNF基因及其与环境的交互作用部分解释了该个体差异。 本质上，青少年抑郁初始水平的影响因素与起始时间点的同时性预测因素相一致。即BDNF基因型(MetMetvs.ValMet)与同伴拒绝的交互作用显著预测抑郁的初始水平，但BDNF基因主效应与抑郁初始水平无关，对此结果不再赘述。然而，与抑郁的初始水平不同，BDNF基因能够预测青少年抑郁的发展速度。这提示，影响抑郁的遗传差异从童年晚期到青少年早期逐渐扩大。这一现象与双生子研究结果相似，从童年期到青少年期，抑郁的遗传率逐渐增加(Scourfieldetal.，2003)。抑郁的遗传差异随着年龄逐渐扩大或显现的原因可能有两种：第一，在青少年期出现了导致BDNF基因遗传效应增大的新兴神经生理过程或遗传基因表达。如前所述，在青少年期BDNF浓度的新变化会改变BDNF基因功能，使得个体的易感基因发生变化，进而影响了抑郁的发展速度。这可能是一些发展敏感的“Clock”基因在青春期开始启动(Picklesetal.，1998)，不仅引起BDNF浓度变化导致BDNF基因效应差异逐渐扩大，并且使得在青少年期观察到新兴的遗传效应，导致抑郁遗传率的增加。第二，在发展过程中，生命早期较小的遗传差异会像滚雪球一样越来越大。根据发展级联模型，在生命早期发展阶段，个体所携带BDNF基因型的差异会导致其对消极环境的敏感性存在差异，遭