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)。




随着测序技术的发展，最初使用全转录组测序对单细胞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]，以此推断该软骨细胞群可能在骨性关节炎软骨的昼夜节律调节中起着重要作用，这为骨性关节炎的治疗提供了新的靶标及方向。

[1] OUYANG HW, CAO T, ZOU XH, et al. Mesenchymal stem cell sheets revitalize nonviable dense grafts: implications for repair of large-bone and tendon defects. Transplantation. 2006;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 repair. Ann 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 era. Sarcoma. 2012; 2012:704872. [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 Invest. 2009;119(4):837-851. [6] HINGORANI P, ZHANG W, GORLICK R, et al. Inhibition of Src Phosphorylation Alters Metastatic Potential of Osteosarcoma In vitro but not In vivo. Clin Cancer Res. 2009;15(10):3416-3422. [7] CAI Y, CAI T, CHEN Y. Wnt pathway in osteosarcoma, from oncogenic to therapeutic. J Cell Biochem. 2014;115(4):625-631. [8] GUO M, CAI C, ZHAO G, et al. Hypoxia promotes migration and induces CXCR4 expression via HIF-1α activation in human osteosarcoma. PLoS ONE. 2014;9:e90518. [9] SAVAGE SA, MIRABELLO L, WANG Z, et al. Genome-wide association study identifies two susceptibility loci for osteosarcoma. Nat Genet. 2013;45:799-803. [10] DUNN SL, SOUL J, ANAND S, et al. Gene expression changes in damaged osteoarthritic cartilage identify a signature of non-chondrogenic and mechanical responses. Osteoarthritis Cartilage. 2016;24(8):1431-1440. [11] WANG J, SONG Y. Single cell sequencing: a distinct new field. Clin 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 sequencing. Nat Immunol. 2016;17(4):451-460. [13] MEREDITH M, ZEMMOUR D, MATHIS D, et al. Aire controls gene expression in the thymic epithelium with ordered stochasticity. Nat Immunol. 2015;16(9):942-949. [14] POULIN JF, GAERTNER Z, MORENO-RAMOS OA, et al. Classification of midbrain dopamine neurons using single-cell gene expression profiling approaches. Trends Neurosci. 2020;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 tendons. BioRxiv. 2019. doi:10.1101/801266. [16] RUBENSTEIN AB,SMITH GR,RAUE U,et al. Single-cell transcriptional profiles in human skeletal muscle.Sci Rep. 2020;10(1):229. [17] JI Q, ZHENG Y, ZHANG G, et al. Single-cell RNA-seq analysis reveals the progression of human osteoarthritis. Ann Rheum Dis. 2019;78(1):100-110. [18] TANG F, BARBACIORU C, WANG Y, et al. mRNA-Seq whole-transcriptome analysis of a single cell. Nature methods. 2009;6(5):377-382. [19] SHAPIRO E, BIEZUNER T, LINNARSSON S, et al. Single-cell sequencing-based technologies will revolutionize whole-organism science. Nat Rev Genet. 2013;14(9):618-630. [20] ZHANG XN, LI TQ, LIU F, et al. Comparative Analysis of Droplet-Based Ultra-High-Throughput Single-Cell RNA-Seq Systems. Molecular Cell. 2019;73(1):130-142.e5. [21] MACOSKO EZ, BASU A, SATIJA R, et al. Highly parallel genome-wide expression profiling of individual cells using nanoliter droplets. Cell. 2015;161(5):1202-1214. [22] WANG Y, NAVIN NE. Advances and applications of single-cell sequencing technologies. Molecular Cell. 2015;58(4):598-609. [23] OZSOLAK F, PLATT AR, JONES DR, et al. Direct RNA sequencing. Nature. 2009;461(7265):814-818. [24] RODRIGUE S, MALMSTROM RR, BERLIN AM, et al. Whole genome amplification and de novo assembly of single bacterial cells. PloS One. 2009;4(9):e6864. [25] ZONG CH, LU SJ, CHAPMAN AR, et al. Genome-wide detection of single-nucleotide and copy-number variations of a single human cell. Science. 2012;338(6114):1622-1626. [26] ISCOVE NN, BARBARA M, GU M, et al. Representation is faithfully preserved in global cDNA amplified exponentially from sub-picogram quantities of mRNA. Nat Biotechnol. 2002;20(9):940-943. [27] PARK J, SHRESTHA R, QIU C, et al. Single-cell transcriptomics of the mouse kidney reveals potential cellular targets of kidney disease.Science. 2018;360(6390):758-763. [28] WU H, HUMPHREYS BD. The promise of single-cell RNA sequencing for kidney disease investigation. Kidney Int. 2017;92(6):1334-1342. [29] CHAN CKF, SEO EY, CHEN JY, et al. Identification and specification of the mouse skeletal stem cell. Cell. 2015;160:285-298. [30] MIZUHASHI K, ONO W, MATSUSHITA Y, et al. Resting zone of the growth plate houses a unique class of skeletal stem cells. Nature. 2018; 563(7730):254-258. [31] CHAN CKF, GULATI GS, SINHA R, et al. Identification of the human skeletal stem cell. Cell. 2018;175(1):43-56.e21. [32] MIZUHASHI K, NAGATA M, MATSUSHITA Y, et al. Growth 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 al. Autocrine regulation of mesenchymal progenitor cell fates orchestrates tooth eruption. Proc Natl Acad Sci U S A. 2019;116(2):575-580. [34] DEBNATH S, YALLOWITZ AR, MCCORMICK J, et al. Discovery of a periosteal stem cell mediating intramembranous bone formation. Nature. 2018;562(7725):133-139. [35] ZHOU B, YUE R, MURPHY MM, et al. Leptin-receptor-expressing mesenchymal stromal cells represent the main source of bone formed by adult bone marrow. Cell Stem Cell. 2014;15(2):154-168. [36] SACCHETTI B, FUNARI A, MICHIENZI S, et al. Self-renewing osteoprogenitors in bone marrow sinusoids can organize a hematopoietic microenvironment. Cell. 2007;131(2):324-336. [37] BI W, DENG JM, ZHANG Z, et al. Sox9 is required for cartilage formation.Nature genetics. 1999;22:85-89. [38] YAMASHIRO T, WANG XP, LI Z, et al. Possible roles of Runx1 and Sox9 in incipient intramembranous ossification. J 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 bones. Nature Cell Biology. 2014;16(12):1157-1167. [40] AKIYAMA H, KIM JE, NAKASHIMA K, et al. Osteo-chondroprogenitor cells are derived from Sox9 expressing precursors. Proc Natl Acad Sci      U S A. 2005;102(41):14665-14670. [41] KELLY NH, HUYNH NPT, GUILAK F. Single cell RNA-sequencing reveals cellular heterogeneity and trajectories of lineage specification during murine embryonic limb development. Matrix Biol. 2020;89:1-10. [42] DECKER RS, KOYAMA E, PACIFICI M. Genesis and morphogenesis of limb synovial joints and articular cartilage. Matrix Biol. 2014;39:5-10. [43] LEE CH, LEE FY, TARAFDER S, et al. Harnessing endogenous stem/progenitor cells for tendon regeneration. J Clin Invest. 2015;125(7): 2690-2701.  [44] PRYCE BA, BRENT AE, MURCHISON ND, et al. Generation 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 R. Molecular regulation of tendon cell fate during development. J Orthop Res. 2015;33(6):800-812.  [46] DISSER NP, GHAHRAMANI GC, SWANSON J, et al. Widespread diversity in the transcriptomes of functionally divergent limb tendons. J Physiol. 2020;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 fibrosis.Nat Cell Biol. 2019;21(12): 1490-1503.  [48] LESCROART F, WANG XN, LIN XH, et al. Defining the earliest step of cardiovascular lineage segregation by single-cell RNA-seq. Science. 2018;359:1177-1181. [49] HAREL I, NATHAN E, TIROSH-FINKEL L, et al. Distinct origins and genetic programs of head muscle satellite cells. Dev Cell. 2009;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 context-dependent manner. Cell Stem Cell. 2013;12:587-601. [51] PENALOZA JS, PAPPAS MP, HAGEN HR, et al. Single-cell RNA-seq analysis of Mesp1-induced skeletal myogenic development. Biochem Biophys Res Commun. 2019;520(2):284-290.  [52] BARRUET E, GARCIA SM, STRIEDINGER K, et al. Functionally heterogeneous human satellite cells identified by single cell RNA sequencing. Elife. 2020;9:e51576.  [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 populations. Skeletal Muscle. 2020;10:19. [54] KANN AP, KRAUSS RS. Multiplexed RNAscope and immunofluorescence on whole-mount skeletal myofibers and their associated stem cells.Development. 2019;146(20):dev179259. [55] GNOCCHI VF, WHITE RB, ONO Y, et al. Further characterisation of the molecular signature of quiescent and activated mouse muscle satellite cells. PLoS One. 2009;4(4):e5205. [56] EINHORNA T, GERSTENFELD LC. Fracture healing: mechanisms and interventions. Nat Rev Rheumatol. 2015;11(1):45-54. [57] SALAZAR VS, GAMER LW, ROSEN V. BMP signalling in skeletal development, disease and repair. Nat Rev Endocrinol. 2016;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 cells. PLoS One. 2014;9(5):e96862. [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 phosphates. Biomaterials. 2016;86: 106-118. [60] BOLANDER J, JI W, GERIS L, et al. The combined mechanism of bone morphogenetic protein- and calcium phosphate-induced skeletal tissue formation by human periosteum derived cells. Eur Cell Mater. 2016;31:11-25. [61] BOLANDER J, HERPELINCK T, CHAKLADER M, et al. Single-cell characterization and metabolic profiling of in vitro cultured human skeletal progenitors with enhanced in vivo bone forming capacity.Stem Cells Transl Med. 2020;9(3):389-402. [62] BOLANDER J, JI W, LEIJTEN J, et al. Healing of a Large Long-Bone Defect through Serum-Free In Vitro Priming of Human Periosteum-Derived Cells. Stem Cell Reports. 2017;8:758-772. [63] PREIN C, WARMBOLD N, FARKAS Z, et al. Structural and mechanical properties of the proliferative zone of the developing murine growth plate cartilage assessed by atomic force microscopy. Matrix Biol. 2016;50:1-15.  [64] JIANG Y, TUANR S. Origin and function of cartilage stem/progenitor cells in osteoarthritis. Nat Rev Rheumatol. 2015;11(4):206-212. [65] BUCKLEY CD, BARONE F, NAYAR S, et al. Stromal cells in chronic inflammation and tertiary lymphoid organ formation. Annu Rev Immunol. 2015;33:715-745.  [66] CROFT AP, CAMPOS J, JANSEN K, et al. Distinct fibroblast subsets drive inflammation and damage in arthritis. Nature. 2019;570:246-251. [67] STEPHENSON W, DONLIN LT, BUTLER A, et al. Single-cell RNA-seq of rheumatoid arthritis synovial tissue using low-cost microfluidic instrumentation. Nat Commun. 2018;9(1):791. [68] CAI SZ, MING BX, YE C, et al. Similar Transition Processes in Synovial Fibroblasts from Rheumatoid Arthritis and Osteoarthritis A Single-Cell Study. J Immunol Res. 2019;2019:4080735. [69] MIZOGUCHI F, SLOWIKOWSKI K, KEVIN W, et al. Functionally distinct disease-associated fibroblast subsets in rheumatoid arthritis. Nat Commun. 2018;9(1):789.  [70] WADE SM, CANAVAN M, MCGARRY T, et al. Association of synovial tissue polyfunctional T-cells with DAPSA in psoriatic arthritis. Ann Rheum Dis. 2019;78(3):350-354.  [71] BERENBAUM F, MENG QJ. The brain-joint axis in osteoarthritis: nerves, circadian clocks and beyond. Nat Rev Rheumatol. 2016;12(9):508-516. [72] SHAH KM,STERN MM,STERN AR,et al.Osteocyte isolation and culture methods. Bonekey Reports. 2016;5:838. [73] TANAKA K-I, XUE Y, NGUYEN‐YAMAMOTO L, et al. FAM210A is a noveldeterminant of bone and muscle structure and strength.Proc Natl Acad Sci U S A. 2018;115(16):E3759-E3768.  [74] ADAM M, POTTER AS, POTTER SS. Psychrophilic proteases dramatically reduce single-cell RNA-seq artifacts: a molecular atlas of kidney development. Development. 2017;144(19):3625-3632. [75] AYTURK U. RNA-seq in Skeletal Biology. Curr Osteoporos Rep. 2019; 17(4):178-185.

[1]	蒲 锐, 陈子扬, 袁凌燕. 不同细胞来源外泌体保护心脏的特点与效应[J]. 中国组织工程研究, 2021, 25(在线): 1-.
[2]	张 超, 吕 欣. 髋臼骨折固定后的异位骨化：危险因素、预防及其治疗进展[J]. 中国组织工程研究, 2021, 25(9): 1434-1439.
[3]	周继辉, 李新志, 周 游, 黄 卫, 陈文瑶. 髌骨骨折修复内植物选择的多重问题[J]. 中国组织工程研究, 2021, 25(9): 1440-1445.
[4]	王德斌, 毕郑刚. 尺骨鹰嘴骨折-脱位解剖力学、损伤特点、固定修复及3D技术应用的相关问题[J]. 中国组织工程研究, 2021, 25(9): 1446-1451.
[5]	姬志祥, 蓝常贡. 尿酸盐转运蛋白在痛风中的多态性和治疗相关性[J]. 中国组织工程研究, 2021, 25(8): 1290-1298.
[6]	袁 美, 张新新, 郭祎莎, 毕 霞. 循环microRNA在血管性认知障碍诊断中的应用[J]. 中国组织工程研究, 2021, 25(8): 1299-1304.
[7]	汪显耀, 关亚琳, 刘忠山. 提高间充质干细胞治疗难愈性创面的策略[J]. 中国组织工程研究, 2021, 25(7): 1081-1087.
[8]	万 然, 史 旭, 刘京松, 王岩松. 间充质干细胞分泌组治疗脊髓损伤的研究进展[J]. 中国组织工程研究, 2021, 25(7): 1088-1095.
[9]	廖成成, 安家兴, 谭张雪, 王 倩, 刘建国. 口腔鳞状细胞癌干细胞的治疗靶点及应用前景[J]. 中国组织工程研究, 2021, 25(7): 1096-1103.
[10]	赵 敏, 冯柳祥, 陈 垚, 顾 霞, 王平义, 李一梅, 李文华. 低氧环境下外泌体可作为疾病的标志物[J]. 中国组织工程研究, 2021, 25(7): 1104-1108.
[11]	谢文佳, 夏天娇, 周卿云, 刘羽佳, 顾小萍. 小胶质细胞介导神经元损伤在神经退行性疾病中的作用[J]. 中国组织工程研究, 2021, 25(7): 1109-1115.
[12]	李珊珊, 郭笑霄, 尤 冉, 杨秀芬, 赵 露, 陈 曦, 王艳玲. 感光细胞替代治疗视网膜变性疾病[J]. 中国组织工程研究, 2021, 25(7): 1116-1121.
[13]	焦 慧, 张一宁, 宋雨晴, 林 宇, 王秀丽. 乳腺癌类器官研究进展及临床应用前景[J]. 中国组织工程研究, 2021, 25(7): 1122-1128.
[14]	王诗琦, 张金生. 中医药调控缺血缺氧微环境对骨髓间充质干细胞增殖、分化及衰老的影响[J]. 中国组织工程研究, 2021, 25(7): 1129-1134.
[15]	曾燕华, 郝延磊. 许旺细胞体外培养及纯化的系统性综述[J]. 中国组织工程研究, 2021, 25(7): 1135-1141.

骨骼系统作为人体八大系统之一，在人体中发挥着重要的作用，也从侧面反映出骨骼系统的复杂性，这给研究者带来了许多挑战；另外随着交通工具的快速发展、生活方式的改变及环境变化，骨骼系统相关性疾病的发生率也越来越高，例如：骨折、肌腱及韧带损伤、骨肿瘤、骨关节炎等。免疫学、流式细胞术及qRT-PCR等传统研究方法因其经济、快速、简单易操作等优点仍被广泛用于探究骨骼系统的发生发育及相关疾病的病理生理、发病机制、发生发展及转归。有研究报道骨髓基质细胞可以促进肌腱修复、加快肌腱与骨的植入愈合并证明间充质干细胞在长骨及肌腱缺损修复中具有重要作用[1-3]；大量研究已发现TP53、PTEN等诸多基因、蛋白及相关信号通路在骨肉瘤中具有重要作用[4-9]；通过RNA测序发现软骨相关基因Sox9、Col11a2、Acan和其他新转录因子在骨关节炎中的表达发生明显改变等等[10]，但上述传统方法无法很好地解释创伤后组织是如何进行修复及愈合？骨肿瘤组织中是否存在细胞异质性？骨骼、软骨、韧带等骨骼组织的起源和形成机制仍不清楚，无法完全了解参与骨骼组织发育的间充质细胞亚群的起源和潜在功能，骨髓基质及间充质细胞在骨骼组织创伤后修复中的作用和潜在功能，以及骨关节性疾病中软骨细胞类型组成和促进疾病进展的细胞异质性等相关问题仍不清楚，因此需要一种新的研究方法来进一步研究。
近些年来，以单细胞RNA测序技术为研究的方法正在得到蓬勃发展，目前单细胞RNA测序技术已广泛使用于肿瘤、免疫学、神经生物学等领域[11-14]，并且已逐渐开始用于骨骼系统生物学相关研究。相当多的研究通过单细胞RNA测序在骨骼系统中发现了以前尚未报道过的新细胞类型及细胞标志物，例如：在肌腱[15]、肌肉的单核细胞及骨关节软骨中都发现了新的细胞亚群[16-17]，这将有助于对该疾病有更深入的认识以及对疾病的精准治疗提供新见解，对未来再生医疗、精准治疗及个体化治疗具有重要意义。文章旨在综述单细胞RNA测序技术的基本原理及其在骨骼系统生物学相关研究中的应用。
中国组织工程研究杂志出版内容重点：组织构建；骨细胞；软骨细胞；细胞培养；成纤维细胞；血管内皮细胞；骨质疏松；组织工程

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。

近年来，单细胞RNA测序技术得到快速发展，其已被广泛用于肿瘤、胚胎发育、免疫学、神经生物学、骨骼系统等学科研究。但是单细胞RNA测序技术仍面临诸多挑战，首先，并不是所有组织都能成功分离成单细胞悬液并获得足够的单细胞数量，例如骨组织，尽管有成功的案例[72-73]；其次，在消化分离成单细胞过程中，应激反应会改变细胞基因表达[74]，无法获得细胞表达的原始状态；再次，单细胞RNA测序成本较高[75]；另外，庞大的数据量使得分析变得更加复杂，需要专业的技术人员进行相关操作等等。总之，相信随着单细胞RNA测序技术的不断发展完善，其会越来越多的用于骨骼系统发育及相关疾病的研究，使人们对其发病机制有了更为深刻和具体的认识。
单细胞RNA测序技术作为一门新技术，具有良好的发展前景，已在诸多领域得到应用，目前在骨骼相关领域的研究已经得到重视并逐渐增多。单细胞RNA测序技术可以帮助更好的了解骨骼发育、创伤后修复、骨骼相关疾病中细胞的异质性，发现新的细胞亚群、细胞标志物以及细胞的潜在功能，加深对疾病发病机制的认识，识别个体间的差异，将推动未来再生医疗、精准治疗及个体化治疗发展。
中国组织工程研究杂志出版内容重点：组织构建；骨细胞；软骨细胞；细胞培养；成纤维细胞；血管内皮细胞；骨质疏松；组织工程

文题释义：

单细胞RNA测序技术：指在单个细胞水平上对mRNA进行扩增与测序，并进行定量分析的一项新技术，能够显示每个细胞的变化，发现新的细胞类型、细胞标志物及潜在功能。#br# 骨骼系统生物学：是指脊椎动物的器官系统之一，包括骨骼、关节、韧带等组织，对骨骼系统的正常生理及相关疾病的病理生理及发病机制的相关研究。

中国组织工程研究杂志出版内容重点：组织构建；骨细胞；软骨细胞；细胞培养；成纤维细胞；血管内皮细胞；骨质疏松；组织工程

目前对骨组织的起源和形成机制、创伤后修复的细胞亚群的分子机制和相互作用及骨关节性疾病、骨肿瘤等#br# 骨相关疾病中细胞异质性及疾病的发生发展等问题仍不清楚。免疫学、流式细胞术及qRT-PCR 等传统方法是从整#br# 体上对骨骼系统进行研究，来探索骨骼系统的发生发育及相关疾病的发生机制，但无法了解单个细胞的类型和状#br# 态，容易忽略单个细胞之间的差异并掩盖稀有样品的特征。文章旨在综述单细胞RNA 测序技术的基本原理#br# 及其在骨骼系统生物学相关研究中的应用。#br# 中国组织工程研究杂志出版内容重点：组织构建；骨细胞；软骨细胞；细胞培养；成纤维细胞；血管内皮细胞；骨质疏松；组织工程

	阅读次数
	全文
	摘要