MTHFR基因1298A>C位點多態性與亞洲人群乳腺癌易感性的Meta分析

葉彥軍 郝福軍 張 玲

(寶雞市人民醫院普通外科，陜西 寶雞 721000)

乳腺癌是城市女性最常見的癌癥，是農村女性第四大常見癌癥，城市地區乳腺癌的發病率是農村地區的2倍〔1，2〕。亞甲基四氫葉酸還原酶(MTHFR)基因是葉酸代謝中的關鍵限速酶之一，與DNA甲基化及DNA合成相關，1298位點是MTHFR基因上研究較為集中的位點，該位點突變與MTHFR活性密切相關。本研究擬采用Meta分析評價MTHFR基因1298A>C位點多態性與乳腺癌發病風險的相關性。

1 資料和方法

1.1文獻檢索 計算機檢索PubMed、EMBASE、web of science、中國生物醫學文獻數據庫、CNKI數據庫、重慶VIP數據庫和萬方數據庫，各數據庫檢索日期均從開始建庫到2016年8月，全面檢索MTHFR基因1298A>C位點多態性與乳腺癌易感性的病例對照研究。英文檢索詞為:MTHFR，methylenetetrahydrofolate reductase和breast cancer，中文檢索詞為:MTHFR、亞甲基四氫葉酸還原酶、乳腺癌。納入標準：①病理學或細胞學檢查確診為乳腺癌的患者；②MTHFR基因1298A>C位點多態性與乳腺癌易感性的病例對照研究；③報道了MTHFR基因1298A>C位點基因多態性的分布頻數；④報道了MTHFR基因1298A>C位點多態性與乳腺癌易感性的OR值及95%CI。

1.2資料提取 由兩位研究者獨立進行文獻篩選，排除明顯不符合納入排除標準的研究，對可能符合納入排除標準的研究通過閱讀全文來確定，資料提取完畢后兩位研究者再交叉核對，對于難以確定是否納入的試驗，通過小組討論解決分歧。

1.3統計分析 采用Stata12.0軟件進行統計分析，OR及95%CI為療效分析統計量。當納入文獻間無明顯的統計學異質性時(I2≤50%,P>0.10)，采用固定效應模型，反之，若納入研究間有顯著的統計學異質性時(I2>50%,P<0.10)，分析異質性可能的來源，采用亞組分析和敏感性分析探討其可能的來源，若無明顯的統計學異質性時采用隨機效應模型進行分析。分別計算以下幾種遺傳模型(AA vs AC、AA vs CC、AC vs CC、AA vs AC+CC、AA+AC vs CC)的基因型分布。采用Stata 12.0軟件進行Begg檢驗以評估納入研究的發表偏倚。

2 結 果

2.1文獻檢索結果 初檢各數據庫共獲得文獻343篇，通過閱讀題目及摘要排除明顯不符合納入排除標準的文獻264篇，剩余79篇文獻通過閱讀全文進一步確定，經過逐層篩選后，最終有17篇文獻〔3～19〕符合納入排除標準。納入研究中2篇研究來自泰國〔12,17〕，1篇來自新加坡和中國的合作研究〔7〕，1篇來自克什米爾地區〔8〕，11篇來自中國大陸〔3,4,6,10,11,13～16,18,19〕，2篇來自中國臺灣〔5,9〕。見表1。

表1 納入研究一般特征(n)

HWE：哈迪-溫伯格平衡定律，HB：對照來源醫院，PB：對照來源社區

2.2Meta分析結果 共有17篇〔3～19〕研究報道了MTHFR基因1298A>C位點多態性與乳腺癌易感性的OR和95%CI，其中，病例組共納入5 680例乳腺癌患者，對照組共納入6 403例志愿者。Meta分析結果結果如下：(1)雜合子模型AA vs AC:CMTHFR基因1298A>C位點多態性與乳腺癌易感性的差異無統計學意義(OR=1.07，95%CI0.99～1.15)，異質性檢驗結果表明納入研究間差異無統計學意義(I2=0.0%,P=0.653)，故采用固定效應模型，見圖1。

(2)雜合子模型AC vs CC:CMTHFR基因1298A>C位點多態性與乳腺癌易感性的差異無統計學意義(OR=0.89，95%CI0.76～1.05)，異質性檢驗結果表明納入研究間差異無統計學意義(I2=0.0%,P=0.891)，故采用固定效應模型，見圖2。

(3)純合子模型AA vs CC:CMTHFR基因1298A>C位點多態性與乳腺癌易感性的差異無統計學意義(OR=0.92，95%CI0.79～1.08)，異質性檢驗結果表明納入研究間差異無統計學意義(I2=0.0%,P=0.906)，故采用固定效應模型，見圖3。

(4)顯性基因模型AA vs AC+CC:CMTHFR基因1298A>C位點多態性與乳腺癌易感性的差異無統計學意義(OR=1.05，95%CI0.97～1.13)，異質性檢驗結果表明納入研究間差異無統計學意義(I2=0.0%,P=0.630)，故采用固定效應模型，見圖4。

(5)隱性基因模型AA+AC vs CC:CMTHFR基因1298A>C位點多態性與乳腺癌易感性的差異無統計學意義(OR=0.92，95%CI0.79～1.07)，異質性檢驗結果表明納入研究間差異無統計學意義(I2=0.0%,P=0.905)，故采用固定效應模型，見圖5。

為了定量評價納入研究發表偏倚的可能性，本研究采用Begg檢驗定量評價發表偏倚，Begg檢驗結果提示差異無統計學意義(P=0.773)，提示納入研究間發表偏倚的可能性小，見圖6。

圖1 MTHFR基因1298A>C位點多態性與乳腺癌相關性的Meta分析(AA vs AC)

圖2 MTHFR基因1298A>C位點多態性與乳腺癌相關性的Meta分析(AC vs CC)

圖3 MTHFR基因1298A>C位點多態性與乳腺癌相關性的Meta分析(AA vs CC)

圖4 MTHFR基因1298A>C位點多態性與乳腺癌相關性的Meta分析(AA vs AC+CC)

圖5 MTHFR基因1298A>C位點多態性與乳腺癌相關性的Meta分析(AA+AC vs CC)

圖6 Begg圖

3 討 論

探討MTHFR基因多態性與乳腺癌的相關性，有助于篩查乳腺癌的易感人群。MTHFR基因位于1號染色體短臂的末端(1p36.6)，參與DNA合成、修復及甲基化過程，是葉酸代謝途徑的關鍵酶〔5〕。在葉酸的代謝過程中，MTHFR可催化5，10-亞甲基四氫葉酸向5-甲基四氫葉酸轉化，MTHFR基因位點突變可導致MTHFR活性降低，導致同型半胱氨酸水平升高，進而影響DNA結構和DNA甲基化水平〔7〕。MTHFR的催化作用決定了DNA甲基化與DNA合成之間的平衡，其基因多態性可能會影響葉酸的代謝，可能會增加了乳腺癌的發生風險〔6〕。MTHFR基因的位點的多態性會使該酶的生物活性及熱穩定性下降，可能會參與乳腺癌的發生及發展，多項研究也顯示MTHFR基因多態性與乳腺癌發生、發展及預后密切相關〔12,16,19〕。MTHFR基因多態性的分布在不同的種族及地區存在明顯差異，且不同地區乳腺癌的發病率存在著明顯的差異，本研究Meta分析結果提示MTHFR基因1298A>C位點多態性與乳腺癌的易感性無相關性。

理論上講，MTHFR基因多態性通過影響葉酸代謝異常，從而誘導DNA甲基化正常模式的畸變，可能會導致DNA修復系統的效率減低，通過影響DNA合成和甲基化導致乳腺癌的發生。也有研究〔12〕提出一種假說，即MTHFR基因多態性與乳腺發病風險受葉酸攝入量的影響，根據這個假說,當葉酸攝入較為充足時,攜帶MTHFR基因變異型者乳腺癌的患病風險降低,因為葉酸足夠可以保證正常的甲基化，從而減少DNA損傷的風險。相反,當葉酸攝入不足時,攜帶MTHFR變異基因型的人DNA甲基化和DNA合成/修復異常成為誘導乳腺癌發生的主要機制。一項Meta分析〔20〕也評價了攝入葉酸和乳腺癌發病率的相關性，該研究共納入14項前瞻性隊列研究共計677 858例個體，Meta分析結果提示葉酸攝入量高的個體和攝入量低者乳腺癌的發病率差異無統計學意義(OR=0.97，95%CI0.90～1.05;P=0.451)，劑量相關分析結果也提示每天增加葉酸攝入量并不能降低乳腺癌的發病風險(OR=0.99，95%CI0.98～1.01;P=0.361)，但是葉酸每天的攝入量大于400 μg時會增加乳腺癌的發病率，葉酸的攝入量與乳腺癌的發病率呈“J”形相關。這也提示乳腺癌的發生受環境因素、激素水平及基因多因素交互作用所致。再者，有研究〔6〕提示不同人群MTHFR基因 1298A>C位點的多態性頻率分布存在較大的差異，這也可能是導致不同地區乳腺癌易感性存在差異的原因之一，Shrubsole 等〔3〕在 研究中，發現對照與病例的 MTHFR基因在1298 位點中，1298位點CC的分布頻率分別為3.3%和3.8%；而Narayanan 等〔21〕在歐洲人群的研究中發現，MTHFR基因1298 位點CC缺高達14.6 %。

本研究尚存在以下局限性：納入研究均為中文和英文，沒有納入亞洲地區其他語種的文獻，雖然Begg檢驗提示發表偏倚的可能性小，仍存在文獻漏檢的可能；有些研究對照人群來源于醫院，有些研究對照人群來源于社區，在年齡和性別方面不完全匹配，可能會影響研究的內在真實性；乳腺癌是環境因素、易感基因及激素水平等多種因素共同作用的結果，由于原始數據的缺乏，無法分析MTHFR基因1 298位點多態性與其他因素的交互作用對亞洲人群乳腺癌易感性的影響，因此，上述結論仍需高質量的研究進一步去證實。