在酵母S.酵母缺氧反应期间, 测量一段时间内的 mRNA 水平

摘要

在这里, 我们提出了一个协议使用 RNA 序列监测 mRNA 水平随着时间的推移在缺氧反应的S. 酿酒酵母细胞。这种方法可以适应于分析基因表达在任何细胞反应。

摘要

基因表达的复杂变化通常会调节细胞反应的很大一部分。每个基因都可能改变表达与独特的动力学, 因为该基因是由特定的时间, 其中一个许多刺激, 信号通路或次要的影响调节。为了在酵母S.酵母中捕获对缺氧的整个基因表达反应, 使用 RNA 序列分析在接触缺氧后的特定时间内监测所有基因的 mRNA 水平。缺氧是由生长在 100% N₂气体中的细胞建立的。重要的是, 不像其他的缺氧研究, 麦角固醇和不饱和脂肪酸不添加到媒体, 因为这些代谢物影响基因表达。时间点选择在 0-4 小时的范围内缺氧, 因为该时期捕获的主要变化的基因表达。在每个时间点, mid-log 缺氧细胞迅速过滤和冻结, 限制暴露于 O₂和伴随变化的基因表达。总 RNA 从细胞被提取并且用于丰富为 mRNA, 然后被转换成 cDNA。从这个 cDNA 中, 创建了多个复合库, 并在下一代排序器的一个泳道中测序了八或更多的样本。描述了 post-sequencing 管道, 包括质量基础修剪、读取映射和确定每个基因的读取数。DESeq2 在 R 统计环境中被用来识别在任何一个缺氧时间点发生显著变化的基因。对三生物复制的分析显示, 高重现性, 不同动力学的基因和大量预期 O₂调控基因。这些方法可用于研究各种生物体的细胞如何随着时间的推移对缺氧反应, 并适应在其他细胞反应过程中研究基因表达。

引言

许多生物体反应缺氧, 或低 O₂, 通过改变基因表达式^1,2,3。这种反应可以帮助细胞应对缺乏对有氧呼吸和一些生物合成反应至关重要的基质, 而且还能改变氧化还原状态⁴。在S. 酿酒酵母中进行的几个微阵列研究表明, 数百个基因的 mRNA 水平因缺氧而发生变化^{5、6、7、8、9,10,11,12}. 最近, RNA seq 被用来表征在缺氧时的基因表达变化¹³。在此, 对实验细节进行了介绍和讨论。

缺氧可以以各种方式实现, 每种方法都产生不同级别的 O₂。在这里, 缺氧是建立在连续流动超高纯度 N₂成烧瓶, 这降低 [O₂] 立即溶解与可再生动力学¹⁰。这是可能的, 有一些 O₂分子存在, 有助于新陈代谢和基因表达, 但这种环境被认为非常接近厌氧。在没有 O₂的情况下, 酵母细胞不能合成血红素、麦角甾醇和不饱和脂肪酸^4、12、14。因此, 以往的研究包括这些代谢物时, 生长酵母没有氧气^5,10,15。然而, 许多缺氧反应是介导的这些代谢物的损耗, 从而补充他们逆转缺氧基因表达反应^12,16。为了模拟自然缺氧, 这些代谢物没有添加到媒体。在短时间内, 细胞暴露在缺氧状态, 没有这些必需的代谢物, 没有明显增加细胞死亡 (数据没有显示), 也没有长期的压力响应¹³。

反应也取决于菌株及其基因型。特别重要的是已知的低氧反应调节器的等位基因²。S288C 应变背景是高度理想的, 以便结果可以与其他基因组研究与此菌株进行比较。然而, S288C 包含了一个部分功能等位基因的HAP1吉恩¹⁷, 转录调节器对缺氧反应至关重要。此等位基因在 S288C 中修复, 使用Σ1278b 应变背景下的野生拷贝¹¹。

基因表达高度依赖于细胞环境。因此, 在进行全基因组的 mRNA 分析时, 保持一个恒定的环境是很重要的, 同时改变另一个参数, 例如时间、刺激或基因型。为了获得高度重现性的结果, 请考虑这三项研究的实践及其所有的生物或技术复制。首先, 同样的实验者应该进行这项研究, 因为技术实践可能在实验者中有所不同。第二, 在培养基中应使用相同批次的配料, 因为每个批次都有一个稍有不同的成分, 可以影响基因表达。第三, 为了最小化细胞周期效应, 每个时间点应该由生长 mid-log 阶段的异步细胞组成 (1-2 x 10⁷细胞/mL)。

当表征一个复杂的反应, 如基因表达反应的缺氧, 一个时间的过程是有利的, 以确定动力学的各种事件。应选择特定的时间点来捕获响应的主要变化。在这项研究中, 观察到0和 4 h 之间的时间点, 因为过去的实验揭示了在这一期间的基因表达的广泛变化¹³。

为了测量全球基因表达, RNA seq 使用了^18,19。这个方法使用下一代测序来确定每个基因转录的相对丰度。与 DNA 微阵列分析相比, RNA 序列表现出更高的灵敏度 (以检测更少的抄本), 更大的动态范围 (测量更大的褶皱变化) 和优异的重现性 (准确地跟随基因表达随着时间的推移)。通常, 大多数细胞 rna 是核糖体 rna, 因此许多方法已经开发, 以丰富的特定 RNA 种类²⁰。在这里, 聚 T 珠被用来纯化多含 A 的 mrna 转录, 虽然各种商业可用的 rRNA 耗竭套件也可以有效的 mrna 丰富。

在这里, 对缺氧的S 酿酒酵母基因表达反应进行了表征。细胞暴露在缺氧, 然后取样在八时间点 (0, 5, 10, 30, 60, 120, 180 和240分钟)。为了确认重复性和识别统计变化的成绩单, 三进行了生物复制。rna 是由机械破坏和柱纯化提取, 然后处理 rna 序列分析。描述了 post-sequencing 管道, 并提供了允许精确复制所执行的分析的编程脚本。具体来说, Trimmomatic ²¹、TopHat2 ²²、HTseq ²³、R 统计环境²⁴和 DESeq2 包²⁵用于处理 RNA 序列数据, 并识别607基因在缺氧.主成分分析 (PCA) 和基因表达的复制表明该技术的重现性。聚类和 heatmaps 揭示了 wide-ranging 表达动力学, 而基因本体 (GO) 分析表明, 许多细胞过程, 如有氧呼吸, 是丰富的一套氧调节基因。

研究方案

1. 诱导缺氧

1. 前一天或更早的缺氧时间课程: 准备孵化器, 细胞过滤系统, 真空, 气罐, 烧瓶, 塞子, 玻璃管, 和油管, 如在材料表。
2. 将 N 个₂容器、孵化器、真空和过滤系统放在接近的位置, 以便能够快速处理电池。
3. 通过将玻璃瓶和灭菌中的成分混合, 制备无菌液体 YPD 培养基 (1% 酵母提取物, 2% 蛋白胨, 2% 葡萄糖)。
4. 规划孵化器的烧瓶布局。
   注: 从 N₂罐, 第一个烧瓶将是一个水陷阱, 第二个烧瓶将是最后的时间点, 第三个烧瓶将是第二的时间点等, 直到最后一个烧瓶是最后的水陷阱。图 1显示了烧瓶之间的顺序和连接。
5. 在时间的前一天, 接种一个酵母菌落到5毫升液体 YPD 在一个无菌试管内的文化。具体地说, 用无菌涂抹棒从盘子里拿出一个完整的菌落。把棍子放入液体 YPD, 摇动棍子直到大部分的细胞悬浮。
6. 旋转管在30° c 过夜 (〜16小时)。
   注:16 小时后, 野生细胞将达到饱和 (〜 2 x 10⁸细胞/毫升), 由一个非常多云的文化表示。其他菌株可能需要更长的时间达到饱和, 并应通过测量细胞浓度的分光光度计, 如步骤1.9.4 所示。
7. 在时间路线的天, 在 16 h 以后孵育, 稀释饱和的文化1:50 通过使用5毫升吸管安置4毫升隔夜文化入196毫升液体 YPD 在无菌500毫升烧瓶之内。在30° c 的4小时在 200 rpm 的震动下成长。
   注意: 4 小时后, 野生的区域性将达到 mid-log 1-2 x 10⁷单元/mL 的浓度。
8. 在4小时孵化, 标签的烧瓶 (缺氧和水陷阱) 和50毫升收集管。如果以上步骤1.7 中的孵化器不用于低氧时间课程, 则打开其他孵化器并设置为30° c。
9. 就在细胞缺氧之前:
   1. 添加〜50毫升水到两个无菌250毫升烧瓶, 将作为水的陷阱。
   2. 填充带有液氮的冰桶的¹/₃ , 并淹没聚苯乙烯泡沫机架以容纳50毫升离心管。
   3. 设置用于收集单元格的筛选系统。使用无菌镊子, 在过滤系统的底部单元上添加无菌滤片, 小心只触及过滤盘的边缘。将过滤器顶部的单位放到过滤器底部, 并与系统提供的夹具一起安全。确保底部和顶部的单位是对齐, 以便有一个紧密的密封和没有泄漏。
   4. 用分光光度计测量细胞浓度。稀释细胞, 使它们可以在分光光度计的线性范围内进行测量--OD₆₀₀介于 0.2-0.6 之间, 这取决于光谱仪。
      注: 一个 OD₆₀₀单位是〜 3 x 10⁷细胞/毫升, 根据细胞形态学和分光光度计。所需浓度为 ~ 1-2 x 10⁷单元/mL, 对应于 ~ 0.33-0.66 的 OD₆₀₀ 。
   5. 快速获取时间0样本。
      1. 通过真空管和1000毫升的烧瓶陷阱将过滤系统连接到一个强 (约10毫巴) 真空。打开真空。将养殖 (通常为20毫升) 倒入顶部, 等待液体通过过滤器 (〜十年代, 取决于真空强度)。
      2. 小心地用干净的镊子将滤片取出, 放入50毫升离心管。立即将离心管放入液体氮气的机架中。重要的是, 在插入过滤器之前不要 pre-chill 管, 因为管子会爆炸。
      3. 三十年代 > 之后, 将管子放入一个-80 ° c 的冰箱中, 以备以后的 RNA 提取。每次过滤后, 清洗并重新组装过滤系统。冲洗系统, 通过拉水通过十年代没有一个过滤器光盘。最后, 用纸巾擦干系统。
   6. 将 4 h 培养液稀释为不同的时间点, 使每个烧瓶达到 mid-log 浓度 (〜 1-2 x 10⁷细胞/毫升) 在指定的时间点缺氧。
      注: 表1显示了用于野生类型 S288C HAP1⁺单倍体菌株的稀释。建议对这些缺氧培养条件下的每株菌株进行测试, 并相应调整稀释。
   7. 一旦细胞和培养基被添加到烧瓶中, 用装有两个玻璃管的塞子取代覆盖烧瓶开口的铝箔。
   8. 将烧瓶放在孵化器中, 并将其放置在先前确定的布局中。
   9. 安全地连接油管, 如图 1所示。
   10. 检查所有的塞子都被紧紧地推入烧瓶开口。
10. 0时, 打开调节阀。然后将流量计设置为3升/分。
    注: 气泡将在两个水瓶中观察到, 表明适当的气体流动。如果不是这样的话, 检查所有的塞子都是紧密的, 油管是正确的连接。
11. 关闭孵化器和设置震动速度〜 200 rpm。启动计时器。
12. 在每个时间点之前, 按照上面的步骤1.9.3 设置过滤系统。
13. 在每个时间点 (例如, 5 分钟, 10 min,等) 中, 有两个人快速处理这些单元格, 如下所示:
    1. 第一人: 从持有人身上取出适当的烧瓶, 取出塞子 (附玻璃管)。
       注意: 这不会将 N₂的流分解为任何剩余的区域性烧瓶, 但会暂时中断流向最后一个水阱的流。
    2. 第二人称: 吸管1毫升的文化, 并免除试管。重新连接现在位于生产线末端的养殖烧瓶中的油管, 将其与最后一个水阱 (一个管和一个塞子将在过程中删除)。然后, 测量试管中的细胞浓度, 如上面的步骤1.9.4 所述。
    3. 第一人称: 将养殖的其余部分倒入真空过滤系统, 然后按照上述步骤1.9.5 进行过滤和冻结。
14. 完成所有时间点后, 关闭 N₂气体。首先关闭稳压器, 等待压力释放, 然后关闭流量计。最后, 拆卸系统, 用清水清洗所有材料, 然后70% 乙醇。

2. RNA 提取

1. 为 RNA 列纯化准备 RLT 缓冲区, 如材料表中所述。
2. 通过将 dnasei 库存解决方案的10µL 添加到每个示例的70µL 缓冲区 RDD, 准备一个 dnasei 的工作解决方案 (请参见材料表)。
3. 从-80 ° c 的冰箱中取出含有过滤器和电池的50毫升管子, 并在冰上放上解冻 (约15分钟)。用冰上的管子执行以下步骤。
4. 标签2毫升螺丝帽管, 并把它们放在冰, 一个管为每个样品。加入〜0.6 毫升酸洗珠每管, 使用1.5 毫升离心管舀和测量的珠子。
5. 在50毫升的管子上加入0.6 毫升的冷 RLT 缓冲器。
6. 吸管向上和向下移除过滤器中的单元格并将其挂起到解决方案中。避免让过滤器留在解决方案中, 因为过滤器会吸收解决方案。通过使用吸管尖端将过滤器挤压到管壁上, 将所有吸收到过滤器中的溶液移除。
7. 将所有的液体从50毫升的管子转移到含有珠子的螺钉盖管。
8. 将螺钉盖管放入一个珠磨均质机中, 并运行一分钟。
   立即把管子放在冰上三分钟。再次, 将管放入均质机, 并运行一分钟。
9. 从均质机上取出管子, 在冰上放置5分钟。珠子会沉淀到管子的底部。
10. 在室温下执行其余步骤。
11. 与吸管, 转移只有裂解 (〜350µL), 避免珠子, 一个新的1.5 毫升离心管。
12. 离心为2分钟的最大速度, 然后转移上清到一个新的1.5 毫升离心管。
13. 添加1量的70% 乙醇 (由200证明分子生物学级乙醇)。混合移。
14. 将溶液和任何沉淀转移到已放置在2毫升收集管内的 RNA 柱上。
15. 离心机为十五年代在≥12,000 x g并且摒弃流动通过 (液体在管子)。
16. 向列中添加350µL 的缓冲区 RW1。离心机十五年代在≥12,000 x g 和放弃流动通过。
17. 添加80µL dnasei I 孵化混合到柱状膜 (不得到两侧的管), 并允许坐15分钟。
18. 向列中添加350µL 的缓冲区 RW1。离心十五年代在≥12,000 x g , 并丢弃流通过。
19. 添加500µL 的缓冲 RPE 到该列。离心十五年代在≥12,000 x g并丢弃流。
20. 添加500µL 的缓冲 RPE 到该列。离心为2分钟, 在≥12,000 x g。
21. 小心地删除的列和放置到一个新的2毫升收集管。离心的速度为1分钟 (完全干燥的膜)。
22. 丢弃收集管, 并将该列放入1.5 毫升离心管中。在柱膜上加入30µL 的无核糖核酸水。
23. 洗在≥12,000 x g的 microcentrifuging 中, 从该列中提取 RNA 1 分钟。当把柱子装入离心时, 要确保收集管的盖子朝向离心机旋转的方向, 以防止盖子折断。
24. 在柱状膜上加入另30µL 的无核糖核酸水, 同时保持柱在同一离心管中。
25. 离心为1分钟, 在≥12,000 x g, 使最终洗卷为60µL。
26. 在-80 ° c 的冷冻库中储存每1.5 毫升含有60µL 的 RNA 的管子。

3. 测定 RNA 浓度和质量

1. 测量 rna 和 dna 在每个 rna 样品中的浓度, 使用 (a) 荧光, (b) dna 和 rna 特异荧光染料和 (c) 化验管, 如在材料表中所列。按照荧光和染料的指示进行操作。
2. 另外, 用紫外分光光度计测量核酸浓度 260 nm。
   注: A₂₆₀读数为 1.0, 相当于 ~ 40 µg/毫升单 RNA。
3. 通过在商业核酸分析仪上运行 rna 来测试 rna 质量 (参见材料表) 或标准甲醛/琼脂糖凝胶²⁶。

4. RNA-seq 分析

1. 从总 rna 的股票, 构成了21µL 100 ng/µL RNA 在核糖核酸自由水。使用1µL 这21µL 测试 RNA 浓度如上所述。将剩余的20µL (2 µg) 总 RNA 提交到外部测序中心, 用于 mRNA 富集, 特定于链路的库准备 (八或更多的条形码用于复用) 和排序 (请参见材料表)。或者, 在本地执行 mRNA 浓缩和库准备, 然后提交库进行排序。
2. 池八或更多的多路复用样本和序列在一个通道的下一代排序器。
3. 下载包含所有序列读取和相应索引读取的 FASTQ 文件。另外, 下载包含多重条码的文本文件。
   注意: 索引读取可能存在于单独的 FASTQ 文件中。将所有这些文件放在一个目录中。
4. 将 shell 脚本放在这里, "fastq_pipeline", 在同一目录中。根据实验和计算机目录编辑此脚本。
   注意: 此脚本包含解释步骤的大量注释。简而言之, 脚本质量修剪读取, 将读取映射到基因组, 并生成一个制表符分隔的文件, 其中包含每个样本的每个基因的读取次数。
5. 在发布²⁷之前, 将 FASTQ 文件和原始读取计数数据存入 NCBI 的基因表达式综合中。按照 "向 GEO 提交高通量序列数据" 的说明: https://www.ncbi.nlm.nih.gov/geo/info/seq.html。
6. 将读取计数数据导入 r 统计环境, 并使用此处提供的 r 脚本进行统计分析, "time_course_script。R "。根据实验和计算机目录编辑此脚本。
   注意: 此脚本包含解释步骤的大量注释。简而言之, 该脚本导入读取数据, 规范化读取, 识别在时间过程中显著变化的基因, 执行 PCA 分析, 并对选定基因的表达进行图表。

结果

三次独立进行缺氧时间过程和 RNA 序列分析。为了检验三复制的重现性, 利用主成分分析 (PCA) 对所有基因的基因表达数据进行了分析。图 2显示了示例在前两个主要组件上的变化情况, 它们一起表示58.9% 的可变性。这一分析表明, 每一次课程都有类似的变化 (如每条曲线的相似形状所示)。另外, 最后两个复制比第一个复制更相似。这与最后两次复制是由...

讨论

在本研究中, 对所有基因的 mRNA 水平进行了测试, 在低氧的酵母S酵母。目的是分析全球基因表达在受控近缺氧环境中的生长变化。采取了若干步骤, 以确保此处所述的方法是经过严格控制和重现的。首先, 细胞暴露在一个精确定义的缺氧环境中: 富媒体 (YPD) 中的 99.999% N₂ 。其他缺氧的研究已经关闭了烧瓶或管空气³², 使用缺氧模拟氯化钴³³, 或雇用缺...

材料

Name	Company	Catalog Number	Comments
Enclosed dry incubator	Thermo Scientific	MaxQ 4000	Set at 30 °C. Modify the door to allow entry of one Tygon tube. Alternatively, use the New Brunswick G25 incubator, which contains a tube port. Do not use an open-air water shaker, as condensation will collect in the tubes between flasks, possibly cross-contaminating cultures.
Micro-analysis Filter Holder	Millipore	XX1002530	25mm diameter, stainless steel support, no. 5 perforated silicone stopper mounts in standard 125 mL filtering flask
Strong vacuum	Edwards	E-LAB 2	The “house” vacuum may be too weak. Alternatively, use an electric-power portable vacuum pump like the one listed here.
1,000 mL flask			To act as vacuum trap.
~2 foot lengths of Heavy Wall Vacuum Tubing, inner diameter 3/8 in, outer diameter 7/8 in	Tygon	38TT	Two pieces: the first connects vacuum to trap, and the second connects trap to filter system.
High-pressure N2 gas tank			99.999% purity, >1,000 psi, with a regulator and gas flow controller
Autoclaved 500 mL flask			Opening covered with aluminum foil. One for each yeast strain.
Autoclaved 250 mL flasks			Openings covered with aluminum foil. One for each time point plus two for water traps.
Flask stoppers (size 6, two holes with 5 mm diameter)			Sterilized with 70% ethanol. One for each flask.
Glass tubing, length 9 cm or 17 cm, inner diameter 2 mm, outer diameter 5 mm			Sterilized with 70% ethanol. Two tubes for each flask. Place into the holes of each stopper. See Figure 1 for placement of 9- vs 17- cm tubes.
~25 cm lengths of plastic tubing, inner diameter 5 mm	Tygon	E-3603	One piece for each flask. Sterilized with 70% ethanol.
Sterile filter discs	Millipore	HAWP02500	25 mm diameter, 0.45 µm pore size, one for each time point
Sterile dH2O (~100 mL)
1 mL cuvettes			For measuring OD600 (i.e., cell concentration)
50 mL sterile centrifuge tubes			One for each time point
Clean and sterile tweezers
liquid nitrogen			For freezing cells
acid-washed beads	Sigma	G8772	Keep at 4 °C for lysing cells
Qiagen RNeasy Mini Kit	Qiagen	74104	For RNA column purification
Qiagen RLT buffer			Prepare by adding 10 µL of β-Mercaptoethanol per 1 mL of RLT buffer, keep at 4 °C.
2 mL collection tubes	Qiagen		included in the Qiagen Rneasy Mini Kit
Buffer RPE	Qiagen		included in the Qiagen Rneasy Mini Kit
Buffer RW1	Qiagen		included in the Qiagen Rneasy Mini Kit
DNase I stock and working solutions	Qiagen	79254	The DNase I enzyme comes as lyophilized powder in a glass vial. Using a sterile needle and syringe, inject 550 µL of RNase-free water (provided in Qiagen kit) into the vial. Mix by gently inverting the bottle. To avoid denaturing the enzyme, do not vortex. Using a pipet, remove this stock solution from the vial and store in freezer (-20 °C) in single-use aliquots (80 µL each). The stock solution should not be thawed and refrozen.
Buffer RDD	Qiagen		included in the Qiagen DNase Kit
Ice cold 2 mL screw-cap tubes			For lysing cells during RNA extraction
bead mill homogenizer	Biospec Mini-Beadbeater-24	112011	Keep in cold room
Bacto Peptone	BD	DF0118	for liquid YPD media
Bacto Yeast Extract	BD	DF0886	for liquid YPD media
glucose	Fisher	D16	for liquid YPD media
Qubit assay tubes	Thermo Fisher	Q32856	for measuring nucleic acid concentration
Quant-iTTM dsDNA BR Assay Kit	Thermo Fisher	Q32853	for measuring nucleic acid concentration
Quant-iTTM RNA Assay Kit	Thermo Fisher	Q32855	for measuring nucleic acid concentration
Qubit Fluorometer	Thermo Fisher	Q33216	for measuring nucleic acid concentration
Commercial electrophoresis system	Agilent	Bioanalyzer 2100	for measuring nucleic acid quality
Next-generation sequencer	Illumina	HiSeq 2500	for sequencing libraries
automated liquid handling system	Wafergen	Apollo 324	for creating sequencing libraries
PrepX PolyA mRNA Isolation Kit	Wafergen	400047	for isolating mRNA from total RNA
PrepX RNA SEQ for Illumina Library Kit	Wafergen	400039	for creating strand-specific sequencing libraries from total RNA
Barcode Splitter			https://toolshed.g2.bx.psu.edu/repository?repository_id=7119c4f7a89efa57&changeset_revision=e7b7cdc1834d
Samtools, which includes the gzip command			http://www.htslib.org/download/
Trimmomatic			http://www.usadellab.org/cms/?page=trimmomatic
Bowtie2 (installed before TopHat)			http://bowtie-bio.sourceforge.net/bowtie2/index.shtml
TopHat			https://ccb.jhu.edu/software/tophat/index.shtml
HTSeq			http://www-huber.embl.de/HTSeq/doc/overview.html
R (installed before R Studio)			https://cran.rstudio.com
R Studio (free version)			https://www.rstudio.com/products/rstudio/download/