作者投稿
专家审稿
编辑办公
邮件订阅
RSS
-
摘要:目的
结合长骨关节端的声像表现，初步探索采用超声定量评估骺软骨的发育情况，并分析其与骨龄的相关性。
方法对2023年3—6月就读于济宁市某体校的青少年进行研究，记录年龄、身高等信息，于1周内行X线骨龄检查，并由1名内分泌科医师使用Greulich-Pyle图谱法评估骨龄。由2名超声医学科医师针对非优势侧手及腕部、膝关节共7个部位(包括第三掌骨头部背侧、尺骨茎突、桡骨茎突、股骨内上髁、股骨外上髁、胫骨背侧、胫骨内侧髁)进行超声检查，测量超声骺软骨厚度(ultrasonic epiphysis cartilage thickness，UECT)，并分析其与骨龄的相关性。采用组内相关系数(intraclass correlation coefficient, ICC)评估UECT测量的信度与可重复性。
结果共纳入141名青少年，其中男性80名，女性61名。中位年龄为13(12，14)岁，男性中位骨龄为15(13，17)岁，女性中位骨龄16(14，17)岁，男、女骨龄与年龄差值＞1岁者分别占比60%、67%，男性身高显著高于女性。男、女7个部位UECT均随骨龄增长而减小，各部位UECT与骨龄均呈负相关。男性单部位UECT与骨龄均呈高度负相关(|r|≥0.80)，女性则表现出中高度负相关(0.65≤|r|≤0.75)。男、女7个部位UECT总和与骨龄的负相关性均达到较高水平(|r|=0.93，|r|=0.80)，单独选取手和腕部3个部位或膝关节4个部位时，女性UECT总和与骨龄呈高度负相关(|r|≥0.76)，而男性则表现出非常高度负相关性(|r|≥0.90)。7个部位UECT的测量一致性均较好(ICC≥0.75)。
结论超声可观察描述长骨关节端的生长发育变化，UECT与骨龄表现出中至高度负相关，具有用于骨成熟度定量评估的潜在价值。
Abstract:ObjectivePreliminary exploration of using ultrasound to quantitatively evaluate the development of epiphyseal cartilage and analyze its correlation with bone age, based on the ultrasound findings of the long bone joint end.
MethodsA study was conducted on adolescents studying at a sports school in Jining from March to June 2023. Age, height and other information were recorded. Bone age assessment by X ray were performed within 1 week with an endocrinologist interpreted the bone age using the Greulich-Pyle atlas. Two sonographers scanned a total of 7 sites in the hand, wrist and knee joint of the non-dominant side (including the dorsal side of the third metacarpal head, the ulnar styloid process, the radial styloid process, the medial and the lateral femoral epicondyle, dorsal tibia, and medial tibial condyle). The ultrasonic epiphyseal cartilage thickness (UECT) was measured and its correlation with bone age was analyzed. Intraclass correlation coefficient (ICC) was used to evaluate the reliability and repeatability of UECT measurement.
ResultsA total of 141 adolescents were included, with 80 males and 61 females. The average age was 13 (12-14) years old. The average bone age was 15(13-17) years in males and 16 (14-17) years in females. The proportion of the males and females whose bone age was 1 year older than chronological age was 60% and 67%, respectively, and the height of the males was significantly higher than that of the females(P < 0.001). The UECT of 7 sites in the males and females decreased with the increase of bone age, and there was a significant negative correlation between UECT and bone age. In males, UECT was highly and negatively correlated with bone age (|r|≥0.80), while in females, it was moderately and negatively correlated (0.65≤|r|≤0.75). The correlation between the sum of UECT at seven sites and bone age reached high levels in both sexes(|r|=0.93, |r|=0.80).When 3 sites of hand and wrist or 4 sites of knee joint were selected separately, the UECT of the females was highly negatively correlated with bone age (|r| ≥0.76), while that of the males showed a very high correlation (|r| ≥0.90). The measurement consistency of UECT in all seven sites was good (ICC≥0.75).
ConclusionUltrasound can observe and describe the developmental change of long bone, and UECT has potential value in quantitative evaluation of bone maturity.
-
Keywords:
- ultrasound /
- bone age /
- epiphysis /
- adolescent /
- growth and development
-
利用测序技术分析肿瘤遗传变异、特异性代谢活动及进化过程，对肿瘤发生及进展机制的研究有重要意义，且为肿瘤临床治疗提供理论基础。相对于整体测序，单细胞测序技术更有利于阐释肿瘤异质性，且能从有限的样品中获得更多的信息。成人最常见的脑恶性肿瘤是脑转移瘤，而原发脑恶性肿瘤以胶质母细胞瘤(glioblastoma, GBM)最为常见，存在较强的异质性。在分子病理指导精准医疗的背景下，脑恶性肿瘤患者或因无法接受颅脑手术而失去分子病理检测和更有效的靶向治疗的机会。同时，由于颅脑手术创伤较大，多数患者无法在病情进展时重复手术治疗，利用外周血、脑脊液等微创活检的方式进行病程监测并指导个体化治疗是脑肿瘤治疗的方向。因此，鉴于单细胞测序技术的优势，其在脑恶性肿瘤临床诊疗中的应用备受期待。本文将综述单细胞测序技术流程及应用该技术进行肿瘤研究的进展，并结合脑恶性肿瘤的特征对单细胞测序技术在其诊疗中的应用进行展望。
1. 单细胞测序技术发展概述
癌细胞具有基因组不稳定性，在肿瘤起始和进展过程中可积累大量遗传突变，如单核苷酸突变、染色体变异等。利用测序技术测定癌细胞突变类型、功能，有助于探究癌症发生发展机制，并有针对性地研制治疗方案。依据DNA复制的分子生物学原理，Sanger等利用双脱氧核苷酸链终止法进行DNA片段测序，开启了基因组测序的新时代。但Sanger测序法操作繁琐、耗时长、成本高，大规模应用受限。随着荧光染色和毛细管微阵列电泳技术应用于基因组测序，第二代测序技术(next- generation sequencing，NGS)基于“边合成边测序”的原则，实现大规模、高通量测序，显著降低了测序成本。NGS法也具有一定局限性，包括测序片段较短，不能反映基因组长片段的变异，如拷贝数变异(copy number alteration，CNA)和染色体结构变异。第三代测序利用单分子实时测序及基于二代测序的合成长读长测序技术，DNA片段读长可至数十Kb，同时可进行完整的mRNA测序, 有利于对特殊基因的外显子剪切进行精确分析[1]。
既往肿瘤多组学研究通常采用NGS和/或Sanger测序技术进行肿瘤组织整体测序，获得一块肿瘤组织中成千上万个细胞的总体特征，在鉴别并筛选肿瘤相关突变基因、信号通路方面取得大量成果，促进了肿瘤临床诊断和治疗的发展。但肿瘤组织存在异质性，且不同亚型的肿瘤细胞之间的相互作用可影响疾病进程。与肿瘤组织整体测序相比，单细胞测序技术在评估肿瘤内异质性及解决特殊临床问题方面具有优势。第一，单细胞测序技术可用于分析肿瘤内异质性。第二，对于一些浸润生长的肿瘤，单细胞测序技术的应用可更高效地分析肿瘤信息。例如胰腺癌活检组织中肿瘤细胞仅占30%，其余为大量基质细胞。此时，肿瘤细胞中基因组变异可因频率过低而不能检出，尤其是CNA等变异可因测序深度不足而被忽略。第三，对于肿瘤组织样本，临床常因取样少等原因样本受限，如细针活检样本、外周血循环肿瘤细胞(circulating tumor cells，CTCs)等，不足以构建测序文库，单细胞测序或可为临床决策提供有价值的信息[2]。
2. 肿瘤单细胞测序技术概述
相对于整体测序，单细胞测序数据的获取依赖于某些技术的发展，包括单个细胞的分离及分选、单细胞全基因组测序和全外显子测序、单细胞转录组扩增及测序、单细胞表观遗传测序，以满足测序底物要求、测序成本降低和足够的测序细胞数量、控制并评估测序结果的偏倚及误差。
2.1 肿瘤单细胞分离及分选
分离及分选用于测序的单个细胞是肿瘤单细胞测序技术的重要步骤，主要包括实体肿瘤细胞、肿瘤干细胞(cancer stem cells，CSCs)、CTCs和骨髓播散肿瘤细胞(disseminated tumor cells，DTCs)[3]。其中，实体瘤可通过机械切割及胶原酶消化分散为单个细胞，而DTCs来源于骨髓穿刺。同时，该技术还可应用于肿瘤微环境中其他类型细胞的分离及分选，如免疫细胞。
CTCs是存在于患者循环系统中的肿瘤细胞，在血液中数量非常少，每毫升血液仅有数十或数百个CTCs[4]，目前分离主要依据其物理或生物学特性。与循环系统中的血细胞相比，CTCs体积更大、密度小且所带电荷与表面蛋白相关，据此可利用梯度离心、微流体惯性聚焦(microfluidic inertial focusing, MIF)、微滤过器(microfabricated filters, MF)法和双向电泳法等进行分离[5]。CTCs的生物学分离方法包括免疫磁珠分离(immunomagnetic separation, IS)、经微流体免疫分离(microfluidics-enabled immunoseparation, MI)、固相酶联免疫斑点技术(epithelial immunospot, EPISPOT)法和入侵实验(invasion assay, IA)等[5]。其中最常见的是IS法，该技术利用携带捕获剂(通常为特异性抗体)的免疫磁珠富集具有特异性膜蛋白的CTCs。基于该技术建立的CellSearch检测系统已获得美国食品药品监督管理局(Food and Drug Administration，FDA)的批准[6]。考虑到测序成本及肿瘤异质性，在将样品分散成单个细胞后，有时需进一步分选用于测序的单细胞，特别是从实体肿瘤中取得的样本。样本细胞筛选通常根据CTCs的特殊分子标志或细胞表型进行，包括流式细胞术(flow cytometry，FACS)、微流体(microfluidics)系统、显微操作(micromanipulation)法和激光捕获显微切割(laser capture microdissection，LCM)[5]。
2.2 单细胞全基因组测序和全外显子测序
从单个细胞中提取的遗传物质不足以进行全基因组测序及分析，因此，通常使用多重置换扩增(multiple displacement amplilication，MDA)、PCR或二者相结合的方法进行基因扩增，如多重退火和成环循环扩增(multiple annealing and looping based amplification cycles, MALBAC)。PCR技术最先用于单细胞基因组扩增，包括选择散布在基因组中的保守序列或随机寡核苷酸序列作为扩增引物。由于基因组中不同位点与引物结合能力不同，在基因组扩增的过程中，存在大量的序列信息丢失。目前选择处于细胞周期中的四倍体细胞、简并寡核苷酸引物PCR(degenerate oligonucleotide primed-PCR，DOP-PCR)进行扩增[7]，同时，选择自我复制的细胞进行测序可能导致结果存在偏倚。MDA技术利用ϕ29 DNA聚合酶和耐核酸外切酶的随机引物进行全基因组扩增，该过程无需重复温度循环，且错配率低于常规PCR[8]。MALBAC技术首先利用特殊的随机引物(包括27个氨基酸固定序列和8个氨基酸的可变序列)进行MDA，再对扩增产物进行PCR扩增[9]。有研究表明，MDA的基因组覆盖率高于DOP-PCR及MALBAC，利用该方法测序对单核苷酸位点变异(single nucleotide variants，SNV)的检出率更高；但DOP-PCR和MALBAC对基因组覆盖的均一性更强，有利于检出长度超过1 Mb的CNA[10]。
单细胞全外显子测序(single-cell whole exome sequencing，scWES)需进行外显子捕获和富集，主要用于编码基因的SNV分析。相对于scWES，单细胞全基因组测序(single-cell whole genome sequencing, scWGC)更利于CNA分析并获得非编码区信息，但由于存在基因组覆盖率的限制，SNV的检出效率较低。在乳腺癌中应用单细胞核基因组测序(nucleus genome sequencing, nuc-Seq)技术对G2/M期多倍体细胞核进行扩增和测序分析，可达到单碱基对分辨率[11]。
2.3 单细胞转录组扩增及测序
单细胞RNA测序(single-cell RNA sequencing, scRNA-Seq)包括cDNA构建、cDNA扩增和测序，其结果可显示细胞特异性表达谱和外显子转录的剪切及修饰过程，有助于分析细胞间的异质性。不同基因的转录水平差异较大，多数mRNA拷贝数较低，因此需通过扩增达到转录组测序所需的富集程度。Smart-Seq/Smart-Seq2、Phi29转录组扩增(phi29-transcrip-tome amplification，PTA)和基于半随机引物PCR的转录组扩增(semi-random primed PCR-based transcriptome amplification, STA)是常见的扩增手段。由于存在扩增偏差，通常将上述技术与特殊分子识别(unique molecular identifiers，UMIs)技术联用。UMIs可进行内参控制并定量检测基因表达。近期，在Smart-Seq2基础上使用Tn5转座酶进行cDNA富集可获得完整的转录起始位点信息，其在B淋巴细胞转录组测序中表现出较好的准确性和可重复性[12]。此外，类似于高通量二代测序技术，Drop-Seq可实现单细胞mRNA的高通量测序，该技术将单个细胞分离至纳升大小的PCR体系中，同一个细胞来源的扩增产物具有相同的序列标记，可显著提高测序效率[13]。
2.4 单细胞表观遗传测序
2.4.1 单细胞甲基化测序
基于传统的重亚硫酸盐甲基化测序过程，单细胞全基因组重亚硫酸盐测序技术的测序效果已在小鼠基因组中得到验证[14]，该技术首先对单细胞DNA进行多次重亚硫酸盐转换处理，然后在DNA片段加上标记序列并进行扩增、测序。相比于传统方法，可有效减少重亚硫酸盐处理导致的DNA降解和信息损失。另一种常见的技术是简化重亚硫酸盐测序(reduced representation bisulfite sequencing, RRBS)，该技术使用核酸内切酶消化DNA并富集特定的片段进行分析，在检测CpG岛方面更有优势[15]。受限于单细胞DNA含量，相较于整体样品甲基化测序，单细胞DNA甲基化测序技术的CpGs位点覆盖率较低。
2.4.2 单细胞染色体结构测序
肿瘤细胞的遗传物质不稳定，染色体结构可发生变化。基于染色质转座酶结合实验，scATAC-Seq技术利用Tn5转座酶对染色质开放区域进行测序[16]，可用于细胞间异质性分析，但该方法仅针对染色质开放结构进行测序。此外，基于染色体构象捕获的全基因组高通量Hi-C可用于检测染色体的三维构象，如染色质环和拓扑结构域[17]。目前，单细胞染色体结构测序普遍具有单次测序数据较少、基因组不能全覆盖的局限性，但上述两种主要的技术目前还无可替代。单细胞核小体结构和甲基化测序技术通过GpC甲基转移酶处理，可同时获得染色质开放区域、DNA甲基化和核小体相位的序列信息，还可筛选转录因子结合位点和脱氧核糖核酸酶高敏感性位点(DNase hypersensitive sites，DHSs)等有重要意义的序列[18]。
2.4.3 单细胞非编码RNA测序
机体细胞中存在大量对细胞功能有重要调节作用的非编码RNA。全RNA测序可获得mRNA及包括长链非编码RNA(long noncoding RNAs，lncRNAs)在内的非多聚腺苷酸RNA[non-poly(A) RNA]序列信息。对于lncRNAs，可根据其序列长度、外显子数量、序列进化保守性等特点从全RNA测序结果中筛选有用信息[19]。但使用寡脱氧胸苷酸(oligo-dT)进行全RNA测序对non-poly(A) RNA的覆盖率仍有待提高。随机置换扩增测序技术(random displacement amplification sequencing，RamDA-Seq)可应用于单细胞RNA全长测序。该技术将逆转录随机置换扩增技术和部分随机引物相结合，对non-poly(A) RNA有较高的灵敏度[20]。同时，另一项研究简化了单细胞RNA测序技术的长度筛选过程，结合信息学分析可获得如microRNA、tRNA片段和小核RNA(small nucleolar RNAs，snoRNAs)等短序列测序结果，且在人类胚胎干细胞中得到验证[21]。
2.4.4 单细胞多组学共测序
单细胞基因组、转录组及表观遗传组学等多组学同时测定，可获得多维度信息，在细胞间异质性、代谢特异性分析等方面有重要意义。目前，已有多个平台可以通过细胞质与细胞核分离、基因组和转录组分离的方法，同时进行单细胞全基因组、转录组和甲基化组测序[22]。
3. 单细胞测序技术在肿瘤中的研究进展
3.1 循环肿瘤细胞
肿瘤患者外周血中存在CTCs和CTCs簇。有关肿瘤转移的研究表明，肿瘤细胞可发生上皮间质转化(epithelial-to-mesenchymal transition, EMT)并侵入血管，在循环系统中存活的部分肿瘤细胞即为CTCs，其与肿瘤转移有关[23]。此外，医源性因素也可导致CTCs。从前列腺癌患者外周血中通过微液滴系统富集CTCs并进行单细胞RNA测序，结果显示患者自身和不同患者之间的CTCs转录组具有异质性，且发现耐药性相关的雄激素受体存在基因突变和剪切变异，该过程受上调的非经典Wnt信号通路调控[24]。相比于单个CTC细胞，CTCs簇具有更强的转移能力。一项研究对结肠癌、乳腺癌、胃癌和前列腺癌原发肿瘤CTCs进行分析，结果表明后者特殊基因的CNA变异显著提高，如Myc和第10号染色体缺失的磷酸酶张力蛋白同源物基因(gene of phosphate and tension homology deleted on chromsome ten, PTEN)，提示潜在的转移机制。且同一来源的CTCs之间显示出相对保守的CNA模式，同种肿瘤不同来源的CTCs之间以及原发肿瘤和CTCs之间SNV谱系存在较大差异，且SNV的发生具有更强的随机性[25]。由于存在扩增偏差和基因组测序覆盖程度的限制，CTCs的SNV研究发展需单细胞全基因组测序技术的进一步优化。
3.2 肿瘤进化及异质性
单细胞水平的研究可用于分析肿瘤内和肿瘤间的异质性，有助于肿瘤诊断、治疗，特别是对治疗过程中出现放化疗抵抗的机制研究有重要意义。利用nuc-Seq技术研究乳腺癌克隆进化可知，单倍体重排发生在肿瘤早期，在肿瘤细胞克隆扩增时染色体相对稳定；点突变则持续发生且与克隆异质性相关。同时，相比于雌激素受体阳性细胞，三阴性乳腺癌肿瘤细胞突变率更高[11]。利用肺腺癌患者肿瘤细胞进行异种移植，并对药物处理前后的异种移植肿瘤细胞进行单细胞mRNA测序, 可根据转录组将细胞分为不同的药物敏感性亚型[26]。依赖于多组学技术的发展，scTrio-Seq应用于人肝细胞肝癌基因组、甲基化组和转录组共测序，并根据SNV、DNA甲基化水平和转录组谱将被测细胞分为2个亚型，不同亚型具有不同的免疫识别潜能[27]。在单细胞基因组学和转录组学研究的基础上，表观遗传学测序结果可对基因组转录和翻译过程提供更多信息。
3.3 分子标志物筛选
既往黑色素瘤、胰腺肿瘤等相关研究提示，特征基因突变和基因重排谱是目前肿瘤诊断和预后评估的标志物之一[28-29]。针对肾透明细胞癌的单细胞外显子测序结果与既往利用整体测序技术得到的突变谱不同，且在被测序的25个细胞中未发现显著的克隆亚型[30]。另一项针对膀胱癌的研究中，对分离的膀胱癌干细胞、膀胱癌非干细胞、膀胱上皮细胞等类型的细胞进行单细胞DNA-Seq，结果显示膀胱癌干细胞存在克隆同源性，并发现了新的与肿瘤起始和细胞增殖相关的突变基因[31]。由于受到肿瘤异质性和测序成本的限制，被测序细胞的选择至关重要，因此，不仅需要发展高通量单细胞测序技术，还应建立用于单细胞检测的细胞分选标准。目前，已有研究利用单细胞RNA测序结果，构建能预测胶质瘤患者生存预后的转录组特征谱[32]。
3.4 肿瘤转移
实体肿瘤可发生远处转移和微转移。肿瘤细胞微转移DTCs通常出现在骨髓中，部分对治疗有抗性，可进展为肿瘤骨转移、内脏转移。利用单细胞测序技术可检出DTCs并定位其原发灶。从乳腺癌患者骨髓样品中分离得到单个细胞并进行全基因组测序，将序列信息中的CNA模式和原发肿瘤、淋巴结转移灶进行比对，可确定DTCs的原发肿瘤来源[33]。另一方面，对肿瘤原发灶和转移灶进行单细胞测序有助于发现肿瘤转移和定植相关的分子标志物及生态位特点。对头颈鳞癌原发灶及淋巴结转移灶进行单细胞全外显子组和转录组测序，结果显示原发灶细胞表达谱具有空间特异性，位于前缘的细胞部分上皮间质转化(partial epithelial-to-mesenchymal transition，p-EMT)水平较高，可作为肿瘤转移的预测标志[34]。
3.5 肿瘤微环境及免疫微环境
利用高通量单细胞测序技术可同时获得肿瘤微环境中的多种信息，有助于肿瘤生态位的研究。一项关于转移性黑色素瘤微环境的研究对患者肿瘤组织中的癌细胞、免疫细胞、间质细胞和内皮细胞进行单细胞RNA测序，结果显示肿瘤微环境中细胞间相互作用的特殊模式。针对肿瘤浸润T淋巴细胞的测序分析显示了该肿瘤组织中T淋巴细胞活化和克隆机制，为靶向治疗和免疫治疗提供了新思路[35]。肿瘤浸润T淋巴细胞在免疫治疗中有重要作用，从肝细胞癌患者血液、癌组织及癌旁组织中分离T淋巴细胞并进行单细胞RNA测序，可根据其分子及功能特性分型，结果显示CD8+T细胞和调节T细胞的特异性富集，测序结果亦提示了与之相关的靶基因[36]。此外，黑色素瘤中免疫细胞的测序结果显示，单核细胞和树突状细胞中干扰素γ诱导的免疫稳态调节分子的表达谱明显变化，导致抗肿瘤免疫反应被抑制，且与肿瘤免疫逃逸相关[37]。而对于脑转移瘤动物模型的单细胞RNA测序结果则提示脑肿瘤免疫微环境中小胶质细胞、骨髓来源巨噬细胞在肿瘤脑转移、脑转移瘤治疗抵抗中起重要作用[38]。
4. 脑恶性肿瘤诊治发展概述
胶质瘤约占原发脑肿瘤的27%，其中，胶质母细胞瘤是最常见的中枢神经系统原发恶性肿瘤。目前，胶质瘤的诊断主要依据影像学和组织病理学，治疗方案主要为手术结合放化疗。异柠檬酸脱氢酶(isocitrate dehydrogenase, IDH1)突变、O6-甲基鸟嘌呤-DNA甲基转移酶(O6-methylguanine-DNA methyltransferase, MGMT)甲基化、TERT启动子突变、1p19q缺失等基因组变异特点的发现是胶质瘤诊治进程中的里程碑，可指导肿瘤治疗，且与患者预后相关。世界卫生组织(World Health Organization, WHO)于2016年发布的中枢神经系统肿瘤分类指南首次将这些分子分型与组织病理整合[39]。此外，肿瘤相关抗原表皮生长因子受体(epidermal growth factorre-ceptor, EGFR)及GBM肿瘤特异性抗原EGFRvⅢ突变与胶质瘤的预后相关[40]，是靶向治疗、免疫治疗(如抗原肽疫苗、过继免疫治疗等)的重要靶点。在胶质瘤的诊疗过程中，除影像学及临床进展评估外，利用血液、脑脊液中细胞学、DNA检测或可为临床进展、肿瘤复发提供辅助信息[41]。目前，胶质瘤治疗方案的优化仍需深入进行基础研究和临床试验。
脑转移瘤是最常见的中枢神经系统恶性肿瘤，发病率约为原发恶性肿瘤的3~10倍，患者预后较差。对于首诊脑转移瘤患者，手术切除及放疗是主要的治疗方法。对于原发肿瘤化疗敏感的患者，也可尝试应用传统化疗方案治疗[42]。部分患者，如EGFR突变或间变性淋巴瘤激酶重排的非小细胞肺癌、人类表皮生长因子受体2阳性的乳腺癌、BRAF突变的黑色素瘤患者，发生脑转移后可能从靶向治疗中获益。近期一项关于酪氨酸激酶抑制剂埃克替尼治疗EGFR突变的非小细胞肺癌多发脑转移的Ⅲ期临床试验(BRAIN, NCT01724801)结果显示，患者无进展生存期显著延长[43]，此研究依赖于对原发肿瘤和脑转移瘤的组学研究。一项对86例对脑转移瘤及其配对原发肿瘤的单细胞测序研究表明，肿瘤发生脑转移后分子标志物发生变化，因此不能完全依赖外周肿瘤的标志物进行治疗决策，该研究亦分析了脑转移瘤在目前临床药物相关作用通路上的变异情况[44]，为脑转移瘤的诊治提供了新的思路。目前关于脑转移瘤的单细胞测序研究非常有限，仍有待进一步探讨。
5. 单细胞测序技术在脑肿瘤中的研究现状
5.1 特殊分子标志及信号通路
转录因子ZEB1在EMT、DNA损伤修复、免疫抑制等多种过程中发挥重要作用，但既往关于胶质瘤中ZEB1的功能研究结果存在争议。利用单细胞测序技术测定胶质瘤组织中不同类型细胞的ZEB1分子表达，结果提示其在肿瘤细胞、反应性胶质细胞和免疫细胞之间呈不同表达模式，且与EGFR倍增、IDH1突变有关[45]。利用单细胞全基因组测序技术检测GBM中EGFR的变异谱，可发现细胞之间EGFR的CNA不同，且存在相关其他类型突变，可能是肿瘤细胞增殖、克隆变异及产生治疗抗性的机制之一[46]。此外，利用单细胞测序技术根据转录组测序结果将肿瘤组织中的细胞进行分类，为进一步分析基因和肿瘤行为的关系奠定了基础。对复发GBM患者原发肿瘤及配对复发肿瘤进行单细胞RNA测序，结果提示RAS信号通路与GBM复发相关，且该结果在一项纳入3000例胶质瘤患者肿瘤组织RNA测序结果的荟萃分析中得到证实[47]。
5.2 胶质瘤的肿瘤异质性
大量实验证实表观遗传学修饰在肿瘤形成过程中发挥重要作用。一项研究对5例原发GBM和2个GBM肿瘤干细胞系的lncRNAs表达谱进行单细胞RNA测序，结果提示肿瘤间广泛存在lncRNAs的差异表达，但存在31个lncRNAs组成的表达谱可作为肿瘤干细胞特性的分子标志物[48]；此外，可根据lncRNAs表达谱将GBM细胞分为4种亚型，此分型方式与按照蛋白编码基因分型之间的一致性约为77%。对原发GBM细胞进行单细胞转录组测序，其结果显示癌基因信号、细胞增殖、补体-免疫应答及缺氧反应等代谢过程具有较强的异质性，可据此对GBM进行重新分类[49]。针对肿瘤内或肿瘤间异质性的研究可辅助治疗及预后评估。
5.3 胶质瘤的肿瘤进化
对IDH1/2突变的少突胶质瘤细胞进行单细胞RNA测序，结果显示其存在具有神经干细胞表达谱的细胞亚群，其中部分细胞高表达增殖相关基因，与全基因组CNA、点突变谱相一致，这部分细胞被认为是起始肿瘤增殖[50]。相比于整体测序，单细胞测序技术在展示肿瘤内细胞间异质性方面有较大优势，通过该技术可以更精确地分析肿瘤进化过程，并筛选有重要意义的突变形式。利用单细胞测序技术探究EGFR或PDGF突变的GBM，可以分析促血管生成因子在肿瘤发生、增殖过程中的作用，且发现PDGF受体编码基因中特定位点的缺失与前神经元性GBM增殖密切相关，提示其可能作为有意义的治疗靶点[51]。
5.4 胶质瘤的循环肿瘤细胞
目前，关于胶质瘤患者循环系统中是否存在CTCs尚存有争议，少量研究报道可在GBM患者外周血中分离得到CTCs。另一项研究针对肿瘤细胞端粒酶设计病毒载体的荧光蛋白以筛选肿瘤CTCs，利用单细胞测序技术可对该体系进行验证，实验结果提示胶质瘤患者外周血中存在CTCs，且高级别胶质瘤患者外周血中CTCs富集程度更强[52]，或可将CTCs富集程度作为一种预后相关的细胞标志。
5.5 胶质瘤肿瘤微环境和免疫应答
胶质瘤存在浸润生长的特点。针对GBM患者肿瘤中心细胞、边缘浸润的肿瘤细胞及不同位置免疫细胞的单细胞测序结果表明，其基因组和转录组具有异质性，位于肿瘤中心和边缘浸润的免疫细胞之间基因表达也具有特异性。此外，不同患者之间的浸润肿瘤细胞具有相对保守的基因突变及表达模式，这些基因可能与肿瘤浸润、侵袭及GBM复发相关[53]。对IDH突变的胶质瘤细胞进行单细胞转录组测序，并结合癌症基因组图谱(the cancer genome atlas, TCGA)数据库RNA谱相关数据进行分析，该研究结果显示胶质瘤免疫微环境中小胶质细胞和巨噬细胞存在富集现象[54]。此外，肿瘤免疫微环境的进展可辅助治疗相关研究。利用单细胞RNA测序技术发现胶质瘤中存在外周血来源的巨噬细胞，高表达具有免疫抑制功能的细胞因子，且在不同级别胶质瘤中的富集程度不同[55]，其可能与不同级别肿瘤恶性程度相关。
6. 未来脑恶性肿瘤诊疗应用展望
6.1 早期诊断及动态监测
目前已有部分研究报道了胶质瘤患者外周血中可分离、获得CTCs，依据胶质瘤进展、复发等过程中肿瘤组织、肿瘤干细胞及CTCs的分子标志物研究，可通过外周血中CTCs相关基因突变、表达及表观遗传学修饰提示疾病进程、治疗效果预估并监测是否复发。作为影像学检查的辅助工具，其时效性、预测的灵敏度及特异度仍有待进一步研究。CTCs富集在微创方面具有一定优势，CTCs单细胞测序或可替代组织活检，用于患者初步病理诊断并指导治疗。近期一项研究从GBM患者外周血中分离并鉴别CTCs，可作为胶质瘤复发的微创检测指标之一[56]。
此外，肺癌、乳腺癌、结直肠癌等可发生脑转移，一般预后不佳。脑转移瘤的早期检出可为治疗提供更多的选择。单细胞测序技术可对转移相关的CTCs及DTCs进行多组学测序。关于肺癌的研究发现，单细胞全基因组测序可检出CTCs中存在的特殊CNA模式，具有指示肿瘤转移的潜能[57]。此外，还可通过单细胞转录组、表观遗传组测序技术靶向CTCs或DTCs原发灶。CTCs可从患者的外周血中分离，创伤较小，且与原发肿瘤联系紧密。因此，可通过CTCs对肿瘤原发灶及转移进展的动态变化进行监测。由于CTCs在血液中富集程度较低，因此，需更高效的CTCs捕获和富集技术。为了更好地诊断和监测肿瘤，不仅需要更多关于肿瘤形成、增殖过程的多组学数据，还需对肿瘤原发灶、CTCs及肿瘤转移灶之间在分子、细胞水平上的关联和变异进行更深入的研究。
6.2 个体化治疗
单细胞测序的优势是发现肿瘤内和肿瘤间的异质性。基因组序列变异和特殊表观遗传修饰将导致细胞表达谱变化，这可能使部分细胞在治疗过程中产生耐药性，导致肿瘤细胞未被完全杀灭。利用单细胞测序技术可获得更全面的肿瘤分子特征，包括药物作用靶点、放化疗抗性标志物等，通过靶向特殊分子设计并完善治疗方案。结合针对CTCs的单细胞测序技术可实时监测肿瘤细胞异质性的变化，在治疗过程中针对其分子特征调节治疗方案。
单细胞测序技术还可对肿瘤微环境及免疫微环境进行分析，可能为临床免疫治疗提供新靶点和思路。
6.3 肿瘤分型和预后评估
相对于肿瘤组织整体测序，单细胞基因组、转录组和表观遗传组测序可获得更多特殊分子信息，结合相关分子的功能实验和临床数据，可筛选与治疗、预后相关的新型分子标志物，或针对现有标志物进行更精确的亚型分类。其次，多组学信息的整合可补充现阶段的肿瘤分子分型。此外，免疫微环境中免疫细胞的类型、数目，以及CTCs富集程度也可作为胶质瘤预后的标志物。对不同级别胶质瘤临床样本中肿瘤相关巨噬细胞(tumor-associated macrophages，TAM)进行单细胞测序分析表明，对于低级别胶质瘤，免疫微环境中高表达外周血细胞相关基因的TAM富集程度与患者预后相关[55]。
作者贡献：赵泽庆负责数据分析、论文初稿撰写；张莉、张怡璇、杨亚梅负责数据收集及论文写作指导；谷怡琳协助数据分析；陈适、王凤丹、潘慧负责论文指导、审核; 杨筱、李建初负责论文选题设计及写作指导。利益冲突：所有作者均声明不存在利益冲突 -
![]()
图 1 次级骨化中心呈从骨干向关节端下降的抛物线形(A)和半圆形(B)强回声带(UECT勾画示意图)
Dia(diaphysis)：骨干；Soc(secondary ossification center)：次级骨化中心；UECT(ultrasonic epiphyseal cartilage thickness)：超声骺软骨厚度
Figure 1. The schematic diagram depicts the way to measure UECT when secondary ossification center shows a parabolic strong echo zone descending from the diaphysis to the joint end (A) and a semicircular strong echo band(B)
![]()
图 2 次级骨化中心呈从骨干向关节端下降的抛物线形(A~G)和半圆形(a~g)强回声带(7个部位UECT超声图示)
A.第三掌骨头部背侧正中矢状面；B.尺骨茎突冠状面；C.桡骨茎突冠状面；D.股骨内上髁冠状面；E.股骨外上髁冠状面；F.胫骨背侧正中矢状面；G.胫骨内侧髁冠状面
Figure 2. Ultrasound diagram of UECT at 7 sites when the secondary ossification center shows a parabolic strong echo zone descending from the diaphysis to the joint end(A-G) and a semicircular strong echo band(a-g)
A.dorsal median sagittal plane of the third metacarpal head(MCP3); B.coronal plane of ulnar styloid process(Ulna); C.coronal plane of radial styloid process(Radius); D.coronal plane of medial epicondyle of femur(FM); E.coronal plane of lateral epicondyle of femur(FL); F.median sagittal plane of dorsal tibia(TD); G.coronal plane of medial condyle of tibia(TM)
UECT: 同图 1![]()
图 3 141名青少年年龄、身高及骨龄分布
A.年龄；B.身高；C.骨龄
Figure 3. Age, height and bone age distribution of 141 adolescents
A.age; B.height; C.bone age
![]()
图 4 不同骨龄男性青少年长骨关节端声像图(以尺骨茎突、桡骨茎突、股骨外上髁、胫骨内上髁为例)
Figure 4. Sonograms of the joint ends of different long bones in male adolescents of different bone ages (taking the ulnar styloid process, radial styloid process, femoral lateral epicondyle, and tibial medial epicondyle as examples)
表 1 男性与女性一般资料比较
Table 1 General information between male and female
项目 男性(n=80) 女性(n=61) P值 年龄[M(P25, P75)，岁] 13(12, 14) 13(12, 14) 0.52 骨龄[M(P25, P75)，岁] 15(13, 17) 16(14, 17) 0.40 骨龄-年龄＞1岁[n(%)] 48(60) 41(67) 0.48 身高[M(P25, P75)，cm] 175.3 (167.8, 181.4) 167.1 (161.4, 173.2) ＜0.001 
下载: 导出CSV
表 2 不同骨龄男性青少年7个部位的UECT测值[M(P25, P75)，mm]
Table 2 Measured values of UECT at 7 sites in male adolescents of different bone ages[M(P25, P75), mm]
骨龄(岁) MCP3 Ulna Radius FM FL TD TM 8(n=1) 1.4 4.2 5.4 5.0 2.6 3.0 4.7 9(n=1) 1.4 3.0 4.8 5.7 3.7 3.4 3.5 11(n=6) 1.7(1.3, 2.2) 3.2(2.7, 3.5) 4.2(4.0, 4.9) 3.5(2.4, 4.2) 2.5(2.2, 2.7) 3.6(2.8, 4.1) 4.0(3.5, 4.8) 12(n=3) 1.3(1.2, 1.9) 4.4(3.6, 5.2) 4.9(3.7, 5.5) 2.4(2.15, 3) 3.1(2.7, 3.8) 1.6(1.5, 2.1) 4.5(4.5, 4.8) 12.5(n=1) 1.0 2.6 3.1 2.1 2.4 2.3 2.8 13(n=11) 1.4(0.9, 1.5) 2.5(2.1, 2.9) 3.0(2.6, 3.3) 2.0(1.3, 2.6) 2.9(2.3, 3.0) 1.8(1.7, 2.0) 3.0(2.6, 3.2) 13.5(n=2) 1.5(1.2, 1.9) 3.0(2.9, 3.0) 3.1(3.0, 3.2) 1.9(1.9, 1.9) 1.7(1.6, 1.9) 3.0(2.4, 3.5) 2.7(2.5, 2.8) 14(n=7) 0.8(0.7, 0.9) 2.0(1.5, 2.1) 2.2(1.9, 2.9) 1.8(1.4, 2.0) 2.0(1.3, 2.3) 1.6(1.2, 2.2) 2.1(2.0, 3.3) 15(n=11) 0.6(0.6, 0.7) 1.9(1.5, 2.2) 2.2(1.6, 2.3) 1.5(1.1, 1.9) 1.7(0.9, 2.7) 1.2(0.7, 1.8) 2.2(1.3, 2.5) 16(n=4) 0.4(0.2, 0.4) 1.1(0.9, 1.4) 1.2(0.8, 1.7) 0.8(0.6, 1.0) 1.2(1.0, 1.3) 0.7(0.5, 0.8) 1.3(1.1, 1.4) 17(n=19) 0.1(0, 0.4) 1.3(0.7, 1.7) 1.3(1.1, 1.5) 0(0, 0.8) 0.9(0.6, 1.4) 0(0, 0.7) 0.6(0, 1.3) 18(n=14) 0(0, 0) 0.6(0.4, 1.0) 0.6(0.4, 0.9) 0(0, 0) 0(0, 0.4) 0(0, 0) 0(0, 0)
下载: 导出CSV
表 3 不同骨龄女性青少年7个部位的UECT测值[M(P25, P75)，mm]
Table 3 Measured values of UECT at 7 sites in female adolescents of different bone ages [M(P25, P75), mm]
骨龄(岁) MCP3 Ulna Radius FM FL TD TM 7(n=1) 1.7 4.7 6.4 3.1 3.8 1.5 2.3 8(n=1) 2.0 5.4 4.9 1.8 4.0 3.8 3.0 11(n=2) 1.1(1.0, 1.2) 2.4(2.0, 2.7) 3.4(2.6, 4.3) 2.5(2.2, 2.7) 2.2(1.7, 2.6) 1.9(1.8, 1.9) 1.7(1.4, 1.9) 11.5(n=1) 0.7 1.5 2.6 2.9 1.9 1.4 2.4 12(n=2) 0.9(0.8, 0.9) 1.9(1.9, 2.0) 2.1(2.0, 2.3) 1.5(1.4, 1.6) 1.7(1.7, 1.8) 0.9(0.9, 1.0) 4.4(4.0, 4.7) 13(n=7) 0.5(0.4, 0.6) 1.7(1.1, 1.8) 1.7(1.5, 1.9) 1.3(1.1, 2.0) 1.1(1.0, 1.3) 1.1(0.8, 1.4) 1.9(1.0, 2.7) 14(n=3) 0.3(0.3, 0.7) 1.3(1.3, 1.8) 1.3(1.2, 1.8) 1.4(0.7, 2.1) 1.0(1.0, 1.3) 0.6(0.3, 0.7) 1.3(1.2, 1.3) 15(n=4) 0.5(0.3, 0.5) 1.9(1.5, 2.1) 1.6(1.5, 1.7) 1.1(0.6, 1.6) 0.9(0.7, 1.4) 0.9(0.6, 1.3) 0.8(0.3, 1.5) 16(n=13) 0(0, 0) 0.9(0.9, 1.1) 0.8(0.6, 1.2) 0(0, 0.6) 0.9(0.5, 1.5) 0(0, 0) 0(0, 0.5) 17(n=19) 0(0, 0.2) 0.8(0.6, 1.1) 1.1(1.0, 1.4) 0(0, 0.7) 0.5(0, 1.0) 0(0, 0.2) 0(0, 0.2) 18(n=8) 0(0, 0) 0.4(0.2, 0.7) 0.6(0.4, 0.8) 0(0, 0) 0(0, 0.1) 0(0, 0) 0(0, 0)
下载: 导出CSV
表 4 单部位UECT与骨龄的相关系数(r)
Table 4 Correlation coefficient between single-site UECT and bone age
性别 MCP3 Ulna Radius FM FL TM TD 男 -0.84 -0.85 -0.90 -0.83 -0.80 -0.87 -0.86 女 -0.74 -0.71 -0.69 -0.67 -0.65 -0.75 -0.71 MCP3、Ulna、Radius、FM、FL、TD、TM：同图 2
下载: 导出CSV
表 5 7个部位UECT总和与骨龄及年龄的相关系数(r)
Table 5 Correlation between the sum of UECT at 7 sites and bone age or chronological age
性别 骨龄 年龄 男 -0.93 -0.79 女 -0.80 -0.64
下载: 导出CSV
表 6 多部位UECT总和与骨龄的相关系数
Table 6 Correlation coefficient between UECT sum of different sites and bone age
性别 FL+TM+TD FM+TM+TD knee joint hand and wrist all 男 -0.91 -0.92 -0.93 -0.90 -0.93 女 -0.80 -0.76 -0.79 -0.77 -0.80 FL、FM、TM、TD：同图 2；knee joint：膝关节4个部位UECT总和；hand and wrist：手和腕部3个部位UECT总和；all：7个部位UECT总和
下载: 导出CSV
表 7 UECT测量的一致性分析
Table 7 Consistency analysis of UECT measurements
部位 ICC 95% CI MCP3 0.91 0.79~0.96 Ulna 0.95 0.84~0.98 Radius 0.91 0.78~0.97 FM 0.98 0.95~0.99 FL 0.93 0.83~0.97 TD 0.93 0.83~0.97 TM 0.95 0.89~0.98 MCP3、Ulna、Radius、FM、FL、TD、TM：同图 2；ICC(intraclass correlation coefficient): 组内相关系数
下载: 导出CSV
-
[1] Gaudino R, De Filippo G, Bozzola E, et al. Current clinical management of constitutional delay of growth and puberty[J]. Ital J Pediatr, 2022, 48(1): 45. DOI: 10.1186/s13052-022-01242-5
[2] 王凤丹, 次旦旺久, 焦洋, 等. 人工智能在骨龄评估中的应用[J]. 基础医学与临床, 2022, 42(11): 1776-1780. Wang F D, Cidanwangjiu, Jiao Y, et al. Use of artificial intelligence in bone age assessment[J]. Basic Clin Med, 2022, 42(11): 1776-1780.
[3] 袁新宇, 潘诗农, 程晓光. 构建儿童骨龄核心数据建立中国儿童骨龄标准体系[J]. 中华放射学杂志, 2023, 57(4): 341-343. DOI: 10.3760/cma.j.cn112149-20230219-00114 Yuan X Y, Pan S N, Cheng X G. Establishing bone age database of contemporary healthy children and building the standard for assessment in China[J]. Chin J Radiol, 2023, 57(4): 341-343. DOI: 10.3760/cma.j.cn112149-20230219-00114
[4] 赵泽庆, 潘慧, 张莉, 等. 超声评估骨龄研究现状及临床应用前景[J]. 协和医学杂志, 2024, 15(2): 400-405. DOI: 10.12290/xhyxzz.2023-0270 Zhao Z Q, Pan H, Zhang L, et al. Research status and application prospect of bone age assessment by ultrasonography[J]. Med J PUMCH, 2024, 15(2): 400-405. DOI: 10.12290/xhyxzz.2023-0270
[5] Pazzaglia U E, Beluffi G, Benetti A, et al. A review of the actual knowledge of the processes governing growth and development of long bones[J]. Fetal Pediatr Pathol, 2011, 30(3): 199-208. DOI: 10.3109/15513815.2010.524693
[6] Rivas R, Shapiro F. Structural stages in the development of the long bones and epiphyses: a study in the New Zealand white rabbit[J]. J Bone Joint Surg Am, 2002, 84(1): 85-100. DOI: 10.2106/00004623-200201000-00013
[7] Mukaka M M. Statistics corner: a guide to appropriate use of correlation coefficient in medical research[J]. Malawi Med J, 2012, 24(3): 69-71.
[8] Ekizoglu O, Er A, Buyuktoka A D, et al. Ultrasonographic assessment of ossification of the distal radial epiphysis for estimating forensic age[J]. Int J Legal Med, 2021, 135(4): 1573-1580. DOI: 10.1007/s00414-021-02521-2
[9] Benito M, Muñoz A, Beltrán I, et al. Assessment of adulthood in the living Spanish population based on ossification of the medial clavicle epiphysis using ultrasound methods[J]. Forensic Sci Int, 2018, 284: 161-166. DOI: 10.1016/j.forsciint.2018.01.007
[10] Torenek Ağirman K, Bilge O M, Miloğlu Ö. Ultrasono-graphy in determining pubertal growth and bone age[J]. Dentomaxillofac Radiol, 2018, 47(7): 20170398. DOI: 10.1259/dmfr.20170398
[11] Castriota-Scanderbeg A, De Micheli V. Ultrasound of femoral head cartilage: a new method of assessing bone age[J]. Skeletal Radiol, 1995, 24(3): 197-200. DOI: 10.1007/BF00228922
[12] Castriota-Scanderbeg A, Sacco M C, Emberti-Gialloreti L, et al. Skeletal age assessment in children and young adults: comparison between a newly developed sonographic method and conventional methods[J]. Skeletal Radiol, 1998, 27(5): 271-277. DOI: 10.1007/s002560050380
[13] Sherif H, Noureldin M, Bakr A F, et al. Sonographic measurement of calcaneal volume for determination of skeletal age in children[J]. J Clin Ultrasound, 2003, 31(9): 457-460. DOI: 10.1002/jcu.10213
[14] Windschall D, Pommerenke M, Haase R. Ultrasound assessment of the skeletal development of the proximal tibial, proximal femoral, and distal femoral epiphyses in premature and mature newborns[J]. Ultrasound Med Biol, 2016, 42(2): 451-458. DOI: 10.1016/j.ultrasmedbio.2015.10.020
[15] Wan J, Zhao Y, Feng Q Q, et al. Statistical confirmation of a method of US determination of bone age[J]. Radiology, 2021, 300(1): 176-183. DOI: 10.1148/radiol.2021204353
[16] Wan J, Zhao Y, Feng Q Q, et al. Potential value of conventional ultrasound in estimation of bone age in patients from birth to near adulthood[J]. Ultrasound Med Biol, 2019, 45(11): 2878-2886. DOI: 10.1016/j.ultrasmedbio.2019.07.681
[17] Wan J, Zhao Y, Feng Q Q, et al. Summation of ossification ratios of radius, ulna and femur: a new parameter to evaluate bone age by ultrasound[J]. Ultrasound Med Biol, 2020, 46(7): 1761-1768. DOI: 10.1016/j.ultrasmedbio.2020.03.021
[18] Roth J, Jousse-Joulin S, Magni-Manzoni S, et al. Definitions for the sonographic features of joints in healthy children[J]. Arthritis Care Res (Hoboken), 2015, 67(1): 136-142. DOI: 10.1002/acr.22410
[19] Benedick A J, Hogue B, Furdock R J, et al. Estimating skeletal maturity using knee radiographs during preadolescence: the epiphyseal: metaphyseal ratio[J]. J Pediatr Orthop, 2021, 41(9): 566-570. DOI: 10.1097/BPO.0000000000001921
[20] Huang L F, Furdock R J, Uli N, et al. Estimating skeletal maturity using wrist radiographs during preadolescence: the epiphyseal: metaphyseal ratio[J]. J Pediatr Orthop, 2022, 42(7): e801-e805. DOI: 10.1097/BPO.0000000000002174
-
期刊类型引用(5)
1. 郭文杰，刘芳，高娜，陈彻，楚惠媛. 肿瘤细胞周期调控研究进展. 中国免疫学杂志. 2024(06): 1315-1319 . 
百度学术
2. 刘远亭，韩淑萍，刘沙，杨证玉，何俊超，张燕琴. 基于网络药理学探讨长春新碱治疗乳腺癌的作用机制. 食品与药品. 2023(06): 567-575 . 百度学术
3. 姜惠心，祁子泰，高修歌，郭大伟，季辉，彭麟，江善祥. 犬乳腺肿瘤化学治疗药物的研究进展. 畜牧与兽医. 2021(12): 116-122 . 百度学术
4. 毛倩倩，逯遥，黄彬，林久茂. 片仔癀对人结肠癌细胞HCT-116裸鼠皮下移植瘤的抑制作用. 福建中医药. 2021(12): 15-19 . 百度学术
5. 焦秋香，陈新艳. miR-198在卵巢癌中的表达及对卵巢癌细胞增殖作用的机制研究. 标记免疫分析与临床. 2021(12): 2103-2108+2151 . 百度学术
其他类型引用(2)
计量
- 文章访问数: 671
- HTML全文浏览量: 2041
- PDF下载量: 41
- 被引次数: 7
