自然。作者手稿;PMC 2012年2月18日提供。
以最终编辑形式发布为:
预防性维修识别码:项目经理3210554
NIHMSID公司:NIHMS313098标准
非霍奇金淋巴瘤组蛋白修饰基因的频繁突变
,*,1 ,*,1 ,1 ,1 ,1 ,1 ,2 ,1 ,1 ,1 ,1 ,1 ,2 ,1 ,1 ,1 ,1 ,2 ,1 ,1 ,1 ,9 ,1 ,1 ,1 ,1 ,4 ,1,7 ,6中,8 ,2 ,2 ,2 ,2 ,1 ,1 ,1 ,1 ,1 ,1 ,1 ,1 ,1 ,1 ,2 ,1 ,1 ,2 ,1 ,2,5和§,1,三
瑞安·D·莫林
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
玛丽亚·门德斯·拉戈
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
安德鲁·蒙格尔
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
罗德里戈·戈雅
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
卡伦·L·蒙格尔
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
理查德·科贝特
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
泰萨·M·塞弗森
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
邱瑞德曼(Readman Chiu)
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
马修·菲尔德
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
肖恩·杰克曼
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
马丁
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
大卫·W·斯科特
2不列颠哥伦比亚省癌症机构淋巴癌中心
黛安娜·L·特里恩
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
杰西卡·塔穆拉-威尔斯
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
萨利
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
马洛公司
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
马拉奇·格里菲斯
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
苏珊娜·陈
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
奥列克桑德·雅科文科
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
埃里克·Y·赵
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
杜安·斯迈卢斯
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
米歇尔·莫克萨
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
苏甘蒂·奇塔兰詹
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
安吉拉·布鲁克斯·威尔逊
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
7西蒙·弗雷泽大学生物医学生理学和运动学系
约翰·斯皮内利
6不列颠哥伦比亚省癌症机构癌症控制研究
8不列颠哥伦比亚大学人口与公共卫生学院
苏珊娜·本·内里亚
2不列颠哥伦比亚省癌症局淋巴癌中心
布鲁斯·沃尔科克
2不列颠哥伦比亚省癌症局淋巴癌中心
海伦·麦克唐纳
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
安吉拉·谭
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
赵永军
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
阿伦·德莱尼
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
托马斯·曾
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
谢霆锋(Kane Tse)
1加拿大不列颠哥伦比亚省癌症机构迈克尔·史密斯基因组科学中心
亚伦·巴特菲尔德
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
Inanc Birol公司
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
罗布·霍尔特
1加拿大不列颠哥伦比亚省癌症机构迈克尔·史密斯基因组科学中心
杰奎琳·谢恩
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
道格拉斯·霍斯曼
2不列颠哥伦比亚省癌症局淋巴癌中心
毛利小五郎
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
史蒂文·琼斯
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
马丁·赫斯特
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
兰迪·D·加斯科因
2不列颠哥伦比亚省癌症局淋巴癌中心
5不列颠哥伦比亚大学病理学系
马可·马拉
1加拿大不列颠哥伦比亚省癌症机构迈克尔·史密斯基因组科学中心
三不列颠哥伦比亚大学医学遗传学系
1加拿大不列颠哥伦比亚省癌症局迈克尔·史密斯基因组科学中心
2不列颠哥伦比亚省癌症机构淋巴癌中心
三不列颠哥伦比亚大学医学遗传学系
4亚利桑那大学病理学系
5不列颠哥伦比亚大学病理学系
6不列颠哥伦比亚省癌症机构癌症控制研究
7西蒙·弗雷泽大学生物医学生理学和运动学系
8不列颠哥伦比亚大学人口与公共卫生学院
9计算机科学系高通量生物学中心
*这些作者贡献均等
- 补充资料
1
GUID:212ED8CC-BA0A-4F2E-B11D-BC99BC1EF6C1
8
GUID:C5625F66-C62F-4684-92BE-42043D5CB9CA
9
GUID:9B7867B2-9E44-492E-B8CD-02F3B0243B5F
10
GUID:6AD146B2-333F-4F60-AC5A-578D591148E1
11
GUID:8AECE0A0-9B30-4F2B-9B34-0C1097EA84CC
2
GUID:A7BEACB6-7588-4417-B95F-B8C2505A7B16
三。
GUID:FB41DE66-A5CF-4549-A7F8-C61AC0430B1D
4
GUID:D8E72045-BBC5-4795-AE2B-22AF48449423
5
GUID:D9365C84-C42F-4D85-A21C-10ED4BAE1CEA
6
GUID:4EC07E1D-CA70-4048-B899-17DB34014AC5
7
GUID:0361AC49-65A4-42FE-93D7-11D68BBA25A7
摘要
滤泡性淋巴瘤(FL)和弥漫性大B细胞淋巴瘤(DLBCL)是两种最常见的非霍奇金淋巴瘤(NHL)。为了确定B细胞NHL突变基因,我们对肿瘤进行了测序,并匹配了13例DLBCL和1例FL患者的正常DNA。我们分析了这些NHL和其他113个NHL的RNA-seq数据,以确定具有候选突变的基因,然后对这些病例中的肿瘤进行重新排序并匹配正常DNA,以确认109个具有多个体细胞突变的基因。组蛋白修饰基因是体细胞突变的常见靶点。例如,32%的DLBCL和89%的FL患者在MLL2级,编码组蛋白甲基转移酶。11.4%的DLBCL和13.4%的FL患者在MEF2B型是一种钙调节基因,在乙酰化组蛋白中与CREBBP和EP300合作。因此,我们的分析表明,淋巴图像病中染色质生物学的破坏是以前未被重视的。
关键词:H3K27、H3K4、HDAC、HAT、癌症基因组学、甲基化、乙酰化、EZH2、驱动因素、癌症测序
介绍
非霍奇金淋巴瘤(NHL)是B、T或自然杀伤淋巴细胞的癌症。在北美每年新诊断的B细胞淋巴瘤中,滤泡性淋巴瘤(FL)和弥漫性大B细胞淋巴瘤(DLBCL)这两种最常见的NHL占60%1FL是一种无痛且典型的不治之症,其特征是临床和遗传异质性。DLBCL具有攻击性,同样具有异质性,至少包含两种不同的亚型,对标准治疗的反应不同。FL和DLBCL的生发中心B细胞(GCB)起源细胞(COO)亚型均来源于生发中心B细胞,而表现出更具侵略性临床病程的活化B细胞(ABC)被认为来源于已经退出或即将退出生发中心的B细胞2目前对导致DLBCL和FL的特定遗传事件的了解仅限于存在一些复发性遗传异常2例如,85-90%的FL和30-40%的GCB DLBCL病例三,4harbort t(14;18)(q32;q21),导致BCL2癌蛋白的表达失控。GCB DLBCL特有的其他遗传异常包括c-继电器miR-17-92 microRNA簇的基因和5与GCB病例相比,24%的ABC DLBCL存在结构改变或失活突变,影响项目风险管理1参与将GCB细胞分化为抗体分泌浆细胞6ABC特异性突变也影响调节NF-κB信号的基因7,8,9,使用TNFAIP3(A20)和88马来西亚第纳尔10突变最丰富的分别为24%和39%。为了加深我们对B细胞NHL遗传结构的理解,我们进行了一项研究,以(1)确定体细胞突变,(2)确定FL和DLBCL突变的流行率、表达和局部复发。使用我们和他人应用于癌症基因组和转录组特征分析的策略和技术11,12,13,我们对117个肿瘤样本和10个细胞系的肿瘤DNA和/或RNA进行了测序(补充表S1和S2)并鉴定出651个基因(补充图S1)有证据表明B细胞NHL发生体细胞突变。经过验证,我们发现109个基因在2个或更多的NHL病例中发生了体细胞突变。我们进一步描述了基因突变的频率和性质MLL2级和MEF2B型这些基因是最常见的突变基因之一,之前在淋巴瘤中没有已知的作用。
复发性突变基因的鉴定
我们对14例非霍奇金淋巴瘤患者的基因组或外显子进行了测序,所有病例的结构DNA测序深度相当(补充表S1和S2). 在筛选单核苷酸变体,然后减去已知多态性并目视检查序列读取比对后,我们确定了717个影响651个基因的非同义(编码单核苷酸变体;cSNV)(补充图S1; 方法)。我们在每个基因组中鉴定了20到135个cSNV。在651个携带cSNV的基因中,只有25个基因参与了癌症基因普查(2010年12月发布)14.
我们对这14例NHL病例和113个扩大样本进行了RNA测序(RNA-seq),其中包括83例DLBCL、12例FL和8例B细胞NHL病例以及其他组织学和10例DLBCL-衍生细胞系(补充表S2). 我们分析了这些数据以确定新的融合转录本(补充表S3)和cSNV(). 我们在基因组/外显子组或RNA-seq“突变热点”(见下文)中鉴定了240个至少含有一个cSNV的基因,总共至少有三个病例含有cSNV(补充表S4). 我们从这240个基因中的每一个中选择cSNV进行重新测序,以确认其体细胞状态。我们没有对先前记录的淋巴瘤突变(例如。CD79B、BCL2). 我们确认了317个基因中543个cSNV的体细胞状态,其中109个基因至少有两个已确认的体细胞突变(补充表S4和S5). 在从基因组预测的成功重测序的cSNV中,171个(94.5%)被确认为体细胞,7个为假呼,3个存在于生殖系中。这109个反复突变的基因显著富集了与淋巴细胞活化相关的基因(P=8.3×10-4; 例如STAT6、BCL10),淋巴细胞分化(P=3.5×10-3; 例如卡片11)和凋亡调节(P=1.9×10-3; 例如BTG1、BTG2). 与转录调控相关的基因也显著富集(P=5.4×10-4; 例如TP53型)和参与甲基化的基因(P=2.2×10-4)和乙酰化(P=1.2×10-2)包括组蛋白甲基转移酶(HMT)和乙酰转移酶(HAT),这些酶以前在淋巴瘤(例如。EZH2型13和CREBBP公司15; 方法)。
NHL体细胞突变靶点的全基因组可视化我们分析中发现的11个DLBCL基因组中的结构重排和拷贝数变异(CNV)以及109个复发突变基因中的cSNV概述。内弧代表在11个基因组中的一个基因组中鉴定的体细胞融合转录物。在11个DLBCL肿瘤/正常配对中检测到的CNV和LOH显示在同心环组上。内11环显示用蓝色绘制的增强纯合子区域(解释为LOH)。外11环显示体细胞CNV。紫色圆圈表示至少有两个已确认体细胞突变的基因的位置,圆圈直径与该基因中检测到的cSNVs病例数成比例。标记了代表具有显著正选择证据的基因的圆圈。反复突变的基因和增益/损耗区域之间的一致性在标签中用彩色编码(绿色=损耗,红色=增益)。例如企业对企业编码β-2-微球蛋白的基因反复突变,有两例被删除。
突变热点可能是由强选择性压力下的位点突变引起的,我们之前已经使用RNA-seq数据确定了此类位点13。我们搜索了RNA-seq数据中有突变热点的基因,并在14个基因组中确定了10个没有突变的基因(PIM1、FOXO1、CCND3、TP53、IRF4、BTG2、CD79B、BCL7A、IKZF3和企业对企业),其中五个(FOXO1、CCND3、BTG2、IKZF3和B2M)不是以前已知的NHL点突变靶点(补充表S6; 方法)。FOXO1、BCL7A和企业对企业显示出影响其起始密码子的热点。The effect of aFOXO1系列使用起始ATG突变为TTG的细胞系,进一步研究了在三例患者中观察到的起始密码子突变。用FOXO1抗体进行的蛋白质印迹显示一条分子量减少的带,表明FOXO1N末端截断(补充图S2)与使用下一帧内ATG启动翻译一致。第二个热点FOXO1系列2例T24位点突变。据报道,在B细胞受体(BCR)刺激后,T24被AKT磷酸化16诱导FOXO1核出口。
我们分析了RNA-seq数据,以确定109个反复突变基因中的任何体细胞突变是否显示出等位基因不平衡,表达偏向于一个等位基因。在380个表达杂合突变等位基因中,我们观察到突变优先表达16.8%(64/380),野生型优先表达27.8%(106/380;补充表S7). 至少在两种情况下,有七个基因显示出突变等位基因的显著优先表达:(BCL2、CARD11、CD79B、EZH2、IRF4、MEF2B和TP53型; 方法)。43例患者中有27例BCL2级cSNV的表达倾向于突变等位基因,这与之前描述的假设一致,即易位(因此,转录解除调控)等位基因BCL2级以体细胞超突变为目标17已知致癌热点位置(如F123I)的突变示例卡片1118在某些情况下,表现出有利于突变等位基因的等位基因失衡。类似地,我们注意到表达有利于两种新的热点突变MEF2B型(Y69和D83)和EZH2型此前未报道淋巴瘤(A682G和A692V)发生突变。
我们试图使用Greenman提出的方法来区分新的癌症相关突变和乘客突变等19我们推断,这将揭示具有强烈选择特征的基因,而这些基因的突变将是很好的候选癌症驱动因素。我们鉴定了26个具有显著阳性选择证据的基因(FDR 0.03,Methods),这些基因要么存在获得非同义点突变的选择压力,要么存在截断/无义突变(Methods;;补充表S8). 包括已知的淋巴瘤致癌基因(BCL2、CD79B9,卡片1118,88马来西亚第纳尔10和EZH2型13),所有这些都表现出指示选择非同义变体的特征。
表1
| 案例 | 总计 | |
---|
基因 | NS公司 | S公司 | T型 | NS公司 | S公司 | T型 | 体细胞cSNV(RNA-seq队列)* | P(生) | q个 | NS SP | T SP(T SP) | 倾斜(M、WT、两者)*** |
---|
MLL2级† | 16 | 8 | 17 | 17 | 8 | 18 | 10 | 6.85×10-8 | 8.50×10-7 | 0.834 | 14.4 | 重量 |
TNFRSF14型G公司† | 7 | 1 | 7 | 8 | 1 | 7 | 11 | 6.85×10-8 | 8.50×10-7 | 7.52 | 118 | 二者都 |
新加坡元1G公司† | 18 | 6 | 6 | 37 | 10 | 6 | 9 | 6.85×10-8 | 8.50×10-7 | 19.5 | 61.7 | - |
BCL10(BCL10)† | 2 | 0 | 4 | 三 | 0 | 4 | 4 | 6.85×10-8 | 8.50×10-7 | 3.62 | 112 | 重量 |
GNA13号机组G公司† | 21 | 1 | 2 | 33 | 1 | 2 | 5 | 6.85×10-8 | 8.50×10-7 | 24.1 | 25.7 | 二者都 |
TP53型G公司† | 20 | 2 | 1 | 23 | 三 | 1 | 22 | 6.85×10-8 | 8.50×10-7 | 15.6 | 14.1 | 二者都 |
EZH2型G公司† | 33 | 0 | 0 | 33 | 0 | 0 | 33 | 6.85×10-8 | 8.50×10-7 | 11.4 | 0 | 二者都 |
基站2† | 12 | 6 | 1 | 14 | 6 | 1 | 2 | 6.85×10-8 | 8.50×10-7 | 23.9 | 35.1 | - |
BCL2级G公司† | 42 | 45 | 0 | 96 | 105 | 0 | 43 | 9.35×10-8 | 8.50×10-7 | 3.78 | 0 | M(M) |
BCL6号机组†** | 11 | 2 | 0 | 12 | 2 | 0 | 2 | 9.35×10-8 | 8.50×10-7 | 0.175 | 0 | M(M) |
CIITA公司†** | 5 | 三 | 0 | 6 | 三 | 0 | 2 | 9.35×10-8 | 8.50×10-7 | 0.086 | 0 | |
船边交货† | 2 | 0 | 4 | 三 | 0 | 4 | 2 | 1.52×10-7 | 1.17×10-6 | 2.54 | 66.5 | 重量 |
英国电信1号† | 11 | 6 | 2 | 11 | 7 | 2 | 10 | 1.52×10-7 | 1.17×10-6 | 17.5 | 52.5 | 二者都 |
MEF2B型G公司† | 20 | 2 | 0 | 20 | 2 | 0 | 10 | 2.05×10-7 | 1.47×10-6 | 14.2 | 0 | M(M) |
IRF8型† | 11 | 5 | 三 | 14 | 5 | 三 | 三 | 4.55×10-7 | 3.03×10-6 | 8.82 | 28.2 | 重量 |
TMEM30A型† | 1 | 0 | 4 | 1 | 0 | 4 | 4 | 6.06×10-7 | 3.79×10-6 | 0.785 | 65 | 重量 |
CD58型† | 2 | 0 | 三 | 2 | 0 | 三 | 2 | 2.42×10-6 | 1.43×10-5 | 2.29 | 69.2 | - |
KLHL6号机组† | 10 | 2 | 2 | 12 | 2 | 2 | 4 | 1.00×10-5 | 5.26×10-5 | 5.42 | 16.4 | - |
88马来西亚第纳尔A类† | 13 | 2 | 0 | 14 | 2 | 0 | 9 | 1.00×10-5 | 5.26×10-5 | 12.4 | 0 | 重量 |
CD70型† | 5 | 0 | 1 | 5 | 0 | 2 | 三 | 1.70×10-5 | 8.48×10-5 | 7.08 | 44 | - |
CD79B型A类† | 7 | 2 | 1 | 9 | 2 | 1 | 5 | 2.00×10-5 | 9.52×10-5 | 10.9 | 18.3 | M(M) |
CCND3号机组† | 7 | 1 | 2 | 7 | 1 | 2 | 6 | 2.80×10-5 | 1.27×10-4 | 6.55 | 36.3 | 重量 |
CREBBP公司† | 20 | 7 | 4 | 24 | 7 | 4 | 9 | 1.00×10-4 | 4.35×10-4 | 2.72 | 6.04 | 二者都 |
HIST1H1C公司† | 9 | 0 | 0 | 10 | 0 | 0 | 6 | 1.80×10-4 | 7.50×10-4 | 11.9 | 0 | 二者都 |
企业对企业† | 7 | 0 | 0 | 7 | 0 | 0 | 4 | 3.90×10-4 | 1.56×10-3 | 16.6 | 0 | 重量 |
ETS1型† | 10 | 1 | 0 | 10 | 1 | 0 | 4 | 4.10×10-4 | 1.58×10-3 | 5.76 | 0 | 重量 |
卡片11† | 14 | 三 | 0 | 14 | 三 | 0 | 三 | 1.90×10-3 | 7.04×10-3 | 3.37 | 0 | 二者都 |
工厂验收试验2†** | 2 | 1 | 0 | 2 | 1 | 0 | 2 | 6.30×10-3 | 2.25×10-2 | 0.128 | 0 | - |
IRF4型†** | 9 | 4 | 0 | 26 | 5 | 0 | 5 | 7.00×10-3 | 2.41×10-2 | 0.569 | 0 | 二者都 |
FOXO1系列† | 8 | 4 | 0 | 10 | 4 | 0 | 4 | 7.60×10三 | 2.53×10-2 | 4.02 | 0 | - |
STAT3(状态3) | 9 | 0 | 0 | 9 | 0 | 0 | 4 | 2.19×10-2 | 6.08×10-2 | - | - | 二者都 |
RAPGEF1型 | 8 | 三 | 0 | 10 | 三 | 0 | 三 | 2.98×10-2 | 7.45×10-2 | - | - | 重量 |
ABCA7公司 | 12 | 三 | 0 | 15 | 三 | 0 | 2 | 7.76×10-2 | 1.67×10-1 | - | - | 重量 |
RNF213号机组 | 10 | 8 | 0 | 10 | 8 | 0 | 2 | 7.87×10-2 | 1.67×10-1 | - | - | - |
MUC16型 | 17 | 12 | 0 | 39 | 25 | 0 | 2 | 8.32×10-2 | 1.73×10-1 | - | - | - |
HDAC7型 | 8 | 4 | 0 | 8 | 4 | 0 | 2 | 8.94×10-2 | 1.82×10-1 | - | - | 重量 |
PRKDC公司 | 7 | 三 | 0 | 7 | 4 | 0 | 2 | 1.06×10-1 | 2.05×10-1 | - | - | - |
SAMD9公司 | 9 | 2 | 0 | 9 | 2 | 0 | 2 | 1.79×10-1 | 3.01×10-1 | - | - | - |
TAF1型 | 10 | 0 | 0 | 10 | 0 | 0 | 2 | 3.03×10-1 | 4.74×10-1 | - | - | - |
PIM1型 | 20 | 19 | 0 | 33 | 34 | 0 | 11 | 3.40×10-1 | 5.23×10-1 | - | - | 重量 |
COL4A2系列 | 8 | 2 | 0 | 8 | 2 | 0 | 2 | 7.64×10-1 | 8.99×10-1 | - | - | - |
EP300型 | 8 | 7 | 1 | 8 | 7 | 1 | 三 | 9.54×10-1 | 1 | - | - | 重量 |
选择失效变更的证据
我们预计肿瘤抑制基因在获得无义突变方面表现出强烈的选择性。在我们的分析中,8个最重要的基因包括7个对无义突变具有强烈选择性压力的基因,包括已知的肿瘤抑制基因TP53型和TNFRSF14型20().CREBBP公司,最近报告为DLBCL中常见的失活15,也显示了获得无义突变和cSNV的一些证据(补充图S3;补充表S9). 我们还观察到无义突变在BCL10(BCL10)一种NF-κB的阳性调节物,在淋巴瘤中描述了致癌截断产物21.其余的极显著基因(BTG1、GNA13、SGK1和MLL2级)在淋巴瘤中没有报道作用。GNA13号机组受22例突变影响,包括多个无义突变。GNA13号机组编码负责调节RhoA活性的异源三聚体G蛋白偶联受体的α亚单位22.一些突变残基对其功能产生负面影响23,24,包括T203A突变,该突变也表现出有利于突变等位基因的等位基因失衡(补充表S7). 在含有无义突变、终止密码子缺失、移码缺失或影响剪接位点变化的细胞系中,GNA13蛋白在Western blot上减少或缺失(方法;补充图S4).
新加坡元1编码PI3K调节激酶,具有调节FOXO转录因子的功能25,通过磷酸化IkB激酶调节NF-κB26和NOTCH信号的负调节27.新加坡元1也存在于DLBCL中通常缺失的6号染色体区域内()5.其机制新加坡元1和GNA13号机组失活可能导致淋巴瘤尚不清楚,但对其失活的明显选择程度很强,且其整体高突变频率(每106例DLBCL中有18例发生突变)表明,其缺失可能导致B细胞NHL。已知某些基因在GCB DLBCL中更常见的突变(例如。TP53型28和EZH2型13). 这里,两个新加坡元1和GNA13号机组仅在GCB病例中发现突变(P=1.93×10-3和2.28×10-4Fisher精确检验;n=15和18)(). 另外两个基因(MEF2B型和TNFRSF14型)在DLBCL中没有先前描述的角色,显示出与GCB病例类似的限制().
非霍奇金淋巴瘤突变和潜在合作相互作用概述该热图显示了体细胞突变和结构重排的共现(红色)和互斥(蓝色)的可能趋势。通过左尾和右尾Fisher精确测试的最小值分配颜色。为了捕捉趋势,使用了0.3的P值阈值,颜色的最深阴影表示符合统计学意义(P<0.05)。左侧显示了ABC(蓝色)、GCB(红色)、不可分类(黑色)DLBCL和FL(黄色)病例的相对突变频率。基因排列为ABC病例(蓝色三角形)突变显著富集(P<0.05,Fisher精确检验)的向上(和左侧)和GCB病例(红色三角形)突变明显富集的向下(和右侧)。顶部显示了每个基因包含cSNV或经证实的体细胞突变的病例总数。蓝色方块簇(右上方)是ABC富集突变相互排斥的结果(例如。MYD88、CD79B)从富含GCB的突变(例如。EZH2、GNA13). 存在涉及两个癌基因的结构重排2006年和BCL2级(表示为BCL6s和BCL2s)通过FISH技术利用断裂部分探针进行测定(方法)。
正在停用MLL2级突变
MLL2级显示了最重要的选择证据和最多的无意义SNV。我们的RNA-seq分析表明,26.0%(33/127)的病例携带至少一种MLL2级cSNV公司。为了解决可变RNA-seq覆盖率MLL2级未能捕获一些突变,我们PCR扩增了整个MLL2级89例患者(35个原发性FLs、17个DLBCL细胞系和37个DLBCLs)的基因座(~36kb)。其中58例属于RNA-seq队列。Illumina扩增子重测序(方法)揭示了78个突变,证实了重叠病例中的RNA-seq突变,并确定了33个额外突变。我们使用Sanger测序法确认了46个变异体的体细胞状态(补充表S10),并显示33个额外突变中有20个是插入或缺失(indels)。剪接位点也检测到三个SNV,还有10个新的cSNV没有被RNA-seq检测到。
体细胞突变分布在MLL2级(). 其中37%(n=29/78)为无义突变,46%(n=36/78)为改变阅读框的indels,8%(n=6/78)是剪接位点的点突变,9%(n=7/78)是非同义氨基酸替换(). 四个体细胞剪接位点突变对MLL2级转录长度和结构。例如,两个杂合剪接位点突变导致使用了一个新的剪接供体位点和一个内含子保留事件。
影响体细胞突变的总结和影响MLL2级和MEF2B型(A) 重新排序MLL2级89个样本的基因座主要表现为无意义(红圈)和移码诱导的indel突变(橙色三角形)。还观察到少量影响剪接位点(黄星)的非同义体细胞突变(绿圈)和点突变或缺失。所有非同义点突变都影响催化SET结构域、FYRC结构域(“富FY-C末端结构域”)或PHD锌指结构域中的残基。这些剪接位点突变对MLL2级还探索了剪接(补充图S7). (B) 发现的cSNV和体细胞突变MEF2B型在所有FL和DLBCL病例中,用相同的符号示出了测序。只有在至少两名患者中含有变异的氨基酸才被标记。cSNV在MEF2B的前两个蛋白编码外显子(外显子2和3)中最常见。与EP300结合的MEF2的晶体结构支持这样的观点,即两个突变位点(L67和Y69)在这些蛋白质之间的相互作用中很重要(补充图S8;补充讨论)50.
表2
样品类型 | 佛罗里达州 | DLBCL(DLBCL) | DLBCL单元线 | 成中心细胞 |
---|
截断 | 18 | 4 | 7 | 0 |
带移帧的索引 | 22 | 8 | 6 | 0 |
拼接现场 | 4 | 2 | 0 | 0 |
SNV公司 | 三 | 2 | 2 | 0 |
|
任何突变(病例数) | 31/35 | 12/37 | 10/17 | 0/8 |
|
百分比 | 89% | 32% | 59% | 0% |
我们测序的NHL病例中大约有一半有两例MLL2级突变(补充表S10). 我们对8例FL患者进行了BAC克隆测序,以表明在所有8例患者中,突变均位于反式,影响两者MLL2级等位基因。这一观察结果与以下概念一致:MLL2级在这些病人的肿瘤细胞中。
除了两个原发性FL病例和两个DLBCL细胞系(Pfeiffer和SU-DHL-9)MLL2级突变似乎是杂合的。对两例具有明显纯合子突变的FL患者的Affymetrix 500k SNP序列数据进行分析,结果表明,这两种肿瘤在12号染色体包含MLL2级(方法)。因此,除了双等位基因突变外,LOH是第二种机制,尽管其不太常见MLL2级功能丢失。
MLL2级是FL中最常见的突变基因,也是DLBCL中最常见突变基因之一(). 我们确认了MLL2级35例FL患者中的31例(89%)、37例DLBCL患者中的12例(32%)、17例DLBCL细胞系中的10例(59%)以及我们测序的8例正常中心体样本中无一例发生突变。我们的分析预测MLL2级正在失活(91%破坏了阅读框或截断了点突变),这向我们表明MLL2级是NHL中重要的肿瘤抑制因子。
复发性点突变MEF2B型
我们的选择性压力分析还显示,与无义突变相比,氨基酸替代的获得压力更大的基因。其中一个基因是MEF2B型之前未与淋巴瘤相关。20例(15.7%)患者MEF2B型cSNV和4例(3.1%)患者MEF2C公司cSNV。RNA-seq检测到的所有cSNV都影响MADS盒或MEF2域。确定MEF2B型我们对261例原发性FL样本的外显子2和3进行了Sanger测序;259例DLBCL原发肿瘤;17个细胞株;各类NHL 35例(IBL、复合FL、PBMCL);和8个非恶性中心粒细胞样本。我们还使用捕获策略(方法)对整个MEF2B型261个FL样本中的编码区,揭示了外显子2和3之外的六个额外变体。因此,我们确定了69例(34例DLBCL;12.67%和35例FL;15.33%)MEF2B型cSNV或indels未能在其他NHL和非恶性样本中观察到新的变异。55(80%)个变体影响MADS盒和外显子2和3编码的MEF2结构域内的残基(补充表S11;). 每个患者通常都有一个MEF2B型我们观察到相对较少(总共8个,10.7%)的截短诱导SNV或indels。到目前为止,非同义SNV是观察到的最常见的变化类型,59.4%的检测到的变异影响K4、Y69、N81或D83。12例MEF2B型突变显示为体细胞突变,包括K4、Y69、N81和D83的代表性突变(补充表S12). 我们没有检测到ABC病例中的突变,这表明MEF2B型在GCB DLBCL和FL的发展中发挥独特的作用().
讨论
在我们对127例B细胞NHL病例的基因组、转录组和外显子组序列的研究中,我们鉴定了109个基因,这些基因在多个个体中有明确的体细胞突变证据。显著的选择似乎作用于其中至少26个,以获得无意义或错义突变。据我们所知,这些基因中的大多数以前没有与任何癌症类型相关。我们观察到影响转录调控基因的体细胞突变增多,更具体地说,是染色质修饰。
MLL2级从我们的分析中发现,它是NHL中的一个主要肿瘤抑制基因座。它是六种人类H3K4特异性甲基转移酶之一MLL公司所有这些都与果蝇属三胸基因29三甲基化H3K4(H3K4me3)是与活性转录基因启动子相关的表观遗传标记。通过设置这个标记,MLL负责发育基因的转录调控,包括同源异型盒(Hox)基因家族30它们共同控制发育中胚胎的片段特异性和细胞命运31,32每个MLL家族成员被认为针对不同的Hox基因亚群33此外,已知MLL2可调节多种基因的转录34。最近,MLL2级小细胞肺癌细胞系发生突变35和肾癌36但无意义突变的频率影响MLL2级这些报告没有证实这些癌症。Parsons及其同事最近报道了MLL2级或MLL3级16%的髓母细胞瘤患者37进一步牵连MLL2级作为一种癌症基因。
我们的数据链接MLL2级B细胞NHL的体细胞突变。报告的突变很可能是失活的,在8例多重突变病例中,我们证实两个等位基因都受到影响,可能导致MLL2功能基本完全丧失。The high prevalence ofMLL2级FL突变(89%)等于t(14;18)(q32;q21)易位的频率,这被认为是FL中最常见的遗传异常三在DLBCL肿瘤样本和细胞系中,MLL2级突变频率分别为32%和59%,也超过了最常见的细胞遗传学异常的流行率,例如涉及3q27的各种易位,这些易位发生在25-30%的DLBCL中,并在ABC病例中富集38重要的是,我们发现MLL2级两种DLBCL亚型均发生突变(). 因此,我们的分析表明MLL2级在FL和两种DLBCL亚型中充当中心肿瘤抑制物。
这个MEF2型该基因家族编码四种相关转录因子,以钙调节的方式招募组蛋白修饰酶,包括组蛋白脱乙酰酶(HDAC)和HAT。虽然在我们的分析中检测到截断变体MEF2型基因家族成员,我们的分析表明,与MLL2、MEF2家族成员倾向于选择性地获得非同义氨基酸取代。在以下情况下MEF2B型,59.4%的cSNVs在蛋白质内的四个位点(K4、Y69、N81和D83)发现,并且这些位点中的所有四个都被证实是体细胞突变的靶点。39%的MEF2B型这些改变影响D83,导致用丙氨酸、甘氨酸或缬氨酸替换带电的天冬氨酸。虽然我们还无法预测这些替代对蛋白质功能的影响,但它们的作用可能会影响MEF2B促进基因表达的能力,从而在促进生发中心B细胞向淋巴瘤的恶性转化中发挥作用(补充讨论).
MEF2B型突变可能与CREBBP公司和EP300型突变和复发性Y641突变EZH2型13CREBBP/EP300 HAT活性的一个靶点是H3K27,它被EZH2甲基化以抑制转录。有证据表明EZH2的作用对抗CREBBP/EP300的作用39MEF2的一个功能是招募HDAC或CREBBP/EP300来靶向基因40有人建议HDAC与CREBBP/EP300竞争MEF2上的相同结合位点41.正常Ca下2+MEF2与IIa型HDAC结合,使组蛋白尾部保持在脱乙酰抑制染色质状态42.胞浆钙增加2+水平诱导HDAC的核输出,使HAT(如CREBBP/EP300)的招募成为可能,促进MEF2靶基因的转录。的突变爬行器,EP300或MEF2B型由于这些基因附近核小体乙酰化减少,可能影响MEF2靶基因的表达(补充图S5;补充讨论). 根据最近的发现,杂合子EZH2型Y641突变增强PCR2的整体H3K27三甲基化活性43,44,可能两者的突变MLL2级和EZH2型可以合作减少一些相同靶基因的表达。我们的数据表明:(1)组蛋白的转录后修饰在生发中心B细胞中至关重要;(2)由于这些突变而引起的组蛋白修饰被解除管制可能导致乙酰化减少和甲基化增强,并在NHL的发展中起核心驱动事件的作用(补充图S5).
方法总结
所有分析的样品中都含有至少50%的肿瘤细胞。使用Illumina GAIIx和HiSeq 2000仪器组合对基因组、外显子和转录组进行测序,以读取36到100个核苷酸的长度。外显子捕获是使用安捷伦SureSelect靶点富集系统协议(版本1.0,2009年9月)执行的。使用BWA完成校准45并使用SNVmix鉴定变异体46对IGV中的变异体进行手动审查,并通过PCR(如适用)进行确认,然后进行Sanger测序或Illumina重新测序。使用ABySS鉴定基因组和转录组中的结构重排47用于亚型分配的基因表达值计算为RPKM值48根据Affymetrix表达阵列数据开发的方法,对亚型进行了调整49使用之前按此标准方法分类的样本进行训练。
致谢
本研究的部分资金来源于国家癌症研究所癌症基因组办公室(合同号HHSN261200800001E)、特里·福克斯基金会(019001号赠款,癌症生物学:淋巴瘤基因组分析的见解)、,Genome Canada/Genome BC Grant Competition III(项目名称:滤泡淋巴瘤基因组的高分辨率分析)授予J.M.C.、R.D.G.和M.A.M.。我们感谢NIH赠款P50CA130805-01“淋巴瘤中的SPORE,组织资源核心(PI Fisher)”的支持和1U01CA114778“改善淋巴瘤诊断和预后的分子特征(PI Chan)”。A.J.M.是白血病和淋巴瘤协会职业发展项目研究员。N.A.J.是特里·福克斯基金会(NCIC 019005奖)和迈克尔·史密斯健康研究基金会(ST-PDF-01793)的研究员。M.A.M.是特里·福克斯(Terry Fox)的年轻调查员和迈克尔·史密斯(Michael Smith)的高级研究学者。R.D.M是Vanier学者(CIHR),持有MSFHR高级研究生学位。M.M.L感谢西班牙教育部博士后奖学金项目“2008-2011年国家I-D+I计划”的支持。D.W.S得到了加拿大卫生研究院Terry Fox基金会癌症研究战略健康研究培训项目(批准号TGT-53912)的支持。JJS获得了加拿大癌症协会和加拿大卫生研究院的认可资金。RG由UBC四年奖学金资助。IMM承认加拿大领导人机会基金创新基金会。本研究的实验室工作由基因组科学中心、不列颠哥伦比亚省癌症研究中心和转化与应用基因组中心承担,该中心是省级卫生服务局实验室的一个项目。作者感谢Chris Greenman博士提供了他的软件,并感谢Daniela Gerhard博士和Samuel Aparicio博士的有益讨论和指导。特别感谢Cecelia Suragh、Robyn Roscoe、Armelle Troussard和Adrienne Drobnies提供的专家项目管理协助,以及基因组科学中心的图书馆建设、测序和生物信息学团队。本出版物的内容不一定反映卫生与公共服务部的政策观点,提及商品名称、商业产品或组织也不意味着美国政府的认可。可以使用DCC门户访问所有原始序列数据(http://cgap.nci.nih.gov/数据访问),研究登录phs000235.v2.p1。
脚注
作者贡献:MAM、RDG、DEH、MH和JMC构思了该研究并领导了实验设计。RDM对序列数据进行了分析,确定了突变,并与MML、AJM和MAM一起制作了数字并撰写了手稿。MML、AJM、DLT、SC、SC、DS、HM、JS、MM、TZ、AD、KT、YB、MF、JTW和TMS设计并进行了扩增、发现和验证突变的实验。RG、MG和IMM参与分析并审阅了手稿。NAJ、MB、BW和BM准备了样本,进行了样本分类和COO分析,并对文本做出了贡献。ABW和JJS收集并制备了宪法DNA样本。生成KLM、RC、SL、MF和SJ从头开始组装和鉴定突变。MK、SR、MG、OY和EYZ编写了软件并为Figures做出了贡献。RC进行拷贝数分析并生成Figure,SBN进行验证性FISH实验。YZ和AT产生了测序库。IB、RH、SJMJ、RM、JS、MH为实验和分析方案的开发做出了贡献。LR提供了材料并审阅了手稿。
作者声明没有竞争性的经济利益。
提交的SRA登录号为SRP001599,与dbGAP研究登录号phs000235.v2.p1链接。
参考文献
1Anderson JR,Armitage JO,Weisenburger DD。非霍奇金淋巴瘤的流行病学:主要亚型的分布因地理位置而异。非霍奇金淋巴瘤分类项目。安·昂科尔。1998;9:717–720.[公共医学][谷歌学者] 2Lenz G,Staudt LM。侵袭性淋巴瘤。N英格兰医学杂志。2010;362:1417–1429. [PMC免费文章][公共医学][谷歌学者] 三。Horsman DE等。缺乏t(14;18)(q32;q21)的滤泡性淋巴瘤:两种疾病亚型的鉴定。英国血液学杂志。2003;120:424–433.[公共医学][谷歌学者] 4Iqbal J等。BCL2易位在生发中心B细胞样弥漫性大B细胞淋巴瘤中定义了一个独特的肿瘤亚群。《美国病理学杂志》。2004;165:159–166. [PMC免费文章][公共医学][谷歌学者] 5Lenz G等。弥漫性大B细胞淋巴瘤的分子亚型由不同的遗传途径引起。美国国家科学院程序。2008;105:13520–13525. [PMC免费文章][公共医学][谷歌学者] 6Pasqualucci L等。弥漫性大B细胞淋巴瘤中PRDM1/BLIMP1基因的失活。《实验医学杂志》。2006;203:311–317. [PMC免费文章][公共医学][谷歌学者] 7Kato M等。B细胞淋巴瘤中A20的频繁失活。自然。2009;459:712–716.[公共医学][谷歌学者] 8Compagno M等人。在弥漫性大B细胞淋巴瘤中,多个基因的突变导致NF-kappaB的失调。自然。2009;459:717–721. [PMC免费文章][公共医学][谷歌学者] 9Davis RE等。弥漫性大B细胞淋巴瘤中的慢性活性B细胞受体信号。自然。2010;463:88–92. [PMC免费文章][公共医学][谷歌学者] 11Mardis ER等。通过对急性髓细胞白血病基因组进行测序发现的复发突变。N英格兰医学杂志。2009;361:1058–1066. [PMC免费文章][公共医学][谷歌学者] 12Shah SP等。单核苷酸分辨率下小叶乳腺肿瘤的突变进化。自然。2009;461:809–813.[公共医学][谷歌学者] 13Morin RD等。生殖中心起源的滤泡性和弥漫性大B细胞淋巴瘤中改变EZH2(Tyr641)的体细胞突变。自然遗传学。2010;42:181–185. [PMC免费文章][公共医学][谷歌学者] 15Pasqualucci L等。B细胞淋巴瘤乙酰转移酶基因失活突变。自然。2011;471:189–195. [PMC免费文章][公共医学][谷歌学者] 16Yusuf I,Zhu X,Kharas MG,Chen J,Fruman DA。最佳的B细胞增殖需要磷酸肌醇3-激酶依赖性的FOXO转录因子失活。鲜血。2004;104:784–787.[公共医学][谷歌学者] 17Saito M等。BCL6通过Miz1抑制BCL2及其在弥漫性大B细胞淋巴瘤中的破坏。美国国家科学院程序。2009;106:11294–11299. [PMC免费文章][公共医学][谷歌学者] 18Lenz G等。人类弥漫性大B细胞淋巴瘤中的致癌CARD11突变。科学。2008;319:1676–1679.[公共医学][谷歌学者] 19Greenman C、Wooster R、Futreal PA、Stratton MR、Easton DF。癌症体细胞突变致病性的统计分析。遗传学。2006;173:2187–2198. [PMC免费文章][公共医学][谷歌学者] 20Cheung KJ等。滤泡性淋巴瘤中获得性TNFRSF14突变与不良预后相关。癌症研究。2010;70:9166–9174.[公共医学][谷歌学者] 21杜MQ等。淋巴瘤中BCL10基因突变。鲜血。2000;95:3885–3890.[公共医学][谷歌学者] 22Kreutz B、Hajicek N、Yau DM、Nakamura S、Kozasa T。Galpha13的不同区域参与其与含有RGS同源结构域的RhoGEFs的调节相互作用。细胞信号。2007;19:1681–1689.[公共医学][谷歌学者] 23Bhattacharyya R,Wedegaertner P.Galpha 13要求棕榈酰化用于质膜定位、Rho依赖性信号传递和促进p115-RhoGEF膜结合。生物化学杂志。2000;275:14992–14999.[公共医学][谷歌学者] 24Manganello JM、Huang J、Kozasa T、Voyno-Yasenetskaya TA、Le Breton GC。蛋白激酶A介导的Galpha13开关I区磷酸化改变Galphabetagamma13-G蛋白偶联受体复合物并抑制Rho活化。生物化学杂志。2003;278:124–130.[公共医学][谷歌学者] 25Brunet A等。蛋白激酶SGK通过磷酸化叉头转录因子FKHRL1(FOXO3a)介导生存信号分子细胞生物学。2001;21:952–965. [PMC免费文章][公共医学][谷歌学者] 26Tai DJC、Su C、Ma Y、Lee EHY。IkappaB激酶α和p300的SGK1磷酸化上调NF-κB活性并增加N-甲基-D-天冬氨酸受体NR2A和NR2B的表达。生物化学杂志。2009;284:4073–4089.[公共医学][谷歌学者] 27Mo J等。血清和糖皮质激素诱导激酶1(SGK1)通过Fbw7泛素连接酶下调蛋白质稳定性来控制Notch1信号。细胞科学杂志。2011;124:100–112.[公共医学][谷歌学者] 28Young KH等。TP53基因突变的结构特征预测弥漫性大B细胞淋巴瘤的临床结局:一项国际合作研究。鲜血。2008;112:3088–3098. [PMC免费文章][公共医学][谷歌学者] 29Shilatifard A.组蛋白H3赖氨酸4(H3K4)甲基化的分子实现和生理作用。细胞生物学的最新观点。2008;20:341–348. [PMC免费文章][公共医学][谷歌学者] 30Milne T等。MLL将SET域甲基转移酶活性靶向Hox基因启动子。分子细胞。2002;10:1107–1117.[公共医学][谷歌学者] 31脊椎动物发育中的Krumlauf R.Hox基因。单元格。1994;78:191–201.[公共医学][谷歌学者] 32Canaani E等,ALL-1//MLL1,果蝇TRITHORAX的同源物,修饰染色质,直接参与婴儿急性白血病。英国癌症杂志。2004;90:756–760. [PMC免费文章][公共医学][谷歌学者] 33Wiedemann L等。H3K4甲基化的全球分析定义了MLL家族成员靶点,并指出了MLL1介导的H3K4-甲基化在RNA聚合酶II转录起始调控中的作用。分子细胞生物学。2009;29:6074–6085。 [PMC免费文章][公共医学][谷歌学者] 34Issaeva I等。敲除ALR(MLL2)揭示ALR靶基因并导致细胞粘附和生长的改变。分子细胞生物学。2007;27:1889–1903. [PMC免费文章][公共医学][谷歌学者] 35Pleasance ED等人。一种具有烟草暴露复杂特征的小细胞肺癌基因组。自然。2010;463:184–190. [PMC免费文章][公共医学][谷歌学者] 36Dalgliesh GL等。肾癌的系统测序揭示了组蛋白修饰基因的失活。自然。2010;463:360–363。 [PMC免费文章][公共医学][谷歌学者] 38Iqbal J等。弥漫性大B细胞淋巴瘤不同亚群中BCL6分子改变的不同模式及其功能后果。白血病。2007;21:2332–2343. [PMC免费文章][公共医学][谷歌学者] 39Pasini D等。Polycomb组靶基因转录调控中组蛋白H3赖氨酸27甲基化和乙酰化之间拮抗开关的特征。核酸研究。2010年doi:10.1093/nar/gkq244。 [PMC免费文章][公共医学] [交叉参考][谷歌学者] 40Giordano A、Avantaggiati M.p300和CBP:生死伴侣。细胞生理学杂志。1999;181:218–230.[公共医学][谷歌学者] 41Han A,He J,Wu Y,Liu JO,Chen L.心肌细胞增强因子-2募集II类组蛋白脱乙酰酶的机制。分子生物学杂志。2005;345:91–102.[公共医学][谷歌学者] 42Youn H,Liu J.Cabin1通过控制组蛋白脱乙酰化酶和乙酰化酶与MEF2的关联,抑制MEF2依赖性Nur77表达和T细胞凋亡。免疫。2000;13:85–94.[公共医学][谷歌学者] 43Yap DB等。EZH2 Y641的体细胞突变主要通过选择性改变PRC2催化活性的机制发挥作用,以增加H3K27三甲基化。鲜血。2011;117:2451–2459. [PMC免费文章][公共医学][谷歌学者] 44Snenger CJ,等。野生型加突变型EZH2的协同活性驱动人类B细胞淋巴瘤组蛋白H3(H3K27)上赖氨酸27的肿瘤相关超三甲基化。美国国家科学院程序。2010;107:20980–20985. [PMC免费文章][公共医学][谷歌学者] 45Li H,Durbin R.使用Burrows-Wheeler变换快速准确地进行短读对齐。生物信息学。2009;25:1754–1760. [PMC免费文章][公共医学][谷歌学者] 46Goya R等人。SNVMix:预测下一代肿瘤测序中的单核苷酸变体。生物信息学。2010;26:730–736. [PMC免费文章][公共医学][谷歌学者] 47Robertson G等人。RNA-seq数据的从头组装和分析。Nat Meth公司。2010;7:909–912.[公共医学][谷歌学者] 48Mortazavi A、Williams BA、Mccue K、Schaeffer L、Wold B.通过RNA-Seq对哺乳动物转录体进行定位和量化。Nat Meth公司。2008;5:621–628.[公共医学][谷歌学者] 49Wright G等。一种基于基因表达的方法,用于诊断弥漫性大B细胞淋巴瘤的临床不同亚群。美国国家科学院程序。2003;100:9991–9996. [PMC免费文章][公共医学][谷歌学者] 50He J,等。DNA上与MEF2结合的p300的结构揭示了增强子组装的机制。核酸研究。2011年doi:10.1093/nar/gkr030。 [PMC免费文章][公共医学] [交叉参考][谷歌学者]