TW202208857A - 評估glp-1類似物療法之反應的方法 - Google Patents

Abstract

本發明關於一種將新穎生物標記用於評估糖尿病患者對類升糖素胜肽1(GLP-1)類似物療法之反應的方法。本發明的新穎生物標記顯示相對於對GLP-1類似物具反應性的糖尿病患者，對類升糖素胜肽1類似物反應較弱的糖尿病患者之某些CpG基因座被調控的胞嘧啶甲基化狀態以及特定的單核苷酸多型性(SNP)之存在。

Description

評估GLP-1類似物療法之反應的方法

本發明關於一種用於評估療程之反應的方法，更具體而言係用於評估糖尿病患者對類升糖素胜肽1(GLP-1)類似物療法之反應的方法。

類升糖素胜肽1(GLP-1)係由腸黏膜的腸內分泌L細胞分泌，並對進食做出反應而經由激素原轉化酶將升糖素原(proglucagon)在其分泌細胞中的轉錄後加工釋放到門脈循環。GLP-1以葡萄糖依賴性方式增強胰島素分泌並抑制升糖素釋放，促進以GLP-1為基礎的療法用於糖尿病之治療的發展。以GLP-1為基礎的療法經由多種機制影響血糖控制，包含減緩胃排空、餐後升糖素的控制、食物攝取的減少以及葡萄糖依賴性胰島素分泌的強化而無低血糖之風險。

然而，在患有第二型糖尿病的患者當中對於GLP-1類似物的臨床反應性不盡相同，意味著遺傳因素可能對這些患者的藥理反應至關重要。由於臨床規範以及這些新藥物的高價格，釐清第二型糖尿病患者的這個關鍵問題是非常重要。

因此，仍有需要評估對GLP-1類似物療法之反應並提供個人化的治療方法。本發明意在解決這些需求與其他需求。

在一實施例中，本發明揭露一種用於評估受試者對GLP-1類似物療法之反應的方法，其包含：(a)測定至少一個選自由cg06536614(SEQ ID NO:1)、cg26328633(SEQ ID NO:2)以及cg25340688(SEQ ID NO: 3)所組成之群組的CpG二核苷酸之胞嘧啶甲基化狀態，或測定來自受試者之樣品的單核苷酸多型性(SNP) rs2346018(風險型等位基因為等位基因A)；以及(b)當在步驟(a)的SNP或在步驟(a)的至少一個CpG二核苷酸之胞嘧啶甲基化相較於對GLP-1類似物療法有反應之受試者低時，鑑定受試者為對GLP-1類似物療法具有低度反應。

本發明將在閱讀下列所附圖式及下文的詳細說明而變得更顯而易見。

用語「受試者(subject)」可指的是具有糖尿病的脊椎動物。受試者包含如哺乳類的溫血動物，如靈長類，且更佳為人類。用語受試者包含如貓、狗等的馴養動物、家畜(例如，牛、馬、豬、綿羊、山羊等)及實驗動物(例如，小鼠、兔、大鼠、沙鼠、天竺鼠等)。

用語「樣品(sample)」指的是含有能夠被偵測的材料之組合物，並包含例如「生物樣品(biological samples)」指的是任何取自生命源的材料，例如人類或其他可能患有糖尿病的哺乳類之動物。生物樣品可以存在任何形式，包含如組織、細胞、細胞沉澱物(cell pellet)、細胞萃取物、手術樣品、活體組織切片(biopsy)或細針抽吸物(fine needle aspirate)的固體材料，或其可以是如尿液、全血、血漿或血清的生物流體之形式，或如唾液(saliva)之由受試者產生的其他任何流體樣品。

當將用語「甲基化程度(methylation level)」應用於基因時指的是CpG環境中存在的一個或多個胞嘧啶殘基是否具有甲基化基團。甲基化程度也可以指的是樣品中的細胞萃段之胞嘧啶C5或N4位置上有或沒有甲基化基團。CpG二核苷酸的甲基化可以駐留在非編碼轉錄控制序列中(例如啟動子、增強子等)或在包含相關基因的內含子和外顯子的編碼序列中。

用語「低甲基化(hypomethylation)」指的是相對於正常對照DNA樣本中相應之CpG二核苷酸處發現的5-mCyt數量的平均甲基化狀態，該狀態對應於在測試DNA樣本的DNA序列中一個或多個CpG二核苷酸處5-mCyt的存在少。

用語CpG基因座也稱為CpG二核苷酸或CpG位點。

本發明提供一種評估受試者對GLP-1類似物療法之反應的方法，其包含：a)測定選自主要由SEQ ID NO:1、SEQ ID NO:2以及SEQ ID NO: 3所組成之群組的CpG二核苷酸之胞嘧啶甲基化狀態，或測定來自受試者之樣品的單核苷酸多型性(SNP) rs2346018；以及(b)當在步驟(a)的SNP或在步驟(a)的至少一個CpG二核苷酸之胞嘧啶甲基化相較於對GLP-1類似物療法有反應之受試者低時，鑑定受試者為對GLP-1類似物療法具有低度反應。

根據本發明之實施例用於評估對GLP-1類似物療法的單核苷酸多型性(SNP)包含rs2346018(此為國家生物技術信息中心數據庫中的GenBank登錄號)之SEQ ID NO:4的多核苷酸序列。所屬技術領域具有通常知識者可以輕易使用GenBank登錄號辨識SNP的序列及位置。在NCBI中註冊的SNP之rs編號對應的特定序列可能隨時間改變。對所屬技術領域具有通常知識者顯而易見的是即使對應的rs編號改變，這些序列也在本發明的範圍內。SEQ ID NO:4的核苷酸序列為多型性序列(polymorphic sequences)(包含代表SNP的多型性位點)。

低甲基化程度或低甲基化為相對的用語，且可以藉由將糖尿病患者之測試樣品的SEQ ID NO:1、SEQ ID NO:2以及SEQ ID NO: 3之CpG基因座的甲基化程度與參考之控制程度(對GLP-1類似物有反應的糖尿病患者)的比較而決定。在一些實施例中，測試樣品中SEQ ID NO:1至3之CpG基因座的甲基化程度為低於對GLP-1類似物有反應的糖尿病患者的1%、5%、10%、15%、20%、25%、30%、35%、40%、45%、50%、55%、60%、65%、70%、75%、80%、85%、90%、95%或100%。表1顯示染色體位置、SEQ ID編號以及DNA序列。

表1

CpG 基因座*	染色體/座標	SEQID 編號	序列
cg06536614		1	AGAAGAGTGACGATGTGGAGAGGGACGGGCTGCATGTGCTCCCCGCCCCCGAGAGGCCTG[CG]TCATGCGGTCTCGCCCGCTCTGCGCCAGGCGTCCTGCTAACGTGTCCTGGAGGGACTCTC
cg26328633		2	TGTGGAGAGGGACGGGCTGCATGTGCTCCCCGCCCCCGAGAGGCCTGCGTCATGCGGTCT[CG]CCCGCTCTGCGCCAGGCGTCCTGCTAACGTGTCCTGGAGGGACTCTCAGTTCCTCCCGCC
cg25340688		3	GAGAGGGACGGGCTGCATGTGCTCCCCGCCCCCGAGAGGCCTGCGTCATGCGGTCTCGCC[CG]CTCTGCGCCAGGCGTCCTGCTAACGTGTCCTGGAGGGACTCTCAGTTCCTCCCGCCCGCA

*「CpG 基因座」一欄為參照Illumina's ® Golden Gate及Infinium ®用於甲基化分析時提供的編號

任何用於評估DNA甲基化程度的方法皆可用於本發明的方法中。評估甲基化程度的方法之非限制性實例包含甲基化敏感性限制酶分析、藉由使用包含亞硫酸鹽、焦亞硫酸鹽、亞硫酸氫鹽或其組合的試劑處理的化學分析或焦磷酸定序。

當亞硫酸氫鈉(sodium bisulfite)接觸DNA時，未甲基化的胞嘧啶轉變成尿嘧啶，同時甲基化的胞嘧啶不會被修飾。如此附加的實施例包含如Illumina® Human Methylation450 BeadChip以及多重PCR分析之基於陣列的分析之使用。

定量擴增法(Quantitative amplification) 可用於定量限制性內切酶消化後側接擴增引物的基因座中完整DNA的數量。定量擴增的方法揭露於例如美國專利號6,180,349；6,033,854以及5,972,602中。

其他甲基化程度的測定方法包含但不限於核酸擴增、聚合酶鏈鎖反應(PCR)、甲基化特異性PCR (MCP)、甲基化CpG島恢復分析(methylated-CpG island recovery assay；MIRA)、結合亞硫酸氫鹽限制分析(combined bisulfite-restriction analysis ；COBRA)、單股結構多形性(single-strand conformation polymorphism ；SSCP)分析、限制性分析及微陣列分析(restriction analysis and microarray analysis) (Lian ZQ et al,使用基於微陣列的甲基化CpG島回收法篩選乳腺癌中明顯高甲基化的基因並鑑定其表達程度，國際腫瘤學雜誌，41: 629-638, 2012)。

GLP-1類似物的非限制性實例包含利拉鲁肽(liraglutide)、艾塞那肽(exenatide)或度拉魯肽(dulaglutide)。

在一些實施例中，用於評估受試者對GLP-1類似物療法之反應的方法包含測定基本上選自SEQ ID NO:1、SEQ ID NO:2以及SEQ ID NO: 3的CpG二核苷酸的胞嘧啶甲基化狀態，以及測定來自受試者之樣品的rs2346018。

測定受試者中選自SEQ ID NO:1、SEQ ID NO:2以及SEQ ID NO: 3中的至少一個CpG基因座的差異甲基化程度及/或rs2346018的存在能使對GLP-1類似物有反應的受試者迅速使用這種類似物治療。

本發明之實施例將藉由下列實例進行說明，但不應以任何方式解釋為對其範圍加以限制。除非另有指明，在描述下列實例的研究中，將遵循習知程序。一些程序將在以下描述以用於說明目的。

實例

研究參與者：在本研究中登錄138位患有第2型糖尿病且接受GLP-1類似物療法的患者。GLP-1類似物選自利拉鲁肽或度拉魯肽。

納入條件為：(a)年齡＞ 20歲；(b)經診斷患有糖尿病＞ 2年；以及(c) 8%至12%的HbA1c量。排除條件為：(a)近期有藥物或酒精成癮之病史；(b)對類似物產品敏感；(c)嚴重心血管疾病；(d)參與其他臨床研究；(e)持續流感、自體免疫疾病或其他代謝病症；以及(f)懷孕或哺乳期婦女。此研究規程由台灣長庚紀念醫院的人體試驗倫理委員會(Institutional Review Board)核准。已取得每個受試者的知情同意書(IRB No.104-9629B)。

在訓練階段，藉由Infinium® HumanMethylation450及 MethylationEPIC BeadChip platform對10位對治療產生GLP-1反應的參與者進行DNA甲基化分析。分析了在GLP-1類似物治療之前及之後，來自10位受試者的基因組DNA的全甲基化圖譜。平均年齡、平均BMI、糖尿病持續期間以及基線A₁ C量在GLP-反應組及無反應組之間並沒有顯著差異。在數據分析之前，根據性別檢查在X及Y染色體基因中的表現差異。結果顯示於第1圖中。可以清楚分辨X及Y染色體特定基因甲基化程度在男性與女性受試者之間的分布。

首先我們分析治療前及治療後的組之間有關甲基化表現的甲基化差異。第2圖中所示的結果顯示在火山圖及熱圖分析中，反應組及無反應組在GLP-1類似物治療之前及之後的甲基話並沒有顯著變化。

參照第3圖的(A)部分，VTRNA2-1 的遺傳性甲基化狀態在GLP-1反應組及無反應組之間最具顯著差異(P＜1x10^-10 )。VTRNA2-1 的遺傳性甲基化狀態在GLP-1類似物治療前後無明顯變化。參照下列表2，除了VTRNA2-1 其他遺傳性DMR區域也顯示了在無反應組及反應組之間不同的甲基化狀態。在候選的DMRs當中，VTRNA2-1 在反應組及無反應組之間顯示出最明顯的差異(P＜1x10^-10 )。

表2：在人類基因組中的差異甲基化區(DMRs)以及在訓練組中與GLP-1治療之血糖反應的相關性

染色體	起始	結束	長度 (bp)	P	基因
chr5	135415693	135416613	920	＜0.0000000001	VTRNA2-1
chr6	29648161	29649084	923	0.000306038	-
chr1	205819179	205819609	430	0.000823387	PM20D1
chr2	113992694	113993313	619	0.00126787	PAX8
chr6	31275551	31276797	1246	0.00180708	HLA-B
chr1	202310824	202311492	668	0.002572174	UBE2T
chr13	23309689	23310675	986	0.002768912	-
chr20	5485144	5485475	331	0.004823737	LINC00654
chr6	32729274	32729876	602	0.005435812	HLA-DQB2
chr10	134600301	134600998	697	0.005996881	-
chr19	58220080	58220818	738	0.006659963	ZNF154
chr1	67600172	67601159	987	0.006776549	Clorf141
chr6	74063982	74064710	728	0.007250176	DPPA5
chr15	91473059	91473569	510	0.007257465	UNC45A
chr19	37825211	37825679	468	0.007315758	HKR1
chr6	149805995	149806652	657	0.009006252	ZC3H12D

為了確認VTRNA2-1 啟動子的甲基化，使用甲基化特異性PCR以及隨後進行DNA定序反應的亞硫酸鹽生物複製(bisulfite cloning following DNA sequencing reactions)分析在VTRNA2-1 啟動子(從-7至-202)中16個CpGs中的每一個之DNA甲基化狀態。在第3圖的(B)部分中，VTRNA2-1 啟動子的差異性甲基化由MS-PCR分析進行確認。在反應組中只有來自對甲基化狀態具有特異性的引子對的PCR產物是明顯的(甲基化狀態)，而在無反應組中則相反(未甲基化狀態)。如第3圖的(C)部分所示，在無反應組中幾乎16個CpGs為未甲基化，而在反應組中16個CpGs中的大部分為過甲基化。總之，在訓練組中VTRNA2-1 啟動子的甲基化與對GLP-1類似物治療的血糖反相關，且VTRNA2-1 啟動子的甲基化狀態經由甲基化特異性PCR以及亞硫酸鹽生物複製進行確認。

在驗證階段，登陸接受結合GLP-1類似物與其他抗糖尿病藥物的128位患者。表3顯示從基線至3個月之A₁ C的變化。A₁ C的變化＜ 0.5%定義為對GLP-1類似物無反應。

表3：登陸本研究之參與者(驗證組)的統計學變數(Demographic variable)

變數	有反應 n=93	無反應 n=35	P
年齡	54.6±13.9	50.7±13.7	0.189
性別 (女性 %)	47.3	58.3	0.326
BMI	27.9±5.3	28.6±5.2	0.563
DM 期間(年)	12.3±1.4	11.6±3.7	0.684
GLP-1治療前的A1 C (%)	10.4±0.3	10.3±0.5	0.875
GLP-1治療後的A1 C (%)	8.6±0.5	11.1±0.2	＜0.0001***

***P＜0.001

如上表3所示，驗證組中登錄了總共為128位患者進行焦磷酸定序分析。在128位患者中，52%每日一次以利拉鲁肽治療，32%每日兩次以艾塞那肽治療，以及16%每周一次以度拉魯肽治療。其中年齡、性別、BMI、糖尿病持續期間以及基線A₁ C量在反應組及無反應組之間沒有顯著差異。在GLP-1治療之後A₁ C量在反應組及無反應組之間具有顯著差異(P＜0.0001)。

DNA萃取

在GLP-1類似物治療之前及之後，使用基因組DNA萃取套組根據製造商(RBC Bioscience, Taiwan)的建議從138位患者的周邊血液的白血球中萃取基因組DNA。

全基因組甲基化分析

在進行分析之前，使用EZ DNA甲基化套組根據製造商(Zymo Research, Orange, A, USA)的指示將1 μg的基因組DNA經過亞硫酸鹽修飾。修飾之後，經亞硫酸鹽轉化的DNA樣品使用Infinium® Human MethylationEPIC BeadChip根據製造商(Illumina Inc, CA, USA)的建議進行評估。由Illumina GenomeStudio software V2011.1萃取每個探針的原始訊號強度(不含背景校正或標準化)。所得的原始數據經過標準化(控制標準化)並藉由製造商軟體進行背景校正以產生甲基化β值。

為分析數據，進行火山圖並以MATLAB (版本 2015a)進行統計分析(t-test)。使用Cluster 3.0以及具有平均鏈接聚類算法的Tree View 軟體將差異甲基化位點(DMPs)可視化。在本發明中，表現最具差異的選擇被用來選擇探針。

甲基化特異性PCR(MS-PCR)

使用MethPrimer software設計針對VTRNA2-1的MS-PCR引子及用於MS-PCR的溫度條件。所述引子列於下表4中。

亞硫酸鹽PCR、生物複製以及DNA定序

使用特定引子藉由PCR將經亞硫酸鹽處理的基因組DNA進行放大(Titanium Taq DNA Polymerase, Takara Bio Inc. USA)。經放大的PCR產物經過純化並選殖到pCRII-TOPO載體中(Invitrogen, CarIsbad, CA)以選擇殖株(clones)。DNA定序在至少10個獨立的殖株上使用Automatic Sequencer (Applied Biosystems, Foster City, CA, US)進行。引子列於下表4中。

亞硫酸鹽焦磷酸定序反應

在VTRNA2-1之CpG位點特異性DNA甲基化的定量分析使用亞硫酸鹽焦磷酸定序(Qiagen PyroMark Q96)進行。焦磷酸定序分析使用由100% 甲基化、0% 未甲基化以及未修飾的人類基因組DNA (EpiTect Control DNA以及Control DNA Set, Qiagen)組成的標準驗證。使用EpiTect Fast BISULFITE Conversino Kit (Qiagen)根據製造商的建議對1 ug的基因組DNA進行修飾。大約10至20 ng的亞硫酸鹽DNA被用於與焦磷酸測序vtRNA2-1正向和反向引子進行PCR反應，引子係使用PyroMark Assay desingn software (Qiagen)設計並列於下表4中。PCR使用經過濾的(0.22 um)無核酸酶水(Invitrogem)在PCR Workstation (AirClean®Systems, Inc)上進行。在進行MS-PCR和焦磷酸測序之前應將胜肽校正及清除。每個生物素化的PCR產物(20 ul)與1 ul的Streptavidin Sepharose^TM (GE Healthcare)以及40 ul 的結合緩衝液(Qiagen)混合，然後搖動10分鐘以進行固定化。使用vacuum workstation (Qiagen)根據製造商的建議將固定化的PCR產物純化程單股DNA、在黏接緩衝液中以0.3 uM的濃度與定序引子混合、加熱至80℃5分鐘，然後冷卻至室溫以使定序引子黏接(annealing)。用於焦磷酸定序反應的模板、酵素、基質及核苷酸在PyroMark Q96 (Qiagen)上設定。

定量反轉錄聚合酶連鎖反應(qRT-PCR)

RNA係使用Tempus Blood RNA Tubes (Thermo Fisher Scientific Inc., US)以及QIAamp RNA Blood Mini Kit (Qiagen, Hilden, Germany)從患者的周邊血液樣品萃取。藉由5 ug的總RNA及具隨機六具體的反轉錄酶(Promega, Madison, WI, USA)合成cDNA，接著使用含vtRNA2-1及5S RNA引子的SYBR Green Master Mix (Applied Biosystems, USA)進行定量PCR。選擇5S RNA作為內參標準化基因。

表4：用於甲基化特異性PCR(MS-PCR)、亞硫酸鹽定序引子(BSP)以及焦磷酸定序反應的引子

引子名稱	正義(sense) 5’-3’	反義(antisense) 5’-3’	大小(bp)	位置 (GRCh37/hg19 組裝)
VTRNA2-1 MSP(M)	TTTTTAAATTATAGAAGAGTGACGA(SEQ ID NO:5)	ACGTTAACAAAACGCCTAACG(SEQ ID NO:6)	117	chr5:135,416,309-135,416,425
VTRNA2-1 MSP(U)	TTTTTAAATTATAGAAGAGTGATGA (SEQ ID NO:7)	ACACATTAACAAAACACCTAACACA (SEQ ID NO:8)	119	chr5:135,416,309-135,416,427
VTRNA2-1 BSP	TTGAAATTTTTAAATTATAGAAGAGTGA (SEQ ID NO:9)	ATAAATAAATTTTACCCCCTTCCAC (SEQ ID NO:10)	197	chr5:135,416,303-135,416,499
焦磷酸定序- VTRNA2-1	AGAGGGAAGGGTTGTATGTGT (SEQ ID NO:11)	AAAAAAACTAAAAATCCCTCCAA (Biotin) (SEQ ID NO:12)	113	chr5:135,416,339-135,416,451
焦磷酸定序- VTRNA2-1 定序引子	TTTAGTTTTAGAGAGGTTTG (SEQ ID NO:13)			chr5:135,416,361-135,416,381

在第4圖的(A)部分中，設計PCR與定序引子以放大表4中chr5: 135416381(cg06536614) (SEQ ID NO:1)、 chr5: 135416388、chr5:135416394 (cg 26328633)(SEQ ID NO:2)以及chr5:135416398 (cg25340688)(SEQ ID NO:3)的區域，用以偵測來自患者之基因組DNA中VTRNA2-1 啟動子的甲基化狀態。使用我們具有甲基化/未甲基化對照基因組(EpiTech control DNA and Control DNA set, Qiagen)的引子或亞硫酸鹽對照基因組DNA的引子放大VTRNA2-1 啟動子，並修飾癌自患者的基因組DNA，但在未修飾的基因組DNA中無法放大(第4圖的(A))。這些PCR產物進一步使用焦磷酸定序進行分析，且那些衍生自控制組DNA的甲基化(100%)及未甲基化控制(0%)顯示出相同的甲基化狀態。此外，參照第4圖的(B)至(D)部分，從未經修飾的基因組DNA衍生的PCR反應中仍然沒有明顯的訊號，表示所設計之引子具有特異性。

藉由焦磷酸定性在反應與非反應性基團之間的VTRNA2-1之啟動子區域上的平均甲基化百分比示於第4圖的(E)部分中。在反應組中的平均甲基化程度高於無反應性組(40.07±2.78 vs 37.48±4.52%, P=0.025)。我們將平均甲基化在40%的甲基化處進行二分化，定義高度甲基化(≥40%)及低度甲基化(＜40%)。如表5所示，反應性的勝算比(odds ratio)為2.757(95% CI 1.240-6.130, P=0.011)，代表具有可遺傳性低甲基化(＜40%)的患者在接受GLP-1類似物治療後傾向於出現較差的血糖反應。

表5：患有糖尿病之患者的GLP-1類似物反應

	GLP-1類似物	勝算比	95% CI	P 值
VTRNA2-1 (所有基因型的 rs2346018)	有反應	無反應
高度甲基化(≥40%)	65	16
低度甲基化(＜40%)	28	19	2.757	1.240-6.130	0.011*
VTRNA2-1 (具CC或CA的基因型)
高度甲基化(≥40%)	31	7
低度甲基化(＜40%)	14	4	1.265	0.318-5.035	0.738
VTRNA2-1 (具CA或 AA的基因型)
高度甲基化(≥40%)	34	9
低度甲基化(＜40%)	14	15	4.048	1.438-11.389	0.007**

P值使用χ2 test測定，*P ＜0.05, **P ＜0.01

表6：在多個邏輯回歸模型中，根據臨床特徵與受試者特徵對GLP-1類似物治療的勝算比

	模型 1	模型2	模型3	模型4
年齡 (歲)	1.023(0.991-1.055)P =0.158	1.019(0.983-1.057)P =0.304	1.027(0.987-1.067)P =0.188	1.023(0.983-1.064)P =0.264
BMI (kg/m2 )		0.936(0.845-1.037)P =0.204	0.938(0.846-1.039)P =0.218	0.959(0.866-1.062)P =0.419
性別 (男性 vs. 女性)	0.438(0.183-1.047)P =0.063		0.523(0.166-1.650)P =0.269	0.545(0.172-1.724)P =0.301
VTRNA2-1啟動子區域的甲基化百分比	1.021(1.005-1.038)P =0.009**	1.026(1.007-1.045)P =0.008**	1.029(1.009-1.050)P =0.004**
VTRNA2-1啟動子區域中的高度-(≥40%) vs. 低度甲基化 (＜40%)				5.264(1.743-15.899)P =0.003**

*統計學上顯著 (**P ＜0.01)

參照表6，在多元邏輯回歸中針對年齡、性別及BMI調整之後，勝算比增加至5.264 (95% CI 1.743-15.899, P=0.003)。進一步地，vtRNA2-1的甲基化狀態藉由CTCF結合位點所接鄰(flanked)。中心粒CTCF結合位點含有A/C多型性(rs2346018)，且已經報導了A等位基因的存在可以增加VTRNA2-1中的DNA甲基化百分比。當VTRNA2-1啟動子區域之甲基化程度對GLP-1類似物反應的影響使用rs2346018的隱性模型(即，主要等位基因C與具有一個或多個次要等位基因A的個體)分析時，與主要等位基因C的純合子相比，次要等位基因A的雜合子或純合子根據啟動子區域的甲基化程度表現出更可預測的GLP-1類似物反應。因此，如表6所示，在具有A等位基因的患者中的勝算比增加至4.408 (95% CI .438-11.389, P=0.007)。總之，在VTRNA2-1啟動子及rs2346018的A等位基因為低度甲基化(＜40%)時傾向為對GLP-1治療無反應。

在VTRNA2-1啟動子的甲基化程度對GLP-1類似物治療的反應性的預測分數

在針對表6中的年齡、性別及BMI調整之後，VTRNA2-1啟動子中的高度遺傳性甲基化程度與第2型糖尿病患者對GLP-1類似物治療有更大反應機會相關。我們進一步使用操作特性(ROC)曲線(AUC)評估了VTRNA2-1的甲基化狀態作為GLP-1治療的預測指標的表現。預測分數的公式計算如下：[0.022×年齡(歲)-0.042×BMI-0.608×性別(男性=1，女性=0)+1.661×甲基化程度 (高甲基化 =1，低甲基化 =0)+1.095]。如第4圖的(E)部分所示，VTRNA2-1啟動子的甲基化狀態用於評估GLP-1反應性具有良好的表現(具有 95% CI 0.628-0.863之AUC, 0.745)。使用Youden’s法之1.542的最佳切點敏感性及特異性分別為53%及86%。藉由VTRNA2-1啟動子區域中每個CpG位點的甲基化狀態的預測模型的其他ROC曲線也藉由年齡、性別及BMI或或rs2346018的已登記多態性調整，且第5圖所示的結果具有類似的預測能力。總之，控制性差的第2型糖尿病患者對GLP-1類似物治療的反應性可能會使用羅吉斯模型(logistic model)及VTRNA2-1的可遺傳甲基化狀態進行預測。

TGF-β1 ELISA分析

取得血液樣品並在離心之後立刻收及血清。直到分析之前將血清樣品儲存於-20℃。使用TGF-β1 ELISA kit (R&D System Inc., Minneapolis, MN)測定血清TGF-β1。

統計分析

使用配對的Student’s t-test對連續變量的組間差異進行檢驗。適當地使用卡方檢驗(chi-square test)或費雪精確檢驗(Fisher’s exact test)評估比例差異。結果以平均值±標準差或百分比表示。同濟學上顯著設定為P值為0.05。使用R package “ChAMP” 進行MethylationEPIC BeadChip 表準化(使用參數“BMIQ”將樣品標準化)、 calling DMPs以及calling DMRs (使用參數“Bumphunter”)。操作特性(ROC)曲線分析及其他統計分析使用SAS (v9.3, SAS Institute, Cary, NC, USA)進行。

VTRNA2-1啟動子的甲基化程基因表現之的相關性

如第6圖所示，甲基化水平與VTRNA2-1的基因表達適度相關(皮爾森值=0.315, P=0.014)，代表VTRNA2-1啟動子中較高的甲基化百分比抑制了其基因表達。此外，VTRNA2-1啟動子位於人類第5染色體並由乙型轉化生長因子(transforming growth factor beta induced,TGFBI )以及SMAD Family Member 5 (SMAD5 )接鄰，兩者都與TGF-β訊號路徑有關。如第7圖的(A)部分所示，VTRNA2-1啟動子的甲基化百分比與血清TGF-β1水平呈正相關，但相關係數為中等(皮爾森值=0.458, P=0.005)。如第7圖的(B)部分所示，使用血清TGF-β1的含量預測對GLP-1類似物組的血糖反應，表示多種因素可調節人體TGF-β1的分泌量。

在表觀遺傳甲基化經晶片分析中，VTRNA2-1啟動子區域的低度甲基化狀態(＜40%)接近於GLP-1類似物的無反應性。VTRNA2-1為RNA聚合酶III轉錄的非編碼RNA，其大小為101 bp，又稱為nc886 。在某些治療藥物耐藥性研究中已顯示出VTRNA2-1受其啟動子甲基化調節，但他可以與被TGF-β1誘導而能夠與Dicer結合以抑制miRNA的成熟並隨後抑制miRNA路徑。在我們近期研究中，我們也發現以腸促素(incretin)為基礎的方案可能涉及某些抗糖尿病治療的不同結果。VTRNA2-1啟動子區域的甲基化狀態在GLP-1類似物反應者與無反應者之間具有明顯不同的模式。因此，在時下精密醫學的時代對患者對GLP-1類似物的反應進行詳細研究至關重要。

本發明中VTRNA2-1啟動子區域上的低甲基化狀態與GLP-1類似物無反應性密切相關。經報導隨後向上調節的VTRNA2-1 與其靶基因TGF-β1相關，可能對β細胞存活產生影響。已知GLP-1用於刺激胰島β細胞產生葡萄糖依賴性胰島素的產生和分泌。初步證據代表GLP-1也可以影響胰腺祖細胞(progenitor)的內分泌分化。此外，TGF-β訊號也可以藉由抑制內分泌祖細胞增殖及增強其分化成成熟的β細胞控制內分泌分化。一些研究代表同功型TGF-β在早期胰腺發育中於朝向內分泌分化和遠離腺泡分化的GLP-1訊號路徑中具有特定的協同作用。TGF-β/SMAD3訊號路徑在調節葡萄糖與能量平衡(energy homeostasis)方面相當重要。肥胖中TGF-β訊號的流失會導致發炎；然而，近期研究中血清TGF-β1含量不是GLP-1反應的預測因子。這可能是由於TGF-β訊號的調控網路複雜。VTRNA2-1與TGF-β1之間對GLP-1反應之相互作用的評估可能不僅僅是藉由患者的TGF-β1量來評估。

DMRs是與其他樣品比較具有不同DNA甲基化模式的DNA片段與生物體基因組。這些樣品可以為相同個體中的不同細胞或組織、不同時間點的相同細胞或來自不同個體的細胞。在多個個體當中已經報導了具有不同甲基化圖譜的個體間DMRs(inter-DMRs)。這種遺傳表觀遺傳標記可以被用於評估疾病及藥物反應。有趣的是，VTRNA2-1 為介穩的表觀等位基因(metastable epiallele)，其對如孕婦營養之類的臨孕期環境敏感，可以指的是在臨孕期階段母體的飲食可以對小孩包含癌症風險、代謝狀態或甚至是藥物反應的終身健康產生影響。

綜上所述，GLP-1類似物的血糖反應與VTRNA2-1啟動子的可遺傳甲基化及rs2346018的多態性有關。篩選鄰近VTRNA2-1 之DMNs的臨床應用可以在開始治療之前識別患者可能對GLP-1類似物無反應。這可以助於在現代以腸促素為基礎的第2型糖尿病治療方案中有效避免此昂貴的藥物被無效地使用。總之，在VTRNA2-1附近的DMR中，對GLP-1類似物療法的可變反應與可遺傳的甲基化狀態相關。可以藉由精確篩選該表觀遺傳標記來製定個人化治療方案。

無。

第1圖為經由(A)火山圖以及(B)熱圖顯示本發明之訓練群(training group)在Infinium^® Human MethylationEPIC BeadChip分析中的性別差異。每個編號代表登錄的患者、編號-1代表GLP-1類似物療法之前萃取的基因組DNA，而編號-2代表GLP-1類似物療法之後萃取的基因組DNA。

第2圖為比較反應組與無反應組在GLP-1類似物治療之前及之後甲基化差異的火山圖。x軸：在GLP-1類似物治療之前及之後反應組與無反應組之β值的差異，y軸：(A) -log₁₀ (P)值以及(B) -log₁₀ (Q) 值。P值為以t test比較在GLP-1類似物治療之前及之後的β值之間的差異。Q值為錯誤發現率調整的P值。在多重測試調整Q值下檢測到的基因甲基化無差異。

第3圖顯示DNA甲基化程度的不同能夠區分以GLP-1類似物治療之患有第二型糖尿病的反應組與無反應組患者。(A) 非監督式(Unsupervised) 階層分群36個甲基化探針。圖譜使用Cluster 3.0以及採用平均連結分群演算法(average-linkage clustering algorithms)的Tree View programs進行可視化。將所有參與者與最重要的差異甲基化區域(DMRs)的熱圖以及基因列表。基於t test選擇最差異表達的差異甲基化點(DMPs)。基於t test(P＜0.0001)以及四倍表達變化選擇在甲基化中最差異表達的β值的火山圖(紅色代表高β值；綠色代表低β值)。每個編號代表登錄的患者、編號-1及編號-2分別代表在GLP-1類似物治療之前及之後萃取的基因組DNA。(B)經由甲基化特異性PCR反應區分VTRNA2-1 的甲基化狀態(NR：無反應；R：有反應；M：甲基化；U：未甲基化)。(C) 在VTRNA2-1 啟動子區域中各CpG二核苷酸之亞硫酸鹽定序分析。每一列表示在每個CpG二核苷酸上甲基化狀態的百分比(+1代表轉錄起始點；NR：無反應；R：有反應)。實體位置係基於GRCh37/hg19染色體。

第4圖顯示VTRNA2-1 啟動子的GLP-1藥物反應與甲基化百分比之間的關聯性。(A)衍生自患者(泳道1-7)的亞硫酸鹽基因組DNA經由經設計的焦磷酸定序PCR引子放大，甲基化控制組(100%甲基化，泳道8)、未甲基化控制組(0%甲基化，泳道9)以及未修飾基因組DNA(泳道10)也使用相同引子組放大。本發明的負控制組為供應有PCR試劑及引子的ddH₂ O (泳道11)。M： 100 bp marker。隨後PCR產物將進行焦磷酸定序反應(pyrosequencing reaction)。(B)至(D)為使用Q96 platform (Qiagen)的PCR產物定序結果之熱解圖(pyrogram)。在各CpG位點的甲基化百分比以灰框表示。由左至右為各樣品中cg06536614、chr5: 135416388、cg26328633、cg 25340688、cg 26896946的甲基化百分比。箭頭代表在焦磷酸定序反應中應為無訊號的非特異性控制。衍生自100% 甲基化、0%甲基化以及 未修飾基因組DNA的PCR焦磷酸定序反應結果分別示於(B)、(C)以及(D)。(E)在VTRNA2-1 之啟動子四個位點：cg06536614、chr5: 135416388、cg 26328633、cg25340688、cg26896946中平均甲基化百分比(Y軸：四個位點的甲基化%之平均)以及對GLP-1治療之藥物反應(無反應、有反應；X軸)之間的關聯性，*P＜0.05。(F)經由年齡、性別及BMI調整之VTRNA2-1 啟動子區域中，甲基化程度的預測模型之AUC。箭頭代表根據Youden’s法之1.542的最佳切點(optimal cutoff point)。

第5圖顯示經由年齡、性別及BMI調整之VTRNA2-1 啟動子區域中，CpG位點甲基化狀態的預測模型之ROC曲線。(A)進一步針對rs2346018 將vtRNA2-1中4個CpG位點的平均甲基化百分比進行調整(1：CC；2：CA或AA)。(B)在vtRNA2-1啟動子中的cg06536614。(C)在vtRNA2-1啟動子中的chr5:135416388。(D) cg26328633。(E) cg25340688。(F) 在vtRNA2-1中的4個CpG位點的平均甲基化百分比。

第6圖顯示根據皮爾森相關係數(Pearson's correlation coefficient)vtRNA2-1 的表達與VTRNA2-1啟動子的甲基化百分比相關。在本發明中用5S RNA的Ct值對vtRNA2-1 的表達進行標準化。

第7圖顯示作為對GLP-1類似物治療的血糖反應之預測子的TGF-β1之血清含量。(A)VTRNA2-1 啟動子甲基化與TGF-β1之血清含量的關聯性。(B)TGF-β1之血清含量可以作為針對GLP-1類似物治療之血糖反應的預測子。

Claims (7)

1. 一種評估受試者對類升糖素胜肽1(GLP-1)類似物療法之反應的方法，其包含： (a)測定至少一個選自由SEQ ID NO:1、SEQ ID NO:2以及SEQ ID NO: 3所組成之群組的CpG二核苷酸之胞嘧啶甲基化狀態，或測定來自受試者之樣品的單核苷酸多型性(SNP) rs2346018；以及 (b)當在步驟(a)的SNP或在步驟(a)的至少一個CpG二核苷酸之胞嘧啶甲基化相較於對GLP-1類似物療法有反應之受試者低時，鑑定受試者為對GLP-1類似物療法具有低度反應。
2. 如請求項1所述之方法，其中該樣品取自來自由組織切片、石蠟包埋組織、活體組織切片、包埋於在石蠟或其切片中的組織、血液、血漿、血清、淋巴液、淋巴組織、腦脊髓液、骨髓其組合所組成之群組的細胞或細胞成分。
3. 如請求項1所述之法，其中測定甲基化狀態包含使用選自由核酸擴增、聚合酶鏈鎖反應(PCR)、甲基化特異性PCR (MCP)、甲基化CpG島恢復分析(MIRA)、結合亞硫酸氫鹽限制分析(COBRA)、單股結構多形性(SSCP)分析、限制性分析及微陣列分析所組成之群組的一個或多個方法。
4. 如請求項1所述之方法，其中GLP-1類似物包含利拉鲁肽(liraglutide)、艾塞那肽(exenatide)或度拉魯肽(dulaglutide)。
5. 如請求項1所述之方法，其中CpG二核苷酸為SEQ ID NO:1。
6. 如請求項1所述之方法，其中CpG二核苷酸為SEQ ID NO:2。
7. 如請求項1所述之方法，其中CpG二核苷酸為SEQ ID NO:3。

Images