为什么有必要在植物中做单细胞测序
如果简单地说，单细胞测序就是获取单个细胞遗传信息的测序技术，似乎没有多大的吸引力。那么为了理解这个问题，我们换个方式提问，在植物中为什么非要用单细胞测序不可？
世界上没有两片相同的叶子。对于多细胞生物来说，细胞与细胞之间是有差异的。当然了，这个差异可大可小。但是当老师的样本存在下图这样一个情况，组织中不同细胞的基因和蛋白表达不同，为了检测组织中蛋白或mRNA的表达水平，老师们可以通过Westernblot和PCR或BulkRNA-seq来实现。
但是用上述方法的话，老师们并没有办法区分这些的情况：目的蛋白只在10%的细胞中强表达，还是在50%的细胞里中等表达，还是在所有细胞中弱表达呢？
举个例子，当老师送的植物样本，一个基因（gene1）可能在分生组织中上调了，在保护组织中下调了，但是bulkRNA测序结果发现没有区别，这样的没有区别其实是个假阴性，老师可以通过单细胞测序找出这种区别，进行更为深入的研究。单细胞转录组测序（scRNA-seq）技术提供了在单细胞水平上观测基因表达的方法，这样可以更好地研究这些组织中不同的细胞类型。
那么我们为什么要在植物中做单细胞测序呢？我们可以看大这里展示的拟南芥以及水稻的根组织，他们所包含的细胞类型是有一个很大的异质性的，而同样拟南芥中，根组织和叶片组织细胞类型也是不同的，同一组织不同细胞类型间同样具有不同的基因表达模式，这些异质性我们是可以借助单细胞测序技术进行深入研究的。
植物中单细胞测序实际应用
随着单细胞测序技术的突破，单细胞测序的时代已然到来。其技术已经大量发表在CELL、Nature、Science等高质量期刊中，其中大多数集中在人和动物方面。近年来，单细胞测序在植物方向上的应用也在逐步扩大，涉及的物种包括拟南芥、水稻，以及玉米，本文就带大家一起看一下由欧易生物助力发表的植物相关单细胞测序文章。
首先是河南大学孙旭武教授课题组在MolecularPlant杂志发表的一篇研究拟南芥气孔谱系细胞发育进程的文章。
MolecularPlant(IF21.)；
PubDate:-06-24；
材料：5日龄拟南芥幼苗子叶中分离的原生质体；
技术：10×GenomicsscRNA-seq、单细胞转录组测序
这篇文章思路首先是一个常规的细胞类型的鉴定，通过文献调研作者最终定义了9种细胞类型，以及另外两种不能使用已知的标记基因来确定的细胞类型。热图和小提琴图分别对top10的基因进行展示。
由于植物marker基因鉴定尚不成熟，所以作者为了验证top10marker基因的可靠性，采取了以下2种验证，第一种方法通过研究top10的分子功能来推测细胞功能。作者构建了bag6和bzip6的突变株，发现bag6突变株表面气孔指数增加，推测BAG6参与调控气孔的发育，进而推测Cluster9是一群参与调控气孔的发育的细胞。
以下是第二种验证top10marker基因的可靠性的方法:
年，Adrian那篇文章是用流式把拟南芥子叶分成了几类细胞，然后回收这些细胞去做测序，公布了测序数据集给大家用。这篇文章就是把Figure2A中top10基因挑出来，看看这些基因在人家的哪个细胞类型的数据集中表达了。简单来说是把Figure2A中top10的marker基因，映射到这两篇文章中，看看这些基因在这两篇文章中是不是也在同一个细胞类型在表达。
世上本无marker基因，top10基因用的人多了，也就变成了marker基因。举个例子，PC（扁平表皮细胞）他没有特定的标志基因，这篇文章他把鉴定到的PC高表达基因和人家数据集比对，比对到了40%，也进一步可以说明细胞类型鉴定的可靠程度。
既然作者都用这么多方法证明了细胞类型鉴定的可靠性了，那么就可以放心开展接下来的研究了。为了发现参与调节气孔谱系细胞早期发育的潜在TF，作者统计了在不同细胞类型中显示高表达的TFs。unknown(u.k.)_1中TFs数量最多，MPC_5中TFs数量最少。不知道大家记不记得，通过文献调研作者最终定义了9种细胞类型，以及另外两种不能使用已知的标记基因来确定的细胞类型。
Cluster9这群不为人知的细胞前面作者已经推断好了，他是参与调控气孔的发育的细胞。这里发现Cluster1是一群参与气孔早期发育的基因。那话不多说我们接着看下去。
MMC，LM，EM和GMC均是气孔的发育早期的重要细胞。因此作者着重研究MMC，LM，EM和GMC四种细胞共表达基因的转录因子调控网络。分析结果表明ATML1,BPC1,BPC6,SCRM,PIF5,andWRKY33可能是调控MMCs、EMs、LMs和GMCs中靶基因的表达的核心转录因子。
通过文献调研，发现BPC蛋白和WRKY33是参与调节幼苗生长发育和胁迫反应的重要TFs。而BPC1，BPC2，BPC4，BPC6和WRKY33的表达在MMC，LM，EM和GMC中表达也相对很高。
综合这些因素，于是作者又研究了BPC蛋白、WRKY33和这次鉴定到的共表达转录因子和已知关于这四种细胞的转录因子之间的调控关系。结果发现BPC1和BPC6可能通过这些相互作用调控早期气孔谱系细胞的发育。
为了研究已鉴定的TFs对气孔谱系细胞发育的影响，作者分析了它们的表达特征图。发现BPC1、BPC2、BPC4和BPC6的表达主要在MMCs、EMs和GMCs中增强，而WRKY33在所有细胞类型中均可检测到表达。于是作者构建了WRKY33突变株，发现WRKY33突变株的M和GMC细胞比例上升，但GC细胞比例下降，表明BPCs蛋白可能参与介导M到GMC的发育过程，从而促进气孔谱系早期发育。作者构建了BPC突变株，也发现了同样的规律。
最后为了分析气孔分化过程中的发育轨迹，作者使用Monocle2软件对scRNA-seq单细胞转录组测序数据进行了拟时序分析，拟时序路径共有三个分支。总体来看，从MMC到GC可以看到气孔谱系细胞的不同发育过程。
这篇文章作者对拟南芥气孔谱系细胞发育过程进行了单细胞转录组测序分析，鉴定了气孔发育过程中各细胞类群显著表达的markergene，并分析了其所参与的GO功能。为了分析气孔细胞早期发育阶段的基因表达调控机制，构建了气孔分生母细胞中基因表达的转录因子调控网络。进一步通过分析所鉴定的标记基因和转录因子的拟南芥突变体的气孔发育表型，确定了这些基因在调控气孔发育中所扮演的重要角色。这项研究为利用单细胞RNA-测序技术在鉴定新的标记基因，以及解析气孔谱系细胞发育的调控机制提供了新的研究思路和技术体系。
接下来这篇文章是福建农林大学等单位的研究人员在水稻叶片原生质体中构建了TF（OsNAC78）的差异表达系统，对构建好的原生质体进行scRNA-seq。
FrontiersinPlantScience(IF6.)；
PubDate:-11-16；
材料：MH63水稻幼苗原生质体；
技术：10×GenomicsscRNA-seq单细胞测序
这篇文章老师重点是利用单细胞测序来筛选转录因子的一个靶基因。
老师为了研究OsNAC78这个转录因子的靶基因，将OsNAC78过表达到MH63水稻幼苗中，对照组为野生型MH63水稻幼苗。
通过对OsNAC78高表达的细胞进行细胞亚型分析，发现cluster1和cluter4一个是表达OsNAC78最多的cluster一个是表达OsNAC78最少的cluster。
于是作者针对cluter1和4做了一个差异分析，将cluster4中上调的19个候选基因基因单出拎了出来。
随后作者通过过表达OsNAC78，看看有哪些基因可以随着OsNAC78的过表达也跟着上调。结果筛选出了两个基因。然后我们通过文献调研等手段构建了一些Os01g和Os01g0900启动子的报告质粒，然后酵母单杂来筛OsNAC78会和哪个启动子结合，找到启动子后。作者又做了EMSA和荧光素酶报告实验，验证了一下Os01g和Os01g0900确实是OsNAC78的靶基因。
年2月26日，河南大学孙旭武教授团队在ThePlantJournal杂志发表了一篇通过对拟南芥幼苗子叶的单细胞测序，揭示了子叶叶脉的早期发育动态的文章。
ThePlantJournal(IF7.)；
PubDate:-02-26；
材料：拟南芥3天龄子叶原生质体；
技术：10×GenomicsscRNA-seq单细胞转录组测序
这篇文章作者重点
