Abstract
随着组学的出现接近我们的慢性疾病，如癌症和代谢综合征的认识有所提高。但是，这是由基因表达微阵列，深度测序实验或代谢分析获得大规模数据集的信息的有效开采是必不可少的发现，然后有效的病变靶细胞表型的关键调节因子。雌激素受体α（ERα）是调节这对于雌激素反应乳腺癌重要的基因程序的主转录因子中的一个。为了理解到的ERα在乳腺癌代谢信令作用我们从雌二醇处理的MCF-7细胞利用转录，cistromic和代谢数据。在这份报告中，我们描述了一代样品进行RNA测序，芯片起和代谢组学实验和获得的数据综合计算分析。这种方法是在描绘小说摩尔有用丘拉尔机制和基因调节电路是由一个特定的转录因子调节的正常或患病细胞的这影响代谢。
Introduction
雌激素是两个女性和男性包括生殖组织，组织代谢，脑和骨1许多生理过程的重要调节。除了在这些组织中的有益作用，雌激素也驱动，从乳腺和生殖组织出现的癌症。雌激素主要是合作，通过减排额来诱导细胞类型特异性作用。通过使用芯片起RNA测序和全基因组ERαDNA结合位点分析ERα调节成绩单深度测序被证明是有用的了解ERα是如何工作的，从他们的出现的不同的组织和癌症。我们和其他人发表不同的受体（ERαVSERβ）相关基因表达谱2,3，不同的配体3-5和不同的共调节因子-2,4,6,7-。
RNA测序的主要方法检测转录，提供更高的精度和效率相比，芯片巴SED基因表达分析8。从细胞系2-4,7，组织或肿瘤样品获得的RNA进行测序，映射到可用的基因组组件和差异调节的基因被鉴定出来。染色质免疫沉淀（ChIP）被用来解剖转录因子和coregulator染色质结合已知调节结合位点。芯片起（芯片后跟高通量测序）提供全球性的结合位点的检测偏见。代谢组学是另一个越来越多地采用系统生物学的方法，它定量措施，生命系统的各种刺激，包括化学和遗传扰动的动态响应多参数。
通过进行全局代谢轮廓，功能性读出可以从细胞，组织和血液中获得。此外，从转录实验信息并不总是反映在有助于生化途径的酶的水平的实际变化。综合转录和代谢分析数据，使我们能够识别并关联与实际代谢物的改变的基因表达的变化。利用从所有这些大型数据集提供的机械细节的信息，了解调控复杂生物系统，特别是那些涉及到人类发展和疾病，如癌症和糖尿病的转录因子的作用。
哺乳动物基因组的复杂性使得它具有挑战性的整合并充分解释从转录，顺反组和代谢组实验中获得的数据。识别功能的结合事件，这将导致的变化，因为一旦功能性结合位点识别的目标基因的表达是重要的;随后的分析，包括转录结合（TF）的基序分析，可以以更高的精度来执行。这导致生物学上有意义的TF级联和机制的鉴定。也可直接comparis上的RNA测序和芯片起实验并不总是可能的，因为从每个实验的数据具有不同的尺度和噪音，在某些情况下，有意义的信号是由人口噪声掩盖。我们不知道，集成这三个独立而又相关的方法信息，了解乳腺癌ERα由直接代谢调控任何研究的。因此，我们在本文的总体目标是生产关系的结合事件的基因表达和代谢物的变化。为了实现这一目标，我们从RNA测序，芯片起和代谢实验的综合数据，并确定那些诱导的ERα的雌激素结合，这将导致在代谢途径基因表达的变化的事件。这是第一次，我们提供了一套完整的生成芯片起，RNA测序和代谢组学分析和数据进行综合分析方案（ 图1）来发现新的基因调节电路乳房的新陈代谢癌细胞系。
Subscription Required. Please recommend JoVE to your librarian.
Protocol
1. RNA测序样品的制备
- 细胞培养和治疗
- 培养的MCF7细胞中改进的MEM（IMEM）加了10％热非活性的FBS，1％青霉素/在湿润的5％CO 2的气氛链霉素，在37℃的酚红培养基。
注意：相同的细胞系和细胞培养条件在整个研究中使用。 - 种子100000 MCF7细胞每孔在6孔板中。对每个治疗一式三份。在第3天，除去培养基，加2ml含或不含10 -8 M的E2向每个孔新鲜培养基中。孵育所述板在湿润的5％CO 2的气氛中在37℃24小时。
- 收获细胞，加入1ml RNA分离试剂（酸性硫氰酸胍 - 苯酚 - 氯仿）到每个孔中，并孵育在室温（RT）下5分钟。然后收集吹打并沉积在细胞裂解液试剂进入1.5毫升管
- 培养的MCF7细胞中改进的MEM（IMEM）加了10％热非活性的FBS，1％青霉素/在湿润的5％CO 2的气氛链霉素，在37℃的酚红培养基。
- RNA分离
- 加入200μlchlorofoRM到1ml细胞溶解物试剂（1.1.3），涡旋10秒。在室温下孵育5分钟。离心样品以16,000 xg离心在4℃下15分钟。离心后，将水相（上部）转移到新的管中，而不会干扰其它相。
- 将500μl100％异丙醇添加到与所传送的水相的新管中，并通过颠倒混合。孵育在-20℃过夜。离心样品以16,000 xg离心在4℃下10分钟。在管的底部观察的颗粒。弃去上清液并通过简短涡旋洗净用750微升不含RNA酶的70％的乙醇沉淀。
- 离心样品在4℃6000×g离心5分钟。除去上清液，并在4℃下以16,000 xg离心离心样品1分钟以收集所有剩余的乙醇。通过移液用凝胶装载尖除去剩余的乙醇，并设置在管上侧向下风干10分钟。溶解在依赖新生20-50微升DEPC处理水沉淀丁对丸粒的大小。
- 使用RNA清理试剂盒按照协议4纯化的RNA。使用基于荧光检测试剂盒按照协议9测量RNA浓度。
注意：OD 260nm处取来确定在260nm的OD的RNA浓度和比率和280nm被用来确定样品的纯度。建议检查使用Bio-分析器测定RNA的质量。 - 取用于测序约0.5至1微克纯化的RNA。
注：生物技术中心准备测序文库，并执行使用以前2详细方法配对末端测序读。
2.芯片起样品的制备
- 细胞培养和治疗
- 种子每10cm盘500000 MCF7细胞在生长培养基中。
注意：对于每个拉下使用了16片。在第3天，除去培养基并加入5毫升或不含10 -8 M的E2到每个新鲜培养基盘子。处理细胞45分钟（37℃，在潮湿的5％CO 2的气氛中）。
- 种子每10cm盘500000 MCF7细胞在生长培养基中。
- 交联和细胞收获
- 加入250微升37％甲醛入5ml培养基交联染色质，在室温下孵育15分钟。停止洗涤细胞用5ml冰冷的1×MEM交联。
- 收获在每板250μl的细胞裂解缓冲液（10mM EDTA，50mM的三盐酸，1％SDS，0.5％N，N-二dimethyldodecan -1-胺氧化物）的细胞。
- 染色质的碎片
- 超声处理的细胞30秒×8以30放大器来达到所需的染色质的大小。 30秒的超声处理后在冰上冷却各样品。
- 离心机在16，000×g离心在4℃下10分钟，并小心地吸取上清液，而不在底部干扰颗粒。
- 取上清的5微升等分试样，并通过1.5％琼脂糖凝胶检查如先前10描述的电泳的DNA大小。理想的大小的DNA应当是200-500碱基对之间。
- 染色质免疫沉淀
- 保存25微升细胞裂解液作为输入。在50ml管中，将4- ml细胞裂解液用20毫升的IP缓冲液（2mM EDTA的，100毫摩尔NaCl，20毫摩尔Tris-HCl和0.5％的Triton X），ERα抗体（500微升F10和500μl的HC- 20），和500微升PROTA和ProtG蛋白涂覆的磁珠。孵育通过旋转该混合物在4℃下过夜，使染色质 - 抗体复合物的形成。
- 清洗和反向交联
- 用磁力架收集关于磁化侧磁珠去除洗涤缓冲液。
- 洗磁珠在25ml洗涤缓冲液I（2毫摩尔EDTA，20毫摩尔Tris盐酸，01％SDS，1％的Triton X，和150毫摩尔NaCl）。彻底清洗磁珠不干扰DNA的蛋白质复合物，加入25 1ml洗涤慢慢缓冲到沿着侧的管然后反转管，直到所有的磁珠混合以及无线TH洗涤缓冲液没有出现一个小球。作为替代，使用旋转器以反转管。不要涡旋样品。用同样的方法在下面的步骤洗涤磁珠。
- 在25ml洗涤缓冲液II（2毫摩尔EDTA，20毫摩尔Tris盐酸，01％SDS，1％的Triton X，和500mM的氯化钠）洗涤磁珠。
- 洗磁珠在25ml洗涤缓冲液III（1.6毫摩尔EDTA，10毫摩尔Tris盐酸，1％NP-40，1％Dexycholate和250mM的氯化锂）。
- 在25ml 1×TE缓冲液（1mM EDTA和10毫摩尔Tris HCl）的两次洗涤磁珠。
- 洗脱的DNA于500μl洗脱缓冲液（1％SDS，和10mM的NaHCO 3），然后分配500微升的DNA洗脱到4个试管。从投入的DNA，洗脱125微升洗脱缓冲液的DNA。
- 在65℃孵育管过夜加热块上。用铝箔覆盖在加热块上的管，以保持热甚至在所有的管子。将一个砝码放在管，以防止他们打开。</ LI>
- DNA提取
- 冷却在RT管5分钟。隔离使用DNA芯片分离试剂盒下拉样本的DNA。请遵循制造商的指示。
- 预热在65℃的洗脱缓冲液。添加625微升结合缓冲液至每个管。如果形成沉淀，再添加入250μl结合缓冲液，直至溶液澄清。
- 加载样品上柱和离心机以16,000×g离心30秒。洗净用250微升洗涤液样品。记得要加乙醇洗涤缓冲液，当第一次使用它。重复洗。
- 洗脱在15微升洗脱缓冲液将DNA。联合收割机一样的DNA，从4管拔下来，实现约55微升的DNA。
- 隔离输入DNA用DNA纯化试剂盒。
- 预热在65℃的洗脱缓冲液。添加625微升结合缓冲液中的每个管中并孵育10分钟。加载样品到柱上在2 ml collecti上管。
- 洗用750微升洗涤缓冲液，离心样品以16,000 xg离心1分钟。重复离心以除去任何剩余的缓冲区。洗脱DNA插入EB的30微升。
- 测量DNA浓度使用荧光染料的商业与试验修改11。
- 稀释20倍TE缓冲到1倍的TE，作为稀释试剂和DNA样品的测定缓冲液。取5微升提取的DNA，并将其稀释成50微升的TE 1X。
- 从2微克/毫升DNA原液进行DNA标准从毛坯到1，10，100，1000微克/毫升。在一个96孔板，分配每个标准的25微升以及在重复DNA稀释。 （推荐标准一式三份）。
- 通过在塑料制管，从双链DNA试剂的200倍稀释于1×TE缓冲液制备工作试剂。用铝箔覆盖工作试剂溶液，以避免光。加入200μl工作试剂到每个样本。在cubate在室温下5分钟。用铝箔覆盖该板，以避免光。
- 测量使用荧光板读数器（激发在480nm，发射520纳米）的样品的荧光。生成花期与DNA浓度的DNA标准曲线。确定基于所产生的DNA的标准曲线对样品的浓度。
- 冷却在RT管5分钟。隔离使用DNA芯片分离试剂盒下拉样本的DNA。请遵循制造商的指示。
- 沉淀实验验证使用qPCR
- 取2微升提取的DNA和稀释成20微升水。通过混合2.5微升DNA稀释液，5微升绿我PCR预混，1微升底漆和1.5微升无核酸酶水设置一式三份的qPCR实验。
- 使用以下程序的快速实时PCR系统上运行的qPCR：步骤1：95℃60秒，第2步：95℃，10秒，60℃60秒，重复执行步骤2 39次。
- 送10纳克DNA样本，以生物技术中心进行测序。
注：生物技术H中央根据指令准备DNA芯片进入图书馆，并执行使用以前2详细方法单读测序。
3.代谢组学分析样品制备
- 细胞培养和治疗
- 种子每10cm平板400000 MCF7细胞在生长培养基中。每个处理制备三个样品。使用每个样品中两个板块获得足够的细胞用于检测。在第3天，除去培养基并加入5毫升或不含10 -8 M的E2到细胞的新鲜培养基。处理细胞24小时（37℃下在湿润的5％CO 2的气氛中）。
- 样品采集
- 加入5毫升冰水预冷的1×PBS的板块，经过几次短暂倾斜板吸出PBS。重复两次，在最后一次洗涤后除去尽可能多的PBS越好。
- 添加750微升预冷的丙酮至各板和刮细胞。结合细胞的丙酮从两个板在冰上2ml的管中。
- 正常化细胞计数
- 对于每个处理，使用细胞量化两个额外板。从板到解附着细胞，加入2毫升，HE缓冲液（1×HBSS，10mM的HEPES，0.075％ 碳酸氢钠，和1mM EDTA）插入每个板并孵育5分钟。
- 收集并通过在HE缓冲吹打细胞拌匀细胞。使用总共20微升细胞溶液到计数使用血球细胞数。
- 发送到代谢中心来识别和使用气相色谱质谱（GC-MS）分析量化代谢物之前储存于-80℃的样品。
4.综合分析
- RNA测序数据分析：
- 执行使用FastQC FASTQ文件质量检查：读QC（http://galaxy.illinois.edu）
- 选择FastQC：NGS的工具下阅读：QC和操纵。载从当前histor原始数据文件年。
注：建议给一个标题为输出文件提醒你什么工作了，因为对于可能有多个输出。
- 选择FastQC：NGS的工具下阅读：QC和操纵。载从当前histor原始数据文件年。
- 执行和输出文件将历史记录下。重复相同的步骤，对每个序列文件。修剪适配器使用NGS下剪辑：FASTX工具包。
- 对齐读取使用Tophat2与RefSeq的基因组参考注释12（默认值）hg19。
- 根据NGS的工具选择Tophat2：RNA分析。指示库配对类型（单端VS末端配对）。
- 上传从目前的历史上输入FASTQ文件。选择参考基因组。选择默认设置。或使用完整的参数表进行修改。
- 执行和它会产生两个输出文件：一个是路口UCSC床轨道;另一种是在BAM格式读出比对的列表。
- 上传BAM文件测序数据分析工具。计算使用定量表达值ification工具。
- 确定使用默认值使用表达分析工具差异上调的基因和倍数变化> 2且p值 <0.05。同样地，确定使用默认值使用表达分析工具差异下调的基因和倍数变化<0.5和p值<0.05。
- 执行使用FastQC FASTQ文件质量检查：读QC（http://galaxy.illinois.edu）
- 芯片起数据分析：
- 对齐FASTQ读取到银河系中的使用Bowtie2 13（默认值）hg19。
- NSC的工具下选择Bowtie2：映射。验证该库是队友配对（单端VS末端配对）。上传从目前的历史FASTQ文件。
- 选择参考基因组。使用默认的参数设置。执行和它会产生输出文件，映射顺序读为BAM文件。重复该过程的每个序列文件。
- 识别使用MACS 14 6.0E-7的p值截止和0.01 FDR峰，如前所述2。
- 代谢组学数据分析：
- 转换代谢产物的化学名称KEGG_IDs。
- 通过比较在车辆检测VS。E2水平鉴定差异调节的代谢物处理的样品（2倍截止在本研究中使用）。
- 综合分析：
- 识别已经使用的测序数据分析工具ERα结合位点的基因。翻译地区使用20 KB作为距离截止基因的功能。
- 使用倍数变化> 2确定从RNA测序数据分析差异上调和下调的基因和倍数变化<0.5，分别与p值 <0.05为切断。比较使用维恩图比较ERα结合位点含有基因的列表中。
- 使用此列表中使用的Cytoscape的Metscape 15模块识别E2调制的代谢途径。
- 通过从应用Metscape模块中，选择建立网络，然后途径-BAS编辑选项。
- 选择生物体（人类）。 E2无管制化合物列表中选择的数据文件（KEGG ID），以及E2上调基因列表。
- 选择化合物-Reaction - 酶 - 基因作为网络类型。选择化合物/基因作为查询。
- 构建网络。重复使用E2下调化合物和基因。
- 对齐FASTQ读取到银河系中的使用Bowtie2 13（默认值）hg19。
Subscription Required. Please recommend JoVE to your librarian.
Representative Results
转录
分析差异表达的基因由E2治疗，我们选择了要执行的RNA测序实验。除了提供关于mRNA水平的信息，RNA测序数据也可被用来监测在非编码RNA（长的非编码RNA，微小RNA）和选择性剪接事件的变化。我们没有提供的非编码RNA或备选转录基因的分析信息，因为我们的研究范围是确定蛋白质编码基因是对于代谢调节很重要的。然而，RNA-Seq中获得的差异表达基因事件的信息的一个很好的方式。
为了获得高质量的读出，它获得高纯度的RNA是重要的。分离RNA后，我们使用纯化RNA一个清理工具包中的RNA。不幸的是近50％的RNA的在此过程中会丢失，起始RNA的量应采取到c中onsideration以获得由纯化步骤结束时所需的量（100-500纳克）。
我们量化使用荧光计的RNA HS测定，这是非常准确的，是低丰度RNA的定量优选的方法中的RNA产量。接着，将RNA样品的完整性和整体质量使用bioanlayzer，这是使用RNA最小量的凝胶电泳一个可行的替代进行了检查。分析表明所有样品都18S和28S个rRNA的锋利和清洁带指示样本（图2B）的完整性和纯度。 RNA完整性（RIN）用作RNA的质量控制标准。 RIN 10表示完整的RNA，RIN 6部分降解和RIN3强烈降低16。
它建议RIN为至少7用于测序实验。我们所有的RNA样品收到9.6-9.9之间RIN值<STRONG>（图2A）。测序文库从高质量的RNA为500ng采用超高通量测序系统构成从每个样品如Illumina公司2500平均1800万的高品质的测序读取每个样品（表1），这是准备用于下行得到分析。要检查高通量序列原始数据的整体素质，FastQC对每个样品的读取运行。一辆车RNA测序数据的代表性FastQC报告显示（ 图3）。最读取序列长度为约100bp的（ 图3A）。对于每个读，质量得分是32-34的前5 bp和36-38来自6-100bp（ 图3B）。其中17990863读取此样品产生，其中大部分存在质量得分高于34（ 图3C）。在序列中每读GC含量也跟着，平均48％GC的正态分布。整体质量的读取成绩是非常高的。
ERα的结合位点的Cistromic分析
为了实现高效的深度测序，在所希望的尺寸获得的DNA片段是其中的重要因素之一。当前测序方法可能对DNA片段有不同的要求，例如300-600碱基对用于超高通量测序系统，100-200碱基对AB /固17。在该实验的大部分DNA片段为300 bp和600 bp的作为凝胶所示（图4）之间。作为替代方案，读者可能会选择使用MNase消化编写自己的染色质，但是我们觉得我们既定的协议提供高品质要求的样品，并提供可重复的结果。从每个处理的DNA和其对应的输入样本的十纳克在EB缓冲液制备。该芯片是DNA单读取使用超高通量测序系统测序。 DNA芯片样品产生约20000000读取而输入的DNA产生多个35.000,000读取（表1）。
代谢组学
以用GCMS量化在哺乳动物细胞中的代谢物，细胞质量的至少10微克是必需的。我们开始40万MCF7细胞与实验。通过收获的时候，大约有500,000个细胞在每块板。由峰表示的代谢物的概述进行控制和雌二醇处理的样品（图5）之间进行比较。约100代谢物进行鉴定，其计数在哺乳动物已知总数的60％。这里，我们提出的顶部10的代谢物，其显示在E2治疗（表2）显著变化的列表。总体而言，E2上调途径，包括果糖和甘露糖代谢，组氨酸代谢，嘌呤代谢紊乱，胆固醇的生物合成和维生素B3的新陈代谢。该下调的途径是丁metabolis男，脂肪酸从头生物合成，戊糖磷酸途径，尿素循环以及精氨酸代谢（ 表3）。
数据分析和整合
从RNA测序，芯片起和GC-MS的数据进行处理和分析，通过一系列步骤（图6）。在RNA测序和芯片起鉴定的基因之间的比较之后，一组的上调和由雌二醇下调的基因的一组中提取。从代谢数据，我们通过雌二醇处理得到两套上调和下调的化合物。下一步，我们使用Metscape，对于Cytoscape的应用程式，以可视化和解释人体代谢15的上下文中的基因表达和代谢的数据。首先，我们选择了包括在网络中的化合物的元素，反应，酶和基因。我们选择的化合物/基因作为查询和E2上调（ 图7A）的两个网络和向下-regulation（ 图7B）被示出。该Metscape还提供了选择性途径构建网络的选项。例如，使用相同的一组数据，作为构建下调整体网络，通过从化合物/基因的具体途径切换查询生成雄激素的雌激素的生物合成和代谢途径的子网。有显著改变的基因或化合物在亮的颜色得到强调。除了选择从所描述的下拉窗口的特定途径，还可以通过施加“创建子网”功能所选区域创建子网络。如在图6A和B，ERα活化和直接结合到基因调控区示出影响在乳腺癌细胞中许多代谢途径。我们的综合分析表明，ERα直接结合并调节FASN基因，它是脂肪酸从头合成途径重要的表达（Figure 7C）。另外E2治疗通过这些途径抑制CYP4Z1，CYP2E1，CYP4B1，CYP2C8，TPO，GSTA4和GSTM2表达下调雄激素雌激素的生物合成途径。
图1：制剂的RNA，DNA样本和细胞裂解物的用于RNA测序，芯片起，和代谢组分析的工作流程。
图2：与安捷伦2100生物分析仪 RNA完整性（RIN），RNA的质量指标分析RNA的代表电泳 ，被分析为9.9给出。 （ 一 ）在RNA的一个虚拟的凝胶图片样张26S和18S的两个频段是雪亮的。 （B）的叔单个跟踪的荧光信号他的RNA显示，28S峰的强度大于18S高峰，无污染看到。 请点击此处查看该图的放大版本。
图3：在RNA测序原始数据的代表性质量检查其中的车辆样本的序列数据，使用FastQC分析。 （A），它是100碱基对的文库的序列长度分布。 （ 二 ）质量分数沿读取的位置。 （ 三 ）每个序列质量评分显示所有序列的通用高品质。 （D）相匹配的理论分布每次读取GC含量。 请点击此处查看拉rger版本这个数字。
图4：超声处理的染色质片段的凝胶电泳将5μl各样品的细胞裂解物在具有1000 bp的加梯在第一车道1.5％琼脂糖凝胶上运行。大多数碎片是300和BP之间600 请点击此处查看该图的放大版本。
图5：表示从由GC-MS代谢物峰的概况的上面板显示，同时从控制下面板E2类处理的样品中确定的峰。封闭单间是在该地区以显示每个P代表的确切代谢产物 EAK。 请点击此处查看该图的放大版本。
图6：综合聚类方法的流程图。
图7：在Metscape综合基因-代谢物数据的可视化 （A）的途径上调E2的网络。 （B）的途径的网络由E2下调。 （℃）E2刺激ERα招募FASN附近的区域，该区域被发现，通过施加子网函数产生的脂肪酸从头合成网络中来调节。3832fig7large.jpg“目标=”_空白“>点击此处查看该图的放大版本。
表1：RNA 测序和芯片起的超高通量测序的总结。
表2： 排名前10位的代谢产物与E2治疗显著的更改列表。
表3：途径上调和下调由E2治疗。
Subscription Required. Please recommend JoVE to your librarian.
Discussion
在本文中，我们描述了从与E2处理MCF-7细胞生成和RNA测序，芯片起和代谢组学数据的综合分析。我们制定了一套运筹帷幄，利用最有效的方法和用户友好的软商品即产生生物相关的发现协议。据我们所知，我们是从三个不同的分析和ERα在乳腺癌细胞直接调节确定的新的代谢途径整合组学数据的第一项研究。
转录
关键是要开始与高品质的RNA的转录成功分析。好的做法将消除RNA降解的可能性，包括创建无RNase长椅，经常改变手套，在-80存储℃并在使用中保持在冰上。在-20ºC过夜培养的加入提高了RNA回收。为了避免冻融RNA反反复复，可以分配给下游再修改的的的RNA在分离和纯化上课。对测序错误的存在，定量RT-PCR，使用该cDNA从样品中被经常执行以验证RNA测序鉴定的差异表达的基因。因此，足够的RNA应保留用于此目的。我们进行配对末端测序，使我们能够分析剪接变体，非编码RNA和微RNA。为了避免高测序成本，单末端测序也可被执行。相比微阵列，RNA测序更灵敏，并不限定于设置在芯片上的探针。我们得到不仅对编码序列，而且还对非编码基因和差异剪接转录信息。目前，我们正在使用这种技术来分析从小鼠实验中获得的配体的治疗以及组织后，各种细胞系的转录。
Cistromics
超高通量测序方法需要图书馆建设最少5纳克DNA芯片。在标准通道IP-SEQ协议18，建议有至少10纳克。定量用荧光是至关重要的。在这个协议中，我们使用的双链DNA检测分析以确定DNA浓度。人们也可以使用荧光计测定精确DNA定量。对于设备和材料的不同，每个实验室建立自己的超声协议，产生理想的DNA片段的大小是非常重要的。过度剪切DNA芯片可能会导致目标DNA的富集的损失，同时大尺寸芯片DNA创造高效的排序问题。该DNA提交测序之前，它始终是重要通过qPCR来验证一些已知区域的富集。该技术仍然被用于免疫沉淀的抗体的效率和特异性的限制。为了获得能够在比较分析可靠地使用结合位点的信息，我们验证了我们以下ENCODE的ChIP-SEQ标准试剂，19和使用STANDARdized协议为我们提供统一的数据集。再加上后期转录因子结合位点富集分析，芯片起仍然是最好的方法之一，以确定直接结合的转录因子染色质。目前，我们正在使用该协议招聘激酶，共调节因子和其他转录因子的研究在其他细胞系以及雌激素反应小鼠组织染色质。
代谢组学
在人类乳腺癌细胞成功代谢轮廓，仔细的样品制备是至关重要的。在这个协议中使用和细胞培养条件的细胞的量是最佳的MCF7细胞。其它细胞系，可能需要检测的代谢物的一个显著数目或多或少的细胞。我们获得围绕在每个实验中100代谢物。该方法只限于检测那些以低浓度存在短暂的代谢物或代谢物。这种方法提供了方便的USI纳克单个样品制备，以识别多个代谢物。我们目前正在使用这种技术，以确定在血清和小鼠不同组织中的雌激素调节的代谢物。
数据集成
众多转录和cistromics研究，完成调查雌激素信号和基因调节在乳腺癌细胞，包括MCF7细胞2,5-7,20。在这项研究中，在第一次，我们增加代谢数据，并试图推断是直接由ERα在配体除了调节代谢网络。这种新方法使我们能够通过针的雌激素影响细胞代谢调控分子点设计电路。总的来说，我们的研究结果支持ERα的作用，乳腺癌生物学的主要调控。
总之，我们的雌激素治疗和INTEGR后提供协议的详细汇编用于产生RNA测序，从乳腺癌细胞芯片起和代谢样品从这些组学的方法所获得的数据的ative分析。这些协议将使研究人员能够对其他转录因子，核受体配体和营养条件做类似的分析。
Subscription Required. Please recommend JoVE to your librarian.
Disclosures
伊利诺伊大学从辉瑞公司ZME收到调查发起的赠款和YCZ从这笔款项收到薪水的支持。
Acknowledgments
这项工作是由粮食和农业研究所，美国农业部，屡获ILLU-698-909（ZME）和辉瑞公司（ZME）资助。其内容完全是作者的责任，并不一定代表美国农业部的官方意见。
Materials
| Name | Company | Catalog Number | Comments |
| Triglyceride Mix | Sigma | 17810-1AMP-S | |
| Oleic acid | Sigma | O1257-10MG | Oleic Acid-Water Soluble powder |
| TRIzol Reagent | Life technologies | 15596-026 | |
| Dynabeads Protein A | Life technologies | 10006D | |
| Dynabeads Protein G | Life technologies | 10007D | |
| QIAquick PCR Purification Kit | QIAGEN | 28106 | DNA isolation of input sample |
| ZYMO ChIP-DNA Isolation kit | Zymo Research | D5205 | ChIP DNA Clean & Concentrator |
| Quant-iTª PicoGreen¨ dsDNA Assay Kit | Invitrogen | P7589 | DNA Quantitation |
| RNase-Free DNase Set | QIAGEN | 79254 | RNA purification |
| DEPC Treated Water | LIFE TECH | 462224 | Dnase/Rnase Free |
| Model 120 Sonic Dismembrator | Fisher Scientific | FB120110 | With 1/8 probe |
| Fisher Scientific Sound Enclosure | Fisher Scientific | FB-432-A | For sound reduction |
| Eppendorf Tubes 5.0 ml | Eppendorf | 0030 119.401 | 200 tubes (2 bags of 100 ea.) |
| Random Primer Mix | NEB | S1330S | |
| DNTP MIX | NEB | N0447S | |
| M-MuLV Reverse Transcriptase | NEB | M0253S | |
| FastStart SYBR Green Master | ROCHE | 4673484001 | |
| Agarose | Fisher | BP160100 | 1.5% agarose gel |
| Qubit RNA HS Assay Kit | Life technologies | Q32852 | RNA quantification (100 assays) |
| Formaldehyde | Sigma | F8775 | |
| Protease Inhibitor tablets | Roche | 4693116001 | Each tablet is sufficient for 10 ml lysis buffer |
| PhosSTOP Phosphatase Inhibitor Cocktail Tablets | Roche | 4906845001 | Each tablet is sufficient for 10 ml lysis buffer |
| Qubit Assay Kits | Life technologies | Q32850 | For 100 assays |
| Automated cell counter | ORFLO | MXZ000 | |
| tube Rotator | VWR | 10136-084 | |
| Victor X5 Multilple plate reader | PerkinElmer | 2030-0050 | |
| Microcentrifuge 5417R | Eppendorf | 22621807 | |
| Magnetic Stand | Diagenode | B04000001 | |
| Fast Real-time PCR system | Applied Science | 4351405 |
References
- Hamilton, K. J., Arao, Y., Korach, K. S. Estrogen hormone physiology: reproductive findings from estrogen receptor mutant mice. Reprod Biol. 14 (1), 3-8 (2014).
- Madak-Erdogan, Z., et al. Integrative genomics of gene and metabolic regulation by estrogen receptors alpha and beta, and their coregulators. Mol Syst Biol. 9, 676 (2013).
- Gong, P., et al. Transcriptomic analysis identifies gene networks regulated by estrogen receptor alpha (ERalpha) and ERbeta that control distinct effects of different botanical estrogens. Nucl Recept Signal. 12, e001 (2014).
- Bergamaschi, A., et al. The forkhead transcription factor FOXM1 promotes endocrine resistance and invasiveness in estrogen receptor-positive breast cancer by expansion of stem-like cancer cells. Breast Cancer Res. 16 (5), 436 (2014).
- Madak-Erdogan, Z., et al. Nuclear and extranuclear pathway inputs in the regulation of global gene expression by estrogen receptors. Mol Endocrinol. 22 (9), 2116-2127 (2008).
- Madak-Erdogan, Z., Lupien, M., Stossi, F., Brown, M., Katzenellenbogen, B. S. Genomic collaboration of estrogen receptor alpha and extracellular signal-regulated kinase 2 in regulating gene and proliferation programs. Mol Cell Biol. 31, 226-236 (2011).
- Madak-Erdogan, Z., Ventrella, R., Petry, L., Katzenellenbogen, B. S. Novel roles for ERK5 and cofilin as critical mediators linking ERalpha-driven transcription, actin reorganization, and invasiveness in breast cancer. Mol Cancer Res. 12 (5), 714-727 (2014).
- Wang, Z., Gerstein, M., Snyder, M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 10 (1), 57-63 (2009).
- Invitrogen. Quant-iT™ PicoGreen ® dsDNA Reagent and Kits. , Available from: http://www.nature.com/protocolexchange/system/uploads/3551/original/Quant-iT_Picogreen_dsDNA.pdf?1423512049 (2008).
- Lee, P. Y., Costumbrado, J., Hsu, C. Y., Kim, Y. H. Agarose gel electrophoresis for the separation of DNA fragments. J Vis Exp. (62), (2012).
- Quant-iT™ Assays for high-throughput quantitation of DNA, RNA, and oligos. , Available from: https://tools.thermofisher.com/content/sfs/brochures/F-065432%20quantit%20htp_FLR.pdf (2006).
- Kim, D., et al. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 14 (4), R36 (2013).
- Langdon, W. B. Performance of genetic programming optimised Bowtie2 on genome comparison and analytic testing (GCAT) benchmarks. BioData Min. 8, (2015).
- Zhang, Y., et al. Model-based analysis of ChIP-Seq(MACS). Genome Biol. 9 (9), R137 (2008).
- Gao, J., et al. Metscape: a Cytoscape plug-in for visualizing and interpreting metabolomic data in the context of human metabolic networks. Bioinformatics. 26, 971-973 (2010).
- Schroeder, A., et al. The RIN: an RNA integrity number for assigning integrity values to RNA measurements. BMC Mol Biol. 7 (3), (2006).
- Mokry, M., et al. Efficient double fragmentation ChIP-seq provides nucleotide resolution protein-DNA binding profiles. PLoS One. 5 (11), e15092 (2010).
- Park, P. J. ChIP-seq: advantages and challenges of a maturing technology. Nat Rev Genet. 10, 669-680 (2009).
- Landt, S. G., et al. ChIP-seq guidelines and practices of the ENCODE and modENCODE consortia. Genome Res. 22 (9), 1813-1831 (2012).
- Hah, N., Kraus, W. L. Hormone-regulated transcriptomes: lessons learned from estrogen signaling pathways in breast cancer cells. Mol Cell Endocrinol. 382, 652-664 (2014).
