Name: 斑马鱼 RNASeq-based 基因表达数据的样品制备与分析 (视频) | JoVE Journal
Uploaded: 2017-10-27T12:30:00.000+00:00
Duration: 11 min 42 s
Description: 学习斑马鱼RNASeq转录组分析方法，用于识别胚胎和幼鱼基因表达谱并进行定量比较。该技术结合RNA提取、通路分析和qRT-PCR验证，为发育生物学和神经科学研究提供可靠框架。

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

本工作由 R01DK102001（N.A.Z.）、P30DK072488（N.A.Z.）和 T32DK098107（T.L.H. 和 J.E.N.）资助。

以下所述的所有动物实验方案均符合马里兰大学机构动物护理与使用委员会（IACUC）的规定，并已获得其批准。
1. 胚胎制备
<ol>
	<li>通过自然交配产生胚胎
	<ol>
		<li>将胚胎培养至3月龄，达到性成熟5,8。</li>
		<li>在收集胚胎的前一晚，将所需品系的成年雄鱼和雌鱼分开放入分隔式交配缸中，每缸放入2尾雄鱼和3尾雌鱼。 
		注：使用转基因insulin2a:mCherry荧光报告品系可实现对胰腺β细胞的分析。</li>
		<li>次日清晨开灯后，将鱼转移至含有新鲜系统水的交配缸中，并立即移除分隔板。</li>
		<li>让鱼自然交配，直至在缸底观察到胚胎。每30分钟收集一次胚胎，直至获得所需数量。将每次收集的胚胎按时间点分别置于不同的培养皿中，使用胚胎培养基在28.5 °C下保存。</li>
		<li>如需进行遗传物质显微注射或转入实验性培养基6，可在此时操作，并将胚胎置于10 cm培养皿中的新鲜Hank's胚胎培养基8中，于28.5 °C下继续培养。 
		注：对于注射吗啉代（MO）或突变体动物的基因表达分析，需注意每种操作可能对基因表达产生非特异性影响。突变体在转录水平可能发生遗传补偿，而基于MO的基因靶向则通常不出现此类现象9。</li>
	</ol>
	</li>
	<li>胚胎分期
	<ol>
		<li>将胚胎以每10 cm培养皿50–75枚胚胎的密度分组培养，以保证所有胚胎发育时序一致。</li>
		<li>在卵裂期、外包期和体节期使用解剖显微镜监测胚胎形态发育，确保正常发育进程10。 
		注：及时移除死亡或畸形胚胎，以防止对同皿其他胚胎造成发育延迟。</li>
		<li>根据胚胎的发育年龄进行分离。在原肠胚形成后（受精后10.33小时，hpf）至约24 hpf期间，通过体节数量测定胚胎年龄；24 hpf之后，依据总体长对胚胎和幼鱼进行分期。 
		注：体节是位于胚胎背部的V形中胚层组织。</li>
		<li>将分选后的胚胎置于28.5 °C培养箱中，继续发育至目标阶段。</li>
	</ol>
	</li>
</ol>
2. 单细胞解离：完整胚胎和分选细胞群体
<ol>
	<li>整个胚胎解离
	<ol>
		<li>将斑马鱼胚胎置于冰上5-15分钟，直至观察不到任何运动，以在所需阶段处死胚胎。</li>
		<li>将20枚胚胎转移至一个已标记的1.5 mL微量离心管中。移除微量离心管中多余的胚胎培养液。</li>
		<li>添加 200 µL 裂解试剂（参见 材料清单) 加入含有胚胎的离心管中。使用研杵在裂解试剂中进行机械匀浆以裂解胚胎。添加 800 µL 裂解试剂，使总体积达到1 mL。储存于-80 °C 直至 RNA 提取。</li>
	</ol>
	</li>
	<li>流式分选细胞分离11
	<ol>
		<li>将胚胎从培养皿转移至1.5 mL微量离心管中，并尽可能吸除胚胎培养液。 
		注意：根据胚胎的发育时期，此步骤可能还需要进行机械解离。对于胚胎 &gt; 48 hpf时，建议使用研杵破坏胚胎完整性后再进行后续操作。</li>
		<li>将胚胎与1 mL解离缓冲液1孵育（参见 材料清单在1.5 mL微量离心管中。室温孵育15–30分钟。孵育期间每3–4分钟用P1000枪头轻轻吹打溶液上下混匀。切勿涡旋。</li>
		<li>通过离心收集胚胎，300 × g 离心 3 分钟，弃去上清液。将胚胎重悬于 1 mL 解离缓冲液 2 中（见 材料清单）。在室温下孵育15–30分钟，期间定期用移液器吹打混匀。</li>
		<li>通过稀释1-评估消化情况2 µL 在显微镜下，用荧光激活细胞分选（FACS）缓冲液重悬细胞。 
		注意：当检测的等分试样显示为单个细胞且极少细胞团及胚胎组织碎片时，表明消化完全。</li>
		<li>以 300 × g 离心 5 分钟收集细胞。移除上清液。将细胞重悬于 1 mL FACS 缓冲液中，并通过重力过滤装置过滤 40 µm 使用细胞滤筛去除样品中未消化的组织。</li>
		<li>使用血细胞计数板对细胞进行计数，并用FACS缓冲液稀释至约1x10⁶个/mL6 细胞/mL，置于FACS管中。 
		注意：对于每胚胎中数量相对较少的细胞类型（占总数不到5%），例如胰腺β细胞，应使用至少1,000个发育阶段匹配的胚胎。</li>
		<li>向FACS分选设施提供含有1 mL细胞悬液样本的FACS管，以及仅含裂解试剂或FACS缓冲液的FACS管。设置FACS流式分选仪的分选门，仅收集表达荧光报告蛋白的单个细胞。</li>
		<li>进行FACS流式分选，直至收集到所需数量的细胞。 
		注意：进行RNA测序需要一定最低量的总RNA。我们发现，3,000–5,000个细胞足以分离出高质量、高浓度的RNA。</li>
		<li>将收集的细胞置于冰上保存，并立即进行RNA制备。</li>
	</ol>
	</li>
</ol>
3. RNA 制备
<ol>
	<li>提取 RNA
	<ol>
		<li>将收集的样品（来自步骤 2）与裂解试剂在室温下孵育 5 分钟。</li>
		<li>每使用 1.0 mL 裂解试剂，加入 0.2 mL 氯仿，用手倒置混匀试管 15 秒。在室温下孵育 2–3 分钟。在 4 °C 条件下以 12,000 ×g 离心 15 分钟。</li>
		<li>将分离后的水相转移至新离心管中，每使用 1.0 mL 裂解试剂加入 0.5 mL 异丙醇。在室温下孵育 10 分钟。在 4 °C 条件下以 12,000 ×g 离心 10 分钟。</li>
		<li>用 75% 乙醇洗涤，在 4 °C 条件下以 7,500 ×g 离心 5 分钟。彻底去除上清液，并在室温下晾干。用 15–30 µL 经焦碳酸二乙酯（DEPC）处理的水重悬 RNA。</li>
	</ol>
	</li>
	<li>纯化高质量 RNA
	<ol>
		<li>将 RNA 悬液与 1/10 体积的 3M 醋酸钠（pH 5.5）和 1 体积的异丙醇混合。在室温下孵育 20 分钟，在 4 °C 条件下以 12,500 ×g 离心 10 分钟。</li>
		<li>用冰预冷的 70% 乙醇（DEPC 处理水配制）洗涤沉淀，在 4 °C 条件下以 10,000 ×g 离心 5 分钟。重复洗涤一次。</li>
		<li>去除上清液，在室温下晾干。用 DEPC 处理水重悬 RNA。</li>
		<li>使用吸收光谱仪测定提取 RNA 的浓度（ng/µL）和纯度（260/230 比值）。在进行 RNA 测序前，确保 260/230 比值约为 2.0。储存于 -80 °C 直至使用（最长 6 个月）。</li>
		<li>将 RNA 样品送至供应商或核心实验室进行 RNA 测序及基于测序读数定量的基因表达变化分析。供应商返回的结果通常以在样品间转录本读数发生显著倍数变化的基因列表形式提供。 
		注：若上述方法获得的 RNA 质量较差，可使用层析柱进行纯化。RNA 样品的用量、浓度和 RNA 完整性数（RIN）应与供应商或核心实验室确认。供应商或核心实验室通常也会评估 RIN。</li>
	</ol>
	</li>
</ol>
4. 通路与 GO 术语分析
注意：有关核心机构或供应商提供的 RNASeq 基因表达分析的代表性结果，请参见图 1。
<ol>
	<li>筛选差异表达基因
	<ol>
		<li>比较单一实验条件与对照组
		<ol>
			<li>在电子表格管理软件中打开数据（结果）（参见 材料表）。</li>
			<li>在包含实验组与对照组之间差异表达基因的电子表格中，点击“排序”按钮旁边的下拉箭头，并选择“自定义排序”。</li>
			<li>在弹出窗口的“列”下方选择框，并选择按“LFC”列进行排序。确保“排序依据”列中选择的是“数值”，在“次序”列中选择“从大到小”，然后点击“确定”。 
			注意：此操作将对差异表达基因进行排序，使表达上调（正LFC）的基因位于列表顶部，表达下调（负LFC）的基因位于列表底部。</li>
		</ol>
		</li>
		<li>按以下步骤将多个实验条件与单一对照组进行比较。
		<ol>
			<li>在“实验1 vs 对照”工作表的空白列中选择第一个单元格，以确定在两个实验条件下均出现的差异表达基因。在该单元格中输入以下公式： 
			<img alt="用于重复检测的Excel公式，方程公式，教育用途。" src="/files/ftp_upload/56187/56187eq1.jpg" data-altupdatedat="1788250015000" data-alt="方程"> 
			注意：该公式结合了IF、ISERROR和MATCH函数。（i）MATCH函数MATCH(A1,Experimental2_vs_Control!A:A,0)将在名为Experimental2_vs_Control的工作表的A列（A:A）中查找A1单元格中的值（ENSEMBL ID）。“0”表示查找完全匹配项，若找到完全匹配项，则函数返回值“1”；若未找到匹配项，则返回错误“N/A”。（ii）MATCH函数的结果将作为输入传递给ISERROR函数ISERROR(MATCH(A1,Experimental2_vs_Control!A:A,0))，若输入为“1”（表示找到匹配项），则返回“FALSE”；若输入为“N/A”（表示未找到匹配项），则返回“TRUE”。（iii）“TRUE”或“FALSE”结果再输入到IF函数IF(ISERROR(MATCH(A1,Experimental2_vs_Control!A:A,0)),"","Duplicate")中，若输入为“TRUE”（表示未找到匹配项），则返回第一对引号之间的值（即空值），因此该单元格为空；若输入为“FALSE”（表示找到匹配项），则返回第二对引号之间的值（即“Duplicate”），该单元格将显示“Duplicate”。</li>
			<li>按“Enter”或“Return”运行公式。选中包含公式的单元格，点击该单元格右下角的小方框，按住鼠标并向下拖动至最后一个“Feature ID”，以将公式复制到该列的每个单元格中。</li>
			<li>再次选择“自定义排序”，通过点击左下角的“+”图标添加第二级排序。在第一级“排序依据”中，“列”选择“Duplicate”，“排序依据”选择“数值”，“次序”选择“Z到A”；在第二级“然后依据”中，“列”选择“LFC”，“排序依据”选择“数值”，“次序”选择“从大到小”，然后点击“确定”。 
			注意：此操作将基因列表根据表达变化的方向性分为四组。检查在两个实验条件下均出现的基因非常重要，以确认其表达变化方向是相同还是相反。</li>
			<li>在继续操作前，从基因符号列中删除所有括号。选中包含基因符号的整列，选择“编辑”菜单中的“替换”功能。在“查找内容：”栏中输入“(*)”，“替换为：”栏留空，点击“全部替换”以删除所有括号实例。 
			注意：“*”代表任意数量的字符，将此字符置于两个括号之间可提示程序查找高亮区域中所有括号实例，并将其替换为空，从而实现删除。</li>
		</ol>
		</li>
	</ol>
	</li>
	<li>确定富集通路
	<ol>
		<li>将希望确定富集通路的基因集合的基因符号复制到剪贴板。访问ConsensusPathDB12，在网页左侧边栏选择“Gene set analysis”，随后选择“Over-representation analysis”。</li>
		<li>在“Paste a list of gene/protein identifiers”框中粘贴基因列表。在“Gene/protein identifier type”框中选择“gene symbol”，然后点击“Proceed”。在“Pathway-based sets”部分下，勾选“Pathways as defined by pathway databases”旁边的复选框。 
		注意：将出现多个分析选项，包括可供搜索的数据库列表。我们建议仅保留KEGG和Reactome数据库，取消选择其他所有数据库，因为这两个数据库最为成熟且全面。最小重叠基因数和p值截断设置也可根据需要调整。我们建议将这些设置保留为默认值：最小重叠基因数为2，p值截断为0.01。</li>
		<li>点击“Find enriched sets”以获取包含输入基因列表中基因的通路列表。 
		注意：输出结果将以表格形式呈现，包括每条通路的名称、“Set size”（表示通路中总基因数）、“Candidates contained”（表示输入列表中属于该通路的基因数）、p值和q值（FDR），以及识别该通路的数据库。</li>
	</ol>
	</li>
	<li>通路网络的构建
	<ol>
		<li>按照上述方法确定富集通路。</li>
		<li>通过勾选各通路名称旁的复选框，或点击选择框列标题上方“Select”下的“All”，选择所有希望在通路网络中可视化的富集通路，然后点击“Visualize selected sets”。 
		注意：如果通路数量少于30条，建议可视化所有通路；如果超过30条，建议选择富集程度最高的前30条通路（即包含输入列表中基因数最多的通路）。</li>
		<li>通过点击页面顶部中央的“relative overlap”或“shared candidates”筛选框，并输入所需的相对重叠百分比或共享候选基因数，以提高筛选严格性，使通路网络中的关键重叠关系更加清晰，然后点击“apply”。 
		注意：我们建议最小相对重叠值设为0.2，即两条通路之间需有20%的基因重叠才能连接，且共享候选基因至少为2个。可通过点击页面左上角的图例查看图形说明。</li>
	</ol>
	</li>
	<li>富集GO术语的确定
	<ol>
		<li>将希望确定富集基因本体（GO）术语的基因组的基因符号复制到剪贴板。访问基因本体联盟（Gene Ontology Consortium）的“GO富集分析工具”13。在页面左侧“Your gene IDs here…”框中粘贴基因符号列表。</li>
		<li>在基因ID框下方选择要使用的GO术语集合：生物过程（biological process）、分子功能（molecular function）或细胞组分（cellular component）。建议选择“生物过程”。在GO术语框下方选择'Danio rerio'，然后点击“Submit”。 
		注意：我们建议使用默认的p值截断值0.05。</li>
	</ol>
	</li>
</ol>
5. 通过qRT-PCR验证
注意：在RNA测序中发现具有显著基因表达变化的单个基因，应通过靶向qRT-PCR在重复实验中进行验证。
<ol>
	<li>cDNA合成
	<ol>
		<li>使用第3.1节“RNA提取”中所述的方法从胚胎中提取RNA。</li>
		<li>转化 1 µg 使用cDNA转化试剂盒将RNA转为cDNA（参见 材料清单). 将RNA、dNTPs和逆转录酶混合。根据制造商说明书进行PCR仪反应。</li>
		<li>将cDNA用经DEPC处理的水按1:3比例稀释。储存于 4 °C 长达 1-2 个月。</li>
	</ol>
	</li>
	<li>qRT-PCR 验证
	<ol>
		<li>从RNA测序结果中选择用于qRT-PCR验证的基因。筛选表达倍数变化较高的基因，并确定其中哪些基因同时具有较高的读段数，因为这表明其表达丰度较高。</li>
		<li>从高倍数变化、高读数的基因列表中选取10-12个目标基因。设计引物以扩增至少跨越1个内含子-外显子边界的靶基因。</li>
		<li>输入基因或基因的靶向区域至qPCR引物设计网站14. 选择“设计 2 对 qPCR 引物”，并提示网站生成引物组合。 
		注意：务必包含内参引物，例如 β-actin，以标准化样本间的比较。β-actin 引物为 F: 5'-TCGAGCTGTCTTCCCATCCA-3' 和 R: 5'-TCACCAACGTAGCTGTCTTTCTG-3'。</li>
		<li>将 cDNA、上游引物、下游引物和聚合酶混合。进行 qRT-PCR 扩增以测定基因的相对表达水平15.</li>
		<li>通过相对定量分析目标基因的相对mRNA表达水平15.</li>
		<li>将qRT-PCR结果与RNA测序结果进行比较，以确保差异表达的方向性一致。</li>
	</ol>
	</li>
</ol>

<ol><li>Nusslein-Volhard, C., Dahm, R. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aNusslein-Volhard%2C+C+%22Zebrafish%22">Zebrafish</a>. , Oxford University Press. (2002).</li>
<li>Detrich, H. W., Zon, L., Westerfield, M. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aDetrich%2C+HW+%22The+Zebrafish%3A+Disease+Models+and+Chemical+Screens%22">The Zebrafish: Disease Models and Chemical Screens</a>. , 4th ed, Academic Press. (2017).</li>
<li>Detrich, H. W., Zon, L. I., Westerfield, M. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aDetrich%2C+HW+%22The+Zebrafish%3A+Genetics%2C+Genomics%2C+and+Transcriptomics%22">The Zebrafish: Genetics, Genomics, and Transcriptomics</a>. , 4th ed, Academic Press. (2016).</li>
<li>Detrich, H. W. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aDetrich%2C+HW+%22The+Zebrafish%3A+Genetics%2C+Genomics+and+Informatics%22">The Zebrafish: Genetics, Genomics and Informatics</a>. , 3rd ed, Academic Press. (2011).</li>
<li>Avdesh, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Regular+care+and+maintenance+of+a+zebrafish+%28Danio+rerio%29+laboratory%3A+an+introduction%5BTitle%5D)+AND+%22J+Vis+Exp%22%5BJournal%5D)">Regular care and maintenance of a zebrafish (Danio rerio) laboratory: an introduction.</a> J Vis Exp. (69), e4196(2012).</li>
<li>Rosen, J. N., Sweeney, M. F., Mably, J. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Microinjection+of+zebrafish+embryos+to+analyze+gene+function%5BTitle%5D)+AND+%22J+Vis+Exp%22%5BJournal%5D)">Microinjection of zebrafish embryos to analyze gene function.</a> J Vis Exp. (25), (2009).</li>
<li>Hwang, W. Y., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Efficient+genome+editing+in+zebrafish+using+a+CRISPR-Cas+system%5BTitle%5D)+AND+%22Nat+Biotechnol%22%5BJournal%5D)">Efficient genome editing in zebrafish using a CRISPR-Cas system.</a> Nat Biotechnol. 31 (3), 227-229 (2013).</li>
<li>Westerfield, M. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aWesterfield%2C+M+%22The+zebrafish+book.+A+guide+for+the+laboratory+use+of+zebrafish+%28Danio+rerio%29%22">The zebrafish book. A guide for the laboratory use of zebrafish (Danio rerio)</a>. , 4th ed, Univ. of Oregon Press. Eugene. (2000).</li>
<li>Rossi, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Genetic+compensation+induced+by+deleterious+mutations+but+not+gene+knockdowns%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">Genetic compensation induced by deleterious mutations but not gene knockdowns.</a> Nature. 524 (7564), 230-233 (2015).</li>
<li>Kimmel, C. B., Ballard, W. W., Kimmel, S. R., Ullmann, B., Schilling, T. F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Stages+of+embryonic+development+of+the+zebrafish%5BTitle%5D)+AND+%22Dev+Dyn%22%5BJournal%5D)">Stages of embryonic development of the zebrafish.</a> Dev Dyn. 203 (3), 253-310 (1995).</li>
<li>Samsa, L. A., Fleming, N., Magness, S., Qian, L., Liu, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Isolation+and+Characterization+of+Single+Cells+from+Zebrafish+Embryos%5BTitle%5D)+AND+%22J+Vis+Exp%22%5BJournal%5D)">Isolation and Characterization of Single Cells from Zebrafish Embryos.</a> J Vis Exp. (109), (2016).</li>
<li>MPIMG. ConsenusPathDB. , Available from:
 <a target="_blank" href="http://cpdb.molgen.mpg.de/">http://cpdb.molgen.mpg.de/</a> (2017).</li>
<li>Consortium, G. O. Enrichment analysis Tool. , Available from:
 <a target="_blank" href="http://www.geneontology.org/">http://www.geneontology.org/</a> (2017).</li>
<li>IDT. IDT Primerquest Tool. , Available from:
 <a target="_blank" href="https://www.idtdna.com/Primerquest/Home/Index">https://www.idtdna.com/Primerquest/Home/Index</a> (2017).</li>
<li>Heid, C. A., Stevens, J., Livak, K. J., Williams, P. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Real+time+quantitative+PCR%5BTitle%5D)+AND+%22Genome+Res%22%5BJournal%5D)">Real time quantitative PCR.</a> Genome Res. 6 (10), 986-994 (1996).</li>
<li>Leitch, C. C., Lodh, S., Prieto-Echague, V., Badano, J. L., Zaghloul, N. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Basal+body+proteins+regulate+Notch+signaling+through+endosomal+trafficking%5BTitle%5D)+AND+%22J+Cell+Sci%22%5BJournal%5D)">Basal body proteins regulate Notch signaling through endosomal trafficking.</a> J Cell Sci. 127 (Pt 11), 2407-2419 (2014).</li>
<li>Lodh, S., Hostelley, T. L., Leitch, C. C., O'Hare, E. A., Zaghloul, N. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Differential+effects+on+beta-cell+mass+by+disruption+of+Bardet-Biedl+syndrome+or+Alstrom+syndrome+genes%5BTitle%5D)+AND+%22Hum+Mol+Genet%22%5BJournal%5D)">Differential effects on beta-cell mass by disruption of Bardet-Biedl syndrome or Alstrom syndrome genes.</a> Hum Mol Genet. 25 (1), 57-68 (2016).</li>
<li>Hostelley, T. L., Lodh, S., Zaghloul, N. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Whole+organism+transcriptome+analysis+of+zebrafish+models+of+Bardet-Biedl+Syndrome+and+Alstrom+Syndrome+provides+mechanistic+insight+into+shared+and+divergent+phenotypes%5BTitle%5D)+AND+%22BMC+Genomics%22%5BJournal%5D)">Whole organism transcriptome analysis of zebrafish models of Bardet-Biedl Syndrome and Alstrom Syndrome provides mechanistic insight into shared and divergent phenotypes.</a> BMC Genomics. 17, 318(2016).</li>
</ol>

作者无任何利益冲突需要披露。

本方案所述方法为全动物或特定分选细胞群体的转录组水平分析提供了一种相对快速且经济有效的策略。斑马鱼是此类研究的理想模型，因其易于快速获得大量起始材料，易于实施遗传或环境实验条件，并且拥有大量转基因报告品系，可用于分离特定细胞类型和特定组织的细胞群体。
尽管此处未详细说明，该方法也可轻松适用于从幼年或成年鱼类中采集的组织切片，作为流式细胞分选和收集的替代方案。我们还开发了一种基本的后续分析策略，仅依赖于基础的电子表格命令以及现有的通路和GO数据库。该方法的主要优势在于易于获取和使用，无需具备高水平的计算机编程或数据库管理专业知识。因此，该策略适用于各种背景的研究人员对大规模基因表达数据集进行信息丰富的分析。
我们的方法结合了基础的斑马鱼胚胎和幼体培养技术，以及标准的RNA提取技术，随后进行RNA测序（RNASeq）。RNA测序需要大量高质量的起始材料；部分供应商要求至少2 µg的总RNA。因此，稀有或罕见细胞类型，以及单个胚胎的分析目前尚难以实现。然而，过去五年已在低产量样本分析、小样本量数据生成及低RNA投入量方面取得了显著进...

全局基因表达的比较分析是识别参与特定表型形成的新基因的有力工具。此类分析通常依赖于对实验组与对照组样本之间转录本丰度的定量评估。针对特定基因的方法，如定量逆转录PCR（qRT-PCR），在单个基因表达变化的研究中具有较快的速度和较高的准确性。而RNA测序（RNASeq）则提供了一种广泛且无需预先假设的方法，可用于识别样本间基因表达的显著差异，现已成为各类实验体系中此类研究的标准方法。
斑马鱼已成为多种疾病研究领域的重要模型。最初因其在发育生物学研究中的实用性而被开发，由于其繁殖力高且维持成本相对较低，斑马鱼的实验应用已扩展至涵盖从胚胎到成体阶段的广泛表型，以及多种分子检测方法。1,2,3事实上，这些优势使得分子机制研究变得快速且经济高效，因为可以轻松获得大量实验材料，同时在生命各个阶段均易于进行遗传和环境操控。此外，斑马鱼胚胎和幼鱼具有透明特性，非常适用于构建细胞和组织特异性的转基因报告品系，从而实现 体内 离散细胞群体的可视化4利用此类细胞系可基于报告基因的表达，对特定分离细胞类型进行全基因表达分析。
本文介绍了一种利用RNA测序（RNASeq）对斑马鱼胚胎进行全基因表达分析的综合实验方案。关于基因操作实验方法，包括基于吗啉环寡核苷酸（MO）的瞬时基因敲低或CRISPR介导的基因组编辑，已在其他文献中有所报道5,6,7。因此，本文重点介绍从整胚或经分选的转基因报告基因表达细胞中提取RNA的详细操作流程，并结合通路分析工具和基因本体（GO）术语对RNASeq结果进行简单的计算分析。最后，本文还提供了通过定量逆转录酶PCR（qRT-PCR）验证基因表达变化的策略。这些实验方案适用于经历多种实验处理的斑马鱼胚胎，包括遗传突变体或不同环境条件之间的比较。

<table><tbody><tr><td>&lt;strong&gt;商用试剂&lt;/strong&gt;</td><td></td><td></td><td></td></tr><tr><td>TriZol</td><td>Thermo Scientific</td><td>15596026</td><td>裂解试剂</td></tr><tr><td>TrypLE</td><td>Gibco</td><td>12604013</td><td>解离缓冲液 1</td></tr><tr><td>FACSMax</td><td>Genlantis</td><td>T200100</td><td>解离缓冲液 2</td></tr><tr><td>DEPC 处理水</td><td>Sigma</td><td>95284</td><td></td></tr><tr><td>FirstStrand cDNA 合成</td><td>Thermo Scientific</td><td>K1621</td><td>cDNA 合成试剂盒</td></tr><tr><td>2X SYBR Green Master Mix</td><td>Roche</td><td>4707516001</td><td>qRT-PCR Master Mix</td></tr><tr><td>FACS 缓冲液</td><td>Fisher Scientific</td><td>50-105-9042</td><td></td></tr><tr><td>氯仿</td><td>Sigma Aldrich</td><td>288306</td><td></td></tr><tr><td>乙酸钠</td><td>Sigma Aldrich</td><td>S2889</td><td></td></tr><tr><td>&lt;strong&gt;名称&lt;/strong&gt;</td><td>&lt;strong&gt;公司&lt;/strong&gt;</td><td>&lt;strong&gt;目录编号&lt;/strong&gt;</td><td>&lt;strong&gt;备注&lt;/strong&gt;</td></tr><tr><td>&lt;strong&gt;斑马鱼品系&lt;/strong&gt;</td><td></td><td></td><td></td></tr><tr><td>Tuebingen</td><td>ZIRC</td><td>ZL57</td><td></td></tr><tr><td>ins2a:mCherry</td><td>ZIRC</td><td>ZL1483</td><td></td></tr><tr><td>&lt;strong&gt;名称&lt;/strong&gt;</td><td>&lt;strong&gt;公司&lt;/strong&gt;</td><td>&lt;strong&gt;目录编号&lt;/strong&gt;</td><td>&lt;strong&gt;备注&lt;/strong&gt;</td></tr><tr><td>&lt;strong&gt;设备&lt;/strong&gt;</td><td></td><td></td><td></td></tr><tr><td>40 微米细胞筛</td><td>Sigma</td><td>CLS431750</td><td></td></tr><tr><td>FACS 管</td><td>BD Falcon</td><td>352063</td><td></td></tr><tr><td>血细胞计数板</td><td>Sigma</td><td>Z359629</td><td></td></tr><tr><td>解剖显微镜</td><td>Zeiss</td><td></td><td></td></tr><tr><td>倒置显微镜</td><td>Zeiss</td><td></td><td></td></tr><tr><td>Nanodrop</td><td>Thermo Scientific</td><td></td><td></td></tr><tr><td>Illumina HiSeq</td><td>Illumina</td><td></td><td></td></tr><tr><td>LightCycler 480</td><td>Roche</td><td></td><td></td></tr><tr><td>交配罐 1.0L 交配罐套装</td><td>Aquaneering</td><td>ZHCT100</td><td></td></tr><tr><td>FACS 管 5 mL 聚丙烯管</td><td>BD Falcon</td><td>352063</td><td></td></tr><tr><td>&lt;strong&gt;名称&lt;/strong&gt;</td><td>&lt;strong&gt;公司&lt;/strong&gt;</td><td>&lt;strong&gt;目录编号&lt;/strong&gt;</td><td>&lt;strong&gt;备注&lt;/strong&gt;</td></tr><tr><td>&lt;strong&gt;软件&lt;/strong&gt;</td><td></td><td></td><td></td></tr><tr><td>Excel</td><td>Microsoft</td><td></td><td></td></tr><tr><td>Consensus Path DB</td><td>http://cpdb.molgen.mpg.de/</td><td></td><td></td></tr><tr><td>GO 富集分析</td><td>http://geneontology.org/page/go-enrichment-analysis</td><td></td><td></td></tr></tbody></table>

斑马鱼rna测序基因表达数据的样本制备与分析

Sample Preparation and Analysis of RNASeq-based Gene Expression Data from Zebrafish

全局基因表达变化的分析是识别与特定表型相关的新型通路的有力工具。由于可从大量个体中便捷地分离RNA，斑马鱼成为快速评估整体转录组（来自完整动物或特定细胞群体）的理想模型。本文介绍了一种利用RNA测序（RNASeq）技术对斑马鱼胚胎进行全局基因表达分析的实验方案。我们描述了从完整胚胎或通过转基因动物中细胞分选获得的特定细胞群体中提取RNA的方法。同时，本文还介绍了一种分析RNASeq数据的策略，用于识别全局基因表达数据集中富集的通路及基因本体（GO）术语。最后，我们提供了一种利用定量反转录PCR（qRT-PCR）验证基因表达变化的实验流程。这些方案可用于对照组与实验组斑马鱼之间的比较分析，以发现新的基因表达变化，并为感兴趣表型提供分子层面的解释。

差异表达基因的筛选：
为了鉴定Alström综合征和Bardet-Biedl综合征（BBS）斑马鱼模型幼虫阶段差异表达的基因，我们通过向野生型斑马鱼胚胎中注射先前已验证的剪接阻断型吗啉代寡核苷酸（MO），靶向alms1或bbs1转录本16,17。受精后5天（dpf），按照上文第3节所述方法，从每种实验条件下分别提取两份RNA样本，同时提取对照MO注射胚胎的RNA。每个样本的RNA测序深度约为1.2亿条reads，其中约1.07亿条reads比对至斑马鱼基因组。比对上的reads中有85–90%定位到基因组的外显子区域。
Alström 模型中总共有 1,348 个基因的表达发生显著变化，...

本方案介绍了一种从斑马鱼胚胎、幼体或分选细胞中进行全转录组分析的方法。我们包括了RNA的分离、RNA测序数据的通路分析，以及基于qRT-PCR的基因表达变化验证。

斑马鱼RNA测序基因表达数据的样本制备与分析

全局基因表达变化的分析是识别潜在表型相关新通路的重要工具。斑马鱼是一种……

本研究对斑马鱼幼体样本进行RNA测序分析的总体目标是鉴定斑马鱼胚胎和幼体的基因表达谱，并对不同样本间的基因表达变化做出定量比较分析。该方法可用于回答多个研究领域中的关键问题，只要这些领域中斑马鱼胚胎或幼体具有信息价值，例如：在与其他条件相比时，特定靶基因的抑制如何引起基因表达的变化或导致某些通路的紊乱。该技术的主要优势在于斑马鱼易于大量繁殖，且可方便地获取整个生物体的转录组数据。此外，利用转基因动物的分选可以轻松实现特定细胞类型的分离。演示该实验操作的是我实验室的博士研究生莱恩·霍斯特利（Lain Hostelley）和博士后杰西卡·内斯米思（Jessica Nesmith）。首先，将胚胎培养至三个月龄，即达到生殖成熟阶段。在收集胚胎的前一晚，将两条雄性成鱼和三条雌性成鱼从目标品系中分离出来，放入装有新鲜系统水的分隔式交配缸中。次日早晨开灯后，移除隔板，让鱼自然交配，直至观察到胚胎出现在缸底。每隔30分钟，在不同的含胚胎培养基的培养皿中收集胚胎，直至获得所需数量。为对胚胎进行分期，将每10厘米的培养皿中放入50至75枚胚胎进行群体培养，以促进所有胚胎发育时间的一致性。随后将培养皿置于28.5°C条件下培养。在原肠胚形成后至受精后约24小时（HPF）期间，通过体节数量来测定胚胎年龄，并根据发育阶段对胚胎进行分类。根据文本方案中所述，在目标阶段处死胚胎后，将约20个胚胎转移至已标记的1.5毫升微量离心管中。随后，去除微量离心管中多余的胚胎培养基。向管中加入200 µL裂解试剂。使用研杵对胚胎进行机械匀浆。然后加入额外的800微升裂解试剂，使总体积达到1毫升。为提取RNA，向收集的样品中加入裂解试剂，并在室温下孵育5分钟。每使用1毫升裂解试剂，加入0.2毫升氯仿，用手倒置离心管15秒。将样品在室温下孵育2至3分钟后，在4摄氏度条件下以12,000 × g离心15分钟。将分离后的水相转移至新的离心管中，每使用1毫升裂解试剂加入0.5毫升异丙醇。室温孵育试管10分钟后，将样品离心10分钟。接着，向RNA中加入75%乙醇，并在7,500 × g、4 °C条件下离心5分钟。彻底去除上清液，将样品置于室温下空气干燥。然后，将RNA重悬于15至30微升经DEPC处理的水中。在根据文本方案提取沉淀并洗涤样品后，为纯化高质量RNA，使用吸收光谱仪检测所提取RNA的浓度和纯度。在将RNA样品送至供应商或公司进行基于测序读长定量的RNA测序和基因表达变化分析之前，确保A260/A230比值约为2.0。为了比较单个实验条件与对照条件，使用电子表格管理软件打开差异表达基因的数据。点击排序按钮旁边的下拉箭头，在电子表格中选择“自定义排序”。在弹出的窗口中，选择“列”下方的框，并选择按 LFC 列进行排序。在“排序依据列”中，确保选择“数值”。在“顺序列”中，选择从大到小排序，然后点击确定。为了确定在两种实验条件下差异表达的基因，在实验组与对照组的电子表格中，选择空白列的第一个单元格。然后，在单元格中输入以下方程。按 Enter 或 Return 键运行该方程。选中包含该方程的单元格，然后单击单元格右下角的方框。按住鼠标并向下拖动所选区域，沿该列一直选到最后一个需要填入公式的特征ID，以将公式复制到该列的每个单元格中。再次选择自定义排序，并通过点击左下角的加号图标添加第二级排序。在第一级排序中，按列选择“duplicate”进行排序。在“排序依据”中选择“数值”，在“次序”中选择“Z 到 A”。在第二级排序中，在“列”下选择“LFC”；在“排序依据”中选择“数值”，在“次序”中选择“从大到小”，然后点击“确定”。在继续操作前，若需去除基因符号列中的所有括号，请选中包含基因符号的整列。在下拉文件菜单中选择“编辑”，然后选择“替换”。在“查找内容”栏中输入左括号、星号、右括号，即“(*)”，并将“替换为”栏留空。选择“全部替换”，以删除所有括号实例。为了确定富集的通路，将目标基因符号复制到剪贴板。前往 ConsensusPathDB 网站，然后在网页左侧边栏选择基因集分析，接着选择过度代表分析。在“粘贴基因和蛋白质标识符列表”框中，粘贴基因列表。在基因/蛋白质标识符类型框中选择基因符号，然后点击继续。在基于通路的集合部分，勾选路径由通路数据库定义旁边的复选框。选择查找富集集合，以获得包含输入列表中基因的通路列表。通过选中通路名称旁的每个复选框，或点击选择框上方列标题“select”中的“all”，来选择所有富集通路，以在通路网络中进行可视化。然后，点击“visualize selected sets”。通过点击页面顶部中央的“relative overlap”或“shared candidates”复选框，输入期望的百分比、相对重叠值或共享候选分子数量，然后点击“apply”以调整相应过滤条件。为了确定富集的基因本体，将相应基因组的基因符号复制到剪贴板。前往基因本体联盟（Gene Ontology Consortium）的基因本体富集分析工具。在页面左侧的“在此处输入您的基因ID”部分，将基因符号列表粘贴到输入框中。在基因 ID 框下方，选择 GO 术语“生物过程”。然后在 GO 术语框下方选择 danio rerio，然后点击提交。根据文本方案进行 QRT PCR，并与 RNASeq 进行比较。对注射了靶向 alms1 或 bbs1 的吗啉代寡核苷酸的幼虫提取RNA进行测序，揭示了在Alström综合征模型和Bardet-Biedl综合征模型中特异性上调和下调的基因，以及在两种模型中均发生显著变化的基因。为了更清晰地阐明Alström模型的分子特征，鉴定了差异表达基因中富集的通路和基因本体。如图所示，共有31条通路被上调。除了代谢这一大类之外，受影响最显著的通路是先天性和适应性免疫系统，分别涉及32个和20个基因。B细胞受体信号通路的下游事件也出现富集。在上调基因中，多个信号通路同样呈现富集，这与Alstrom综合征和初级纤毛功能障碍之间的关联一致。此外，与胰岛素分泌相关的三条通路被上调，分别为与GPR40结合的脂肪酸、游离脂肪酸和乙酰胆碱。最后，在Alstrom模型中上调基因富集了六个GO条目，包括红细胞分化、红细胞稳态、髓样细胞稳态以及稳态过程。一旦掌握，通路和GO条目分析可在一小时之内完成。在进行该实验流程时，重要的是要牢记，通路参数和阈值的选择是解读RNASeq表达比较结果时最关键的考虑因素。完成该流程后，还可采用其他方法，如利用突变体品系或对特定细胞类型进行转录组分析，以回答更多科学问题，例如细胞类型的转录组特征分析，以及特定基因调控下的基因表达变化。观看本视频后，您应能够充分理解如何利用RNASeq生成的斑马鱼转录组数据，来鉴定实验样本间显著的基因表达差异。

sample-preparation-analysis-rnaseq-based-gene-expression-data-from

Research

科研

JoVE Journal

Genetics

遗传学

斑马鱼 RNASeq-based 基因表达数据的样品制备与分析

高质量完整RNA的分离是许多实验室实验方案中的关键步骤。在此，我们展示了从完整的斑马鱼胚胎中提取RNA并合成cDNA，以用于后续多种实验操作，包括基因表达微阵列分析。...

从斑马鱼胚胎中提取RNA及用于基因表达分析的cDNA合成

基因微阵列是在全基因组水平上进行基因表达谱分析的强大工具。该技术广泛应用于包括发育生物学和毒理学在内的多种生物学领域。在本视频中，我们详细介绍了一种利用全面的寡核苷酸微阵列平台对斑马鱼进行全基因组表达分析的实验方案。...

利用斑马鱼寡核苷酸微阵列平台进行全基因组表达分析

本文介绍了一种通过自然交配产生斑马鱼胚胎并对其发育阶段进行分选的实验方案。该方案进一步用于利用转基因鱼分析胰腺β细胞。...

交配与卵子分期：一种按发育阶段生成并分选胚胎的方法

为了分析斑马鱼心脏发育过程中的心脏基因表达谱，必须从分离的心脏中提取总RNA。本文介绍了一种通过快速手工解剖斑马鱼胚胎以获取功能性/跳动心脏的方案，从而获得心脏特异性的mRNA。...

大规模斑马鱼胚胎心脏解剖用于转录分析

本方案描述了一种从斑马鱼胚胎中分离单个细胞、富集目标细胞、在基于微流控的单细胞多重系统中捕获斑马鱼细胞，以及检测单个细胞基因表达的方法。...

斑马鱼胚胎单细胞的分离与表征

使用斑马鱼进行DNA复制时间分析 在体 内 模型系统

本文介绍了一种利用自然病原体Mycobacterium...

在感染Mycobacterium marinum的成年斑马鱼中建立结核病模型

骨骼肌再生由组织驻留的肌肉干细胞驱动，这些细胞在多种肌肉疾病（如肌营养不良症）中功能受损，导致肌肉无法再生。本文介绍了一种可在斑马鱼肌肉疾病模型中研究肌肉再生的实验方案。...

研究肌肉疾病斑马鱼模型中的肌肉再生

斑马鱼靶向反应性亲电体和氧化剂（Z-REX）是一种基于化学生物学的方法，用于研究反应性小分子信号传导。该技术可应用于不同发育阶段的活体斑马鱼。在此，我们将斑马鱼中的标准检测方法与Z-REX相结合，用于信号通路分析。...

监测斑马鱼幼虫中的靶向信号响应——Z-REX揭示亲电性药物与代谢物的精确作用机制

在此，我们介绍一种在斑马鱼幼鱼中用荧光蛋白标签标记蛋白质的实验方案，该方案为新近开发并经过改进的方法 体内...

荧光标记法用于可视化斑马鱼中低表达蛋白

RAS病是一类由RAS-MAPK信号通路过度激活引起的多系统遗传综合征。随着大量潜在致病变异不断被发现，亟需对其进行功能验证，然而临床前证据的缺乏限制了治疗策略的开发。本文中，我们描述了一种在Teen-报告基因斑马鱼中通过活体FRET成像检测RAS病相关ERK活化水平及其在胚胎发育过程中药物调控作用的体内实验方案。...

一种通过活细胞FRET显微镜和免疫荧光高效评估斑马鱼努南综合征早期模型中ERK活性调控的实验方案

学习斑马鱼RNASeq转录组分析方法，用于识别胚胎和幼鱼基因表达谱并进行定量比较。该技术结合RNA提取、通路分析和qRT-PCR验证，为发育生物学和神经科学研究提供可靠框架。

斑马鱼RNA测序基因表达数据的样本制备与分析

探索更多视频