骨成熟度的评估对于生长发育异常患儿的诊断、治疗及随访具有重要作用^[1]。目前临床主要通过拍摄左手及腕部X线片评估骨龄以了解骨成熟度，但其具有辐射风险，且存在来源人群陈旧、一致性低、复杂费时等问题^[2-3]。而超声可清晰显示发育期低回声骺软骨及其内强回声的次级骨化中心并进行相应测量，且无辐射，可用于动态随诊^[4]。基于发育过程中长骨关节端的大体及微观特点^[5-6]，本研究结合长骨关节端声像表现，初步探索采用超声骺软骨厚度（ UECT）定量评估骺软骨发育情况并分析其与骨龄的相关性，以期为骨成熟度评估提供新思路。

1 资料与方法

1.1 研究对象

本研究为前瞻性观察性研究，以2023年3—6月就读于济宁市某体校的青少年为研究对象，记录年龄、身高等信息。

纳入标准：（1）年龄为7~18岁；（2）外形及营养状况良好。

排除标准：（1）患有可能导致生长发育异常的慢性肝肾疾病、骨骼系统疾病、内分泌遗传综合征等；（2）超声或骨龄X线图像质量差、无法识别；（3）手和腕部及膝关节组成骨存在近期骨折史。

本研究已通过北京协和医院伦理审查委员会审批（审批号：K4393），研究对象或其亲属均签署知情同意书。

1.2 研究方法

1.2.1 X线检查

于1周内行X线骨龄检查，并由1名具有15年阅读骨龄片经验的内分泌科医师使用Greulich-Pyle（GP）图谱法评估骨龄。

1.2.2 超声检查

分别由具有20年及8年超声工作经验的2名医师针对非优势侧手、腕部及膝关节7个组成部位进行超声检查（GE Logiq E20彩色多普勒超声诊断仪），探头选用L6-24靴型探头及L2-9线阵探头。

首先使患者手背朝上，采用L6-24靴型探头顺序扫查第三掌骨头部背侧正中矢状面、尺骨茎突冠状面、桡骨茎突冠状面。

随后使患者采取俯卧位，暴露膝关节背侧，采用L2-9线阵探头扫查股骨内上髁冠状面、股骨外上髁冠状面、胫骨背侧正中矢状面、胫骨内侧髁冠状面。所有切面均清晰显示干骺端、次级骨化中心及骺软骨。

1.2.3 UECT测量

首先沿骨干强回声带外缘内侧作骨干方向延长线，其与骨干强回声带外缘产生交点，随后根据次级骨化中心所示强回声带的形状选择相应测量方法：

（1）若呈从骨干向关节端下降的抛物线形，则以次级骨化中心骨干端外缘顶点与前述骨干外缘交点的距离作为UECT；

（2）若呈半圆形，作与前述骨干方向延长线的平行线，其中与次级骨化中心强回声带相切的切点与前述骨干外缘交点距离为UECT（图1、2）。

对于单个研究对象的多部位UECT联合分析，则将各部位UECT之和作为总UECT并分析其与骨龄的相关性。

图1 次级骨化中心呈从骨干向关节端下降的抛物线形(A)和半圆形（B)强回声带（UECT勾画示意图）

Dia（diaphysis）：骨干；Soc（secondary ossification center）：次级骨化中心；UECT（ultrasonic epiphyseal cartilage thickness）：超声骺软骨厚度

图2 次级骨化中心呈从骨干向关节端下降的抛物线形(A~G)和半圆形（a~g)强回声带（7个部位UECT超声图示）

A.第三掌骨头部背侧正中矢状面；B.尺骨茎突冠状面；C.桡骨茎突冠状面；D.股骨内上髁冠状面；E.股骨外上髁冠状面；F.胫骨背侧正中矢状面；G.胫骨内侧髁冠状面

1.2.4 重复测量信度检验

以组内相关系数（ICC）评估UECT测量的信度与可重复性。研究按照性别与骨龄采用分层抽样法随机抽取50人超声图像，由另一操作者再次测量UECT值，计算ICC并评估一致性。

1.3 样本量估算

采用G*power 3.1.9.7软件分析纳入的样本量，参数选择效应量f2=0.15、α=0.05，β=0.80，计算所需最小样本量为55例。

1.4 统计学处理

采用R Studio 4.3.0软件进行统计学分析。正态性分析采用Shapiro-Wilk检验，年龄、骨龄及身高均不符合正态分布，以中位数（四分位数）表示。组间比较采用秩和检验。计数资料以频数（百分数）表示，组间比较采用卡方检验。采用Spearman法分析UECT与骨龄的相关性，高低分级参考前述研究^[7]，|r|为0.5~0.7表示中度相关；|r|为0.7~0.9表示高度相关；|r|为0.9~1.0表示非常高度相关。采用SPSS 24.0软件进行ICC分析，计算参数选择single measurement、2-way random effects、absolute agreement，ICC≥0.75为一致性较好，ICC＜0.4为一致性较差。以P＜0.05为差异具有统计学意义。

2 结果

2.1 一般资料

共纳入141名青少年，其中男性80名，女性61名。中位年龄为13（12，14）岁。男性身高显著高于女性（表1）。

表1 男性与女性一般资料比较

男性中位骨龄15(13，17)岁，女性中位骨龄16（14，17）岁。男女骨龄与年龄差值＞1岁者分别占比60%、67%。男女年龄、身高及骨龄分布见图3。

图3 141名青少年年龄、身高及骨龄分布

A.年龄；B.身高；C.骨龄

2.2 发育过程长骨关节端的声像表现

利用超声观察不同骨龄阶段研究对象的长骨关节端，其声像表现存在以下规律：生长发育初期，次级骨化中心首先在低回声骨骺中央萌出，声像图显示为小半圆形强回声带。次级骨化中心进一步扩展增大，但仍然保持半圆形的形状。随着骨成熟度提高，次级骨化中心所示强回声带呈从骨干向关节端下降的似抛物线形，强回声带的骨干端从低于骨干外缘逐渐上升至与其平齐。生长晚期至终止时，次级骨化中心与骨干不断接近，最终形成一条完整、平滑的强回声带，关节面被覆薄层低回声关节软骨（图4）。

图4 不同骨龄男性青少年长骨关节端声像图（以尺骨茎突、桡骨茎突、股骨外上髁、胫骨内上髁为例）

2.3 UECT与骨龄相关性

男、女7个部位的UECT均随骨龄增长而减小，至骨龄18岁时测量值趋近于0（表2、3）。

表2 不同骨龄男性青少年7个部位的UECT测值 ［M(P₂₅, P₇₅)，mm］

表3 不同骨龄女性青少年7个部位的UECT测值 ［M(P₂₅, P₇₅)，mm］

各部位UECT与骨龄均呈负相关，男性单部位UECT与骨龄均呈高度负相关(|r|≥0.83)，而女性单部位UECT与骨龄呈中高度负相关（0.65≤|r|≤0.75）（表4）。

表4 单部位UECT与骨龄的相关系数(r)

男、女7个部位UECT总和与骨龄的负相关性均达到较高水平（男：|r|=0.93；女：|r|=0.80）（表5）, 单独选取手和腕部3个部位或膝关节4个部位时，女性UECT总和与骨龄高度负相关（|r|≥0.76），而男性则表现出非常高度负相关性（|r|≥0.90）(表6)。

表5 7个部位UECT总和与骨龄及年龄的相关系数（r)

表6 多部位UECT总和与骨龄的相关系数

2.4 不同操作者测量UECT的可重复性

7个部位UECT测量的一致性分析均具有显著统计学意义（P均＜0.001），ICC值均≥0.75提示各部位UECT测量一致性均较好（表7）。

表7 UECT测量的一致性分析

3 讨论

本研究显示，UECT与骨龄呈中至高度负相关，提示UECT可一定程度上反映发育过程中骺软骨的变化，具有较好的可重复性，可用于骨成熟度水平定量评估。声像图显示发育过程中长骨次级骨化中心形态变化及其与骨干的位置关系改变，与大体及微观研究结果相符^[5-6]。

与其他利用骨骺生长发育过程中声像改变划分阶段推测法医学年龄的研究不同^[8-9]，本研究以UECT连续定量评估骨成熟度，而Torenek等^[10]根据手和腕部的所有指骨和桡骨远端骨骺以及第一掌骨头部的籽骨共13个部位的声像表现构建半定量评估系统，结果发现除第一、二指骨近节指骨外，超声与X线的评分结果无统计学差异（P＜0.05）。

本研究发现，UECT与骨龄相关性在不同性别中存在差异，男性单个及多个部位均表现为高度负相关，而女性多为中度负相关，分析原因可能是大骨龄女性比例较高（16~18岁骨龄的女性占比67%），且女性的发育起点和终点均早于男性，使女性UECT值多数较低，未来可通过控制研究群体的分布开展进一步研究。

多项研究通过测量发育过程中长骨关节端结构进行骨成熟度评估，Castriota-Scanderbeg等^[11-12]研究发现，无骨龄提前或延后的研究对象股骨头软骨超声测量厚度与年龄、骨龄、身高等均高度相关，但对于发育异常对象而言，该方法与放射学评估的一致性较低。本研究显示，男女UECT与骨龄均呈负相关性，且明显高于其与年龄的相关性（表5），说明UECT可反映研究对象的发育提前特点。

另有多个研究聚焦于骨关节端的次级骨化中心研究。Sherif等^[13]利用超声测量正常儿童的根骨体积并以此作为标准准确预测了验证组健康儿童的年龄，但从出生后第5年开始根骨的体积增长率仅为12%左右，因此其应用范围较有限，也缺乏在发育异常患儿中的验证。

Windschall等^[14]利用超声测量早产儿与足月儿的股骨远端及胫骨近端次级骨化中心大小并评估其与胎龄的相关性，发现存在中度相关（R²=0.35~0.50）。

Wan等^[15-17]将骨化中心最大直径与无回声骨骺最大直径比值定义为骨化率，发现桡骨茎突冠状面、尺骨茎突冠状面、股骨内上髁冠状面三者的骨化率总和与骨龄呈正相关（r=0.97），对异常骨龄具有非常高的诊断性与特异性（灵敏度≥93%，特异度≥98%）。

骨化率反映了骺软骨骨化不断扩大并最终发育成骨关节末端的过程，而UECT则反映骨化中心不断接近骨干端至最后愈合的过程。前者是比值型指标，受身高影响可能更小，但关节端的骺软骨与发育成熟后保留于骨表面的薄层关节软骨均表现为同质性的低回声^[18]，常难以获得准确的骨骺最大直径。

有研究采用手和腕部以及膝关节X线片测量数个长骨次级骨化中心以及骨干的横径比来评估骨成熟度，结果均表现为高度相关^[19-20]，而本研究声像图显示，发育成熟过程中长骨的次级骨化中心往往先上升至与骨干外缘同一高度，随后表现为二者间骺软骨不断变薄至消失，因此前述方法可能对生长发育晚期的变化描述有限。

本研究局限性：（1）单中心研究；（2）研究对象多处于发育晚期，骨成熟度呈偏态分布；（3）对次级骨化中心的研究不足，可确定次级骨化中心形成似抛物线形及其骨干端与骨干外缘相平齐所对应的成熟度水平，进一步划分发育阶段；（4）缺乏与放射学方法评估骨龄的一致性评价。

综上所述，超声可观察描述长骨关节端的生长发育变化，UECT与骨龄表现出中至高度负相关，具有用于骨成熟度定量评估的潜在价值。

参考文献

［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.

［2］王凤丹, 次旦旺久, 焦洋, 等. 人工智能在骨龄评估中的应用［J］. 基础医学与临床, 2022, 42(11): 1776-1780.

［3］袁新宇, 潘诗农, 程晓光. 构建儿童骨龄核心数据 建立中国儿童骨龄标准体系［J］. 中华放射学杂志, 2023, 57(4): 341-343.

［4］赵泽庆, 潘慧, 张莉, 等. 超声评估骨龄研究现状及临床应用前景［J］. 协和医学杂志, 2024, 15(2): 400-405.

［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.

［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.

［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.

［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.

［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.

［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.

［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.

［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.

［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.

［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.

［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.

［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.

［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.

［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.

［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.