0 引言 Introduction
骨骼系统作为人体八大系统之一，在人体中发挥着重要的作用，也从侧面反映出骨骼系统的复杂性，这给研究者带来了许多挑战；另外随着交通工具的快速发展、生活方式的改变及环境变化，骨骼系统相关性疾病的发生率也越来越高，例如：骨折、肌腱及韧带损伤、骨肿瘤、骨关节炎等。免疫学、流式细胞术及qRT-PCR 等传统研究方法因其经济、快速、简单易操作等优点仍被广泛用于探究骨骼系统的发生发育及相关疾病的病理生理、发病机制、发生发展及转归。有研究报道骨髓基质细胞可以促进肌腱修复、加快肌腱与骨的植入愈合并证明间充质干细胞在长骨及肌腱缺损修复中具有重要作用[1-3]；大量研究已发现TP53、PTEN 等诸多基因、蛋白及相关信号通路在骨肉瘤中具有重要作用[4-9]；通过RNA 测序发现软骨相关基因Sox9、Col11a2、Acan 和其他新转录因子在骨关节炎中的表达发生明显改变等等[10]，但上述传统方法无法很好地解释创伤后组织是如何进行修复及愈合？骨肿瘤组织中是否存在细胞异质性？骨骼、软骨、韧带等骨骼组织的起源和形成机制仍不清楚，无法完全了解参与骨骼组织发育的间充质细胞亚群的起源和潜在功能，骨髓基质及间充质细胞在骨骼组织创伤后修复中的作用和潜在功能，以及骨关节性疾病中软骨细胞类型组成和促进疾病进展的细胞异质性等相关问题仍不清楚，因此需要一种新的研究方法来进一步研究。
近些年来，以单细胞RNA 测序技术为研究的方法正在得到蓬勃发展，目前单细胞RNA 测序技术已广泛使用于肿瘤、免疫学、神经生物学等领域[11-14]，并且已逐渐开始用于骨骼系统生物学相关研究。相当多的研究通过单细胞RNA 测序在骨骼系统中发现了以前尚未报道过的新细胞类型及细胞标志物，例如：在肌腱[15]、肌肉的单核细胞及骨关节软骨中都发现了新的细胞亚群[16-17]，这将有助于对该疾病有更深入的认识以及对疾病的精准治疗提供新见解，对未来再生医疗、精准治疗及个体化治疗具有重要意义。文章旨在综述单细胞RNA 测序技术的基本原理及其在骨骼系统生物学相关研究中的应用。
1 资料和方法 Data and methods
1.1 资料来源第一作者通过 PubMed、CNKI、万方数据库检索文献，英文检索词为“single cell RNA sequencing technology，bone/skeletal，tendons，ligaments，bone/skeletal diseases”；中文检索词为“单细胞RNA 测序技术，骨骼，韧带，肌腱，骨关节疾病”。
1.2 文献入选标准
1.2.1 纳入标准 与骨骼、肌肉、韧带等骨组织的发生发育及相关疾病、单细胞RNA 测序密切相关的文献。
1.2.2 排除标准 ①重复多次的文章；②与此次综述内容相关性低的文献；③对此文参考价值不高的文章及会议文献。
1.3 数据的提取共检索出2 907 篇相关文献，排除重复性文章143 篇后，初步通过阅读标题、摘要排除与研究目的相关性差的文献并进行二次筛选，最后共纳入文献 75 篇。文献筛选流程见图 1。
图1｜文献筛选流程图
2 结果 Results
2.1 单细胞RNA 测序原理自TANG 等[18]于2009 年首次使用单细胞RNA 测序分析基因表达谱以来，目前该技术已广泛用于肿瘤、免疫学、神经生物学等领域[11-14]，完整的单细胞RNA测序实验包括单细胞分离、文库制备和数据分析3 个部分。
2.1.1 单细胞分离 单细胞分离是单细胞RNA 测序技术的关键步骤，决定了研究的成败。目前有多种方法可进行单细胞分离，包括显微操作法、微流控法、荧光激活细胞分选法(FACS)，激光捕获显微切割法(LCM)等，其中荧光激活细胞分选法和微流控法是最常用的方法。荧光激活细胞分选法(比如Smart-seq、CEL-seq/CEL-seq2 等分离技术)是基于细胞表面特异性分子标记或利用预先用荧光标记的细胞光散射特性并通过流式细胞术来进行单细胞分离，然后将得到的单细胞分配到96 或384 孔板中进行RNA 测序并建立文库；基于荧光激活细胞分选法分类的优势之一是能够在特异性和非特异性分离之间进行选择，具有高准确度[19]，但其一次检测到的细胞数量有限。此外该方法也有几个明显的缺点：第一，其需要大量的细胞悬浮液，这降低了获得低丰度细胞的速度；第二，检测细胞数量有限导致每个细胞的测序成本明显升高。相反，最近开发的微流体技术一次可以捕获并处理大量细胞[20-21]；该技术使用独特的标记磁珠与单个细胞结合形成单个液滴；当细胞裂解后，磁珠上的标记与释放的mRNA 杂交，并反转录成带有独特标记的cDNA，然后构建测序文库。该方法具有样品消耗少、分析速度快、自动化程度高等显着优势，而且还可以将反应、预处理、检测等过程集成到微流体系统集中完成，由于其低成本和高通量，目前已成为主流的单细胞分离方法。
2.1.2 文库制备 单细胞RNA 测序文库制备首先需要将单个细胞中的RNA 通过PCR 反转录成cDNA，由于单个细胞中的总RNA 含量非常少，仅约10 pg[22]，故直接对单个细胞的RNA 进行测序不太可能[23]，因此首先需要对RNA 进行扩增，然后再进行测序、建立文库。既往进行全基因组扩增的方法主要包括3 种：PCR 技术、多重置换扩增技术及多次退火环状循环扩增技术(表1)。
表1｜既往进行全基因组扩增的主要方法
随着测序技术的发展，最初使用全转录组测序对单细胞RNA 进行扩增以获得基因表达谱[19-20]，但其鉴定的基因及未知的剪接点有限；TANG 等[18]改进了单细胞全转录组扩增方法，使用新一代测序技术对单细胞全转录组进行扩增，该技术使用寡脱氧胸腺苷酸作为引物来进行反转录，然后通过PCR 技术选择性扩增聚腺苷酸化mRNA。SMART-seq 技术则使用莫洛尼鼠类白血病病毒(moloney murine leukemia virus，MMLV)反转录酶来构建全长cDNA 文库[26]，其最大的优点是能将单细胞中的RNA 扩增为全长cDNA，从而发现可变剪切的外显子，灵敏度低是其最大的缺点；而改良的SMART-seq2技术可以同时对数百个细胞的转录组进行测序，明显提高了测序的效率及灵敏度。近几年来被广泛使用的基于液滴的微流控技术可以对成千上万个细胞同时进行测序，并且可以通过消除扩增偏差来提高全基因转录组测序的准确性；独特的分子标识符标记技术的应用可以对每个基因的来源进行追踪，以此来区分细胞类型及进行复杂的分析。
2.1.3 数据分析 最后，将测序得到的数据进行统计分析，其中无监督聚类分析是最常使用的分析方法之一，可以根据细胞基因表达的不同来揭示细胞间的异质性并对其进行分类，以寻找出感兴趣的基因及细胞进行研究，对探索相关疾病的发生发展、预防及指导临床治疗具有重要意义；此外，可以根据特征基因的表达来推断出细胞种类及其来源，甚至发现新的细胞类型，例如：PARK 等[27]使用单细胞RNA 测序技术对肾脏组织进行测序发现3 个新的细胞种类；还可以根据基因的表达对细胞的发育轨迹进行追踪、筛选差异基因以及进行拟时序分析等[28]。
2.2 单细胞RNA 测序在骨骼系统中的应用目前单细胞RNA测序在肿瘤、免疫学、神经生物学等领域研究中取得较大的进展[11-14]，并逐渐开始用于骨骼系统生物学相关研究。
2.2.1 单细胞RNA 测序在骨骼干细胞中的应用 骨骼组织正常的发生发育及创伤后的自我修复对于维持人体的完整及肢体的功能具有重要作用，但对参与其中的具体细胞类型及是否有骨骼干细胞的参与目前尚不清楚。CHAN 及其同事[29]对围产期小鼠生长板及周围细胞进行单细胞RNA 测序，发现被分离出来的AlphaV(CD51)+Thy1(CD90)-CD105-CD200+细胞具有多能骨骼干细胞的特性，能进行自我更新，并将其称为小鼠骨骼干细胞(mSSCs)，随后他们通过谱系追踪实验发现骨骼干细胞确实存在于生长板的静息区[30]；最近他们通过对人类胎儿生长板进行单细胞RNA 测序发现在生长板肥大前区和肥大区中的细胞表达小鼠骨骼干细胞特异基因的人类直系同源基因，故将该细胞称为人类骨骼干细胞(hSSCs)，并推断人类骨骼干细胞也存在于生长板中[31]。通过单细胞RNA 测序证实了骨骼干细胞的存在，发现其在生长板及骨骼完整性的维持中起着重要作用，这将对治疗创伤引起的骨骼组织损伤提供新的临床见解，但骨骼干细胞是通过何种细胞分子信号或通路参与骨骼组织的维持及修复尚未能阐明，仍需进一步的研究。
生长板在骨的生长与维持中起着重要作用，但存在于靠近软骨膜生长板边缘的“交界性”软骨细胞在体内的具体作用及命运却不得而知，有研究报道称在软骨膜附近的特殊微环境中，生长板周围的“边缘”软骨细胞可以在培养条件下产生成骨样细胞[30]，为了探究生长板周围“边缘”软骨细胞的性质及潜在的功能，MIZUHASHI 等[32]对新生小鼠生长板进行单细胞RNA 测序发现了一种新的软骨细胞群，将其称为“交界性”软骨细胞，后续的血统追踪实验证实其表现为短暂的间充质前体细胞，发现“交界性”软骨细胞早期可以在干骺端分化产生成骨细胞及骨髓基质细胞，到成年期后其细胞数量明显减少，表现为短暂的分化过程；而TAKAHASHI 等[33]和MIZUHASHI 等[32]利用单细胞RNA 测序技术也分别标记和描述了牙囊和长骨生长板中的间充质细胞群，这有助于理解骨的生长发育及停止的全过程，以及进一步了解这些细胞群的起源、命运及细胞间的异质性。
骨膜是附着在骨表面的坚固结缔组织包膜，有丰富的血管与神经，对骨的营养及新陈代谢具有重要意义，另外其膜内成骨的能力对骨折迅速愈合尤为重要，虽然对骨膜进行了大量的研究，证实了其基本的细胞类型及相关功能，但是否存在其他细胞亚型及其作用尚无法得知。单细胞RNA 测序可以从细胞分子层面去解决上述问题。DEBNATH 等[34]对骨膜组织进行单细胞RNA 测序，发现了一种被组织蛋白酶K(CTSK)标记的新的骨膜干细胞群体，该细胞群与存在于骨膜腔的间充质干细胞不同[35-36]，它们不表达有造血功能的间充质细胞相关标记物，如LEPR 或CD146，且只在基线时介导膜内骨化。他们将骨膜细胞聚类为4 组：第1 组高表达Sox9 和Col2a1；第2 组表达成骨细胞标记物，如Bglap 和Alpl；第3组表达Ly6a；第4 组以高表达α-平滑肌肌动蛋白(Acta2)为特征。他们发现被鉴定为骨膜干细胞的细胞几乎都高表达软骨细胞的特异性转录物Sox9 和Col2a1[37]。另外有研究也发现Sox9 和Col2a1 在长骨及颅骨的骨膜细胞中高表达[38-39]，并且Sox9 缺失将严重影响骨的形成[40]，因此DEBNATH 等[34]推测Sox9 和Col2a1 可以作为骨膜干细胞的标记物；同时第2 组成骨细胞标记物(如Bglap 和Alpl)的高表达也确定了骨膜干细胞可以产生成骨细胞，并发现了Ifitm5 基因的高表达，这表明Ifitm5 可能成为单细胞RNA 测序研究中用于鉴定成骨细胞的重要标志物；另外他们还推断Acta2+细胞可能是骨膜干细胞与成熟成骨细胞过渡的中间祖细胞，但需进一步的体外移植研究来验证。因此可以通过单细胞RNA 测序去鉴定新的细胞类型及细胞标志物，从而分析该细胞类型的潜在功能。
2.2.2 单细胞RNA 测序在肢体发育中的应用 哺乳动物肢体的正常发育，涉及到一系列复杂且相互协调的分子过程，例如以高度协同及时间依赖的方式进行各种基质蛋白和信号分子的合成与降解，虽然在肢体发育的某些方面已经有相关阐述，但参与其中的相关细胞亚群的发育轨迹、在肢体发育过程中的潜在分化能力及相关功能暂不清楚。单细胞RNA 测序技术可以采用拟时序分析来阐明在肢体发育过程中细胞亚群与时间、空间的变化关系，更加全面的了解肢体细胞的发育轨迹及细胞成熟的连续性。KELLY 等[41]对胚胎小鼠后肢行单细胞RNA 测序，发现在不同时间点包含不同的细胞类型，包括：软骨前体(Hoxd13 和Col2a1)、骨/肌腱前体(Tbx13和Col1a1)、皮肤(Krt14)、脉管系统(Cdh5)、血液(Lyz2)、生长板(Col10a1)和肌肉(Myog)，并且这些细胞类型在不同时间点所占百分比也不尽相同，还确定了小鼠后肢中发育的5个主要分支。KELLY 等在胚胎的第13.5 和15.5 天发现了表达Col10a1 的细胞簇，但在胚胎第18.5 天时未发现，可能是由于细胞捕获受限而导致的细胞类型脱落，也可能是肥大软骨细胞在最成熟的时间点更类似于骨骼或软骨细胞；另外他们并未发现与关节形成相关的重要标记基因Gdf5，有可能是由于细胞捕获的局限性所致，但也有研究表明Gdf5 本身仅在少量细胞中表达[42]。虽然单细胞RNA 测序在单细胞分离过程中可能存在捕获细胞类型不全的现象，但可以为已捕获的不同类型细胞之间的相互作用提供新的见解，以透彻了解组织发育过程，促进组织修复或最终促进组织或肢体再生。
2.2.3 单细胞RNA 测序在肌腱和韧带中的应用 肌腱及韧带的正常发育对关节活动的维持具有重要意义，但对参与肌腱及韧带的细胞类型组成及创伤后组织修复的内在细胞异质性仍不清楚，为了探讨及解决上述相关问题，SWANSON 等[15]对小鼠肌腱进行单细胞RNA 测序，他们共鉴定出11 种不同的细胞类型，其中包括以前尚未报道过的3 种成纤维细胞类型，将其分别称为肌腱成纤维细胞1 和2，以及连接成纤维细胞；令人惊奇的是他们发现这3 种成纤维细胞都高表达Ⅰ型胶原蛋白，另外连接成纤维细胞还高表达Col22a1，并且支持表达跨膜蛋白CD146 的周细胞是成人肌腱成纤维细胞祖细胞群的观点[43]。肌腱组织细胞亚群之间和内部显著的异质性，可能有助于加深对肌腱细胞外基质的组装和维持的理解；另外成纤维细胞专一性地产生肌腱细胞外基质的不同成分，这可能给肌腱病变的治疗提供新靶点。此外，多项研究发现Scx、Tnmd 和Mkx 是肌腱发生发育过程中重要的标记[44-45]，但SWANSON 等[15]和KELLY 等[41]的研究结果却与之相反，他们通过单细胞RNA 测序和RNA ISH 检测并未发现Scx、Tnmd 和Mkx 在肌腱细胞中广泛表达；SWANSON 等[15]与DISSER 等[46]在髌骨肌腱中未发现成骨细胞标记物骨唾液蛋白(Bsp)、Msx2 或osterix，这与之前的报道称在髌骨肌腱中发现了大量成骨细胞的结果相悖[47]。故单细胞RNA 测序分析可以从细胞分子层面去验证或者推翻先前研究的结果，也可以寻找出新的细胞分子标志物来发现细胞间的异质性，并进一步阐明其相关功能。
2.2.4 单细胞RNA 测序在骨骼肌中的应用 骨骼肌是一种复杂的异质组织，由多种细胞类型组成，除了参与运动外，还影响人体的免疫、代谢和认知等方面，了解骨骼肌的细胞组成类型及其潜在的功能将有助于加深肌肉对运动的各种影响以及肌肉相关疾病的认识。多项研究表明Mesp1 细胞在心脏及骨骼肌的发生发育中起着重要作用[48-50]；而有趣的是，HAREL 等[49]和CHAN 等[50]发现来源于Mesp1+祖细胞的肌肉干细胞主要存在于颅面部的骨骼肌中，几乎不存在于躯干及后肢，这打破了传统的认知观念，为了深入了解头面部和躯干/四肢骨骼肌肉发育过程的差异，PENALOZA等[51]对Mesp1 诱导的骨骼肌细胞进行单细胞RNA 测序，对骨骼肌源性祖细胞标志物PDGFRα 和VCAM 的动态观察发现Mesp1 诱导的骨骼肌细胞分化是从PDGFRα+VCAM-(Tcf15+) 到PDGFRα+VCAM+(Myf5+Tbx1+Pitx2+)，然 后 再 到PDGFRα-VCAM+(Myod1+Myog+)的过程，同时他们还观察到在Mesp1 诱导的骨骼肌细胞分化时出现了不同的细胞亚群和非肌源性细胞，包括一些未分化或分化早期的中胚层及外胚层细胞，以及表达于胚胎/胎儿成肌细胞中的蛋白Acta2 和Tagln，这说明Mesp1 诱导的骨骼肌成肌祖细胞可能是胚胎性质的。通过单细胞RNA 测序可以重建不同细胞亚群的单细胞轨迹来预测细胞的命运，这将有助于破译组织发生的分子机制，从而产生这种独特的细胞类型为未来的再生医学提供见解。
骨骼肌组织中不同细胞亚群内的细胞组成是否会随着年龄、性别、运动等变化尚不清楚，单细胞RNA 测序可以通过对不同年龄人群及各种干预后的骨骼组织细胞进行深入分析，来剖析其细胞类型比例变化并推测其潜在的功能。RUBENSTEIN 等[16]利用单细胞RNA 测序技术在人的股外侧肌中发现了11 种人类骨骼肌单核细胞亚群，包括两种新的成纤维脂肪祖细胞(FAP)亚群：纤维蛋白1+(FBN1+)和lumican+(LUM+)FAP 亚型，并且揭示了特定单核细胞类型比例的显著变化与年龄、性别、急性运动和训练有关，这将有助于对运动的各种影响以及肌肉疾病的病理生理学的 了 解；相 反，BARRUET 等[52]和DE MICHELI 等[53]的 研究发现肌纤维类型组成和基因表达并不会根据供者的年龄而改变，这一冲突可能需要进一步的研究去探索。另外，RUBENSTEIN 等在小鼠的骨骼肌中也发现了同样的成纤维脂肪祖细胞亚型，并且在小鼠的FBN1+FAP 亚型中发现了在人类中没有的特异性标记物TEK。此外，小鼠体内发现的具有软骨生成/成骨潜能的Tie2+FAPs，可能是骨骼肌异位骨化的主要促因，而Tie2 蛋白的基因前体TEK 在小鼠的FBN1+FAP 亚型中表达，因此RUBENSTEIN 等推测FBN1+FAP 亚型可能在人类的骨骼肌异位骨化起作用。BARRUET 等[52]在人骨骼肌卫星细胞中发现了一种新的亚群：CAV1+卫星细胞。CAV1 先前已被证明由小鼠卫星细胞表达[54-55]，BARRUET 等通过对CAV1+卫星细胞中静止基因标记物的表达发现CAV1+卫星细胞与人类卫星细胞的静止和细胞黏附特性有关，可以抵抗细胞的活化、维持卫星细胞生态位及表型。
2.2.5 单细胞RNA 测序在创伤后骨组织修复中的应用 骨折后促进骨形成及愈合过程中细胞内分子改变及细胞间相互作用的具体机制尚不清楚，在组织修复过程中需要细胞增殖与凋亡达到一个新的平衡，有研究表明临床中发现的骨折不愈合很大程度上是由于缺乏足够的祖细胞所造成的[56]；而骨形态发生蛋白2 能够在体外扩增骨骼干细胞并在软组织中诱导骨骼干细胞再生，是骨形成和骨折愈合所必需的[57-58]。有研究报道称人骨膜来源的细胞(hPDC)可以通过激活骨形态发生蛋白信号来增强异位体内骨形成[59]，而激活的骨形态发生蛋白可进一步增强骨愈合能力[60]，并且人骨膜来源的细胞的成骨潜能可能与细胞内骨形态发生蛋白1 型和2 型受体的mRNA 转录升高有关[61]。但骨形态发生蛋白在骨折愈合中的具体分子机制、参与的信号通路及激活下游何种靶分子目前尚不清楚，单细胞RNA 测序可以在单个细胞水平获得基因组全转录信息，以此来探究参与骨折愈合过程中的相关通路及分子机制。BOLANDER 等[62]对经无血清处理的人骨膜来源细胞进行单细胞RNA 测序及体外验证实验发现，激活的骨形态发生蛋白可进一步引起下游骨软骨祖细胞、CD34+祖细胞、软骨细胞和BMPR2+骨软骨祖细胞簇中的骨形态发生蛋白2型受体表达明显上调，然后进一步诱导磷酸核糖焦磷酸合成酶1 基因表达、增加谷氨酰胺的消耗及谷氨酰胺酶的表达来增强骨形成及骨愈合能力，这将为临床中治疗骨折不愈合或延迟愈合提供新的靶点。
2.2.6 单细胞RNA 测序在骨关节性疾病中的应用 骨性关节炎、类风湿性关节炎及银屑病关节炎等是常见的骨关节性疾病，其主要病变在关节的软骨及其滑膜组织，先前的研究已经报道在关节软骨中存在几种不同亚型的软骨细胞，包括增殖性软骨细胞(ProCs)、肥大前软骨细胞(preHTCs)、肥大软骨细胞(HTCs)、纤维软骨细胞(FC)、衰老细胞(SNCs)[63-64]，它们分布在软骨中的不同区域并具有不同的功能；另外成纤维细胞在炎症性疾病生物学中具有重要意义，但对其潜在的细胞异质性来源了解甚少[65]。虽然发现了不同的细胞亚群与某些骨关节疾病相关，但对于某个特定疾病中存在的细胞亚群的了解尚不完全清楚，同时上述细胞群在疾病的发病机制中的作用也尚不清楚。单细胞RNA 测序可以从单个细胞基因表达的差异去剖析、了解不同细胞亚群在疾病的发生发展全过程中的作用机制。最近多项研究利用单细胞RNA 测序技术对骨关节炎滑膜组织进行测序后确定了组织中主要的一些细胞亚群，并发现了以前尚未报道过的成纤维细胞亚群，它们分布在滑膜组织的不同区域，其发挥的作用也不一样，揭示了疾病的特定病理改变与不同细胞亚群之间的关系[66-69]。另外不同的疾病其滑膜组织内细胞亚群组成也不尽相同，且同一类型细胞亚群在不同的疾病中具有不同的病理改变及作用。STEPHENSON 等[67]对类风湿性关节炎患者滑膜组织进行单细胞RNA 测序，将获得的细胞归类为13 个亚群，鉴定出多种以前尚未报道过的新的细胞亚群，同时推测出其在疾病过程中可能发挥的作用。WADE 等[70]应用单细胞RNA 测序技术在银屑病关节炎滑膜组织中分离出来一种多功能性T 细胞的特定亚群的富集即：CD4+CD161+T 细胞，这些细胞表达多种细胞因子：包括白细胞介素17A、干扰素γ、肿瘤坏死因子、集落刺激因子，并且这些细胞因子与银屑病关节炎的疾病活动密切相关。单细胞RNA 测序在滑膜组织中的运用及发现可以为未来滑膜性疾病的诊断及治疗提供新的靶点。
最近JI 等[17]对骨性关节炎患者关节软骨进行单细胞RNA 测序，鉴定了以前尚未发现的3 种新的软骨细胞簇：即效应软骨细胞(EC)、调节性软骨细胞(RegC)和体内稳态软骨细胞(HomCs)，以及软骨细胞簇的新标记物，其中效应软骨细胞、调节性软骨细胞和体内稳态软骨细胞可以保护软骨延缓骨性关节炎的发展，而之前报道过的软骨细胞亚群可以促进骨性关节炎的发生；此外，他们发现体内稳态软骨细胞高表达人类昼夜节律标志物，如PER1 和SIRT1[71]，以此推断该软骨细胞群可能在骨性关节炎软骨的昼夜节律调节中起着重要作用，这为骨性关节炎的治疗提供了新的靶标及方向。
3 讨论 Discussion
近年来，单细胞RNA 测序技术得到快速发展，其已被广泛用于肿瘤、胚胎发育、免疫学、神经生物学、骨骼系统等学科研究。但是单细胞RNA 测序技术仍面临诸多挑战，首先，并不是所有组织都能成功分离成单细胞悬液并获得足够的单细胞数量，例如骨组织，尽管有成功的案例[72-73]；其次，在消化分离成单细胞过程中，应激反应会改变细胞基因表达[74]，无法获得细胞表达的原始状态；再次，单细胞RNA 测序成本较高[75]；另外，庞大的数据量使得分析变得更加复杂，需要专业的技术人员进行相关操作等等。总之，相信随着单细胞RNA 测序技术的不断发展完善，其会越来越多的用于骨骼系统发育及相关疾病的研究，使人们对其发病机制有了更为深刻和具体的认识。
单细胞RNA 测序技术作为一门新技术，具有良好的发展前景，已在诸多领域得到应用，目前在骨骼相关领域的研究已经得到重视并逐渐增多。单细胞RNA 测序技术可以帮助更好的了解骨骼发育、创伤后修复、骨骼相关疾病中细胞的异质性，发现新的细胞亚群、细胞标志物以及细胞的潜在功能，加深对疾病发病机制的认识，识别个体间的差异，将推动未来再生医疗、精准治疗及个体化治疗发展。
作者贡献：文章设计为晏怡果、朱斌，资料收集、分析及论文撰写为朱斌，晏怡果负责文章审校。
经费支持：该文章接受了“湖南省卫生计生委科研计划课题重点项目(A)、湖南省技术创新引导计划临床医疗技术创新引导项目(2017SK)、中国博士后科学基金面上项目(2017M)及衡阳市骨科内植物研发重点实验室(2018KJ115)”的资助。所有作者声明，经费支持没有影响文章观点和对研究数据客观结果的统计分析及其报道。
利益冲突：文章的全部作者声明，在课题研究和文章撰写过程中不存在利益冲突。
写作指南：该研究遵守《系统综述和荟萃分析报告规范》(PRISMA指南)。
文章查重：文章出版前已经过专业反剽窃文献检测系统进行3 次查重。
文章外审：文章经小同行外审专家双盲外审，同行评议认为文章符合期刊发稿宗旨。
文章版权：文章出版前杂志已与全体作者授权人签署了版权相关协议。
开放获取声明：这是一篇开放获取文章，根据《知识共享许可协议》“署名-非商业性使用-相同方式共享4.0”条款，在合理引用的情况下，允许他人以非商业性目的基于原文内容编辑、调整和扩展，同时允许任何用户阅读、下载、拷贝、传递、打印、检索、超级链接该文献，并为之建立索引，用作软件的输入数据或其它任何合法用途。
[1]OUYANG HW,CAO T,ZOU XH,et stem cell sheets revitalize nonviable dense grafts:implications for repair of large-bone and tendon ;82(2):170-174.
[2]OUYANG HW,GOH JC,LEE EH.Use of bone marrow stromal cells for tendon graft-to-bone healing:histological and immunohistochemical studies in a rabbit model.Am J Sports Med.2004;32(2):321-327.
[3]OUYANG HW,GOH JC,MO XM,et al.The efficacy of bone marrow stromal cell-seeded knitted PLGA fiber scaffold for Achilles tendon N Y Acad Sci.2002;961:126-129.
[4]ALLISON DC,CARNEY SC,AHLMANN ER,et al.A meta-analysis of osteosarcoma outcomes in the modern medical ;2012:.
[5]KANSARA M,TSANG M,KODJABACHIAN L,et al.Wnt inhibitory factor 1 is epigenetically silenced in human osteosarcoma,and targeted disruption accelerates osteosarcomagenesis in mice.J Clin ;119(4):837-851.
[6]HINGORANI P,ZHANG W,GORLICK R,et of Src Phosphorylation Alters metastatic Potential of Osteosarcoma In vitro but not In Cancer Res.2009;15(10):3416-3422.
[7]CAI Y,CAI T,CHEN Y.Wnt pathway in osteosarcoma,from oncogenic to Cell ;115(4):625-631.
[8]GUO M,CAI C,ZHAO G,et promotes migration and induces CXCR4 expression via HIF-1α activation in human ONE.2014;9:e.
[9]SAVAGE SA,MIRABELLO L,WANG Z,et association study identifies two susceptibility loci for ;45:799-803.
[10] DUNN SL,SOUL J,ANAND S,et al.Gene expression changes in damaged osteoarthritic cartilage identify a signature of nonchondrogenic and mechanical ;24(8):1431-1440.
[11] WANG J,SONG Y.Single cell sequencing:a distinct new Transl Med.2017;6(1):10.
[12] BJORKLUND AK,FORKEL M,PICELLI S,et al.The heterogeneity of human CD127(+) innate lymphoid cells revealed by single-cell RNA ;17(4):451-460.
[13] MEREDITH M,ZEMMOUR D,MATHIS D,et al.Aire controls gene expression in the thymic epithelium with ordered ;16(9):942-949.
[14] POULIN JF,GAERTNER Z,MORENO-RAMOS OA,et of midbrain dopamine neurons using single-cell gene expression profiling ;43(3):155-169.
[15] SWANSON JB,DE MICHELI AJ,DISSER NP,et al.A single-cell transcriptional atlas identifies extensive heterogeneity in the cellular composition of :10.1101/.
[16] RUBENSTEIN AB,SMITH GR,RAUE U,et transcriptional profiles in human skeletal Rep.2020;10(1):229.
[17] JI Q,ZHENG Y,ZHANG G,et RNA-seq analysis reveals the progression of human Rheum Dis.2019;78(1):100-110.
[18] TANG F,BARBACIORU C,WANG Y,et whole-transcriptome analysis of a single ;6(5):377-382.
[19] SHAPIRO E,BIEZUNER T,LINNARSSON S,et sequencingbased technologies will revolutionize whole-organism Rev ;14(9):618-630.
[20] ZHANG XN,LI TQ,LIU F,et Analysis of Droplet-based Ultra-High-Throughput Single-Cell RNA-Seq ;73(1):
[21] MACOSKO EZ,BASU A,SATIJA R,et parallel genome-wide expression profiling of individual cells using nanoliter ;161(5):1202-1214.
[22] WANG Y,NAVIN and applications of single-cell sequencing ;58(4):598-609.
[23] OZSOLAK F,PLATT AR,JonES DR,et RNA ure.2009;461(7265):814-818.
[24] RODRIGUE S,MALMSTROM RR,BERLIN AM,et al.Whole genome amplification and de novo assembly of single bacterial One.2009;4(9):e6864.
[25] ZONG CH,LU SJ,CHAPMAN AR,et detection of singlenucleotide and copy-number variations of a single human ;338(6114):1622-1626.
[26] ISCOVE NN,BARBARA M,GU M,et is faithfully preserved in global cDNA amplified exponentially from sub-picogram quantities of mRNA.Nat ;20(9):940-943.
[27] PARK J,SHRESTHA R,QIU C,et transcriptomics of the mouse kidney reveals potential cellular targets of kidney ;360(6390):758-763.
[28] WU H,HUMPHREYS BD.The promise of single-cell RNA sequencing for kidney disease Int.2017;92(6):1334-1342.
[29] CHAN CKF,SEO EY,CHEN JY,et and specification of the mouse skeletal stem cell.;160:285-298.
[30] MIZUHASHI K,ONO W,MATSUSHITA Y,et zone of the growth plate houses a unique class of skeletal stem ;563(7730):254-258.
[31] CHAN CKF,GULATI GS,SINHA R,et of the human skeletal stem ;175(1):
[32] MIZUHASHI K,NAGATA M,MATSUSHITA Y,et plate borderline chondrocytes behave as transient mesenchymal precursor cells.J Bone Miner Res.2019;34(8):1387-1392.
[33] TAKAHASHI A,NAGATA M,GUPTA A,et regulation of mesenchymal progenitor cell fates orchestrates tooth Natl Acad Sci U S A.2019;116(2):575-580.
[34] DEBNATH S,YALLOWITZ AR,MCCORMICK J,et of a periosteal stem cell mediating intramembranous bone ;562(7725):133-139.
[35] ZHOU B,YUE R,MURPHY MM,et mesenchymal stromal cells represent the main source of bone formed by adult bone Stem ;15(2):154-168.
[36] SACCHETTI B,FUNARI A,MICHIENZI S,et osteoprogenitors in bone marrow sinusoids can organize a hematopoietic ;131(2):324-336.
[37] BI W,DENG JM,ZHANG Z,et al.Sox9 is required for cartilage ;22:85-89.
[38] YAMASHIRO T,WANG XP,LI Z,et roles of Runx1 and Sox9 in incipient intramembranous Bone Miner Res.2004;19(10):1671-1617.
[39] ONO N,ONO W,NAGASAWA T,et al.A subset of chondrogenic cells provides early mesenchymal progenitors in growing Cell ;16(12):1157-1167.
[40] AKIYAMA H,KIM JE,NAKASHIMA K,et cells are derived from Sox9 expressing Natl Acad Sci U S A.2005;102(41):-.
[41] KELLY NH,HUYNH NPT,GUILAK F.Single cell RNA-sequencing reveals cellular heterogeneity and trajectories of lineage specification during murine embryonic limb ;89:1-10.
[42] DECKER RS,KOYAMA E,PACIFICI and morphogenesis of limb synovial joints and articular ;39:5-10.
[43] LEE CH,LEE FY,TARAFDER S,et endogenous stem/progenitor cells for tendon Clin ;125(7):2690-2701.
[44] PRYCE BA,BRENT AE,MURCHISON ND,et of transgenic tendon reporters,ScxGFP and ScxAP,using regulatory elements of the scleraxis gene.Dev Dyn.2007;236(6):1677-1182.
[45] HUANG AH,LU HH,SCHWEITZER regulation of tendon cell fate during Orthop Res.2015;33(6):800-812.
[46] DISSER NP,GHAHRAMANI GC,SWANSON J,et diversity in the transcriptomes of functionally divergent limb ;598(8):1537-1550.
[47] HARVEY T,FLAMENCO S,FAN CM.A Tppp3Pdgfra tendon stem cell population contributes to regeneration and reveals a shared role for PDGF signalling in regeneration and Cell ;21(12):1490-1503.
[48] LESCROART F,WANG XN,LIN XH,et the earliest step of cardiovascular lineage segregation by single-cell ;359:1177-1181.
[49] HAREL I,NATHAN E,TIROSH-FINKEL L,et origins and genetic programs of head muscle satellite ;16(6):822-832.
[50] CHAN S,SHI XZ,TOYAMA A,et al.Mesp1 patterns mesoderm into cardiac,hematopoietic,or skeletal myogenic progenitors in a contextdependent Stem ;12:587-601.
[51] PENALOZA JS,PAPPAS MP,HAGEN HR,et RNA-seq analysis of Mesp1-induced skeletal myogenic Biophys Res ;520(2):284-290.
[52] BARRUET E,GARCIA SM,STRIEDINGER K,et heterogeneous human satellite cells identified by single cell RNA ;9:e.
[53] DE MICHELI AJ,SPECTOR JA,ELEMENTO O,et al.A reference single-cell transcriptomic atlas of human skeletal muscle tissue reveals bifurcated muscle stem cell ;10:19.
[54] KANN AP,KRAUSS RNAscope and immunofluorescence on whole-mount skeletal myofibers and their associated stem elopment.2019;146(20):dev.
[55] GNOCCHI VF,WHITE RB,ONO Y,et characterisation of the molecular signature of quiescent and activated mouse muscle satellite One.2009;4(4):e5205.
[56] EINHORNA T,GERSTENFELD healing:mechanisms and Rev ;11(1):45-54.
[57] SALAZAR VS,GAMER LW,ROSEN V.BMP signalling in skeletal development,disease and Rev ;12(4):203-221.
[58] MCBRIDE SH,MCKENZIE JA,BEDRICK BS,et al.Long bone structure and strength depend on BMP2 from osteoblasts and osteocytes,but not vascular endothelial One.2014;9(5):e.
[59] BOLANDER J,CHAI YC,GERIS L,et al.Early BMP,Wnt and Ca(2+)/PKC pathway activation predicts the bone forming capacity of periosteal cells in combination with calcium ;86:106-118.
[60] BOLANDER J,JI W,GERIS L,et al.The combined mechanism of bone morphogenetic protein-and calcium phosphate-induced skeletaltissue formation by human periosteum derived Cell ;31:11-25.
[61] BOLANDER J,HERPELINCK T,CHAKLADER M,et characterization and metabolic profiling of in vitro cultured human skeletal progenitors with enhanced in vivo bone forming Cells Transl Med.2020;9(3):389-402.
[62] BOLANDER J,JI W,LEIJTEN J,et of a Large Long-Bone Defect through Serum-Free In Vitro Priming of Human Periosteum-Derived Cell ;8:758-772.
[63] PREIN C,WARMBOLD N,FARKAS Z,et and mechanical properties of the proliferative zone of the developing murine growth plate cartilage assessed by atomic force ;50:1-15.
[64] JIANG Y,TUANR S.Origin and function of cartilage stem/progenitor cells in Rev ;11(4):206-212.
[65] BUCKLEY CD,BARONE F,NAYAR S,et cells in chronic inflammation and tertiary lymphoid organ Rev ;33:715-745.
[66] CROFT AP,CAMPOS J,JANSEN K,et fibroblast subsets drive inflammation and damage in ;570:246-251.
[67] STEPHENSON W,DonLIN LT,BUTLER A,et RNA-seq of rheumatoid arthritis synovial tissue using low-cost microfluidic ;9(1):791.
[68] CAI SZ,MING BX,YE C,et Transition Processes in Synovial Fibroblasts from Rheumatoid Arthritis and Osteoarthritis A Single-Cell Study.J Immunol Res.2019;2019:.
[69] MIZOGUCHI F,SLOWIKOWSKI K,KEVIN W,et distinct disease-associated fibroblast subsets in rheumatoid ;9(1):789.
[70] WADE SM,CANAVAN M,MCGARRY T,et of synovialtissue polyfunctional T-cells with DAPSA in psoriatic Rheum Dis.2019;78(3):350-354.
[71] BERENBAUM F,MENG QJ.The brain-joint axis in osteoarthritis:nerves,circadian clocks and Rev ;12(9):508-516.
[72] SHAH KM,STERN MM,STERN AR,et isolation and culture ;5:838.
[73] TANAKA K-I,XUE Y,NGUYEN-YAMAMOTO L,et is a noveldeterminant of bone and muscle structure and Natl Acad Sci U S A.2018;115(16):E3759-E3768.
[74] ADAM M,POTTER AS,POTTER proteases dramatically reduce single-cell RNA-seq artifacts:a molecular atlas of kidney ;144(19):3625-3632.
[75] AYTURK in Skeletal Osteoporos Rep.2019;17(4):178-185.
文章来源：中国免疫学 网址: http://zgmyx.400nongye.com/lunwen/itemid-55530.shtml
上一篇： 以案例为导向的感染与免疫学实验教学探究
下一篇： 刑法论文_浅析医疗事故罪之因果关系的认定
