2017年发表在Nature Methods杂志上的 SCENIC 算法,利用单细胞RNA-seq数据,同时进行基因调控网络重建和细胞状态鉴定,应用于肿瘤和小鼠大脑单细胞图谱数据,提出并证明了 顺式调控网络分析能够用于指导转录因子和细胞状态的鉴定 。SCENIC通过使用生物学驱动的features自动清除肿瘤样本特异性等批次效应。
我们要深刻了解SCENIC并应用在我们的项目上就要对 基因调控网络(GRN, gene regulatory network) 的背景有一些认识。细胞的转录状态来自潜在的基因调控网络,GRN由数量有限的转录因子(TFs)和辅因子相互调节及调节下游靶基因构成。
单细胞转录组分析的最新进展为高分辨率识别转录状态和状态之间的过渡提供了令人兴奋的机会。例如,在分化过程中针对单细胞RNA-seq进行优化的统计和生物信息学方法带来了新的生物学见解。但是, 稳定细胞状态下的特异性和强健的(robust)GRN是否能够被确定还尚不清楚 。
鉴于在单细胞水平来做这件事情,这可能确实具有挑战性。由于转录爆发和其他来源的基因表达的随机变化,基因表达可能会与TF输入的动力学部分断开(不相关)。目前已经有一些从单细胞RNA-seq数据推断共表达网络的方法发表,但是,这些方法未使用调控序列(regulatory sequence)分析来预测TF与靶基因之间的相互作用。
作者认为,将顺式调控序列与单细胞基因表达连接可以克服缺失(基因表达丰度检测不到)和技术多样性,从而优化细胞状态的发现和表征。鉴于此,作者开发了单细胞调控网络推断和聚类(SCENIC)来绘制GRN,通过评估每个细胞中GRN的活性,来识别稳定的细胞状态。
SCENIC workflow 包含3个主要步骤:
使用RcisTarget是SCENIC不同于大多共表达算法的重要区别 。由于GENIE3模块仅基于共表达,因此结果可能包含许多误报和间接target,为了鉴定推断的直接结合的靶标基因,使用RcisTarget对每个共表达模块进行顺式调控基序(motif)分析。 仅保留具有正确基因上游调节子且显着富集TF motif的模块,并对它们进行修剪以除去缺乏基序支持的间接靶标,这些处理后的模块才称为regulon。
作为SCENIC的一部分,作者开发了AUCell算法来对每个细胞中的每个regulon的活性进行评分。
对于一个给定的regulon,通过比较所有细胞间的AUCell打分值,我们可以识别哪些细胞具有更显著高的regulon活性。通过卡阈值得到的二元活性矩阵使矩阵维数减少(可理解为只有 0|1,on|off),对于下游分析很有用。 例如,基于regulon二元活性矩阵的聚类,可以根据某个调控子网络(regulon)的活性来识别细胞群类型和细胞状态。由于regulon是整体评分的,而不是使用单个基因的表达,因此这种方法对于个别基因的dropouts很有效。
作者对SCENIC进行了 准确性 和 鲁棒性 评估。
作者拿成年小鼠大脑的已知细胞类型的scRNA-seq数据集来测试SCENIC,来评估其准确性。
分析了1,046个初始共表达模块中的151个Regulons,这些Regulons包含对应TF的显着富集的motif(占初始TFs的7%)。每个细胞的Regulon活性评分确定了预期的细胞类型( 上图d,e )以及每种细胞类型的潜在主调控因子的列表(例如下图中的小胶质细胞网络),细胞聚类的结果比一些专用的单细胞聚类方法( SC3 )要准确。
作者用以下方式分析了小鼠大脑单细胞RNA-seq图谱数据:
SCENIC鉴定出了仅由少量细胞代表的细胞类型(例如,来自小胶质细胞,星形胶质细胞或中间神经元的2至6个细胞;下图)。此外,预测的TFs与细胞类型的关联也与先前已知的一致,并且此准确性优于标准分析流程。
为了验证鉴定小鼠interneurons的Dlx1/2调控网络,作者分析了人脑的sNuc-Seq(Single nuclei RNA-Seq)数据集。
在人类脑sNuc-Seq数据集上,SCENIC也鉴定出由DLX1/2强烈驱动的interneurons细胞群,该群具有与小鼠相同的motif,并且识别出一组保守的靶标,包括DLX1本身。
接下来,作者将这种跨物种分析扩展到其他细胞类型。 与基于归一化表达的标准聚类(产生强大的物种驱动聚类)不同,SCENIC分析有效地按细胞类型对细胞进行了分组(下图)。 这表明网络活性的评分是可靠的(robust),可以用来克服批处理或测序方法等技术影响。
由于肿瘤特异性突变和复杂的基因组畸变,癌细胞状态的鉴定比正常细胞状态更具挑战性。一般的标准聚类会通过表达矩阵将细胞按其肿瘤起源分组(各个样品聚成一类),但SCENIC的结果揭示了不同的图景。
以下是作者将SCENIC用在 少突胶质细胞瘤 (来自 6个肿瘤的4,043个细胞 )和 黑素瘤 (来自 14个样本的1,252个细胞 )的scRNA-seq数据集上的结果。
对于少突胶质细胞瘤(oligodendroglioma),在所以肿瘤细胞中鉴定出三种癌细胞状态(上图c-e),每种状态均由预期的TF驱动,包括:
此外,作者将扩散映射(Diffusion Maps)应用于二元SCENIC矩阵(上上图)重建了从stem-like到oligodendrocyte-like 和 astrocyte-like分支的分化轨迹。值得注意的是,与正常oligodendrocyte分化相比,此路径代表不同的“轨迹”。
在黑素瘤(melanoma)数据上观察到了类似的肿瘤效应校正,其中SCENIC识别了跨肿瘤的细胞群。包括一群与少突胶质细胞瘤中类似的TF驱动的周期细胞(例如,E2F1/2/8 和 MYBL2)
黑色素瘤细胞大致分为两组:
SCENIC在MITFlow state的细胞群下鉴定到了两个新的TFs:
NFATC2是JNK/MAPK途径中的转录阻遏物,参与黑色素瘤去分化和免疫逃逸;
NFIB与毛囊和黑素细胞干细胞的干细胞行为有关,它在小细胞肺癌转移过程中起着重要作用。
为了进一步探索NFATC2和NFIB在MITFlow状态下的潜在作用,作者对25个具有不同肿瘤进展的黑色素瘤标本进行了免疫组织化学分析。
作者发现NFIB和NFATC2在前哨淋巴结(sentinel lymph nodes)中表达最高,这与ZEB1表达共定位,这与ZEB1表达共定位,这表明这些markers的表达与最早的转移事件之间存在关系。
在A375黑色素瘤细胞系,NFATC2和NFIB的表达很高,当使用siRNA敲除NFATC2时,作者发现NFATC2调节子中的基因被显着上调。这与先前证实的NFATC2作为阻遏物是一致的。此外参与细胞粘附、细胞外基质和几个先前发表的代表黑色素瘤浸润状态调控的基因也被上调。这表明NFATC2确实可能在疾病进展中起重要作用。作为黑色素瘤regulons的第二次验证,作者使用ChIP-seq数据确定了MITF和STAT的预测靶标。
随着单细胞数据集大小的增加,作者建议两种补充方法来扩展网络推断:
SCENIC是一种普遍适用的分析scRNA-seq数据的方法,利用TF和顺式调控序列来指导细胞状态的发现。文章的结果表明,GRNs是确定细胞状态的可靠方法,并且scRNA-seq数据非常适合跟踪基因调控过程,在基因调控过程中特定组合的TFs驱动细胞特异性的转录组。
作者开发的pyscenic使用教程请参考: https://pyscenicreadthedocsio/en/latest/ 或者查看 pyscenic github
最近发表的 单细胞转录组测序绘制小鼠内皮细胞精细图谱 分析研究转录因子在不同ECs中的调控作用,发现了转录因子调控网络呈现了组织特异性。
RNA聚合酶进入DNA非编码区的酶切位点,解旋DNA使成为单链,核糖核苷酸由碱基互补配对法则形成RNA链(信使RNA)。
RNA的合成即DNA的转录:
DNA复制最重要的特征是半保留复制,即DNA复制时,DNA双链中的互补碱基之间的氢键断裂,解为两条单链,各以一条单链为模板,按碱基互补原则合成新的互补链,新合成的两个DNA分子和亲代DNA分子是完全一致。
在子代DNA分子中一条单链来自亲代,另一条是新合成的,这种复制方式称半保留复制(semiconservative replication)。遗传信息按这种方式忠实地从亲代传给子代。
目前有关复制的知识主要来自于原核生物实验,所以主要以原核生物来叙述。
Mpower留学贷款的审批有一套独创的评估体系,和许多贷款机构不同的是,Mpower不需要抵押和担保,而是根据申请人自身的学术背景、申请院校专业,以及未来职业规划进行审核。这为许多因经济条件有限无法留学的同学们提供了非常难得的机会,但同时也对申请人的发展潜力提出了要求,毕竟申请人日后拥有良好的职业发展,才会使Mpower的贷款风险降到最低。STEM和MBA专业因自身学术含金量高、职业发展前景良好,确实相对更容易申请到贷款,但除此之外,其他专业也不是完全没有可能。终究一句话,申请人自身的综合发展潜力是最重要的,如果你对自己有信心,不妨就去申请看看,留学贷款没有你想象得那么“高不可攀”。
是olige dT,成分上来说,属于DNA。
你需要把RNA逆转录成cDNA,olige dT是和mRNA的polyA尾巴配对的。
当然,反转录mRNA成DNA不一定用olige dT,还有random primer,stem loop 等。
看看这个吧,就很清晰啦!
http://imagebaiducom/ict=503316480&z=0&tn=baiduimagedetail&word=t%2DRNA&in=23033&cl=2&cm=1&sc=0&lm=-1&pn=57&rn=1&di=4220140201&ln=734&fr=&ic=0&s=0&se=1
这个链接氨基酸的结构,叫做氨基酸臂(amino arm): 也称之接纳茎(acceptor stem)。tRNA分子中靠近3ˊ端的核苷酸序列和5ˊ端的序列碱基配对,形成的可接收氨基酸的臂(茎)。
欢迎分享,转载请注明来源:品搜搜测评网