Name: 分析肿瘤基因表达因子与CorExplorer门户网站 (视频) | JoVE Journal
Uploaded: 2019-10-11T15:00:00.000+00:00
Duration: 8 min
Description: 探索CorExplorer门户网站，利用CorEx机器学习算法分析肿瘤RNA测序因子，通过交互式可视化研究基因表达与生存结局的关系。该平台整合蛋白质相互作用数据，助力癌症生物学研究并识别潜在治疗靶点。

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

GV 获得了美国国防高级研究计划局（DARPA）奖项 W911NF-16-0575 的支持。

1. 探索包含目的基因的因素
<ol>
	<li>打开网页浏览器，访问 http://corex.isi.edu，进入 CorExplorer 主页。</li>
	<li>在右侧的 快速链接 下方，点击 Ovarian (TCGA-OV) 旁边的 + 展开按钮，查看基于 TCGA 卵巢癌数据训练得到的 CorEx 因子图的摘要（如 图 1 所示）。也可点击其他项目进行比较。</li>
	<li>完成因子图查看后，点击 Lung (TCGA-LUAD)，进入肺癌 RNA-seq 的 CorExplorer 页面。
	<ol>
		<li>使用 CorExplorer 的“因子图”窗口，探索目标基因的 CorEx 因子图。
		<ol>
			<li>将鼠标指针移至因子图显示窗口上，使用鼠标滚轮或触控板缩放因子图，以查看图中细节，例如每个因子中最重要的基因以及不同层级节点之间的连接。或者，点击并拖动以移动视图区域或任意节点。</li>
			<li>要查找目标基因（此处以 BRCA1 为例），点击因子图窗口顶部的 基因 下拉菜单，输入“BRCA1”以在下拉列表中选择该基因，然后按 Return 键，视图将自动缩放至因子 26，即与 BRCA1 相关性最强的因子。</li>
			<li>将鼠标重新移至图显示区域，滚动以缩小视图，查看二级节点 L2_8 及其与因子 26 相邻的关联因子。请注意，仅显示权重高于 最小连接权重 滑块所设阈值的基因。</li>
			<li>要查看与该因子关联的所有基因，点击 L1_26 节点，并在弹出窗口中选择 加载更多基因。当显示“Done”（完成）时，关闭弹出窗口。</li>
			<li>返回因子图窗口上方的标题区域，抓取并拖动 最小连接权重 调节器。当滑块向下移动至 0.05 时，因子 L1_26 中的其他基因（包括 BRCA2）将按权重顺序显示。可选操作：通过抓取并拖动节点以优化布局。</li>
		</ol>
		</li>
		<li>通过生存分析窗口查询，确定患者按因子分层对生存结果的影响。
		<ol>
			<li>在生存分析窗口中，取消勾选 按 p 值排序，然后在 单因子 下拉菜单中选择因子 26，以显示因子 26 的生存曲线。</li>
			<li>向下滚动生存图，以在横轴上显示处于风险中的患者数量。</li>
		</ol>
		</li>
		<li>通过注释窗口查询，寻找与生物学功能的关联。
		<ol>
			<li>在注释窗口中，若要按因子编号而非错误发现率（FDR）对 因子 下拉菜单进行排序，请取消勾选 FDR 排序。</li>
			<li>滚动并点击注释窗口的下拉菜单，选择因子 26，以显示该因子的功能富集注释。</li>
			<li>向下滚动注释列表，直到看到 DNA 修复，点击该项，即可在图中立即看到相关基因以黄色高亮显示。参见 图 2 中间面板。</li>
			<li>请注意，当选择不同的 GO 术语（例如“DNA 损伤响应中的内在凋亡信号通路”）时，因子会根据是否富集了具有所选注释的基因而出现或消失。</li>
		</ol>
		</li>
		<li>通过添加具有不同功能的窗口，进一步探索因子。
		<ol>
			<li>从顶部菜单栏中，选择 添加窗口 下拉菜单中的 PPI，然后点击 添加 按钮，在显示区域中添加一个蛋白质-蛋白质相互作用网络（PPI）窗口。在 PPI 图窗口中，选择因子“Layer1: 26”以显示蛋白质-蛋白质相互作用。注意连接的密度。</li>
			<li>从顶部菜单栏中，选择 添加窗口 下拉菜单中的 热图 而非 PPI，然后点击 添加 按钮，在显示区域中添加一个热图窗口。在热图窗口中，选择因子“Layer1: 26”以显示基因表达模式。</li>
			<li>抓取并重新定位热图窗口，使生存分析窗口也可见。观察热图顶部的橙色/蓝色/灰色条带，其对应于生存图中的患者风险分层。结果如 图 2 底部所示。</li>
		</ol>
		</li>
	</ol>
	</li>
</ol>
2. 使用基因权重、生存和注释数据过滤和解释 CorEx 因子
<ol>
	<li>使用生存数据和聚类质量筛选感兴趣的因子。
	<ol>
		<li>在页面顶部的“Dataset”下拉菜单中，选择 TCGA_OVCA，进入TCGA卵巢癌RNA-seq的CorExplorer页面。</li>
		<li>页面加载完成后，从生存分析窗口中注意到，不同分层之间生存差异最大的因子是114。</li>
		<li>在因子图谱窗口顶部的Factor下拉菜单中，选择“Layer1: 114”。</li>
		<li>用鼠标拖动链接权重滑块至0.5。注意因子114中基因数量较多（共1609个），但没有任何基因的权重&#62;0.35，表明其聚类相对较弱。</li>
		<li>接着，在生存分析窗口中展开因子列表，并在下拉菜单中选择生存差异次优的因子39，以查看其对应的生存曲线。</li>
		<li>在注释窗口中点击选择因子39，系统将显示其显著相关的GO和KEGG注释。</li>
	</ol>
	</li>
	<li>为了更深入理解因子39中基因的生物学功能，可按以下步骤利用邻近注释信息对因子进行解读。
	<ol>
		<li>在因子图谱窗口顶部的因子下拉菜单中选择“Layer1: 39”。然后将鼠标移至因子图谱窗口并缩小视图，以完整显示包含6个因子的L2_14聚类：14、32、39、42、52和82（如图3所示）。</li>
		<li>为理解与L2_14节点相连各因子的相对重要性，首先查看每个L2_14因子的生存差异。在生存分析窗口中取消勾选Sort by p-val，然后依次点击各个因子编号。通过此操作可注意到，仅有因子14、32和39显示出与生存相关的显著性。</li>
		<li>现在从顶部菜单栏再次选择Add Window下拉菜单中的PPI选项，点击Add，在显示区域添加一个PPI图谱窗口。在PPI图谱窗口中选择因子“Layer1: 52”，以显示显著的蛋白质-蛋白质相互作用。此时窗口布局的一个示例如图3所示。</li>
		<li>点击PPI窗口底部的View at StringDB链接，跳转至StringDB在线数据库。在首个页面点击Continue，然后像之前一样选择网络图下方的Analysis标签页，以获取PPI网络基因的在线GO分析结果。其中最显著的细胞组分为“MHC II类蛋白复合物”。</li>
		<li>返回CorExplorer标签页及PPI窗口，从因子下拉菜单中选择因子32，再次点击View at StringDB链接跳转至StringDB分析页面。此时最显著的细胞组分为“MHC I类蛋白复合物”，与上一步中因子52的II类结果形成对比！</li>
		<li>最后，回到PPI窗口，从顶部的因子下拉菜单中选择“Layer1: 39”。</li>
		<li>点击View at StringDB链接跳转至StringDB分析页面。在首个页面点击Continue，然后选择网络图下方的Analysis标签页，获取PPI网络基因的在线GO分析结果。可观察到最显著的分子功能为“CXCR3趋化因子受体结合”。</li>
	</ol>
	</li>
</ol>
3. 利用生存数据和数据库注释寻找有前景的治疗组合
<ol>
	<li>从Dataset下拉菜单中选择TCGA_SKCM，切换至TCGA黑色素瘤CorExplorer界面。</li>
	<li>注意，具有最大生存差异的因子是因子171。通过滚动查看因子171的注释，可注意到“免疫应答”和“细胞因子介导的信号通路”位于注释前列（与之前卵巢癌顶部因子的情况一致）。</li>
	<li>为寻找一个互补因子，请查看与生存显著相关的顶部因子及其对应的顶级注释术语。为此，请点击顶部菜单栏中的Dataset overview链接，打开一个包含数据集处理详情的独立标签页，其中还列出了按生存差异p值排序的顶部因子摘要。请注意，第一个非免疫相关因子为因子88。</li>
	<li>返回至TCGA_SKCM浏览器标签页。</li>
	<li>在生存、注释和图形窗口中选择因子88。其前几项GO术语与“rRNA加工”和“线粒体组织”相关，证实该因子与免疫相关因子明显不同。</li>
	<li>在生存窗口的成对因子下拉菜单中，选择“88_171”，以观察同时具有中等表达水平的因子171和因子88的患者群体其生存情况的改善。注释与生存比较结果如图4所示。</li>
</ol>
4. 使用搜索页面查找不同肿瘤类型间基因表达变异的共性与差异
<ol>
	<li>点击 CorExplorer 返回首页</li>
	<li>点击 搜索 在顶部菜单栏上点击，进入可搜索 CorExplorer 网站所有数据集的页面。</li>
	<li>在 基因 搜索框中输入“FLT1”（VEGFR1）并点击 返回 或按压 搜索FLT1 在以下因素中权重相对较高：OVCA - 76、LUAD - 162、SKCM - 195 和 SKCM - 184，以及 COAD - 112 和 COAD - 74。</li>
	<li>或者，在所有数据集中搜索相关的 GO 术语。您可以在“GO 搜索”框中输入“angiogenesis”并按下回车进行尝试 返回 或按压 搜索所有FLT1因子（除SKCM-195外）均被列为在“血管生成”基因中具有统计学富集意义；因子195实际上也具有该注释，但低于默认的10-8阈值。本步骤及前一步骤的搜索结果如图所示。 图5.</li>
	<li>作为进一步的示例，在GO搜索框中首先输入“epidermal growth factor receptor”。仅有LUAD对该术语富集，这是肺癌中一个众所周知的分层因子。接着，在搜索框中输入“mesenchymal”。该术语在OVCA的基因表达组中富集，而在OVCA中，它是一个被广泛研究的分层因子。</li>
</ol>

<ol><li>Petryszak, R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+RNASeq-er+API-a+gateway+to+systematically+updated+analysis+of+public+RNA-seq+data%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">The RNASeq-er API-a gateway to systematically updated analysis of public RNA-seq data.</a> Bioinformatics. 33, 2218-2220 (2017).</li>
<li>Steeg, G. V., Galstyan, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Maximally+Informative+Hierarchical+Representations+of+High-Dimensional+Data%5BTitle%5D)+AND+%22Proceedings+of+the+Eighteenth+International+Conference+on+Artificial+Intelligence+and+Statistics+%28AISTATS%29%22%5BJournal%5D)">Maximally Informative Hierarchical Representations of High-Dimensional Data.</a> Proceedings of the Eighteenth International Conference on Artificial Intelligence and Statistics (AISTATS). , San Diego, CA. (2015).</li>
<li>Ver Steeg, G., Galstyan, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Discovering+structure+in+high-dimensional+data+through+correlation+explanation%5BTitle%5D)+AND+%22Advances+in+Neural+Information+Processing+Systems%22%5BJournal%5D)">Discovering structure in high-dimensional data through correlation explanation.</a> Advances in Neural Information Processing Systems. , Montreal, Canada. (2014).</li>
<li>Pepke, S., Ver Steeg, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Comprehensive+discovery+of+subsample+gene+expression+components+by+information+explanation%3A+therapeutic+implications+in+cancer%5BTitle%5D)+AND+%22BMC+medical+Genomics%22%5BJournal%5D)">Comprehensive discovery of subsample gene expression components by information explanation: therapeutic implications in cancer.</a> BMC medical Genomics. 10, 12(2017).</li>
<li>Byron, S. A., Van Keuren-Jensen, K. R., Engelthaler, D. M., Carpten, J. D., Craig, D. W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Translating+RNA+sequencing+into+clinical+diagnostics%3A+opportunities+and+challenges%5BTitle%5D)+AND+%22Nature+Reviews+Genetics%22%5BJournal%5D)">Translating RNA sequencing into clinical diagnostics: opportunities and challenges.</a> Nature Reviews Genetics. 17, 257(2016).</li>
<li>Cancer Genome Atlas Research Network. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Comprehensive+molecular+profiling+of+lung+adenocarcinoma%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">Comprehensive molecular profiling of lung adenocarcinoma.</a> Nature. 511, 543(2014).</li>
<li>Cancer Genome Atlas Network. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Comprehensive+molecular+characterization+of+human+colon+and+rectal+cancer%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">Comprehensive molecular characterization of human colon and rectal cancer.</a> Nature. 487, 330(2012).</li>
<li>Akbani, R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Genomic+classification+of+cutaneous+melanoma%5BTitle%5D)+AND+%22Cell%22%5BJournal%5D)">Genomic classification of cutaneous melanoma.</a> Cell. 161, 1681-1696 (2015).</li>
<li>Cancer Genome Atlas Research Network. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Integrated+genomic+analyses+of+ovarian+carcinoma%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">Integrated genomic analyses of ovarian carcinoma.</a> Nature. 474, 609(2011).</li>
<li>Grossman, R. L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Toward+a+shared+vision+for+cancer+genomic+data%5BTitle%5D)+AND+%22New+England+Journal+of+Medicine%22%5BJournal%5D)">Toward a shared vision for cancer genomic data.</a> New England Journal of Medicine. 375, 1109-1112 (2016).</li>
<li>Moynahan, M. E., Chiu, J. W., Koller, B. H., Jasin, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Brca1+controls+homology-directed+DNA+repair%5BTitle%5D)+AND+%22Molecular+Cell%22%5BJournal%5D)">Brca1 controls homology-directed DNA repair.</a> Molecular Cell. 4, 511-518 (1999).</li>
<li>Szklarczyk, D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((STRING+v11%3A+protein%E2%80%93protein+association+networks+with+increased+coverage%2C+supporting+functional+discovery+in+genome-wide+experimental+datasets%5BTitle%5D)+AND+%22Nucleic+Acids+Research%22%5BJournal%5D)">STRING v11: protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets.</a> Nucleic Acids Research. 47, 607-613 (2018).</li>
<li>Durgeau, A., Virk, Y., Corgnac, S., Mami-Chouaib, F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Recent+advances+in+targeting+CD8+T-cell+immunity+for+more+effective+cancer+immunotherapy%5BTitle%5D)+AND+%22Frontiers+in+Immunology%22%5BJournal%5D)">Recent advances in targeting CD8 T-cell immunity for more effective cancer immunotherapy.</a> Frontiers in Immunology. 9, 14(2018).</li>
<li>Sato, E., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Intraepithelial+CD8%2B+tumor-infiltrating+lymphocytes+and+a+high+CD8%2B%2Fregulatory+T+cell+ratio+are+associated+with+favorable+prognosis+in+ovarian+cancer%5BTitle%5D)+AND+%22Proceedings+of+the+National+Academy+of+Sciences+of+the+United+States+of+America%22%5BJournal%5D)">Intraepithelial CD8+ tumor-infiltrating lymphocytes and a high CD8+/regulatory T cell ratio are associated with favorable prognosis in ovarian cancer.</a> Proceedings of the National Academy of Sciences of the United States of America. 102, 18538-18543 (2005).</li>
<li>De Moura, M. B., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Mitochondrial+respiration-an+important+therapeutic+target+in+melanoma%5BTitle%5D)+AND+%22PLoS+One%22%5BJournal%5D)">Mitochondrial respiration-an important therapeutic target in melanoma.</a> PLoS One. 7, 40690(2012).</li>
<li>Folkman, J., Merler, E., Abernathy, C., Williams, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Isolation+of+a+tumor+factor+responsible+for+angiogenesis%5BTitle%5D)+AND+%22Journal+of+Experimental+Medicine%22%5BJournal%5D)">Isolation of a tumor factor responsible for angiogenesis.</a> Journal of Experimental Medicine. 133, 275-288 (1971).</li>
<li>Takahashi, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Vascular+endothelial+growth+factor+%28VEGF%29%2C+VEGF+receptors+and+their+inhibitors+for+antiangiogenic+tumor+therapy%5BTitle%5D)+AND+%22Biological+and+Pharmaceutical+Bulletin%22%5BJournal%5D)">Vascular endothelial growth factor (VEGF), VEGF receptors and their inhibitors for antiangiogenic tumor therapy.</a> Biological and Pharmaceutical Bulletin. 34, 1785-1788 (2011).</li>
<li>Subramanian, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Gene+set+enrichment+analysis%3A+a+knowledge-based+approach+for+interpreting+genome-wide+expression+profiles%5BTitle%5D)+AND+%22Proceedings+of+the+National+Academy+of+Sciences+of+the+United+States+of+America%22%5BJournal%5D)">Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles.</a> Proceedings of the National Academy of Sciences of the United States of America. 102, 15545-15550 (2005).</li>
<li>Cerami, E., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+cBio+Cancer+Genomics+Portal%3A+An+Open+Platform+for+Exploring+Multidimensional+Cancer+Genomics+Data%5BTitle%5D)+AND+%22Cancer+Discovery%22%5BJournal%5D)">The cBio Cancer Genomics Portal: An Open Platform for Exploring Multidimensional Cancer Genomics Data.</a> Cancer Discovery. 2, 401-404 (2012).</li>
<li>Gao, J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Integrative+Analysis+of+Complex+Cancer+Genomics+and+Clinical+Profiles+Using+the+cBioPortal%5BTitle%5D)+AND+%22Science+Signalling%22%5BJournal%5D)">Integrative Analysis of Complex Cancer Genomics and Clinical Profiles Using the cBioPortal.</a> Science Signalling. 6, 1(2013).</li>
<li>Reich, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((GenePattern+2.0%5BTitle%5D)+AND+%22Nature+Genetics%22%5BJournal%5D)">GenePattern 2.0.</a> Nature Genetics. 38, 500(2006).</li>
<li>Wang, Y. E., Kutnetsov, L., Partensky, A., Farid, J., Quackenbush, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((WebMeV%3A+A+Cloud+Platform+for+Analyzing+and+Visualizing+Cancer+Genomic+Data%5BTitle%5D)+AND+%22Cancer+Research%22%5BJournal%5D)">WebMeV: A Cloud Platform for Analyzing and Visualizing Cancer Genomic Data.</a> Cancer Research. 77, 11-14 (2017).</li>
<li>Morpheus. , Available from: 
 <a target="_blank" href="https://software.broadinstitute.org/morpheus">https://software.broadinstitute.org/morpheus</a> (2019).</li>
<li>Weitschek, E., Lauro, S. D., Cappelli, E., Bertolazzi, P., Felici, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((CamurWeb%3A+a+classification+software+and+a+large+knowledge+base+for+gene+expression+data+of+cancer%5BTitle%5D)+AND+%22BMC+Bioinformatics%22%5BJournal%5D)">CamurWeb: a classification software and a large knowledge base for gene expression data of cancer.</a> BMC Bioinformatics. 19, 354(2018).</li>
<li>Chou, P. -H., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((tACCo%2C+a+Database+Connecting+transcriptome+Alterations%2C+pathway+Alterations+and+Clinical+outcomes+in+Cancers%5BTitle%5D)+AND+%22Scientific+Reports%22%5BJournal%5D)">tACCo, a Database Connecting transcriptome Alterations, pathway Alterations and Clinical outcomes in Cancers.</a> Scientific Reports. 9, 3877(2019).</li>
</ol>

作者声明不存在任何竞争性经济利益。

我们介绍了CorExplorer网站，这是一个公开可访问的网络服务器，用于交互式探索由CorEx算法从肿瘤RNA-seq数据中学习到的最大相关基因表达因子。我们展示了如何利用该网站根据肿瘤基因表达对患者进行分层，并说明了这种分层与生物学功能及生存预后之间的对应关系。
已有其他用于RNA-seq分析的网络服务器。在cbioPortal中，可对肿瘤进行差异表达和共表达分析，并将其与其他类型的数据进行整合19,20。GenePattern21、Mev22和Morpheus23等服务器集成了成熟的聚类技术，例如主成分分析（PCA）、k均值（kmeans）或自组织映射（SOMs）。更具创新性的尝试包括CamurWeb24，其基于一种自动规则生成的分类器；以及TACCO25，该工具实现了随机森林分类器和套索回归（lassos）。本文所采用的CorEx算法通过优化多变量信息，以发现能够解释数据模式的因子层级结构。相较于通过PCA获得的线性全局因子4，这种非线性且具有层级结构的因子学习方法在可解释性方面表现更优。此外，该技术对样本信号的细粒度解析能力，使得相较于常用的宽泛亚型分类，能够实现更为精确的肿瘤间比较。这种重叠性与层级性因子分析的结合，使CorExplorer区别于大多数其他方法，同时也需要开发新的可视化与总结工具。
CorExplorer 因子分析的一个关键部分是能够探索的因子不仅限于几个，而是超过 100 个具有信息性基因模式的因子，这些因子被置于一个重叠的层级结构中。CorExplorer 有助于挖掘这些大量因子与生物学和临床特征之间的关联，并实现对单个肿瘤的异常精细的表征。由于如此大量的因子是通过无监督学习获得的，因此并非所有因子都与疾病生物学相关。在这种情况下，必须利用注释信息或已知基因来提取感兴趣的因子，或搜索与生存等临床数据相关的因子。因此，CorExplorer 允许用户执行这一非常重要的筛选步骤。肿瘤中存在特定因子的基因模式，甚至可能提示个性化肿瘤治疗的策略。此外，每个肿瘤具有多个因子评分，这有助于发现潜在的有用治疗组合。
有时，与生存高度相关的因子可能未显示出显著的 GO 注释。这可能是由于数据噪声或采样不足所致，但也可能存在其他原因，例如聚类尺寸过小而无法获得显著的富集分数，或该组别是由来自不同通路且无明确生物学关联的单个基因组成的“混合类别”。此外，可能需要采用不同于 KEGG 和 GO 生物过程的注释类别，例如细胞组分。这些分析可通过本方案中所述链接至 StringDB 获取。目前 CorExplorer 网站上的基因本体富集分析尚未考虑因子中基因的权重，但这一问题有望在不久的将来得到解决。请注意，在“添加窗口”下提供了一个基因列表选项，可用于下载完整的因子基因列表，以便使用外部工具进行进一步分析。
为网站目的，CorEx 对每个数据集运行了五次，并保留产生最高总体总相关性的那次运行结果。对多次运行结果进行统计学表征可能更具信息量，这也是未来工作的目标。此外，服务器上可用的肿瘤类型集合相对较小，但预计将根据用户需求随时间逐步扩展。
如上所述，CorExplorer 可视化 CorEx RNA-seq 因子之间的关系，并整合临床和数据库信息，从而支持多种不同的分析模式。我们希望该工具能够推动进一步的研究，以利用 RNA-seq 分析的强大能力，在肿瘤学领域实现发现和临床应用。

自十余年前问世以来，RNA-seq 已成为测量基因表达的普遍工具1。这得益于其能够快速且低成本地对样本的整个转录组进行从头分析。然而，RNA-seq 肿瘤数据反映的是内在复杂且常被采样不足的生物学背景，而数据本身又具有高维度和噪声大的特点。这为提取可靠信号带来了重大挑战。CorEx 算法利用多变量互信息，在此类情况下识别细微模式2,3。该技术此前已被用于分析癌症基因组图谱（The Cancer Genome Atlas, TCGA）中的卵巢肿瘤 RNA-seq 样本，在此背景下，其表现明显优于更常用的分析方法4。
尽管RNA测序（RNA-seq）在包括肿瘤学在内的各类研究领域中已得到极为广泛的应用，但这些努力尚未推动其在临床干预中的广泛应用5。造成这一现状的部分原因在于，目前尚缺乏针对这些特定问题而设计的用户友好型算法和软件。为弥合这一差距，我们开发了CorExplorer网络平台，使来自不同背景的研究人员能够利用CorEx机器学习算法所发现的结果，研究肿瘤RNA-seq样本中的基因表达因子。CorExplorer平台支持对多种不同类型肿瘤（包括肺癌、结肠癌、黑色素瘤和卵巢癌）的因子进行交互式可视化和查询6,7,8,9,10，旨在帮助研究人员梳理数据中的相关性，识别潜在通路，从而为治疗目的对患者进行分层。
我们预计 CorExplorer 门户网站可能对多种类型的用户具有帮助。该门户在设计时考虑了用户需求，旨在帮助用户理解公共数据库中驱动肿瘤基因表达差异的广泛因素，并可能将单个基因表达谱置于具有相似特征的肿瘤背景下进行分析。除了本文概述的代表性实验方案外，CorExplorer 的研究还可作为提出假设的起点，用于进一步验证；也可用于比较和对照 CorExplorer 以外数据集上的 CorEx 分析结果，以及将单个肿瘤中一个或少数几个基因的病理表达特征与可能受到协同影响的更大基因群组相联系。最后，对于刚进入该领域的研究人员而言，它还可作为将机器学习应用于 RNA-seq 分析的用户友好型入门工具。

<table><tbody><tr><td>CorExplorer 网站的公共服务器</td><td>USC</td><td>http://corex.isi.edu</td><td>Intel Xeon E5-2690 4 核 2.6 GHz，8GB 内存。后端架构为 LAMP：Linux、Apache、MySQL、PHP。</td></tr><tr><td>网页浏览器</td><td>Google/Apple</td><td>Chrome/Safari</td><td>已验证的网页浏览器。</td></tr></tbody></table>

利用 corexplorer 网络平台分析肿瘤基因表达因子

Analyzing Tumor Gene Expression Factors with the CorExplorer Web Portal

差异基因表达分析是理解疾病状态的重要技术。机器学习算法CorEx在分析肿瘤RNA-seq数据中基因群组的差异表达方面已显示出应用价值，可能有助于推动精准肿瘤学的发展。然而，CorEx会产生大量因子，这些因子在分析时可能难以与现有知识建立联系。为了促进此类关联，我们构建了一个名为CorExplorer的网站，使用户能够交互式地探索数据，并回答与其分析相关的常见问题。我们使用四种肿瘤类型的RNA-seq基因表达数据训练了CorEx，分别为卵巢癌、肺癌、黑色素瘤和结直肠癌。随后，我们将相应的生存数据、蛋白质-蛋白质相互作用、基因本体（GO）和京都基因与基因组百科全书（KEGG）通路富集结果以及热图整合到网站中，以关联因子图谱的可视化结果。本文通过示例方案说明如何利用该数据库结合这些外部数据，深入理解所学习到的肿瘤因子的生物学意义。

在肺癌数据集中搜索基因“BRCA1”发现，该基因与CorEx因子26的关联性最强（图2）。该因子的GO术语富集程度极高，其中DNA修复的FDR仅为1 × 10-19。该选择还突出了第二层聚类L2_8，其下包含六个密切相关的子因子。在GO术语注释或因子图的GO富集下拉菜单中选择“DNA修复”，可高亮显示各因子中相关的基因，其中因子26所关联的基因数量远超其他因子，符合预期11。蛋白质-蛋白质相互作用网络显示出强连接性，进一步支持因子26中基因功能的高度关联性。相应的生存曲线图提示其可能与患者生存相关，但需在更大规模的数据集中进一步验证。
从生存分析入手，可以剖析特定基因表达组别与生存改善之间的关联原因。例如，在卵巢癌中影响生存的首要因素是编号39的因子，该因子显著富集了与免疫系统相关的基因（图3）。另有五个与同一二级节点（level 2 node）相关的因子也被识别为与免疫相关，但它们对生存的影响差异显著，其中因子39的影响最强，而因子52的影响最弱。为某一因子添加蛋白质-蛋白质相互作用（PPI）窗口，可显示其直接的相互作用网络，并支持链接至StringDB12网站，以查询该PPI网络基因的各类富集情况。依次对L2_14的每个因子进行此操作后发现，StringDB对PPI网络基因的富集分析提示了以下可能的生存关联解释：因子32包含编码主要组织相容性复合体（MHC）I类蛋白复合物的基因，该复合物可被细胞毒性T淋巴细胞识别；因子39对应于细胞因子信号通路及CXCR3受体结合功能，与CD8+ T淋巴细胞相关。这两个因子在相应基因表达水平较高的患者中均表现出显著的生存优势。细胞毒性CD8+ T淋巴细胞主要负责抗肿瘤免疫。相反，因子52由编码MHC II类复合物蛋白的基因组成，该复合物主要被CD4+辅助性T细胞识别，而非直接被细胞毒性T淋巴细胞识别。其余的L2_14因子则反映了非特异性的免疫系统整体激活，无法区分上述两类淋巴细胞群体。生存关联特异性地体现在细胞毒性T淋巴细胞对MHC I类细胞抗原的识别，这一结果与我们对肿瘤免疫（包括黑色素瘤等其他癌症）的一般认知一致13,14。
该网络平台支持发现具有互补功能的因子对，可能提示有效的肿瘤特异性联合治疗方案。可通过浏览数据集概览，筛选出与生存相关但具有不同基因本体（GO）富集特征的因子。以黑色素瘤（TCGA_SKCM；图4）为例，排名最高的生存相关因子171与免疫相关，而列表中靠后的因子88则在与线粒体组织相关的基因中表现出富集。事实上，已有研究提出将后者作为黑色素瘤的潜在治疗靶点15。在CorExplorer页面中添加生存时间窗口，可比较使用因子对与单独使用每个因子进行分层的效果，结果显示，同时具备两组有利基因表达模式的患者群体，其生存趋势优于仅具备单一因子有利模式的患者。然而，最高风险分层的生存情况并未明显改善，提示对于部分患者，单独使用免疫治疗可能是最佳选择。
通过跨数据集搜索基因或 GO 术语，可以观察到肿瘤之间的共性与差异图5）。例如，FLT1（又称VEGFR1）是一个被广泛研究的促血管生成标志物16,17当将其输入搜索栏时，所有肿瘤均显示出FLT1发挥重要作用的因素。相反，当在搜索页面输入GO术语“angiogenesis”（血管生成）时，6个FLT1组中有5个显示出该富集结果。除SKCM-195外，所有FLT1因素均被列为在“血管生成”基因上具有统计学意义的富集。第六个因素实际上也具有该注释，但低于默认的10-8阈值。当在替代的富集分析工具（例如基因集富集分析，Gene Set Enrichment Analysis, GSEA）中利用因子列表内的权重时18第六个因子也被发现显著富集了“血管生成”相关基因。
检查热图以确保基因表达模式具有足够的质量来支持生物学解释非常重要。显示强烈且清晰变异的热图可能表现出因子基因从低到高协调表达的模式，或者更复杂的模式，其中一些基因的低表达与另一些基因的高表达相关（图6）。高质量聚类的一个关键标志是存在多个基因，其表达随因子评分呈平滑变化。因子热图中的样本按因子评分排序，因此从左到右应呈现平滑的梯度变化。然而，这种情况至少可能在两种不同情况下无法实现。最常见的是，相关性可能极度嘈杂（图5C），从而质疑关于生存和/或生物学功能推论的稳健性和实用性。此外，仅在极少数样本中出现的模式可能不符合CorEx算法所假设的三种表达状态模型，导致对样本的分类产生误导（图5D右侧）。
<img alt="卵巢癌基因表达分析网页界面，含CorEx机器学习示意图。" class="xfigimg" src="/files/ftp_upload/60431/60431fig01.jpg" data-altupdatedat="1788282455000" data-alt="图1"> 
图1：CorExplorer首页界面。在快速链接下的卵巢癌旁点击+后，将显示因子图的详细信息。CorEx分层模型由底层的输入变量（本例中为基因表达）和上层推断出的潜在因子构成。<a href="https://www.jove.com/files/ftp_upload/60431/60431fig01large.jpg" target="_blank">请点击此处查看该图的放大版本。</a>
<img alt="通过 CoR Explorer 进行生物数据中的图论与网络分析；交互式流程步骤。" class="xfigimg" src="/files/ftp_upload/60431/60431fig02.jpg" data-altupdatedat="1788282455000" data-alt="图 2"> 
图 2：使用基因名称引导探索。 该图展示了一系列截图，说明如何探索与 BRCA1 高度相关的 CorEx 肺癌因子。首先，在因子图的 基因 下拉框中选择“BRCA1”，会使图形视图聚焦于 BRCA1 权重最大的因子。稍微缩小视图可显示连接该因子与其他相关因子的第二层节点 L2_8。可比较生存信息和注释：点击 GO 术语 DNA 修复 可高亮显示已注释的基因。添加一个蛋白质相互作用（PPI）窗口，以展示该因子中基因的网络互作关系。使用 添加窗口 按钮添加热图，可显示表达模式与生存之间的关联，提示 DNA 修复基因的表达升高可能与生存率降低相关。<a href="https://www.jove.com/files/ftp_upload/60431/60431fig02large.jpg" target="_blank">请点击此处查看该图的放大版本。</a>
<img alt="结合生存分析的蛋白质相互作用网络，展示基因关系及生存曲线的示意图。" class="xfigimg" src="/files/ftp_upload/60431/60431fig03.jpg" data-altupdatedat="1788282455000" data-alt="图3"> 
图3：利用临床数据（生存信息）指导探索性分析。 对与卵巢癌生存最相关的主成分因子（因子39）进行分析，可揭示其邻近因子之间有趣的关系。在因子关系图中选中因子39并稍作缩小后，可见第二层中与因子39相连的因子还与其他五个因子相关联。额外的生存分析窗口可直接比较这些因子对应的生存差异。因子39和因子32均显示与生存呈正相关，而因子52则无此相关性。所有蛋白质-蛋白质相互作用网络均结构清晰。通过链接至StringDB可进一步比较基因本体（GO）注释信息（图中未显示）：因子39与细胞因子信号通路相关，涉及细胞毒性CD8+ T淋巴细胞的激活；因子32则主要由MHC I类抗原呈递蛋白构成，可触发此类淋巴细胞的识别；然而，其邻近因子则主要由其他免疫系统组分（如CD4+辅助性T细胞）主导，且未显示出与生存的相关性。<a href="https://www.jove.com/files/ftp_upload/60431/60431fig03large.jpg" target="_blank">请点击此处查看该图的放大版本。</a>
<img alt="基因表达分析图表，包含生存曲线、注释和数据可视化。" class="xfigimg" src="/files/ftp_upload/60431/60431fig04v3.jpg" data-altupdatedat="1788282455000" data-alt="Figure 4"> 
图4：探索主要生存因子提示潜在的治疗组合。 在主页菜单栏上的“数据集”链接可导向一个按p值排序的生存因子简明表格，并附有最主要的GO注释（未显示）。利用该信息针对黑色素瘤分析，发现因子171（与免疫功能相关）与因子88（与线粒体组织相关）的组合具有互补性。本图并列展示了两个因子各自的注释窗口，以便进行对比。根据这两个因子单独或联合对患者进行分层所绘制的生存曲线表明，联合使用这两个因子比单独使用任一因子更能增大生存差异。<a href="https://www.jove.com/files/ftp_upload/60431/60431fig04large.jpg" target="_blank">请点击此处查看该图的放大版本。</a>
<img alt="基因和GO搜索结果图表；用于研究分析的数据集、因子、权重、MI。" class="xfigimg" src="/files/ftp_upload/60431/60431fig05.jpg" data-altupdatedat="1788282455000" data-alt="图5"> 
图5：搜索页面支持泛癌分析。可通过主页上的搜索链接，在所有数据集中查找基因或GO生物学过程术语。该图显示了基因FLT1和GO术语“angiogenesis”的搜索结果。结果显示，在多种癌症中，被注释为“angiogenesis”术语的因子均包含FLT1。<a href="https://www.jove.com/files/ftp_upload/60431/60431fig05large.jpg" target="_blank">请点击此处查看此图的放大版本。</a>
<img alt="使用 CoRExplorer 进行基因表达热图分析；TCGA-LUAD 数据集的可视化表示。" class="xfigimg" src="/files/ftp_upload/60431/60431fig06v2.jpg" data-altupdatedat="1788282455000" data-alt="图 6"> 
图 6：热图可用于根据因子得分定性评估基因与样本之间的相关性。 当患者按热图中的因子得分排序时，高质量的基因表达关系表现为平滑的渐变。最左侧的因子18热图即为一个示例。这些模式也可能包含复杂表达特征，例如中间因子11的大热图所示的上下调表达模式。低质量的模式有时会显示部分患者群体的表达出现突变，如右侧因子9的热图所示，或表现为高度噪声的简单相关性，如右下角因子161的热图所示。 <a href="https://www.jove.com/files/ftp_upload/60431/60431fig06large.jpg" target="_blank">请点击此处查看该图的放大版本。</a>

我们介绍 CorExplorer 网络门户，这是一个用于探索由机器学习算法 CorEx（相关性解释）发现的肿瘤 RNA 测序因子的资源平台，并展示如何将这些因子与生存数据、数据库注释、蛋白质-蛋白质相互作用以及彼此之间的关系进行分析，从而深入理解肿瘤生物学机制及治疗干预策略。

利用 CorExplorer 网络平台分析肿瘤基因表达因子

差异基因表达分析是理解疾病状态的重要技术。机器学习算法CorEx已显示出在该领域的应用价值……

CorExplorer 采用数学上严谨的方法，结合交互式可视化以及临床和数据库信息，发现与肿瘤相关的基因表达因子，用于生物学探索。CorEx 在高维数据中寻找聚类时，能够以相对较少的样本实现信息最大化。其所生成的潜在因子层级结构有助于深入理解生物学机制。该方法在应用于癌症生物学时已显示出其重要价值。由于其基于信息论的理论基础，因此可适用于任何具有大量变量且对变量间关系先验知识较少的系统。首先访问 CorExplorer 主页。在“快速链接”下方，点击加号展开按钮，查看基于目标癌症数据训练得到的 CorEx 因子图的摘要。例如，此处展示了基于 TCGA 卵巢癌数据生成的因子图。检查因子图后，点击 lung TCGA-LUAD，进入肺癌 RNA 测序的 CorExplorer 页面，并利用 CorExplorer 因子图窗口探索目标基因的 CorEx 因子图。在因子图显示窗口中移动鼠标指针时，使用鼠标触控板缩放因子图，以查看图中的细节，例如每个因子中最重要的基因以及不同层级节点之间的连接关系。要定位目标基因，请点击“Gene”菜单，输入基因名称以在下拉列表中选择该基因。按下键盘上的 Return 键，视图将自动缩放至与该目标基因相关性最强的因子。将鼠标重新置于图形显示区域，滚动以缩小视图，查看该节点及其邻近的、与最密切相关基因因子关联的因素水平。请注意，仅显示权重超过“最小连接权重”滑块所设定阈值的基因。如需查看与该因子相关联的所有基因，可点击相应节点，并在弹出窗口中选择“加载更多基因”。当显示“完成”时，关闭弹出窗口。在页眉部分，单击并拖动最小链接权重调节器至0.05，以使基因按权重顺序显示。为识别与生物学功能相关的关联，在注释窗口中取消勾选“错误发现率排序”，以按因子编号而非错误发现率对“因子”下拉菜单进行排序。滚动并点击“注释”窗口的下拉菜单，选择感兴趣的因子，以显示该因子的富集注释。然后点击一个富集因子，即可在图表显示中立即查看与之相关的基因，这些基因将以黄色高亮显示。请注意，当选择不同的 GO 术语时，某些因子可能会消失或出现，这取决于它们是否在具有所选注释的基因中发生富集。针对感兴趣的滤过因子，结合生存数据和聚类质量，从“数据集”下拉菜单中选择 TCGA_OVCA，进入 TCGA 卵巢癌 RNA 测序的 CorExplorer 页面。从“生存”窗口中观察具有最大生存差异的因子，并在“因子图”窗口的“因子”下拉菜单中选择该因子。将“连接权重”滑块点击并拖动至 0.5，记录该因子所包含的基因数量。展开生存窗口中的因子列表，然后点击生存窗口下拉菜单中下一个较优的因子，以查看其对应的生存曲线。显著的 GO 和 KEGG 注释将被显示。为了更好地理解该因子中基因的生物学作用，请在因子图谱窗口顶部选择“因子”层，并在缩小视图的同时将鼠标移过窗口，以显示整个簇及其相关因子。为了理解与聚类节点相关的各因素的相对重要性，在生存分析窗口中取消勾选“按 p 值排序”，然后依次点击各个因素编号，查看其对应的生存关联情况，并记录显示出生存关联的因素。在“添加窗口”菜单中选择 PPI，然后点击“添加”，将 PPI 图谱窗口添加到显示区域。在 PPI 图谱窗口中，选择感兴趣的因素层，以展示具有显著性的蛋白质-蛋白质相互作用。点击“在StringDB中查看”链接，连接至STRINGdb在线数据库，然后点击“继续”。接着打开“分析”标签页，以获取PPI网络基因的在线GO分析结果。最显著的细胞组分将被显示出来。返回 CorExplorer 选项卡和 PPI 窗口，选择另一个因子。再次点击 View at StringDB 链接，将显示不同的主要细胞组分。然后，在PPI窗口中为STRING数据库分析选择另一个因子。为了探究不同肿瘤类型间基因表达变异的共性与差异，请点击CorExplorer标题返回首页，然后点击“搜索”，进入可检索CorExplorer网站所有数据集的页面。在基因搜索框中输入感兴趣的基因名称，然后点击搜索。例如，如 FLT1 所示，该基因以相对较高的权重出现在多个不同因子中。在肺癌数据集中搜索 BRCA1 基因，发现该基因与 CorEx 因子 26 的关联性最强。该因子的 GO 术语富集程度极高，其中 DNA 修复的错误发现率仅为 1×10⁻¹⁹。该选择还突显了第二个层级的聚类L2_8，其下有六个密切相关的因子作为子节点。DNA修复相关的蛋白质-蛋白质相互作用网络具有强连接性，进一步支持了因子26中基因功能的高度关联性。相应的生存曲线图提示这些基因可能与患者生存相关，但需在更大规模的数据集中加以验证。从生存评估开始，可以分析与特定基因表达组相关的生存率改善所对应的因素。依次为每个因素添加蛋白质-蛋白质相互作用图谱，有助于确定这些因素与生存率相关联的可能解释。重要的是要检查每个因素的热图，以确认基因表达模式具有足够的质量，能够支持生物学解释。显示基因表达模式显著且清晰变化的热图可能表现出因子基因的协调表达，其表达水平从高到低连续变化，也可能呈现更复杂的模式，其中某些基因的低表达与另一些基因的高表达相关联。本实验方法的总体目标是通过将样本外肿瘤映射到CorExplorer因子上，从而确定潜在的肿瘤特异性治疗方案，实现个体化治疗。

analyzing-tumor-gene-expression-factors-with-corexplorer-web

Research

科研

JoVE Journal

Cancer Research

癌症研究

分析肿瘤基因表达因子与CorExplorer门户网站

核酸在存档组织中的降解、肿瘤异质性以及缺乏新鲜冷冻组织样本，可能对全球病理学实验室的癌症诊断服务产生不利影响。本文描述了利用多重磁珠检测技术优化一组生物标志物，用于乳腺癌的分类。...

利用乳腺癌存档样本优化基于RNA的多重表达检测方法

本文介绍了一种基于对来自多种公开数据库的合并临床数据集进行综合分析，以逐步开展表达、相关性及生存分析策略，探索乳腺癌生物标志物和生存预测因子的实验方案。...

利用多个公开数据库对乳腺癌生物标志物进行数据挖掘与整合分析

生物信息学是一种处理大规模数据集的有效方法。通过应用生物信息学手段，研究人员能够快速、可靠且高效地获得具有洞察力的应用成果和科学发现。本文展示了生物信息学在卵巢癌研究中的应用，并通过实验成功验证了生物信息学分析的结果。...

整合生物信息学方法与实验验证以解析Notch信号通路在卵巢癌中的作用

本文介绍了一种用于发现维持结直肠HT29细胞来源的已建立癌干细胞样细胞的过表达驱动基因的实验方案。通过RNA测序结合现有的生物信息学分析，研究并筛选基因表达网络，以阐明与靶向肿瘤细胞存活相关的潜在机制。...

HT29 来源的结直肠癌干细胞样肿瘤球体中驱动基因的发现

Tropomodulin 3 过表达作为卵巢癌铂类耐药性与免疫浸润的标志物

本方案概述了培养、操作和分析小鼠CAFs的方法，包括 体外/体内...

小鼠乳腺肿瘤来源的癌症相关成纤维细胞作为分子与计算研究的工具

本实验方案演示了如何在果蝇（Drosophila）眼/触角成虫盘（EAD）中生成带有荧光标记且遗传背景明确的克隆肿瘤。方案详细描述了如何从三龄幼虫中解剖出EAD和脑组织，并对其进行处理，以观察和定量基因表达变化及肿瘤侵袭性。...

Drosophila 幼虫盘肿瘤模型：利用遗传镶嵌体可视化与定量分析基因表达及肿瘤侵袭性

通过基因转移在新型RIP-Tag; RIP-tva小鼠模型中鉴定和表征转移因子

本方案旨在介绍如何在癌细胞系中建立功能性Tet-ON系统，并进一步应用该系统，特别是用于研究肿瘤细胞来源的蛋白质在单核细胞/巨噬细胞向肿瘤微环境募集过程中的作用。...

在癌细胞系中条件性敲低基因表达以研究单核细胞/巨噬细胞向肿瘤微环境的募集

我们描述了一种用于临床前体内追踪癌症转移的实验方案。该方法基于一种放射性核素-荧光双模报告系统，结合了可通过无创[18F]四氟硼酸盐-PET检测的钠碘同向转运体以及用于简化体外验证的荧光蛋白。该方法不仅适用于肿瘤生物学研究，还可广泛应用于其他领域的临床前体内细胞示踪。...

放射性核素-荧光报告基因成像技术用于追踪啮齿类动物肿瘤模型中的肿瘤进展

探索CorExplorer门户网站，利用CorEx机器学习算法分析肿瘤RNA测序因子，通过交互式可视化研究基因表达与生存结局的关系。该平台整合蛋白质相互作用数据，助力癌症生物学研究并识别潜在治疗靶点。

利用 CorExplorer 网络平台分析肿瘤基因表达因子

探索更多视频