雌激素受体(ER)是一个关键分子乳腺癌的象征。分子ER配合物与DNA之间的相互作用调节基因的表达负责癌症细胞表型。然而,ER绑定的位置和机制与下游基因的目标还远未被完全理解的。ChIP-Seq的全基因组研究是一个重要的试验protein-DNA交互。在本文中,我们探讨了全基因组染色质的本地化ER-DNA绑定地区通过分析ChIP-Seq MCF-7乳腺癌细胞系的数据。通过整合三个峰值检测算法和两个数据集,我们本地化933 ER结合位点,其中92%位于远离推动者,表明远程控制ER。此外,489个基因附近的ER结合位点被确定为雌激素反应元件与表达数据进行比较。此外,836个单核苷酸多态性(snp)在或接近157 ER-regulated基因被发现附近的结合位点。此外,我们最近邻基因的功能注释这些结合位点使用基因本体论(去),KEGG, GeneGo通路数据库。结果显示小说ER-regulated基因通路进行进一步的实验验证。 ER was found to affect every developed stage of breast cancer by regulating genes related to the development, progression, and metastasis. This study provides a deeper understanding of the regulatory mechanisms of ER and its associated genes.
乳腺癌是一种复杂的疾病发生。它涉及范围广泛的和多样化的临床病理实体课程。基因和蛋白表达都进行了广泛的描述在不同亚型乳腺癌[
一些先前的研究已经致力于ER-regulated基因和它们的功能在乳腺癌细胞系(
ChIP-Seq数据分析管道。
ChIP-Seq相关的乳腺癌基因表达数据集提取综合(GEO): GSE19013 [
CHIP-Seq数据集。
| 数据集 | 平台 | 细胞系 | 样品信息 |
|---|---|---|---|
| GSE19013 | Illumina公司 | MCF-7 | 乙醇处理 |
| E2-treated | |||
| GSE14664 | Illumina公司 | MCF-7 | ER_minus_ligand |
| ER_E2 |
领结(
峰值检测算法对ChIP-Seq数据集的分析是至关重要的。目前,有几个工具可以用来识别全基因组转录因子的结合位点,如FindPeaks [
不同的特点概述Chip-Seq峰值检测算法。
| 算法 | 配置文件 | 背景模型 | 控制样品 | 使用控制计算罗斯福 |
|
|
||||
| F-Seq | 核密度估计(KDE) |
|
||
| FindPeaks | 聚合重叠的标签 | 蒙特卡罗 | ||
| SISSRs | 窗口扫描 | 泊松 |
|
|
| 追求 | 核密度估计(KDE) |
|
|
|
| 苹果电脑 | 标记转移然后窗口扫描 | 动态泊松 |
|
|
| CisGenome | Strand-specific窗口扫描 | 负二项 |
|
|
此外,MEME程序(
表达分析使用相同的包(
使用表SNP (131) (dbSNP建造131)[
三个功能注释系统,基因本体论(去)类别
浓缩的类别确定的基因本体树机(GOTM) [
使用ChIP-Seq数据集,我们确定了全球ER结合位点。人类基因组序列标签被首先对齐组装(UCSC hg19)使用领结。三个ChIP-Seq峰调用程序,CisGenome mac电脑,和追求,选择确定丰富绑定的山峰。使用错误发现率为0.01,933 ER绑定的山峰被揭示了所有的数据集(表三种工具
ER的结合位点被三个ChIP-Seq峰调用程序(罗斯福< 0.01)。
| ER结合位点的数量 | |||||
| 数据集 | 重叠的站点数量 | ||||
| CisGenome | 苹果电脑 | 追求 | |||
|
|
|||||
| GSE19013 | 8137年 | 5583年 | 5418年 | 2019年 | 933年 |
| GSE14664 | 6773年 | 7765年 | 9280年 | 5061年 | |
比较追求、CisGenome和mac预测结果。(一)罗斯福GSE19013数据集的价值。(b)的罗斯福价值GSE14664的数据集。
ER的基因组结合位点。(a)中标识的共识主题之前绑定的网站。新创主题搜索使用MEME程序执行。(b)的百分比出现之前图案的结合位点。(c)比较之前出现的主题发表和新发现的结合位点之间。
此外,我们研究了ER浓缩网站相对的位置最近邻的基因。结果是图所示
位置ER结合位点的分析。(a)的位置相对于加权的基因。(b)基因的位置ER ChIP-Seq峰值。(c)基因的位置ER ChIP-Seq山峰内+ 10 ~ + 100 kb。
为了确定相对应的特定基因反应在MCF-7细胞ER,我们比较最近邻的ER基因结合位点之间的差异表达基因发表研究ER +和ER−乳腺肿瘤。我们使用在表3的研究
乳腺癌的基因表达数据集和不同表达基因数(
| 作者 | 杂志 | 数组类型 | 样本 |
样本 |
不同表达基因 | |
|---|---|---|---|---|---|---|
| 调节 | 表达下调 | |||||
| 格雷厄姆et al。 |
中国癌症Res | Affy | 15 | 15 | 709年 | 333年 |
| 王等人。 |
《柳叶刀》 | Affy | 209年 | 77年 | 2081年 | 2537年 |
| 陆et al。 |
Res乳腺癌治疗 | Affy | 76年 | 53 | 5136年 | 5445年 |
|
|
||||||
| 所有 | 5692年 | 6101年 | ||||
我们的分析发现,更多的结合位点与ER(60%)相比,下调基因差异基因(40%),表明ER更频繁的参与直接调控基因表达的上调。我们还研究了ER结合位点的位置上调和下调的基因。如图
基因差异表达基因的位置附近的结合位点。
目前的研究表明,乳腺癌风险相关的常见的单核苷酸多态性(snp) [
与基因表达的不同设置附近的结合位点,157年或接近836个snp ER-regulated基因被确定(见附加文件2)。大部分的单核苷酸多态性(94.5%)位于内含子和翻译区域。只有5.5%的地区位于near-gene, coding-synon,错义和转移。这些单核苷酸多态性可能与乳腺癌的关系密切。
确定生物过程和途径改变,呃,我们使用三个功能注释系统,基因本体论(去)类别
获得的生物过程的概述近邻ER结合位点的基因存在,我们首先进行了基因集富集分析使用基因本体数据库。统计学意义(超几何测试,
上丰富的比较之间的类别不同的ER基因表达和其他加权结合位点(≥100)基因的数量。
| 基因集 | 生物过程 | 分子功能 | 蜂窝组件 |
|---|---|---|---|
| 不同的表达 | 生物调控、代谢过程, |
蛋白质绑定, |
膜, |
| 其他人 | 生物调控、代谢过程 | 蛋白结合 | 膜 |
有向无环图(熟练的技艺)显著富集的(基因本体)类别(
KEGG通路数据库(2011年5月23日发布)被用来识别功能模块由ER。17大大丰富通路(
KEGG通路富含最近邻ER结合位点的基因(
| KEGG ID | 路径名称 |
|
数量的基因 | 许多不同的基因表达 |
|---|---|---|---|---|
| hsa05200 | 通路在癌症 |
|
22 | 16 |
| hsa04510 | 粘着斑 | 0.0002 | 15 | 14 |
| hsa04360 | 轴突的指导 | 0.0009 | 11 | 8 |
| hsa04810 | 调节肌动蛋白细胞骨架 | 0.0012 | 14 | 11 |
| hsa04010 | MAPK信号通路 | 0.0022 | 15 | 12 |
| hsa04114 | 卵母细胞减数分裂 | 0.0024 | 9 | 8 |
| hsa04144 | 内吞作用 | 0.0024 | 12 | 11 |
| hsa04115 | p53信号通路 | 0.0024 | 7 | 7 |
| hsa05216 | 甲状腺癌 | 0.0024 | 5 | 4 |
| hsa05218 | 黑素瘤 | 0.0033 | 7 | 3 |
| hsa04020 | 钙信号通路 | 0.004 | 11 | 4 |
| hsa04062 | 趋化因子信号通路 | 0.0064 | 11 | 9 |
| hsa04914 | Progesterone-mediated卵母细胞成熟 | 0.0085 | 7 | 7 |
| hsa01100 | 代谢途径 | 0.0086 | 35 | 28 |
| hsa00450 | Selenoamino酸代谢 | 0.0088 | 4 | 3 |
| hsa05414 | 扩张型心肌病 | 0.0096 | 7 | 7 |
| hsa03440 | 同源重组 | 0.0097 | 4 | 3 |
GeneGo也用于执行路径分析。十个通路被发现明显富集
丰富GeneGo通路的地图(
| GeneGo途径方面 |
|
|---|---|
| 细胞凋亡和survival_APRIL高飞球的一击信号 |
|
| Development_prolactin受体信号 |
|
| Development_glucocorticoid受体信号 |
|
| Development_ligand-independent激活ESR1 ESR2 | 0.000295251 |
| 免疫response_IL-22信号通路 | 0.000381484 |
| Development_EPO-induced Jak-STAT通路 | 0.000531744 |
| 通过统计和PLC / IP3 Development_growth激素信号 | 0.000531744 |
| 细胞骨架remodeling_keratin细丝 | 0.000622315 |
| Development_GM-CSF信号 | 0.000660576 |
| Transcription_transcription胺基酸代谢的调节 | 0.000752764 |
呃是一个重要的乳腺癌的分子标志。充分理解ER的分子机制将对研究有用的预测和治疗乳腺癌。ChIP-Seq技术是有用的研究蛋白质和DNA之间的相互作用在全基因组范围内。ChIP-Seq数据可以有效地分析转录因子的调控机制在全基因组规模。在这项研究中,我们使用ChIP-Seq数据识别全球网站受ER MCF-7乳腺癌细胞系。为了得到更可靠的结果,三种不同的工具被用来分析两个数据集。933结合位点被发现,和之前的主题是精制。
全球ER-regulated基因的基因组入住率的分析显示,92%的总933 ER-binding网站位于远离推动者。这表明ER因子函数的规范模式涉及的远程控制。之前的研究已经报道,ER -
我们比较了ER结合位点之间的差异表达基因发现和发表在这项研究ER +和ER−乳腺肿瘤。一组487个基因被发现在歧视呃重要地位在乳腺肿瘤。这似乎表明,这些基因影响ER的回应。只有9%(44)的基因发现了林等。
在这项研究中,我们发现2694个单核苷酸多态性位点位于或接近ER结合位点。在这些snp, 836个snp的157个基因之间差异表达ER +和ER−乳腺肿瘤。它表示,这组单核苷酸多态性可能与ER在乳腺关系密切。
的功能注释提供了一个深入了解ER和ER-associated基因。浓缩的分析去了基因功能的概述。如图
作者宣称没有利益冲突。
这项工作得到了国家自然科学基金资助(没有。91230117和91230117),专门研究中国高等教育的博士项目基金(20113201110015)、苏州(SH201120),国际科技合作项目和国家高技术研究发展计划(863计划,批准号2012 aa02a601)。