導(dǎo)語

乳腺癌是全球女性最常見的癌癥，發(fā)病率和死亡率較高。全基因組和外顯子組遺傳關(guān)聯(lián)研究已成功鑒定出300多個(gè)乳腺癌易感位點(diǎn)，但大多數(shù)基因座和特定基因變異的機(jī)制仍未得到表征，極大的限制了遺傳易感性基因座向新療法和精準(zhǔn)醫(yī)療工具的轉(zhuǎn)化。孟德爾隨機(jī)化（Mendelian Randomization，MR）是一種基于全基因組測(cè)序數(shù)據(jù)（GWAS數(shù)據(jù)），利用單核苷酸多態(tài)（SNP）作為變量工具（Instrumental Variable，IV），用于揭示因果關(guān)系的新型流行病學(xué)方法。簡(jiǎn)而言之，孟德爾隨機(jī)化是利用自然界中隨機(jī)分配的基因?qū)Ρ硇偷挠绊�，來推斷表型�?duì)疾病的影響。與GWAS相比，MR利用了更狹窄的測(cè)試空間，增加了統(tǒng)計(jì)能力，并且本質(zhì)上支持因果基因鑒定。MR可以通過暴露和結(jié)果的遺傳共定位分析得到進(jìn)一步支持。MR的相關(guān)性通過隨機(jī)對(duì)照試驗(yàn)（RCT）中已證實(shí)在疾病中具有病因?qū)W或因果作用的藥物靶點(diǎn)的回顧性分析結(jié)果可靠。

循環(huán)蛋白具有許多適合使用MR發(fā)現(xiàn)乳腺癌生物學(xué)的特征。首先，血漿蛋白質(zhì)組已被證明可以反映癌癥的正常生理和致病生物學(xué)過程。其次，高通量和高精度測(cè)量手段可以檢測(cè)到循環(huán)蛋白。第三，最近的研究表明，大多數(shù)循環(huán)蛋白與順式作用蛋白數(shù)量性狀位點(diǎn) （pQTL） 有關(guān)，即位于距離蛋白質(zhì)編碼基因 1 Mbp 以內(nèi)。第四，單個(gè)順式-pQTL可以解釋蛋白質(zhì)中相對(duì)較大比例的方差，使其成為使用MR進(jìn)行因果推斷的統(tǒng)計(jì)學(xué)上強(qiáng)大的工具。但到目前為止，還沒有關(guān)于女性群體pQTL的報(bào)告。

Olink采用PEA（Proximity Extension Assay）檢測(cè)技術(shù)，將具有特定核苷酸序列探針的一對(duì)抗體與被檢測(cè)蛋白特異性結(jié)合，探針通過末端5bp配對(duì)堿基互補(bǔ)結(jié)合，在延伸酶的作用下形成雙鏈模板，利用qPCR或NGS進(jìn)行檢測(cè)。根據(jù)特異性的核苷酸序列信號(hào)確定常見疾病和生理過程的典型蛋白biomarker含量，用于預(yù)測(cè)疾病和治療結(jié)果、確定新的藥物靶點(diǎn)、驗(yàn)證藥物安全性和有效性，加速蛋白組向精準(zhǔn)醫(yī)療應(yīng)用轉(zhuǎn)化。在這里，本文使用 Olink PEA Explore 檢測(cè)在采樣時(shí)未診斷為乳腺癌的598名女性的血漿樣本，一共測(cè)量了2929種獨(dú)特的蛋白質(zhì)。對(duì)蛋白質(zhì)水平進(jìn)行了遺傳關(guān)聯(lián)分析以鑒定cis-pQTL，在BCAC乳腺癌風(fēng)險(xiǎn)的病例對(duì)照薈萃分析中使用cis-pQTL作為乳腺癌MR分析的工具變量，以及在FinnGen的第二次乳腺癌病例對(duì)照薈萃分析中復(fù)制了MR結(jié)果。最后，通過可視化和評(píng)估蛋白質(zhì)和乳腺癌遺傳關(guān)聯(lián)的共定位來跟蹤MR分析中鑒定的重要蛋白質(zhì)，并評(píng)估與已建立和新出現(xiàn)的乳腺癌風(fēng)險(xiǎn)因素的潛在因果關(guān)系。 
 

主要技術(shù)

olink蛋白質(zhì)組學(xué)；孟德爾隨機(jī)化

研究結(jié)果

1. 樣本特性
KARMA研究招募了瑞典兩個(gè)地區(qū)70877名女性進(jìn)行乳房x光檢查，該項(xiàng)目的目的是確定乳腺癌的危險(xiǎn)因素。該樣本選擇的目的是評(píng)估血漿蛋白生物標(biāo)志物與2年內(nèi)血液采樣中乳腺癌發(fā)病率的關(guān)系；來自同一地區(qū)的299名隨機(jī)對(duì)照女性，截至2021年，這些女性仍然沒有患乳腺癌（表1）。
 

表1
 

2. 蛋白質(zhì)分析、檢測(cè)和質(zhì)量控制
使用親和蛋白質(zhì)組學(xué)方法分析血漿樣本。Olink PEA本質(zhì)上偏向于被測(cè)量的蛋白質(zhì)子集，需要測(cè)量盡可能多的蛋白質(zhì)來最大限度地提高發(fā)現(xiàn)標(biāo)志物的可能性。因此，使用Olink公司新推的Explore I和II panel，包含2949種蛋白質(zhì)。在這組樣本中，歸一化蛋白表達(dá)水平高于檢測(cè)限時(shí)，50%的樣品可以檢測(cè)到2213個(gè)（75%）。Explore II中的一組平均蛋白不如Explore I中的蛋白豐富。

3.血漿蛋白水平與乳腺癌發(fā)病率之間的關(guān)系
為了評(píng)估蛋白質(zhì)與乳腺癌風(fēng)險(xiǎn)的關(guān)聯(lián)，對(duì)每種超過QC的Olink蛋白進(jìn)行回歸模型擬合，該模型調(diào)整了抽血年齡、體重指數(shù)和樣品儲(chǔ)存時(shí)間，使用5%的錯(cuò)誤發(fā)現(xiàn)率來確定統(tǒng)計(jì)顯著性意義，所有蛋白質(zhì)均未超過統(tǒng)計(jì)顯著性閾值，因此，聯(lián)合分析了事件病例和對(duì)照組的蛋白質(zhì)水平。

4. 血漿蛋白水平與臨床特征的關(guān)系
為了檢驗(yàn)觀察性蛋白質(zhì)水平與KARMA女性臨床特征之間的關(guān)系，將每個(gè)測(cè)量到的蛋白質(zhì)與七個(gè)因素（年齡、飲酒量、出生次數(shù)、BMI、激素替代療法(HRT)、絕經(jīng)前后和當(dāng)前吸煙情況）進(jìn)行了回歸分析，顯示共有684種蛋白質(zhì)與BMI相關(guān)，459種蛋白質(zhì)與年齡相關(guān)（圖2）。已有文章報(bào)道，隨著 BMI 的增加，瘦素和脂肪酸結(jié)合蛋白 4 （FABP4） 的血漿水平升高；絕經(jīng)后婦女的卵泡刺激素 （FSHB）較高，吸煙者的胎盤堿性磷酸酶（PLAP）水平較高。同時(shí)還包括一些較少描述的相關(guān)性，包括在絕經(jīng)后和圍絕經(jīng)期婦女中糖氨酸（PAEP）和CHRDL2水平較低，糖蛋白激素α多肽（CGA）水平較高，以及在使用（HRT）的女性中骨調(diào)節(jié)素（OMD）水平較低。已知性狀與蛋白質(zhì)關(guān)聯(lián)表明，數(shù)據(jù)質(zhì)量令人滿意，并且通過增加可檢測(cè)蛋白質(zhì)的數(shù)量，可以實(shí)現(xiàn)額外的性狀與蛋白質(zhì)的關(guān)聯(lián)。
 

圖2

5. 順式pQTL的鑒定和復(fù)制分析
為了確定使用MR進(jìn)行下游因果關(guān)系測(cè)試的遺傳工具，測(cè)試了編碼2929種蛋白質(zhì)中每種蛋白基因上下游1 Mbp范圍內(nèi)的基因變異和蛋白水平的關(guān)聯(lián)，共觀察到超過737個(gè)蛋白質(zhì)的812個(gè)獨(dú)立變異。在Explore I和II panel上很容易檢測(cè)到蛋白質(zhì)的pQTL，為737種蛋白質(zhì)提供了潛在的MR工具。

將KARMA中鑒定的順式pQTL在先前的報(bào)道中進(jìn)行了復(fù)制分析。使用 Olink PEA 測(cè)量的 90 種蛋白質(zhì)子集。KARMA研究中提供了90種蛋白質(zhì)的測(cè)量結(jié)果。在這90種蛋白質(zhì)中，33種蛋白質(zhì)的順式 pQTL 與 KARMA 相關(guān)。先前發(fā)表的基于Somascan蛋白質(zhì)組學(xué)平臺(tái)的cis-pQTL報(bào)告中查找了變異或處于高連鎖不平衡變異。在KARMA研究中質(zhì)量控制后可用的Olink蛋白與先前發(fā)表的基于Somascan平臺(tái)的研究中測(cè)量的蛋白質(zhì)的重疊為569種蛋白質(zhì)。在KARMA中觀察到的重疊蛋白質(zhì)子集的 603 個(gè)顯著順式pQTL中，我們觀察到374個(gè)蛋白質(zhì)復(fù)制的證據(jù)。

6. 孟德爾隨機(jī)化分析
使用來自KARMA順式pQTL的蛋白質(zhì)暴露進(jìn)行了雙樣本MR分析，利用來自BCAC和FinnGen R8 - UK - biobank薈萃分析的結(jié)果數(shù)據(jù)調(diào)查乳腺癌風(fēng)險(xiǎn)的潛在因果效應(yīng)。沒有確定在KARMA中具有順式pQTL的7個(gè)蛋白質(zhì)的遺傳代理，對(duì)730個(gè)蛋白質(zhì)暴露進(jìn)行了測(cè)試。有7個(gè)蛋白在發(fā)現(xiàn)研究中超過了顯著性統(tǒng)計(jì)閾值（圖3），5個(gè)蛋白在FinnGen和UK Biobank的獨(dú)立乳腺癌病例/對(duì)照研究中復(fù)制，效果大小和方向一致（表4）。復(fù)制的蛋白按編碼基因的名稱顯示為CD160、DNPH1、LAYN、LRRC37A2和TLR1。
 

圖3
 

表4
 

在ER+和ER-乳腺癌風(fēng)險(xiǎn)中，效應(yīng)大小相似，表明這五種蛋白與ER+和ER-乳腺癌風(fēng)險(xiǎn)同等相關(guān)（表5）。
 

表5

為了驗(yàn)證蛋白質(zhì)在乳腺癌的病因?qū)W上的作用可能會(huì)通過乳腺癌風(fēng)險(xiǎn)因素影響，MR分析使用潛在乳腺癌危險(xiǎn)因素的GWAS作為結(jié)果，包括初潮年齡、絕經(jīng)年齡、腰臀比、乳房x線攝影密度、性激素結(jié)合球蛋白和IGF-1水平。在兩個(gè)獨(dú)立的數(shù)據(jù)集中，LRRC37A2顯示出較晚的月經(jīng)初潮年齡和較早的絕經(jīng)年齡，以及較高的IGF-1水平的MR證據(jù)。CD160顯示出較低的月經(jīng)初潮年齡在病因?qū)W上的作用。綜上所述，MR分析顯示CD160、DNPH1、LAYN、LRRC37A2和TLR1的基因升高與乳腺癌風(fēng)險(xiǎn)相關(guān)，并且對(duì)ER+和ER -癌癥具有相似的影響。

7. 共定位分析
所有順式pQTL附近有顯著MR證據(jù)的蛋白質(zhì)的變異體都通過鏡像圖與相應(yīng)的基因組區(qū)域一起進(jìn)行了乳腺癌風(fēng)險(xiǎn)的共定位分析。DNPH1和LRRC37A2周圍的順式區(qū)域顯示了蛋白質(zhì)水平和乳腺癌風(fēng)險(xiǎn)之間的最強(qiáng)一致性。CD160、LAYN和TLR1順式區(qū)域的鉛pQTL不是乳腺癌風(fēng)險(xiǎn)p值最低的變異，但定位在相同的大小受限的基因組區(qū)域。我們認(rèn)為順式pQTL與乳腺癌風(fēng)險(xiǎn)共定位。

為了調(diào)查在本研究中發(fā)現(xiàn)的五種蛋白質(zhì)中是否有任何一種曾被作為藥物靶點(diǎn)進(jìn)行過探索，在NIH Pharos Consortium, IUPHAR/BPS Guide To Pharmacology, DrugBank和http://ClinicalTrials.gov數(shù)據(jù)庫中進(jìn)行系統(tǒng)搜索。除了透明質(zhì)酸靶向的LAYN蛋白外，其他都沒有被注冊(cè)為已知的藥物靶標(biāo)。

結(jié)論
在使用MR評(píng)估的730 種血漿蛋白中，五種蛋白的遺傳水平升高與乳腺癌風(fēng)險(xiǎn)相關(guān)，即 CD160、DNPH1、LAYN、LRRC37A2 和 TLR1。這表明這五種蛋白在乳腺癌中起著病因或因果作用，為進(jìn)一步評(píng)估其作為藥物靶點(diǎn)的潛力提供了基礎(chǔ)。

參考文獻(xiàn)：
Mälarstig A, Grassmann F, Dahl L, Dimitriou M, McLeod D, Gabrielson M, Smith-Byrne K, Thomas CE, Huang TH, Forsberg SKG, Eriksson P, Ulfstedt M, Johansson M, Sokolov AV, Schiöth HB, Hall P, Schwenk JM, Czene K, Hedman ÅK. Evaluation of circulating plasma proteins in breast cancer using Mendelian randomisation. Nat Commun. 2023 Nov 24;14(1):7680. doi: 10.1038/s41467-023-43485-8. PMID: 37996402; PMCID: PMC10667261.