基因-環(huán)境交互作用理論模型及其檢驗方法*

王美萍

(山東師范大學心理學院, 濟南 250014)

揭示心理與行為的發(fā)生機制是心理學的經(jīng)典和重要研究任務。20世紀60年代以來, 伴隨著行為遺傳學的興起, 遺傳與環(huán)境的影響及其作用機制問題再次引起心理學研究者的關注。近十年來分子遺傳學技術的發(fā)展, 則使基因(Gene, G)×環(huán)境(Environment, E)交互作用的探究成為可能, 并日益發(fā)展成為心理學、精神病理學、行為遺傳學等領域的前沿和熱點研究問題之一。時至今日,主導基因-環(huán)境交互作用研究的理論模型主要包括素質(zhì)-壓力模型(diathesis-stress model)、不同易感性模型(differential susceptibility model)和優(yōu)勢敏感性模型(vantage sensitivity model)三種。目前,心理學界尚沒有關于這些模型的系統(tǒng)介紹。基于此, 本文擬對這三種模型的主要觀點、檢驗方法進行回顧與梳理, 并在此基礎上展望未來研究方向, 以期為后續(xù)相關研究提供參考。

1 主要理論模型的基本觀點

1.1 素質(zhì)-壓力模型(diathesis-stress model)

“素質(zhì)(diathesis)”源于希臘語中的傾向(disposition)一詞, 通常指個體先天的生物或遺傳特征(如氣質(zhì)、神經(jīng)系統(tǒng)的敏感性、基因等)。素質(zhì)-壓力模型最初由心理學家Bleuler (1963)和Rosenthal (1963)于20世紀60年代在有關精神分裂癥的理論中提出。其核心觀點是(Monroe & Simons, 1991;Zuckerman, 1999), 攜帶某種“風險” (risk)或者“易受” (vulnerability)基因的個體更容易受到不良環(huán)境因素的影響而產(chǎn)生心理或行問題, 即“問題”的產(chǎn)生是由個體的風險素質(zhì)和消極環(huán)境經(jīng)歷共同作用的結果, 因此該模型又被稱為雙重風險模型(dual-risk model, Sameroff, 1983)。在該理論模型的影響下, 早期的基因-環(huán)境交互作用研究主要關注童年期受虐經(jīng)歷(Caspi et al., 2002)、壓力性生活事件(Caspi et al., 2003)、早期創(chuàng)傷(Frazzetto et al., 2007)、母親孕期吸煙(Neuman et al., 2007)、感知到的種族歧視(Enoch, Steer, Newman, Gbison,& Goldman, 2009)等消極環(huán)境因素與基因對精神分裂癥(Walker & Diforio, 1997)、多動癥(Neuman et al., 2007)、反社會行為(Caspi et al., 2002)、攻擊(Frazzetto et al., 2007)、抑郁(Caspi et al., 2003;Taylor et al., 2006)等心理或行為問題的影響。譬如, Caspi等人(2002)的首創(chuàng)性研究考察了單胺氧化酶A基因上游可變串聯(lián)重復序列(monoamine oxidase A gene-upstream variable number of tandem repeat,MAOA-uVNTR)多態(tài)性與兒童期虐待(消極環(huán)境)對反社會行為(不良行為)的交互作用, 結果發(fā)現(xiàn)在那些攜帶低活性基因型的被試中, 童年期遭受過虐待的個體比那些沒有受虐經(jīng)歷的個體更可能表現(xiàn)出反社會行為, 但在那些攜帶高活性基因型的被試中, 受虐經(jīng)歷對反社會行為無顯著影響。在素質(zhì)-壓力模型的理論框架下, 個體通常被區(qū)分為易受性個體(vulnerability individual)和有韌性或者耐挫的個體(resilient individual), 而在積極環(huán)境條件下, 攜帶不同基因型個體的心理與行為表現(xiàn)的差異問題受到了忽略, 這與心理學長期以來居主導地位的消極研究取向(與積極心理學相對)是一脈相承的。

1.2 不同易感性模型(differential susce ptibility model)

不同易感性模型由Belsky基于進化論的思想提出(Belsky, 1997, 2005; Belsky, Bakermans-Kranenburg, & van IJzendoorn, 2007)。在他看來,個體的不同易感性是進化的產(chǎn)物, 自然選擇機制使得親代生育出對環(huán)境具有不同易感性的子代,因為在充滿不確定性、無法預測未來的情況下,這是最有利于后代繁衍的方式。此外, 不同易感性模型是在倡導研究人的積極心理品質(zhì)的積極心理學思潮興起的歷史背景下(20世紀末21世紀初)提出的, 因此二者之間應該有著必然的聯(lián)系。與素質(zhì)-壓力模型不同, 該模型認為, 由于某遺傳特征或者基因的影響, 那些容易受到消極環(huán)境的影響而表現(xiàn)出心理或行為問題的個體, 也容易受到積極環(huán)境的影響而表現(xiàn)優(yōu)秀或良好。即同一遺傳特征或基因型同時具有使個體變得更好(for better)或者更壞(for worse)的效應。因而該理論模型常使用易感性(susceptibility)、可塑性(plasticity或者malleability), 而非易受性(vulnerability)來描述基因。該模型觀點也得到了諸多研究的支持。譬如有研究表明, 當童年期所經(jīng)歷的負性事件(adversity)水平較高時, 那些攜帶低活性單胺氧化酶(monoamine oxidase A,MAOA)等位基因的男孩比攜帶高活性等位基因的男孩更可能表現(xiàn)出行為問題, 但是當經(jīng)歷的負性事件水平較低時, 前者表現(xiàn)出問題行為的可能性則低于后者(Foley et al.,2004)。與攜帶其他基因型的嬰兒相比, 那些攜帶多巴胺受體基因4(dopamine receptor D4,DRD4)7R(7-repeat)等位基因的嬰兒表現(xiàn)出較多(當母親未解決的失業(yè)或者創(chuàng)傷問題水平較高時)或者較少(當母親未解決的失業(yè)或者創(chuàng)傷問題水平較低時)的依戀問題(van IJzendoorn & Bakermans-Kranenburg,2006)。上述早期研究雖然也獲得了支持不同易感性模型的結果, 但令人遺憾的是, 此類研究仍然考察的是基因與消極環(huán)境變量對消極心理與行為發(fā)展變量的影響, 其中的積極發(fā)展環(huán)境實際上代表的是消極環(huán)境得分較低或者缺失的情況。其不足之處是顯而易見的?；诖? 近年來, 有些研究者(Belsky & Beaver, 2011; Bakermans-Kranenburg &van IJzendoorn, 2011)開始采用積極環(huán)境指標, 揭示基因與環(huán)境對個體積極發(fā)展結果的交互作用模式。譬如, Belsky和Beaver的研究(2011)顯示, 在支持性父母教養(yǎng)的條件下, 那些攜帶更多的可塑性等位基因(plasticity alleles)的男青少年的自我調(diào)控得分較高, 而在消極教養(yǎng)條件下, 這些男青少年的自我調(diào)控得分較低。簡言之, 不同易感性理論模型正受到越來越多的基因-環(huán)境交互作用研究的青睞。

1.3 優(yōu)勢反應敏感性模型(vantage se nsitivity model)

該模型是由Manuck和Sweitzer (Manuck, 2011;Sweitzer et al., 2013)新近提出的?！皟?yōu)勢敏感性”是與“不同易感性”相關但又不同的概念, 旨在描述個體對積極環(huán)境因素的不同反應性。Pluess和Belsky (2013)進一步將其詮釋為某些個體更容易受積極環(huán)境因素影響而表現(xiàn)良好的傾向, 而那些對積極環(huán)境因素反應不敏感的傾向則被稱為優(yōu)勢阻抗性(vantage resistance)。該模型也得到了幾項研究的支持。譬如Kegel等(Kegel, Bus, & van IJzendoorn, 2011)的研究顯示, 那些攜帶DRD4基因7R (7-repeat)等位基因且得到積極反饋的兒童表現(xiàn)出較高的讀寫能力, 但是在沒有得到積極反饋的條件下, 攜帶不同基因型的兒童的讀寫能力無顯著差異。目前, 驗證優(yōu)勢反應敏感性模型的研究多數(shù)只考察了積極環(huán)境因素的影響。

1.4 三種模型的比較

圖1 三種模型的基因-環(huán)境交互作用圖

素質(zhì)-壓力模型一直較為盛行, 不同易感性模型是新近興起的且當前較受推崇, 優(yōu)勢敏感性模型也是新興模型但尚未獲得關注。反應的三種模型都主張個體對環(huán)境的反應性存在差別, 有些個體容易受到環(huán)境因素的影響, 有些個體則對環(huán)境影響不敏感(見圖1)。其區(qū)別在于素質(zhì)-壓力模型認為攜帶某基因型的個體更容易受消極環(huán)境條件的影響而表現(xiàn)出更多的心理或行為問題, 但該模型忽視了積極環(huán)境的影響, 或者可以推定該模型認為在積極教養(yǎng)環(huán)境條件下, 不同基因型個體的心理或行為表現(xiàn)無顯著差異(Pluess & Belsky,2013); 優(yōu)勢反應敏感性模型則主張在積極環(huán)境條件下攜帶某基因型的個體比其他個體表現(xiàn)出更少的心理與行為問題, 或者發(fā)展更為良好, 但不關注消極環(huán)境因素的作用; 不同易感性模型則同時關注不同性質(zhì)環(huán)境因素的影響, 認為那些在消極環(huán)境條件下表現(xiàn)出更多問題的個體, 在積極環(huán)境條件下則會表現(xiàn)出更少的問題或者發(fā)展更為良好。簡言之, 素質(zhì)-壓力模型認為“近墨者黑”或“出淤泥而不染”, 優(yōu)勢反應敏感性模型認為“近朱者赤”, 而不同易感性模型則兼收并蓄, 認為某些個體近墨則黑, 近朱則赤。

由上述組圖很容易推測, 素質(zhì)-壓力模型與優(yōu)勢反應敏感性模型分別構成了不同易感性模型的左右兩部分, 因而很可能本質(zhì)上是屬于不同易感性模型, 但是這三者也可能確實是獨立存在的、不同的基因-環(huán)境交互作用模型。因為現(xiàn)實生活中個體對環(huán)境的反應模式可能確實存在差異,即有些個體容易受到消極環(huán)境的影響, 而對積極環(huán)境效應不敏感; 有些個體更容易受到積極環(huán)境的影響, 但是對消極環(huán)境反應不敏感; 而有些個體則可塑性較高, 對積極與消極環(huán)境影響均較為敏感。

2 基因-環(huán)境交互作用模型的檢驗方法

2.1 傳統(tǒng)檢驗方法

目前常用的、傳統(tǒng)的檢驗方法主要包括兩種：分組回歸和包含交互項的分層回歸。其中分層回歸法是應用最為廣泛的方法。在分組回歸中, 研究者(王美萍, 張文新, 2010, 2014)通常是根據(jù)一定的標準將基因型或者環(huán)境變量分為不同的組別,然后分別考察不同組別中, 環(huán)境或者基因效應是否存在差別, 如果存在則表明基因與環(huán)境存在交互效應, 反之, 則不然。而采用分層回歸方法的研究(Caspi et al., 2002, Belsky & Beaver, 2011;Benjet, Thompson, & Gotlib, 2010; van IJzendoorn,Bakermans-Kranenburg & Mesman, 2008)通常是構建基因型(一般依據(jù)攜帶“風險”或者“敏感性”基因的數(shù)量而分為三類或者兩類)與環(huán)境變量的交互項(即乘積項), 然后依次將環(huán)境(X1)與基因型(D2①注：D2為二分虛擬變量, 包括三個類別的虛擬變量的方程見Widaman et al., 2012。)的主效應、交互項(X1D2)納入方程(Y = B0+ B1X1+B2D2+ B3X1D2+ E), 若交互項能夠顯著預測因變量(即所考察的心理或行為變量), 則進行進一步的檢驗(follow-up analysis), 通常是簡單斜率檢驗(simple slope analysis), 以探明交互作用的具體表現(xiàn)形式。該檢驗方法一般以距離平均數(shù)正負一個標準差為標準, 將連續(xù)計分的環(huán)境變量分為高(高于平均數(shù)一個標準差)、低(低于平均數(shù)一個標準差)兩組, 然后考察不同環(huán)境條件下攜帶不同基因型個體的心理或行為變量的差異。

此外, Belsky等人(2007)還專門闡述了對不同易感性模型進行檢驗的步驟, 除了進行傳統(tǒng)的基因-環(huán)境交互效應的統(tǒng)計檢驗(回歸分析)外, 還包括排除基因-環(huán)境相關, 檢驗基因與因變量之間是否存在顯著關聯(lián), 比較基因-環(huán)境交互圖等四個環(huán)節(jié)。篇幅所限, 在此不再贅述。

2.2 新近檢驗方法

鑒于當所測量的環(huán)境指標的變異有限時, 采用傳統(tǒng)方法所發(fā)現(xiàn)的素質(zhì)-壓力交互作用模式可能實際上是屬于不同易感性交互作用模型的一部分, 為了更科學地區(qū)分所考察的基因-環(huán)境交互作用究竟支持二者中的哪一個, 近期有研究者(Roisman et al.,2012; Widaman et al., 2012)頗具創(chuàng)意地提出了檢驗素質(zhì)-壓力模型與不同易感性模型的新方法。

2.2.1 顯著性區(qū)域分析法

顯著性區(qū)域(region of significance, RoS)指的是自變量(基因)顯著預測因變量(心理或行為指標)時調(diào)節(jié)變量的取值范圍(Roisman et al., 2012)。與簡單斜率法相比, 顯著性區(qū)域分析方法考察的是自變量與因變量關聯(lián)顯著時調(diào)節(jié)變量所有的取值,而不是只分析了兩個點上的取值(M±1SD), 其優(yōu)勢顯而易見。當某基因與因變量在環(huán)境變量的低分端區(qū)域(消極環(huán)境), 而不是高分端區(qū)域(消極環(huán)境缺失或者積極環(huán)境)顯著關聯(lián)時, 則該基因-環(huán)境交互作用支持素質(zhì)-壓力模型, 如果在兩端區(qū)域均顯著關聯(lián), 則支持不同易感性模型。該方法的具體測量指標包括(1)交互作用比例指標(the proportion of interaction index, PoI)：該指標意味著在交互作用的區(qū)域中與不同易感性模型相符合的比例, 其計算公式為b/(b+w), w指在消極環(huán)境影響下的消極發(fā)展結果的比例, b則指在積極環(huán)境條件下的積極發(fā)展結果的比例。通常若0.40≤PoI≤0.60, 則說明基因-環(huán)境交互作用符合不同易感性模型, 典型的不同易感性模型中PoI值等于0.50,若PoI趨近0, 則意味著完全符合素質(zhì)-壓力模型,雖然Roisman等(2012)沒有明確說明當趨近于1時,交互作用模式支持哪一種理論模型, 但根據(jù)b與w的含義, 本文推論應符合前面所提及的優(yōu)勢敏感性模型。(2)受影響比例指標(the proportion affected index, PA)：該指標表示大于兩條回歸線的交叉點所對應的X值的比例, 即交叉點右側區(qū)域所占的比例。通常PA的取值范圍為0～0.50, 若PA近似為0.50 (0.40～0.60), 則交互作用模式支持不同易感性模型, 若近似為0, 則支持素質(zhì)-壓力模型, 若小于2%則可以確定該模型不支持不同易感性模型。鑒于PoI會隨環(huán)境變量的取值范圍的變化而變化,因此Roisman等人推薦使用M± 2SD, 而不是M± 1SD的區(qū)間。因為前者幾乎涵蓋了95%的取值區(qū)域,而后者僅涵蓋了68%, 所以前者能更準確地探測出基因-環(huán)境交互作用模式。該方法應用的前提條件是環(huán)境變量與因變量的關系符合線性假設, 因而在模型檢驗前需要考察X2和ZX2是否對因變量具有顯著預測作用, 如有, 則需要在控制其影響的前提下考察交互項是否依然顯著。

2.2.2 新參數(shù)回歸模型(Re-parameterized r egression model)法

該方法由Widaman等(2012)提出。為了克服共線性的問題, 傳統(tǒng)回歸分析法是以平均數(shù)(M)為標準進行中心化[Y = B0+ B1(X1- M) + B2D2+B3(X1- M) D2+ E], 新參數(shù)回歸模型法則是以回歸線的交叉點(crosspoint)為標準進行中心化[Y =A0+ B1(X1- C) + B3(X1- C) D2+ E] (Widaman et al.,2012)。該方程與前面提及的兩個方程是等價的,均擁有相同數(shù)目的自由參數(shù)和R2值, 其優(yōu)勢在于不僅能夠直接計算交叉點的值, 而且能夠估計出其區(qū)間范圍。與傳統(tǒng)檢驗方法以及顯著性區(qū)域方法不同的是, 該方法可以根據(jù)交叉點的值或者位置直接判斷交互作用的性質(zhì), 而不必對基因與環(huán)境的交互項進行顯著性檢驗及其后續(xù)分析。若G×E交互作用呈按序型(ordinal)交互作用, 且兩條回歸線的交叉點位于或接近于積極環(huán)境的一端(橫坐標的右端), 則該交互作用符合素質(zhì)-壓力模型。如果G×E交互作用呈非按序型(disordinal)交互作用, 且兩條回歸線的交叉點位于或接近于橫坐標的中間, 則該交互作用支持不同易感性模型。因此, 在新參數(shù)回歸模型檢驗方法中, 交叉點是關鍵參數(shù)之一。此外, Widaman等(Widaman et al.,2012; Belsky, Pluess, & Widaman, 2013)提出該方法可以同時對以下四種模型進行檢驗(1)弱的不同易感性模型; (2)強的不同易感性模型; (3)弱的素質(zhì)-壓力模型; (4)強的素質(zhì)-壓力模型。強的模型假設某些個體容易受環(huán)境的影響, 而另一些個體則完全不受環(huán)境的影響, 弱的模型則假設所有個體都受環(huán)境的影響, 但是一些個體比另一些個體受的影響更大。模型擬合度通過AIC和BIC (Akaike and Bayesian information criteria)來判斷。此外, 該方法也要求環(huán)境變量與因變量之間符合線性假設。

3 未來探究方向

3.1 理論問題的探討

目前無論是有關素質(zhì)-壓力模型的研究, 還是有關不同易感性模型或者優(yōu)勢反應敏感性模型的研究, 考察的幾乎均是某基因與某單一環(huán)境的交互作用, 而且環(huán)境指標主要集中于父母教養(yǎng)、生活事件、受虐經(jīng)歷等微系統(tǒng)或者近端環(huán)境變量,因而攜帶某種基因型的個體是對所有的環(huán)境都表現(xiàn)出某種反應模式(三種模型中的一種), 還是僅對某特殊環(huán)境表現(xiàn)出該反應模式的問題尚不清楚。此外, 不同基因與某一環(huán)境, 以及多個基因與多種環(huán)境的交互作用模式也有待探索(張衛(wèi), 甄霜菊, 2013)。未來研究可同時考察近端環(huán)境和遠端環(huán)境(例如居住地的經(jīng)濟發(fā)展水平、文化價值觀念等)等多種環(huán)境指標, 亦可采用多基因-環(huán)境設計,以探明不同的基因-環(huán)境交互作用模型是否具有領域特殊性, 以及其確切表現(xiàn)的問題。

既有研究顯示, 基因與環(huán)境的交互作用存在種族差異。近期的一項元分析進一步表明(van IJzendoorn, Belsky, & Bakermans-Kranenburg, 2012),在高加索人群中, 5-羥色胺轉運體啟動子區(qū)(serotonintransporter-linked polymorphic region, 5-HTTLPR)基因型與環(huán)境的交互作用符合不同易感性模型,在包括多種族的美籍人群中, 該交互作用則支持素質(zhì)-壓力的觀點。我國是個多民族的國家, 然而有關我國不同民族的兒童青少年對環(huán)境的反應模式是否存在差異, 與西方文化背景下的兒童青少年的反應模式是否存在差異的問題尚有待考察。

基因的表達不是獨立于個體的發(fā)育成熟或者環(huán)境的靜態(tài)變量, 而是一個動態(tài)系統(tǒng)(Kaffman &Meaney, 2007; van IJzendoorn, Caspers, Bakermans-Kranenburg, Beach, & Philibert, 2010)。然而已有基因-環(huán)境交互作用的相關研究大多是橫段研究設計, 即使有縱向研究, 也主要是考察短期內(nèi)的變化, 幾乎沒有涉及不同年齡階段的基因-環(huán)境交互作用的發(fā)展變化問題。因而基因-環(huán)境的交互作用模式是否隨年齡的增長而發(fā)展變化, 在某些階段表現(xiàn)為素質(zhì)-壓力模式, 而在另一些階段則表現(xiàn)為不同易感性模型的問題也是未來研究可探討的內(nèi)容之一。此外, 基因-環(huán)境交互作用只說明基因與環(huán)境效應的相互依賴性, 雖然諸多研究指出(Belsky, 1997; Bakermans-Kranenburg & van IJzendoorn, 2011), 個體基因對環(huán)境效應的差異性與神經(jīng)生物系統(tǒng)的敏感性有關, 但具體的作用機制仍有待揭示。

3.2 測量與檢驗方法

某基因與環(huán)境的交互作用究竟是符合素質(zhì)-壓力模型、優(yōu)勢敏感性模型還是不同易感性模型,不僅與觀測到的具體交互作用模式有關, 而且還與所測量的環(huán)境變量與結果變量的性質(zhì)和取值范圍, 以及所采用的統(tǒng)計檢驗方法有關。

概覽現(xiàn)有關于三種模型的研究可以發(fā)現(xiàn), 支持素質(zhì)-壓力模型的研究多局限于對消極環(huán)境,以及問題心理與行為的考察; 旨在驗證不同易感性模型的研究雖然開始采用積極環(huán)境指標, 譬如積極父母教養(yǎng)訓練(Bakermans-Kranenburg, van IJzendoorn, Mesman, Alink, & Juffer, 2008)、敏感性教養(yǎng)(Bradley & Corwyn, 2008)考察基因-環(huán)境交互作用形式, 但是因變量指標的選擇仍然主要基于問題視角, 針對積極心理與行為特征, 譬如自我調(diào)控、親社會行為等的研究尚不多見; 而那些支持優(yōu)勢敏感性模型的研究多屬于干預研究,即創(chuàng)設了積極環(huán)境條件。未來研究需要同時測量個體成長的積極與消極環(huán)境, 積極與消極發(fā)展結果, 增大樣本的代表性(譬如所測量的青少年樣本不僅來自于普通學校, 而且來還來自職業(yè)學校等),以最大程度地提高所揭示的G×E模型的可靠性。

目前而言, 綜合運用傳統(tǒng)回歸法與新興研究方法獲得的檢驗結果更為可靠和令人信服, Roisman等(2012)和Hartman等(Hartman, Widaman, & Belsky,2014)在這方面已經(jīng)作了較好的示范, 但是迄今采用該方法的研究仍寥寥可數(shù), 而且新近的檢驗方法主要是針對素質(zhì)-壓力模型和不同易感性模型的, 有關優(yōu)勢反應敏感性模型的檢驗方法仍局限于傳統(tǒng)的回歸法。因而如何采用一種科學、有效地方法同時對三種模型進行檢驗是未來需要繼續(xù)探究的重要方法問題。

4 小結

素質(zhì)-壓力模型是有關基因-環(huán)境交互作用的傳統(tǒng)理論模型, 不同易感性模型和優(yōu)勢敏感性模型則是新近興起的模型, 運用實證數(shù)據(jù)檢驗三種模型的適切性是當前基因-環(huán)境交互作用研究領域的重要和熱點研究內(nèi)容之一, 且愈來愈多的研究支持后兩種模型的合理性。未來研究不僅需要進一步探索三種理論模型的領域特殊性、種族差異等問題, 而且檢驗這些模型的研究方法也仍需繼續(xù)改善。

王美萍, 張文新. (2010). CHRM2基因rs1824024多態(tài)性與青少年早期抑郁的關系.心理學報, 42(8), 853-861.

王美萍, 張文新. (2014). COMT基因rs6267多態(tài)性與青少年期親子親合與沖突的關系: 性別與父母教養(yǎng)行為的調(diào)節(jié)作用分析.心理學報, 46(7), 931-941.

張衛(wèi), 甄霜菊. (2013). 發(fā)展可塑性的個體差異——基因與環(huán)境的交互作用.華南師范大學學報(社會科學版),(5),106-111.

Bakermans-Kranenburg, M. J., & van IJzendoorn, M. H.(2011). Differential susceptibility to rearing environment depending on dopamine-related genes: New evidence and a meta-analysis.Development and Psychopathology, 23(1),39-52.

Bakermans-Kranenburg, M. J., van IJzendoorn, M. H.,Mesman, J., Alink, L. R. A., & Juffer, F. (2008). Effects of an attachment-based intervention on daily cortisol moderated by dopamine receptor D4: A randomized control trial on 1-to 3-year-olds screened for externalizing behavior.Development and Psychopathology, 20(3), 805-820.

Belsky, J. (1997). Variation in susceptibility to environmental influence: An evolutionary argument.Psychological Inquiry, 8(3), 182-186.

Belsky, J. (2005). Differential susceptibility to rearing influences: An evolutionary hypothesis and some evidence.In B. Ellis & D. Bjorklund (Eds.),Origins of the social mind：Evolutionary psychology and child development(pp.139-163). New York: Guilford.

Belsky, J., Bakermans-Kranenburg, M. J., & van IJzendoorn,M. H. (2007). For better and for worse: Differential susceptibility to environmental influences.Current Directions in Psychological Science, 16(6), 300-304.

Belsky, J., & Beaver, K. M. (2011). Cumulative-genetic plasticity, parenting and adolescent self-regulation.Journal of Child Psychology and Psychiatry, 52(5), 619-626

Belsky, J., Pluess, M., & Widaman, K. F. (2013).Confirmatory and competitive evaluation of alternative gene-environment interaction hypotheses.Journal of C hild Psychology and Psychiatry, 54(10), 1135-1143.

Benjet, C., Thompson, R. J., & Gotlib, I. H. (2010).5-HTTLPR moderates the effect of relational peer victimization on depressive symptoms in adolescent girls.Journal of Child Psychology and Psychiatry, 51(2), 173-179.

Bleuler, M. (1963). Conception of schizophrenia within the last fifty years and today.Proceedings of the Royal Society of Medicine, 56(10), 945-952.

Bradley, R. H. & Corwyn, R. F. (2008). Infant temperament,parenting, and externalizing behavior in first grade: A test of the differential susceptibility hypothesis.Journal of Child Psychology and Psychiatry, 49(2), 124-131.

Caspi, A., McClay, J., Moffitt, T. E., Mill, J., Martin, J.,Craig, I. W., …Poulton, R. (2002). Role of genotype in the cycle of violence in maltreated children.Science, 297(5582),851-854.

Caspi, A., Sugden K., Moffitt T. E., Taylor, A., Craig, I. W.,Harrington, H., … Poulton, R. (2003). Influence of life stress on depression: Moderation by a polymorphism in the 5-HTT gene.Science, 301(5631), 386-389.

Enoch, M. -A., Steer, C. D., Newman, T. K., Gbison, N., &Goldman, D. (2009). Early life stress,MAOA, and geneenvironment interactions predict behavioral disinhibition in children.Genes, Brain and Behavior, 9(1), 65-74,

Foley, D. L., Eaves, L. J., Wormley, B., Silberg, J. L., Maes,H. H., Kuhn, J., & Riley, B. (2004). Childhood adversity,monoamine oxidase a genotype, and risk for conduct disorder.Archives of General Psychiatry, 61(7), 738-744.

Frazzetto, G., Di Lorenzo, G., Carola, V., Proietti, L.,Sokolowska, E., Siracusano, A.,... Troisi, A. (2007). Early trauma and increased risk for physical aggression during adulthood: The moderating role of MAOA genotype.PLoS One, 2(5), e486.

Hartman, S., Widaman, K. F., & Belsky, J. (2014). Genetic moderation of effects of maternal sensitivity on girl’s age of menarche: Replication of Manuck et al. study.Development and Psychopathology, doi: 10.1017/S0954579414000856.

Kaffman, A., & Meaney, M. J. (2007). Neurodevelopmental sequelae of postnatal maternal care in rodents: Clinical,and research implications of molecular insights.Journal of Child Psychology and Psychiatry, 48(3-4), 224-244.

Kegel, C. A. T., Bus, A. G., & van IJzendoorn, M. H. (2011).Differential susceptibility in early literacy instruction through computer games: The role of the dopamine D4 receptor gene (DRD4).Mind, Brain, and Education, 5(2),71-78.Manuck, S. B. (2011).Delay discounting covaries with childhood soc ioeconomic status as a funct ion of genet ic variation in the dopamine D4 r eceptor (DRD4). Paper presented at the Society for Research in Child Development,Montreal, Quebec, Canada.Monroe, S. M., & Simons, A. D. (1991). Diathesis-stress theories in the context of life stress research: Implications for the depressive disorders.Psychological B ulletin,110(3), 406-425.

Neuman, R. J., Lobos, E., Reich, W., Henderson, C. A., Sun,L. W., & Todd, R. D. (2007). Prenatal smoking exposure and dopaminergic genotypes interact to cause a severe ADHD subtype.Biological Psychiatry, 61(12), 1320-1328.

Pluess, M., & Belsky, J. (2010). Children’s differential susceptibility to effects of parenting.Family Science, 1(1),14-25.

Pluess, M., & Belsky, J. (2013). Vantage sensitivity:Individual differences in response to positive experiences.Psychological Bulletin, 139(4), 901-916.

Roisman, G. I., Newman, D. A., Fraley, R. C., Haltigan, J. D.,Groh, A. M., & Haydon, K. C. (2012). Distinguishing differential susceptibility from diathesis-stress:Recommendations for evaluating interaction effects.Development and Psychopathology, 24(2), 389-409.

Rosenthal, D. (1963). Theoretical overview: A suggested conceptual framework. In D. Rosenthal (Ed.),The Genain quadruplets：A case study and theoretical ana lysis of heredity and environment in schizophrenia(pp. 505-516).New York: Basic Books Inc.

Sameroff, A. J. (1983). Developmental systems: Contexts and evolution. In P. Mussen (Ed.),Handbook o f child psychology(Vol. 1, pp. 237-294). New York, NY: Wiley.

Sweitzer, M. M., Halder, I., Flory, J. D., Craig, A. E.,Gianaros, P. J., Ferrell, R. E., & Manuck, S. B. (2013).Polymorphic variation in the dopamine D4 receptor predicts delay discounting as a function of childhood socioeconomic status: Evidence for differential susceptibility.Social Cognitive and Affective Neuroscience, 8(5), 499-508.

Taylor, S. E., Way, B. M., Welch, W. T., Hilmert, C. J.,Lehman, B. J., & Eisenberger, N. I. (2006). Early family environment, current adversity, the serotonin transporter promoter polymorphism, and depressive symptomatology.Biological Psychiatry, 60(7), 671-676.

van IJzendoorn, M. H., & Bakermans-Kranenburg, M. J.(2006). DRD4 7-repeat polymorphism moderates the association between maternal unresolved loss or trauma and infant disorganization.Attachment & Human Development, 8(4), 291-307.

van IJzendoorn, M. H., Bakermans-Kranenburg, M. J., &Mesman, J. (2008). Dopamine system genes associated with parenting in the context of daily hassles.Genes,Brain and Behavior, 7(4), 403-410.

van IJzendoorn, M. H., Belsky, J., & Bakermans-Kranenburg,M. J., (2012). Serotonin transporter genotype 5HTTLPR as a marker of differential susceptibility? A meta-analysis of child and adolescent gene-by-environment studies.Translational Psychiatry, 2(8), e147.

van IJzendoorn, M. H., Caspers, K., Bakermans-Kranenburg,M. J., Beach, S. R. H., & Philibert, R. (2010). Methylation matters: Interaction between methylation density and serotonin transporter genotype predicts unresolved loss or trauma.Biological Psychiatry, 68(5), 405-407.

Walker, E., & Diforio, D. (1997). Schizophrenia: A neural diathesis stress model.Psychological R eview, 104(4),667-685.

Widaman, K. F., Helm, J. L., Castro-Schilo, L., Pluess, M.,Stallings, M. C., & Belsky, J. (2012). Distinguishing ordinal and disordinal interactions.Psychological Methods,17(4), 615-622.

Zuckerman, M. (1999).Vulnerability to psychopathology：A biosocial m odel. Washington: American Psychological Association.