一种植物生长发育动态监控的方法技术

本发明专利技术公开了一种植物生长动态监控的方法，包括：获取两组单细胞转录组数据，其中第一组数据为对照组，第二组数据为变化组，分别对应不同的发育阶段；通过奇异值分析，构建第二组数据在对照组对应的发育阶段的虚拟细胞状态矩阵；得到第二组数据对应的单细胞在不同发育阶段之间的状态变化响应特征；利用状态变化响应特征对单细胞进行聚类。本发明专利技术解决了单细胞技术销毁测序细胞因此无法对植物生长发育过程中长时间跨度的每个细胞状态变化追踪的问题，使得本发明专利技术可以挖掘出已有方法无法发现的细胞在植物生长过程中的动态分化特征。的细胞在植物生长过程中的动态分化特征。的细胞在植物生长过程中的动态分化特征。

【技术实现步骤摘要】
一种植物生长发育动态监控的方法


[0001]本专利技术涉及植物
，尤其涉及一种植物生长发育动态监控的方法。

技术介绍

[0002]近年来，单细胞测序技术(single cell sequencing)的发展使得其成为炙手可热的热点技术，越发成熟的产业链也使得相关研究项目和成果如雨后春笋般迸发出来。相比传统的bulk测序，单细胞转录组(scRNA
‑
seq)测序技术可以为研究者提供单个细胞的基因表达情况。2019年，研究者首次将scRNA
‑
seq技术应用于拟南芥的根组织发育研究。该研究鉴定出韧皮部、柱细胞、QC等细胞类型的主要表达特征，并进一步通过拟时间序列分析揭示了拟南芥根部细胞分化的复杂过程以及关键转录因子如LAR3、ATHB
‑
20和GATA4等。随之往后，国内外出现更多植物单细胞相关的研究。这些研究目前主要集中于拟南芥、番茄、水稻等模式植物的细胞分化、组织发育形成等研究领域。
[0003]植物单细胞相关的研究正在为研究者提供全新的研究思路和角度。例如在一项利用单细胞技术对拟南芥根尖发育的研究中，研究者绘制了拟南芥根尖的细胞图谱。研究者鉴定出根尖干细胞、形成层细胞、木质部细胞等多种细胞群，并且通过伪时间序列分析描绘了它们的发育轨迹和对应的部分分化相关基因表达变化。结果表明，拟南芥根尖包含高异质性的细胞群体，即使在同种细胞类型中，基因模态也有很大的差异。
[0004]对于植物生长发育的研究，通常包含多个生长阶段时间点下获取的样本，然后通过对比不同条件下的细胞中基因表达谱特征推测条件改变对细胞的影响。目前植物生长发育单细胞数据分析主要通过对细胞的静态基因表达谱特征，即测量值，来推测其动态过程。目前单细胞技术测序技术在测序时裂解破坏细胞的特性，不同生长阶段的植物单细胞样本间的细胞无法配对，由于细胞间随机的异质性干扰可能掩盖了大部分于植物生长发育相关信息，而基于细胞静态特征的传统方法无法挖掘出这部分信号。这导致了大量隐藏的与植物生长发育相关信息被忽略。
[0005]由于单细胞测序技术在测序获取数据时裂解破坏细胞，为不可逆过程，通过生物学实验无法获取同一细胞在不同时期的单细胞RNA
‑
seq数据，进而不能进行同一样本不同时期的生物学数据分析比对。
[0006]因此，本领域的技术人员致力于开发一种植物生长发育动态监控的方法，通过创建细胞映射的真实
‑
虚拟细胞对，解决了无法比较同一细胞不同时期RNA
‑
seq数据的难点，并成功将其应用于植物细胞生长发育研究。

技术实现思路

[0007]为实现上述目的，本专利技术提供了一种植物生长动态监控的方法，包括以下步骤：
[0008]获取两组单细胞转录组数据，其中第一组数据X为对照组的细胞状态矩阵，第二组数据Y为变化组的细胞状态矩阵，所述对照组和所述变化组分别对应不同的发育阶段；
[0009]通过奇异值分析，构建所述第二组数据Y在所述对照组对应的所述发育阶段的虚
拟细胞状态矩阵；
[0010]得到所述第二组数据Y对应的单细胞在所述对照组对应的所述发育阶段与所述变化组对应的所述发育阶段之间的状态变化响应特征；
[0011]利用所述状态变化响应特征对所述第二组数据Y对应的所述单细胞进行聚类。
[0012]进一步地，所述发育阶段包括：韧皮部组织生长期、植物干细胞生长期、根冠生长期和木质部生长期。
[0013]进一步地，所述对照组对应所述韧皮部组织生长期，所述变化组对应所述植物干细胞生长期。
[0014]进一步地，所述状态变化响应特征是利用所述虚拟细胞状态矩阵与所述第二组数据Y的差值得到的。
[0015]进一步地，使用奇异值分解的方法，将所述不同的发育阶段的单细胞数据的协方差矩阵进行分解，得到使所述单细胞数据间协方差最大化的分量，并基于该分量，构建所述虚拟细胞状态矩阵。
[0016]进一步地，所述第二组数据Y和所述虚拟细胞状态矩阵的关系通过以下公式一表示：
[0017][0018]其中，表示所述虚拟细胞状态矩阵；V表示所述状态改变响应特征矩阵。
[0019]进一步地，所述虚拟细胞状态矩阵通过以下公式二表示：
[0020][0021]其中，P表示载荷矩阵，T表示转置矩阵，S
′
表示奇异值分解矩阵。
[0022]进一步地，所述载荷矩阵通过用NIPALS方法最大化以下公式三获得：
[0023][0024]进一步地，所述状态改变响应特征矩阵的求解过程为：
[0025]由于可以得到由于V和不相关，P
T
(P
T
P)
‑1PV≈0，可近似为P
T
(P
T
P)
‑1PU，V可计算为
[0026]V＝Y
‑
P
T
(P
T
P)
‑1PY
[0027]在给定Y时，通过P估计得到V。
[0028]进一步地，所述状态改变响应特征矩阵的求解过程为：
[0029]假设和X是来自同一特征空间的不同细胞，并使用X替换所述公式三中的第一个用C替换P以表示和X之间可能的不匹配，则所述公式三近似为以下公式四：
[0030][0031]由于X和V不相关，对于任意P和C，Cov(XC,VP)≈0，则所述公式四表示为以下公式五：
[0032][0033]在给定X和Y时，基于所述公式五计算得到P，进而得到并基于所述公式一计算得
到V。
[0034]本专利技术具有以下技术效果：本专利技术适用于植物生长发育不同阶段间样本两两比较的单细胞RNA
‑
seq研究。对于每一个细胞状态随着植物生长发育不同阶段(例如，从韧皮部组织生长期，植物干细胞生长期，根冠生长期，木质部生长期，等)的动态监控和对比。本专利技术解决了单细胞技术销毁测序细胞因此无法对植物生长发育过程中长时间跨度的每个细胞状态变化追踪的问题，使得本专利技术可以挖掘出已有方法无法发现的细胞在植物生长过程中的动态分化特征。
[0035]以下将结合附图对本专利技术的构思、具体结构及产生的技术效果作进一步说明，以充分地了解本专利技术的目的、特征和效果。
附图说明
[0036]图1是本专利技术的一个较佳实施例的方法流程图；
[0037]图2是本专利技术的一个较佳实施例的实验结果图。
具体实施方式
[0038]以下参考说明书附图介绍本专利技术的多个优选实施例，使其
技术实现思路
更加清楚和便于理解。本专利技术可以通过许多不同形式的实施例来得以体现，本专利技术的保护范围并非仅限于文中提到的实施例。
[0039]在附图中，结构相同的部件以相同数字标号表示，各处结构或功能相似的组件以相似数字标号表示。附图所示的每一组件的尺寸和厚度是任意示出的，本专利技术并没有限定每个组件的尺寸和厚度。为了使图示更清晰，附图中有些地方适本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种植物生长动态监控的方法，其特征在于，包括以下步骤：获取两组单细胞转录组数据，其中第一组数据X为对照组的细胞状态矩阵，第二组数据Y为变化组的细胞状态矩阵，所述对照组和所述变化组分别对应不同的发育阶段；通过奇异值分析，构建所述第二组数据Y在所述对照组对应的所述发育阶段的虚拟细胞状态矩阵；得到所述第二组数据Y对应的单细胞在所述对照组对应的所述发育阶段与所述变化组对应的所述发育阶段之间的状态变化响应特征；利用所述状态变化响应特征对所述第二组数据Y对应的所述单细胞进行聚类。2.如权利要求1所述的植物生长发育动态监控的方法，其特征在于，所述发育阶段包括：韧皮部组织生长期、植物干细胞生长期、根冠生长期和木质部生长期。3.如权利要求2所述的植物生长发育动态监控的方法，其特征在于，所述对照组对应所述韧皮部组织生长期，所述变化组对应所述植物干细胞生长期。4.如权利要求1所述的植物生长发育动态监控的方法，其特征在于，所述状态变化响应特征是利用所述虚拟细胞状态矩阵与所述第二组数据Y的差值得到的。5.如权利要求4所述的植物生长发育动态监控的方法，其特征在于，使用奇异值分解的方法，将所述不同的发育阶段的单细胞数据的协方差矩阵进行分解，得到使所述单细胞数据间协方差最大化的分量，并基于该分量，构建所述虚拟细胞状态矩阵。6.如权利要求5所述的植物生长发育动态监控的方法，其特征在于，所述第二组数据Y和所述虚拟细胞状态矩阵的关系通过以下公式一表示：其中，表示所述虚拟细...

【专利技术属性】
技术研发人员：郝洁，邹欣，陈可，
申请(专利权)人：上海辰山植物园，
类型：发明
国别省市：