天津文章成稿指导数据科学活动

    GSEA基本原理从方法上来讲，GSEA主要分为基因集进行排序、计算富集分数（EnrichmentScore，ES）、估计富集分数的***性水平并进行多重假设检验三个步骤。**步对输入的所有基因集L进行排序，通常来说初始输入的基因数据为表达矩阵，排序的过程相当于特定两组中（case-control、upper-lower等等）基因差异表达分析的过程。根据所有基因在两组样本的差异度量不同（共有六种差异度量，默认是signal2noise，GSEA官网有提供公式，也可以选择较为普遍的foldchange)，对基因进行排序，并且Z-score标准化。第二步是GSEA的**步骤，通过分析预先定义基因集S在**步获得的基因序列上的分布计算富集指数EnrichmentScore，并绘制分布趋势图Enrichmentplot。每个基因在基因集S的EnrichmentScore取决于这个基因是否属于基因集S及其差异度量（如foldchange）。差异度量越大基因的EnrichmentScore权重越大，如果基因在基因集S中则EnrichmentScore取正，反则取负。将基因集L在基因集S里的所有基因的EnrichmentScore一个个加起来，就是Enrichmentplot上的EnrichmentScore趋势，直到EnrichmentScore达到**值，就是基因集S**终的EnrichmentScore。第三步是为了检验第二部获得结果的统计学意义。 蛋白组代谢组个性化分析。天津文章成稿指导数据科学活动

Bubbles可以同时展示pvalue和表达量。例如展示motif的pvalue和motif对应的转录因子的表达量，方便快速看出转录因子富集且高表达所在的group，预示着该分组对细胞状态的改变（例如细胞分化、转移、应激）起关键调控作用；例如做基因功能富集分析时，展示富集的通路qvalue和基因数量或geneRatio。

基本原理：

Bubbles的实质是分组数据下基因表达量或通路内基因数量的可视化，同时可以展示pvalue。

数据要求：

表达矩阵，分组 天津组学实验数据科学欢迎咨询自有服务器机房，可随时调用各计算平台算力，且团队成员有多年科研经历。

    pancancer泛**图谱泛*研究是通过整合不同**类型、不同组织起源的**表达数据，查找**之间的共性或者差异的过程。通常使用**数据信息较为***的TCGA数据，通过分裂小提琴图展示某个基因在TCGA**和正常组织中的表达差异。分裂小提琴图(ViolinPlot)结合了箱形图和密度图的特征，主要用来显示数据的分布形状，它一般应用于对比某一基因在TCGA**组织和正常组织基因表达量TPM值或其它表达量数据。基本原理：小提琴图(ViolinPlot)使用一组数据中的最小值、**四分位数、中位数、第三四分位数和**值来反映数据分布的中心位置和散布范围，将多组数据的小提琴图画在同一坐标上，可以清晰地显示各组数据的分布差异。分裂小提琴图在小提琴图的基础上又加入了分组对比项，便于观察多**类型在某一基因上的表达分布情况，或者某一基因在某一**上，其疾病与正常的对比表达差异情况。

    STEM基因表达趋势分析基因调控网络是一个连续且复杂的动态系统。当生物体按照一定顺序发生变化或者受到外界环境刺激（如受到不同浓度的化学药物诱导）时，基因表达变化也会呈现趋势特征。趋势分析就是发现基因表达的趋势特征，将相同变化特征的基因集中在一种变化趋势中，从而找到实验变化过程中相当有有代表性的基因群。STEM（ShortTime-seriesExpressionMiner），中文名短时间序列表达挖掘器。该软件主要用于分析短时间实验数据，也可用于多组小样本数据。推荐3至8组数据。一般可应用的研究方向有：多个时间点的时间序列数据，例如多个发育时期、处理后多个时间点取样。基本原理STEM采用了一种新的聚类算法来分析时间序列基因表达趋势。聚类算法首先选择一组不同的、有代表性的时间表达模式（temporalexpressionprofiles）作为模型（modelprofiles）。模型是**于数据选择的，并从理论上保证了所选择的模型剖面具有代表性。然后，根据每个标准化过后的基因表达模式，分配给模型中相关系数比较高的时间表达模式。由于模型的选择是**于数据的，因此该算法可以通过排列测试，确定哪些时间表达模式在统计意义上***富集基因。对每一个基因都分配时间表达模式完成后。 检测服务及数据分析助力取得2020年国自然面上十项、青年基金十八项。

    STEM基因表达趋势分析数据要求表达谱芯片或测序数据（已经过预处理）下游分析得到***富集的时间表达模式之后的分析有：1.时间表达模式中基因的功能富集2.时间表达模式中基因表达与性状之间的相关性挖掘模块的关键信息：1.找到时间表达模式中的**基因2.利用关系预测该时间表达模式功能文献1：DynamicEBF1occupancydirectssequentialepigeneticandtranscriptionaleventsinB-cellprogramming（于2018年1月发表在GenesDev.，影响因子）EBF1动态占据在B细胞中对序列表观遗传和转录过程的影响该文献采用基因表达趋势分析，探寻了EBF1诱导前后25kb转录起始位点内基因转录水平的差异，来寻找EBF1对特定功能基因的影响以及造成影响的时间节点。文献2：ComprehensivetranscriptionalprofilingofNaCl-stressedArabidopsisrootsrevealsnovelclassesofresponsivegenes（于2016年10月发表在BMCPlantBiol.，影响因子）该文献采用基因表达趋势分析，研究了高浓度盐水作用不同时间下拟南芥根的基因表达差异，来探寻在遇到高浓度盐水时拟南芥在基因层面上的应对方式。 乳腺类疾病预后相关信性基因突变研究数据包。湖北数据科学方案

协助构建各类科研、临床数据库。天津文章成稿指导数据科学活动

    蛋白质主要由碳、氢、氧、氮等化学元素组成，是一类重要的生物大分子。蛋白质的功能由蛋白质的三维结构决定。蛋白质三维结构绘图，可以直观地展示蛋白质三维功能结构，广泛应用于单核苷酸突变功能分析、药物蛋白分子相互作用分析等研究领域。基本原理蛋白质三维结构绘图主要分为蛋白质三维结构预测以及对结构进行可视化两步。蛋白质三维结构预测是基于蛋白质中氨基酸序列预测蛋白质折叠结构的步骤，**常用的预测方法为同源建模，同源建模的原理是序列相似的蛋白质具有相似的蛋白质结构，要推测一个未知结构蛋白的三维结构，只需要找到与之序列高度相似的已知结构模板。在无法进行同源建模（找不到模型）的情况下，还有折叠识别及从头建模法，但是计算量大运行缓慢且建模准确度不如同源建模。获得蛋白质三维结构预测的pbd文件后还需要通过分子三维结构软件绘制可视化的三维图，并分析特殊位点（分子对接或突变位点分析），常用的有pymol和DeepView等。数据要求目标蛋白的氨基酸序列或者编码蛋白的基因序列，突变数据等。下游分析突变位点靶向药物分析等。 天津文章成稿指导数据科学活动