16q12 rs3803662 基因多態性與乳腺癌易感性的關聯性分析

范艷玲　郭子健　楊曉丹　吳小紅　余　波　李莉華

?

16q12 rs3803662 基因多態性與乳腺癌易感性的關聯性分析

范艷玲郭子健楊曉丹吳小紅余波李莉華

【摘要】目的分析16q12 rs3803662位點基因多態性與乳腺癌罹患風險的相關性。方法檢索PubMed、 Ovid、 Medline 和Web of Science數據庫的文獻，應用Stata/SE 12.0 對全部關于16q12 rs3803662基因多態性與乳腺癌易感性的研究結果進行異質性檢驗和效應值合并。結果檢索到30篇相關文獻，25篇納入研究的文獻中，累計病例69 127例，對照95 954例。總研究人群分析結果顯示，與C等位基因相比，16q12 rs3803662位點T等位基因顯著增加高加索人群、亞洲人群和混合人群乳腺癌的罹患風險。分層分析結果顯示，與雌激素受體陰性乳腺癌患者相比，rs3803662 T等位基因顯著增加雌激素受體陽性乳腺癌患者的罹患風險。與攜帶BRCA1基因突變的乳腺癌患者相比，rs3803662 T等位基因顯著增加攜帶BRCA2基因突變乳腺癌患者的易感性。結論排除非洲人群，16q12 rs3803662位點T等位基因顯著增加個體罹患乳腺癌風險。

【關鍵詞】rs3803662；雌激素受體；孕激素受體；乳腺癌；Meta分析

(ThePracticalJournalofCancer,2015,30:0046～0050)

乳腺癌是女性最常見的惡性腫瘤之一，是女性癌癥相關死亡的主要原因[1]。近日出版的《2012中國腫瘤登記年報》顯示，乳腺癌的發病率已躍居我國女性惡性腫瘤發病率的首位，明顯高于其它癌癥的發病率，并以每年2%~3%的增長速率遞增和呈現年輕化趨勢。乳腺癌發病是遺傳與環境因素共同作用的結果。與家族性乳腺癌相比，大多數散發性乳腺癌的發生主要與低外顯率基因的多態性有關。近年，全基因組關聯性研究發現一些新的乳腺癌易感基因位點。Easton等[2]發現位于16q12的rs3803662是一個乳腺癌風險因子。然而隨后在不同人群中進行的獨立研究的結論不一致。因此，我們對目前已發表的數據進行Meta分析，以更精確地評估16q12 rs3803662位點基因多態性與乳腺癌易感性的關聯性。

1材料與方法

1.1文獻檢索

在PubMed、Ovid、Medline和 Web of Science數據庫中，以“16q12 or rs3803662”、“polymorphism or polymorphisms”、“breast cancer”為關鍵詞聯合檢索2013年10月以前發表的相關研究。

1.2納入標準

①關于16q12 rs3803662 C>T位點多態性與乳腺癌罹患風險的研究；②病例對照或隊列研究；③有足夠的原始數據以提供比值比(OR)及95%置信區間(CI)或可以轉化OR值及95%CI的研究；④相同研究重復發表的僅納入樣本量較大的文獻。

1.3排除標準

①無對照人群；②無可利用的基因型頻率分布數據。

1.4資料提取

資料提取由2名研究人員分別獨立完成，不一致部分經討論解決。資料包括：第一作者姓名、文獻發表年限、研究對象種族、樣本量、病例組和對照組基因型、各組等位基因頻率、雌激素受體和孕激素受體狀態、BRCA基因狀態和是否進行Hardy-Weinberg平衡(HWE)檢驗。

1.5統計學分析

通過計算下列基因型模型獲得OR值及其95%CI，用以評估16q12 rs3803662位點多態性與乳腺癌易感性的關聯強度：等位基因對比模型(T vs.C)、雜合子模型(TC vs.CC)、純合子模型(TT vs.CC)以及顯性和隱性模型(TT/TC vs.CC，TT vs.TC/CC)，并按照種族、雌激素受體和孕激素受體狀態以及BRCA基因狀態分層分析。各研究之間的異質性判斷采用基于卡方檢驗的Q檢驗[3]，有異質性者應用隨機效應模型(P<0.10和(或)I2>50%)[4]，否則應用固定效應模型(P≥0.10和(或)I2≤ 50%)[5]。采用卡方檢驗計算對照組中Hardy-Weinberg平衡，P值>0.01視為滿足Hardy-Weinberg平衡。應用Egger’s線性回歸和繪制倒漏斗圖的方法評估潛在發表性偏倚，P<0.05表示存在統計學意義的發表性偏差。用敏感性分析評價結果的穩定性。用卡方檢驗比較低頻等位基因在人群中的頻數分布情況。所有的統計學分析均應用Stata/SE version 12.0軟件(Stata Corporation，College Station，TX，USA)。

2結果

2.1基本情況

共檢索到相關文獻資料30篇，根據上述入選及排除標準，最終納入統計文獻共計25篇。研究對象涉及14個國家及地區，其中12篇文獻涉及高加索人群，6篇文獻涉及亞洲人群，1篇文獻涉及非洲人群，6篇文獻涉及混合人種，累計病例69127例，對照95954例。納入文獻的基本資料見表1，所有文獻的基因頻率分布見表2。

2.2文獻合并分析

總研究人群分析顯示，與C等位基因相比，16q12 rs3803662位點T等位基因顯著增加乳腺癌的罹患風險(等位基因模型T vs.C：OR 1.20，95%CI 1.16-1.24；純合子模型TT vs.CC：OR 1.38，95%CI 1.27-1.50；雜合子模型TC vs.CC：OR 1.16，95%CI 1.10-1.21；顯性模型TT/TC vs.CC：OR 1.20，95%CI 1.14-1.27； 隱性模型TT vs.TC/CC：OR 1.29，95%CI 1.21-1.37)，見圖1。

按人種分層分析發現，在高加索人群中，T等位基因顯著增加乳腺癌的罹患風險(T vs.C：OR 1.24，95% 1.17-1.32； TT vs.CC：OR 1.53，95% 1.30-1.79；TC vs.CC：OR 1.18，95% 1.06-1.31； TT/TC vs.CC：OR 1.24，95% 1.11-1.38；TT vs.TC/CC：OR 1.40，95% 1.24-1.59)，對于亞洲人群，在等位基因模型中T等位基因與增加乳腺癌的罹患風險相關(T vs.C OR 1.19，95%CI 1.12-1.26)。在非洲人群中未發現16q12 rs3803662基因多態性與乳腺癌易感性相關。

分析結果顯示，與雌激素受體陰性乳腺癌(OR 1.08，95%CI 0.97-1.20) 相比，rs3803662 T等位基因顯著增加雌激素受體陽性患者乳腺癌(OR1.23，95%CI 1.13-1.32)的罹患風險。對于孕激素受體狀態進行分層分析，得出相似的結論。與攜帶BRCA1基因突變的乳腺癌患者 (OR 1.09，95%CI 1.04-1.15) 相比，rs3803662 T等位基因顯著增加攜帶BRCA2基因突變患者的乳腺癌(OR 1.23，95%CI 1.05-1.44)易感性。

2.3敏感性分析

從總體分析中逐一剔除單個研究之后進行Mata分析顯示，所有基因模型的合并比值比沒有發生實質性的改變，表明盡管在合并分析過程中存在研究間的異質性，但統計結果是仍然是可靠的。

2.4偏倚性分析

發表性偏倚 Begger’s漏斗圖未顯有明顯的不對稱性。Egger’s檢驗顯示在所有基因模型中未發現有統計學意義的發表性偏倚(P>0.05)，見圖2。

3討論

16q12 rs3803662是一個功能未知的基因位點。盡管如此，近期的全基因組關聯性研究已經發現其基因多態性與乳腺癌的罹患風險相關，但是后來的各獨立研究對于其是否增加乳腺癌的罹患風險的結論尚不一致。

總人群的合并分析結果表明16q12 rs3803662 C>T位點基因多態性與乳腺癌的罹患風險相關。在人種的分層分析中，高加索人群、亞洲人群、混合人群中T等位基因相較于C等位基因顯著增加乳腺癌的罹患風險，而在非洲人群中并未發現兩者有相關性。T 等位基因在不同人種分布不同，在Tapper 等[6]的研究中，歐裔美國人中T等位基因的頻率是 0.341，在Barnholtz-Sloan等[7]的研究中，非洲人群中T等位基因的頻率是0.530，因此推測不同人種間腫瘤易感性的差異性結果可能主要與不同人種的遺傳差異有關。此外，本文納入的25篇文獻中高加索人研究病例數最多，且得出rs3803662 C>T 位點多態性增加乳腺癌發病風險的結論，根據所占權重比，本研究得出rs3803662 C>T 位點多態性增加乳腺癌患病風險結論可能性大。

表1　25組病例的一般情況比較

注：其中NA為未能獲取，PB為基于人群的，HB為基于醫院的，Nested為巢式病例對照，HWE為Hardy-Weinberg平衡。

圖1　rs3803662 C>T 基因多態性與人罹患乳腺癌的關系

第一作者報道年限病例組CCCTTTCT對照組CCCTTTCTStacey200761842924256989456Tapper2008127652244301530Mcinerney2009486382821354546532396581460512Barnholtz-Sloan20105855121331682778589440891618618Tamimi201033330054966408415273501103373Latif201042239584123956321713719571175Gorodnova201074501619882778215236112Hemminki201018969342708952Slattery20117787552042311116397886220228181266Teraoka2011306309889214856406061431886892Butt201235327864984406780512952072702Ottini201314319574481343352323701027463Li20091181413237720512312840374208Long20104061862926564934Liang201012641348666513851274644557181374Seuta201156083413101478Han201136914351481217343975161617136126494339Kim20121408310614772627Barnholtz-Sloan2010166378196710770142333182617697Garcia-Closas2008775971321848226501082813295970520263629513757Antoniou200824222173497701731672244183138263192595Campa2011370635281071109405670572147241150161667024Muligan201131092652585887038222899219742679953049Harlid201217941420330500820802768189835274342602Rinella2013219187326200

圖2　納入文獻的發表偏倚性分析

另外，本篇meta分析結果表明：在亞洲人群中，T等位基因顯著增加乳腺癌的罹患風險，這與先前Li等[8]和Liang等[9]的研究不一致。Liang等[9]報道來自于中國的一個小的病例對照研究，而我們的大樣本量的綜合分析發現16q12 rs3803662 C>T 基因多態性與亞洲人群乳腺癌易感性之間有相關性。這種結果的不同可能與納入樣本量不同有關系，因此以后還需要更多更大樣本含量的獨立研究來證實。

分層分析顯示，與雌、孕激素受體陰性相比，T等位基因顯著增加雌激素受體陽性和孕激素受體陽性乳腺癌的罹患風險，表明不同腫瘤亞型的變異風險不同。內分泌系統產生的性激素對乳腺癌發生有促進作用，大部分乳腺組織表達ER、PR、AR，性激素通過這些受體發揮作用[10]。由于雌、孕激素受體狀態是乳腺癌分型的重要標記，因此測定雌、孕激素受體表達情況，不僅可以對已罹患乳腺癌的患者制訂個體化的內分泌治療方案，提高術后生存率，還可在一定程度上預測乳腺癌發病風險，進而預防、早期篩查。

乳腺癌易感基因BRCA1和BRCA2發生突變直接影響基因組的穩定性，使細胞對環境變化的敏感性增加，同時損傷修復功能降低，由此導致多種突變包括原癌基因突變在細胞內累積，最終發生惡性轉化[11]。我們的結果顯示，與BRCA1基因突變攜帶者相比，T等位基因顯著增加BRCA2基因突變攜帶者乳腺癌的易感性。

異質性的存在可能對結果的解釋產生潛在的影響。我們的研究表明異質性并非來自于某單個研究，基因多態性在種族之間頻率分布的差異以及其他未知的因素可能是真正的異質性來源。由于乳腺癌發生受多因素影響，未來還需結合飲食、生活習慣、環境暴露、乳腺相關疾病史等因素進行綜合研究分析。

參考文獻

[1]Shulman LN，Willett W，Sievers A，et al.Breast cancer in developing countries：opportunities for improved survival〔J〕.J Oncol，2010，2010：595167.

[2]Easton DF，Pooley KA，Dunning AM，et al.Genome-wide association study identifies novel breast cancer susceptibility loci〔J〕.Nature，2007，447(7148)：1087-1093.

[3]Higgins JP，Thompson SG.Quantifying heterogeneity in a -

meta-analysis〔J〕.Stat Med，2002，21(11)：1539-1558.

[4]DerSimonian R ，Laird N . Meta-analysis in clinical trials 〔J〕.Control Clin Trials，1986，7(3)：177-188.

[5]Mantel N，Haenszel W.Statistical aspects of the analysis of data from retrospective studies of disease〔J〕.J Natl Cancer Inst，1959，22(4)：719-748.

[6]Tapper W，Hammond V，Gerty S，et al.The influence of genetic variation in 30 selected genes on the clinical characteristics of early onset breast cancer〔J〕.Breast Cancer Res，2008，10(6)：R108.

[7]Barnholtz-Sloan JS，Shetty PB，Guan X，et al.FGFR2 and -

other loci identified in genome-wide association studies are associated with breast cancer in African-American and younger women〔J〕.Carcinogenesis，2010，31(8)：1417-1423.

[8]Li L，Zhou X，Huang Z，et al.TNRC9/LOC643714 polymorphisms are not associated with breast cancer risk in Chinese women〔J〕.Eur J Cancer Prev，2009，18(4)：285-290.

[9]Liang J，Chen P，Hu Z，et al.Genetic variants in trinucleotide repeat-containing 9 (TNRC9) are associated with risk of estrogen receptor positive breast cancer in a Chinese population〔J〕.Breast Cancer Res Treat，2010，124(1)：237-241.

[10]張晟，劉鵬，只向成，等.ER、PR及C-erbB-2表達與乳腺癌臨床病理特征及其預后的關系〔J〕.實用癌癥雜志，2009，24(5)：447-449.

[11]張銀華，趙峰，楊麗麗，等.BRCA1和BRCA2在基底細胞樣癌中的表達及意義〔J〕.實用醫學雜志，2011，27(3)：375-377.

(編輯：吳小紅)

Correlation Analysis of 16q12 rs3803662 Polymorphism and Breast Cancer Susceptibility

FANYanling，GUOZijian，YANGXiaodan，etal.DepartmentofOncology，TheFourthAffiliatedHospitalofSoochowUniversity，Wuxi，214062

【Abstract】ObjectiveTo study the relationship between 16q12 rs3803662 polymorphisms and breast cancer susceptibility.MethodsPubMed，Ovid，Medline and Web of Science database were searched.Crude odds ratio (OR) with 95% confidence interval (CI) was used to assess the strength of association.ResultsA total of 25 studies including 69，127 cases and 95，954 controls were identified.Breast cancer risk was associated with 16q12 rs3803662 polymorphism when all studies were pooled into the meta-analysis (OR 1.20，95%CI 1.16-1.24).For rs3803662，stronger association was observed with ER+ tumors (OR 1.23，95%CI 1.13-1.32) compared to ER- tumors (OR 1.08，95%CI 0.97-1.20) and a same condition occurred for the polymorphism with PR+ tumors (OR 1.26，95%CI 1.02-1.55) versus PR-tumors (OR 1.15，95%CI 0.90-1.46).When stratified by BRCA mutation status，stronger association was observed with BRCA2 carriers (OR 1.23，95%CI 1.05-1.44) than BRCA1 carriers (OR 1.09，95%CI 1.04-1.15).ConclusionThe T-allele of 16q12 rs3803662 is risk factor associated with increased breast cancer susceptibility.

【Key words】rs3803662；Estrogen receptor；Progesterone receptor；Breast cancer；Meta-analysis

(收稿日期2014-03-03修回日期 2014-05-13)

中圖分類號：R737.9

文獻標識碼：A

文章編號：1001-5930(2015)01-0046-05

DOI：10.3969/j.issn.1001-5930.2015.01.013

通訊作者：李莉華

基金項目：無錫市科技發展計劃(CAE00801-08)
作者單位：214062 蘇州大學附屬第四醫院，無錫市第四人民醫院