固有动力学可视化工具是一个交互式可视化软件包,可连接至基因调控网络推断工具,以增强并简化功能网络模型的生成。该可视化工具可用于更合理地确定推断工具的参数,从而提高所得模型的可靠性。
需要JoVE订阅才能观看此内容。 请登录或开始免费试用
固有动力学可视化工具是一个交互式可视化软件包,可连接至基因调控网络推断工具,以增强并简化功能网络模型的生成。该可视化工具可用于更合理地确定推断工具的参数,从而提高所得模型的可靠性。
构建基因调控网络模型是系统生物学中的一个重大挑战。为应对这一挑战,已开发出多种计算工具和分析流程,其中包括新近推出的固有动力学流程(Inherent Dynamics Pipeline)。该流程由多个先前发表的工具协同组成,以线性方式连接,前一个工具的输出作为下一个工具的输入。与大多数计算技术一样,固有动力学流程的每一步都需要用户对缺乏明确生物学定义的参数做出选择,而这些选择会显著影响分析所生成的基因调控网络模型。因此,能够在每一步可视化并探索不同参数选择所带来的影响,有助于增强对参数设定及其结果的信心。固有动力学可视化工具(Inherent Dynamics Visualizer)是一个综合性的可视化软件包,通过网页浏览器中的交互式界面,简化了评估参数选择的过程。用户可分别检查流程中每一步的输出,依据可视化信息进行直观调整,并自动获得固有动力学流程所需的输入文件。固有动力学可视化工具为从时间序列转录组数据中发现基因调控网络这一高度复杂的工具提供了前所未有的可访问性。
许多重要的生物学过程,例如细胞分化和环境响应,均由基因调控网络(GRN)中相互作用的基因集合所控制。这些基因调控网络能够产生激活和维持其所控制表型所需的转录动力学,因此,确定基因调控网络的组成成分及其拓扑结构是理解众多生物学过程与功能的关键。基因调控网络可被建模为一组相互作用的基因和/或基因产物,其网络中的节点代表基因,边则描述相互作用的方向和形式(例如转录激活/抑制、翻译后修饰等)1。随后,这些相互作用可表示为参数化的数学模型,用以描述调控基因对其靶基因表达产生的影响2,3,4。推断基因调控网络模型需要同时推断出相互作用网络的结构以及估计其底层的相互作用参数。目前已开发出多种计算推断方法,可输入时间序列的基因表达数据并输出基因调控网络模型5。最近,一种名为“固有动力学流程”(Inherent Dynamics Pipeline, IDP)的新型基因调控网络推断方法被提出,该方法利用时间序列基因表达数据生成具有标注的调控-靶向关系的基因调控网络模型,且所生成的模型能够再现基因表达数据中观察到的动力学特征6。IDP 是一系列线性连接的工具组合,可分为三个步骤:第一步为“节点发现”(Node Finding),根据已知或推测与基因调控网络功能相关的基因表达特征对基因进行排序7,8;第二步为“边发现”(Edge Finding),对成对的调控关系进行排序8,9;第三步为“网络发现”(Network Finding),生成能够再现观测到的动力学行为的基因调控网络模型10,11,12,13,14,15。
与大多数计算方法类似,IDP 需要一组用户指定的参数,以决定输入数据的分析方式,而不同的参数组合可能在相同数据上产生不同的结果。例如,包括 IDP 在内的多种方法均包含对数据施加某种阈值的参数,若在连续运行中提高或降低该阈值,可能导致各次运行结果不一致(参见补充说明10:5中的网络推断方法)。理解每个参数如何影响分析过程及后续结果,对于提高结果的可信度至关重要。与大多数基因调控网络(GRN)推断方法不同,IDP 由多个计算工具组成,每个工具都有用户必须指定的独立参数集,并产生各自的输出结果。尽管 IDP 提供了关于如何设置每个工具参数的详尽文档,但由于各工具的输出依赖于前一步骤的结果,因此在不进行中间分析的情况下对整个流程进行参数设置具有挑战性。例如,边(Edge)和网络发现(Network Finding)步骤中的参数通常需要依据先验的生物学知识来确定,因而会依赖于具体的数据集和/或生物体。为了深入探究中间结果,用户需要具备基本的编程知识,以及对 IDP 生成的所有结果文件及其内容的全面理解。
固有动力学可视化工具(IDV)是一个可在用户浏览器窗口中运行的交互式可视化软件包,为IDP用户提供了一种评估其参数选择对IDP任意步骤结果影响的方法。IDV能够导航由IDP生成的复杂目录结构,收集每个步骤所需的数据,并以直观且可交互的图表和表格形式呈现,供用户探索。在浏览这些交互式显示内容后,用户可根据更充分的判断,从IDP的某一步骤生成新的数据。这些新数据可立即用于IDP的下一个相应步骤。此外,对数据的探索还有助于判断某个IDP步骤是否应以调整后的参数重新运行。IDV能够增强IDP的使用效果,使其使用过程更加直观和易于理解,这一点已在对酵母细胞周期核心振荡基因调控网络(GRN)的研究中得到验证。以下实验方案比较了完全参数化运行IDP所获得的结果,与在IDP各步骤(即节点、边和网络发现)运行后结合使用IDV的方法所获得的结果。
访问受限。请登录或开始试用以查看此内容。
1. 安装 IDP 和 IDV
注意:本节假设您已安装 docker、conda、pip 和 git(材料表)。
2. 节点识别
3. 边缘检测
4. 网络识别
访问受限。请登录或开始试用以查看此内容。
上述以文字描述并以图1图示的步骤被应用于酵母细胞周期的核心振荡基因调控网络(GRN),以检验是否能够发现具有功能性的GRN模型,这些模型能够重现从酵母细胞周期研究中收集的时间序列基因表达数据所观察到的动力学特征16。为了说明IDV如何澄清并改进IDP的输出结果,本文比较了两种分析方式的结果:1)一次性完整运行IDP的所有步骤而不使用IDV;2)在IDV的辅助下逐步执行IDP流程,从而能够根据先验生物学知识以及IDP输出结果进行中间结果的调整和更精细的选择。作为示例所使用的酵母细胞周期GRN已被广泛研究,其多数调控关系已通过实验验证。如果研究对象为其他物种或注释较少的生物过程,则对中间结果或参数进行调整的方式可能有所不同。为了展示可用于评估网络的一类查询方法,本文测量了每个网络的鲁棒性,以支持稳定的振荡行为,并使其节点在模型参数变化下与观察到的转录动力学相匹配。
来自 Orlando 2008 的两组重复时间序列基因表达数据16经...
访问受限。请登录或开始试用以查看此内容。
基因调控网络(GRN)的推断是系统生物学中的一个重要挑战。IDP 利用一系列工具,以逐步复杂化的方式处理基因表达数据,从而生成模型化的 GRN。每一步都需要决策如何处理数据,以及哪些元素(基因、功能相互作用)将被传递到 IDP 的下一层。这些决策对 IDP 结果的影响并不明显。为此,IDV 提供了 IDP 中 GRN 推断工具各步骤输出结果的有用交互式可视化功能。IDV 简化并促进了对这些计算推断方法结果的评估过程,加快实验进程并指导分析选择,从而有助于加速生成高置信度的网络模型和科学假设。IDV 还实现了扩展 IDP 功能的特性,包括根据 LEM ODE 参数选择对边进行过滤、按基因表达时间进行分箱,以及基于与模体或网络的相似性对网络进行聚类。重要的是,IDV 允许在 IDP 的每一步之间进行人工干预,使用户能够便捷地整合来自文献的人类知识和先验信息,而这些操作难以完全自动化。IDP 的默认运行方式无法原生整合此类信息,因此在实验特异性信息可用时,使用 IDV 将提高结果的可信度。总体而言,将 IDV 与 IDP 结合使用,可使用户在对真实 GRN 知之甚少甚至一无所知的情况下,更有信心地构建针对生物过程的网络假设。
访问受限。请登录或开始试用以查看此内容。
作者无任何利益冲突需要披露。
本工作由美国国立卫生研究院(NIH)资助项目 R01 GM126555-01 和美国国家科学基金会(NSF)资助项目 DMS-1839299 支持。
访问受限。请登录或开始试用以查看此内容。
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| Docker | https://docs.docker.com/get-docker/ | ||
| Git | https://git-scm.com/ | ||
| 固有动力学分析流程 | https://gitlab.com/biochron/inherent_dynamics_pipeline | ||
| 固有动力学可视化工具 | https://gitlab.com/bertfordley/inherent_dynamics_visualizer | ||
| Miniconda | https://docs.conda.io/en/latest/miniconda.html | ||
| Pip | https://pip.pypa.io/en/stable/ |
访问受限。请登录或开始试用以查看此内容。