使用口腔颌面锥形束CT图像通过第二恒磨牙髓室增龄性变化进行年龄推断的研究

人类活体或尸体的年龄推断是法医学中一项重要的内容，精确的年龄推断方法在实践中有着广泛的应用需求[1-2]。牙齿发育完成后，随着年龄的增长，继发性牙本质沉积导致牙齿髓腔逐渐变小[3-5]。本研究通过口腔颌面锥形束CT(CBCT)图像，研究通过第二恒磨牙髓室的增龄性变化进行年龄推断的方法，为法医学中通过牙齿进行年龄推断提供参考依据。

1.新生到校报名后，即开始接受入学第一课——入学教育。包括在全校、学院开展的各种层次的、丰富多彩的活动，如军训、专题讲座、主题班会等。可以邀请优秀毕业生（小学的骨干教师）回校畅谈感受，上一堂公开课等，可以邀请实践基地的校领导做专题讲座，还可以选派本院教师介绍本专业培养目标、课程设置和学习方法，等等。学生通过入学教育，间接体验，对教师教育的培养有了初步了解，对本专业的学习有了远景目标。

P2P网络借贷是指通过互联网将资金从宽裕的借款人手中借贷给有资金需求的人，从而产生的无抵押小额贷款模式。这是一种来自欧美的新兴模式。P2P网络借贷平台作为一种中介，将资金聚集起来的同时还要负责对借款人的借款资质、贷款使用情况、还款条件等方面进行评定。随着“互联网+”时代的到来，P2P企业迅猛发展，得到了很多创业者的相应并愿意将其作为投资理财产品，这无疑为大学生创业融资提供帮助。P2P网络借贷平台逐步开始推出大学生创业贷款方式，渠道更加便捷。通过限制贷款最高限额，降低投资风险的同时满足创业资金需求。

3.科学设置服务内容。梳理出涉及民政、人力社保、国土资源、建设、规划、农办、农业、卫生、人口计生等部门70项服务事项，由村级代理村民去诸暨市公共服务中心或乡镇公共服务中心代办。村便民服务中心还专门在窗口设立了服务指南，把需要由村民自行提供证件和资料的10多项内容挂在墙上。按照当即可办的事项即受即办、既定代办的事项实行全程代办、不能办理的事项说明理由的原则，向群众提供信息服务、法律保障、政策引导、帮扶弱势群体等服务内容。

1 材料和方法

1.1 材料来源

收集青岛市口腔医院放射科数据库内15～69岁的110例男性和110例女性的220颗上颌第二恒磨牙和220颗下颌第二恒磨牙的CBCT图像。全部CBCT图像均因病人临床方面需求而拍摄。病人的出生日期由医院信息系统确认，图像拍摄日期和出生日期之差记为病人的真实年龄。纳入病人按年龄分为6组：15～20岁、21～30岁、31～40岁、41～50岁、51～60岁和61～69岁组，每一岁年龄均包括2例男性和2例女性病人，每例病人需包含至少一颗符合纳入标准的上颌第二恒磨牙以及下颌第二恒磨牙。

许多学者利用医学影像学的手段对通过牙齿髓腔大小的变化进行了年龄推断的研究[8-12]，但对于通过第二恒磨牙的髓腔变化进行年龄推断的研究很少[13]。GE等[13]通过分割上、下颌第二恒磨牙髓室的CBCT图像计算得出了上、下颌第二恒磨牙的髓室体积与年龄相关性的判定系数R2分别为0.498和0.487，高于本研究所得上、下颌第二恒磨牙髓室高度、髓室颊舌向宽度、牙本质厚度分别与年龄的相关性判定系数。这可能是由于GE等[13]所测指标为髓室体积值，而本研究所测量的是牙齿髓室在某一断层图像上的高度、宽度、厚度，并不能完全反应牙齿髓室的真实改变。刘彦等[14]研究了汉族人下颌第一恒磨牙髓室高度及牙本质厚度与年龄的关系，得出下颌第一恒磨牙髓室高度与年龄之间的R2为0.384，牙本质厚度与年龄的R2为0.026，与本研究中所得的上颌第二恒磨牙髓室高度与年龄的R2为0.296、上颌第二恒磨牙牙本质厚度与年龄的R2为0.053、下颌第二恒磨牙髓室高度与年龄的R2为0.363、下颌第二恒磨牙牙本质厚度与年龄的R2为0.085较为接近。

1.2 图像获取和数据测量

非线性回归分析显示，上颌第二恒磨牙髓室高度、髓室颊舌向宽度、牙本质厚度与年龄的对数回归模型判定系数R2=0.055～0.296；下颌第二恒磨牙髓室高度、髓室颊舌向宽度、牙本质厚度与年龄的对数回归模型判定系数R2=0.050～0.362。所有对数回归方程均具有统计学意义(P<0.05)。

各年龄组上颌第二恒磨牙及下颌第二恒磨牙髓室高度、髓室颊舌向宽度及牙本质厚度测量结果见表1。线性回归分析显示，上颌第二恒磨牙髓室高度与年龄线性回归模型判定系数R2=0.296(P<0.05),直线回归方程y=64.575-11.823x，标准估计误差值为13.38岁；上颌第二恒磨牙髓室颊舌向宽度与年龄的线性回归模型判定系数R2=0.132(P<0.05),直线回归方程y=78.680-7.973x，标准估计误差值为14.85岁；上颌第二恒磨牙牙本质厚度与年龄的线性回归模型判定系数R2=0.053(P<0.05),直线回归方程y=19.633+7.679x，标准估计误差值为15.52岁；下颌第二恒磨牙髓室高度与年龄的线性回归模型判定系数R2=0.363(P<0.05),直线回归方程y=66.388-15.61x，标准估计误差值为12.73岁；下颌第二恒磨牙髓室颊舌向宽度与年龄的线性回归模型判定系数R2=0.044(P<0.05),直线回归方程y=60.708-6.085x，标准估计误差值为15.59岁；下颌第二恒磨牙牙本质的厚度与年龄的线性回归模型判定系数R2=0.085(P<0.05),直线回归方程y=13.461+8.936x，标准估计误差值为15.26岁。

1.3 统计分析

采用SPSS 22.0(Chicago, IL, USA)软件进行统计分析。计量资料以表示；髓室高度、髓室颊舌向宽度、牙本质厚度与年龄之间的相关性采用线性回归分析和对数回归分析；通过髓室高度、髓室颊舌向宽度及牙本质厚度进行年龄推断数学模型的建立采用多元线性回归分析，以上、下颌第二恒磨牙的髓室高度、髓室颊舌向宽度、牙本质厚度作为自变量，以年龄值作为因变量。为了得到最优化及最显著的预测年龄的指标组合，在多元回归分析中对上颌及下颌的第二恒磨牙的髓室高度、髓室颊舌向宽度、牙本质厚度这6个指标采用后退法逐步筛选，将获得最大P值(当P>0.1时)的指标逐一进行排除，直到剩余指标的最大P值小于0.1时即意味着回归方程具有统计学显著意义；组内一致性及组间一致性检验采用配对t检验，以P<0.05为差异有统计学意义。

2 结 果

对2名研究人员进行测量统一性训练，达到标准化之后，由1名研究人员进行测量采集，每个数据测量2次，取平均值。在全部数据测量完成后2周，随机抽取20例病人的20颗上颌第二恒磨牙和20颗下颌第二恒磨牙图像，由2名研究人员再次进行测量，以验证组间一致性和组内一致性。

所有CBCT图像均由NewTom VG(Quantitative Radiology, Verona, Italy)CBCT机拍摄。拍摄时曝光参数为110 kVp，3.69～84.10 mAs，扫描视野15 cm×15 cm大小，体素值0.15 mm。图像采用DICOM格式导出以后导入到RadiAnt DICOM Viewer软件中进行测量。如图1所示，分别在轴位、矢状位和冠状位视窗调整图像十字定位光标与牙体长轴相一致，在冠状位视窗中将十字定位光标与第二恒磨牙髓室冠根向和颊舌向的中点相重叠，定位线与髓室底交于A点，与髓室顶交于B点，与釉牙本质界交于C点，与髓室颊侧壁交于D点，与髓室舌侧壁交于E点。将AB点间距离定义为髓室高度，BC点间距离定义为牙本质厚度，DE点间距离定义为髓室颊舌向宽度。

 

①为矢状位视窗；②为轴位视窗；③为冠状位视窗。

 

图1 测量示意图

 

表1 各年龄组上下颌第二恒磨牙髓室高度、髓室颊舌向宽度及牙本质厚度比较

  

组别n上颌第二恒磨牙髓室高度(h/mm)髓室颊舌向宽度(b/mm)牙本质厚度(d/mm)下颌第二恒磨牙髓室高度(h/mm)髓室颊舌向宽度(b/mm)牙本质厚度(d/mm)15～20岁242.810±0.6545.100±0.4422.785±0.4732.162±0.4283.339±0.4043.029±0.44121～30岁402.124±0.5764.873±0.6502.917±0.5201.866±0.4853.094±0.3673.176±0.52431～40岁401.848±0.6744.475±0.6882.912±0.4501.599±0.6752.980±0.4923.114±0.45541～50岁401.821±0.6194.402±0.6842.919±0.4421.509±0.5942.865±0.6243.201±0.37751～60岁401.552±0.5784.314±0.6133.161±0.4451.250±0.4382.862±0.6623.324±0.60861～69岁361.378±0.6024.263±0.8503.104±0.4610.997±0.2692.955±0.5653.589±0.507总计2201.868±0.7324.538±0.7252.977±0.4761.531±0.6142.993±0.5513.249±0.518

花五奇轻功高妙，武林中人尽皆知，看来乔十二郎毫无退路。老太医情急之下，来不及思考，只想扑上去护住乔十二郎。他刚来得及像护崽的老母鸡那样张开两翼，却见乔十二郎伏低身子往前一窜，居然从敌方胯下钻过，将敌人从背上掀过去。花五奇轻功好，不至于摔跤，等他站稳，已被老太医紧紧拉住。此时，跌伏于花五奇身后的乔十二郎迅速爬起，跑向场中，与赶来的秦铁崖汇合。

本研究结果值得注意的是，上颌第二恒磨牙髓室高度与年龄的相关性低于下颌第二恒磨牙髓室高度与年龄的相关性，而上颌第二恒磨牙髓室颊舌侧宽度与年龄的相关性则高于下颌第二恒磨牙髓室颊舌侧宽度与年龄的相关性，这与我们在研究过程中观察到的随着年龄的增长，上颌第二恒磨牙髓室颊舌向宽度变化更为明显，而下颌第二恒磨牙髓室高度变化更为明显相一致，这可能提示了上、下颌第二恒磨牙髓室继发性牙本质沉积规律的不同。

3 讨 论

第二恒磨牙的纳入标准：无龋坏、无过度磨耗、无修复性充填或治疗、无金属伪影干扰、无牙髓腔钙化。其中，将“无过度磨耗”定义为50岁前牙齿磨耗指数(TWI)≤2，50岁后TWI≤3。具体标准参见相关文献[6-7]。

20颗上颌第二恒磨牙和20颗下颌第二恒磨牙的髓室高度、髓室颊舌向宽度、牙本质厚度的组间一致性分析和组内一致性分析均显示差异无显著性(P>0.05)。

多元线性回归分析显示，上颌第二恒磨牙髓室高度、上颌第二恒磨牙髓室颊舌向宽度、下颌第二恒磨牙髓室高度、下颌第二恒磨牙牙本质厚度四个指标是进行年龄推断的最优组合，通过这4个指标进行年龄推断的公式为：预测年龄=80.646-6.008×上颌第二恒磨牙髓室高度-5.961×上颌第二恒磨牙髓室颊舌向宽度-9.408×下颌第二恒磨牙髓室高度+4.469×下颌第二恒磨牙牙本质厚度，模型的拟合度R2=0.51，模型标准估计误差值为11.24岁。

有研究显示，老年人继发性牙本质沉积较慢，而青年人则较快[15-17]，提示对数回归模型可能更加符合人类髓腔大小增龄性变化的规律。GE等[15]对上、下颌第一恒磨牙髓室体积与年龄的相关性进行了研究，得出了在通过第一恒磨牙髓室体积值进行年龄推断时，对数回归模型优于线性回归模型的结论。不同于GE等[15]的研究结果，本研究结果显示，在上、下颌第二恒磨牙髓室高度、髓室颊舌向宽度、牙本质厚度与年龄的的回归分析中，对数模型的拟合度均不优于线性模型。

本研究结果表明，通过多元线性回归分析，使用多个指标共同预测年龄比仅通过一个指标预测年龄，准确性更高。这与GE等[13]的研究结果一致。

另外，了解髓室的增龄性变化对口腔医师进行牙体牙髓病的治疗具有一定的指导作用。髓室底的医源性穿孔是牙体牙髓病临床治疗中较为常见的并发症之一。了解髓室高度、宽度及牙本质厚度的增龄性变化规律对提高临床治疗效率、降低并发症具有重要意义。

本研究为法医学中利用第二恒磨牙髓室增龄性变化进行年龄推断提供了基础数据。相比于以往研究，本研究每一岁年龄均纳入了相同数目、相同性别比例的病人，这确保了样本分布的均质性。但研究样本量仍相对较小，需进一步扩大样本量进行研究。

无根萍虽然主要依靠出芽法繁殖，但它也会开花结果。开花之前，叶状体上方会形成一个凹陷下去的坑，称为“花腔”，里面长出一枚雌蕊和一枚雄蕊，但是没有花瓣和花萼。雌蕊先成熟，并伸出花腔开口，等雌蕊受粉或枯萎后，雄蕊才发育成熟并伸出花腔口，释放出橘红色的花粉。每一朵花结一个果实。无根萍的花朵大小不足半毫米，果实就更小了，只有在显微镜下才能观察到。

[参考文献]

[1] CIAPPARELLI L. The chronology of dental development and age assessment[M]. Oxford: Wright Butterworth-Heinemann Ltd, 1992:22-42.

[2] MASTHAN K M K. Age and sex. Textbook of forensic odontology[M]. New Delhi: Jaypee Brothers Medical Publishers (P) Ltd, 2009:59-65.

[3] 于世凤. 口腔组织病理学[M]. 北京:人民卫生出版社, 2012:23-80.

[4] MORSE D R. Age-related changes of the dental pulp complex and their relationship to systemic aging[J]. Oral Surg Oral Med Oral Pathol, 1991,72(6):721-745.

[5] SOLHEIM T. Amount of secondary dentin as an indicator of age[J]. Scand J Dent Res, 1992,100(4):193-199.

[6] 陈丽萍,张丁华,阙珍妮. 1033名成年人牙齿磨损情况的调查分析[J]. 临床口腔医学杂志, 2007(3):170-172.

[7] 刘东秀,邵金陵,万毅,等. 416例人群牙齿的平均磨耗程度调查[J]. 第四军医大学学报, 2007(6):550-552.

[8] JAGANNATHAN N, NEELAKANTAN P, THIRUVENGADAM C, et al. Age estimation in an Indian population using pulp/tooth volume ratio of mandibular canines obtained from cone beam computed tomography[J]. J Forensic Odontosto-matol, 2011,29(1):1-6.

[9] PINCHI V, PRADELLA F, BUTI J, et al. A new age estimation procedure based on the 3D CBCT study of the pulp cavity and hard tissues of the teeth for forensic purposes: A pilot study[J]. J Forensic Leg Med, 2015,36:150-157.

[10] PORTO L V, CELESTINO D S, ANJOS P A, et al. Evaluation of volumetric changes of teeth in a Brazilian population by using cone beam computed tomography[J]. J Forensic Leg Med, 2015,36:4-9.

[11] SOMEDA H, SAKA H, MATSUNAGA S, et al. Age estimation based on three-dimensional measurement of mandibular central incisors in Japanese[J]. Forensic Sci Int, 2009,185(1-3):110-114.

[12] STAR H, THEVISSEN P, JACOBS R, et al. Human dental age estimation by calculation of pulp-tooth volume ratios yielded on clinically acquired cone beam computed tomography images of monoradicular teeth[J]. J Forensic Sci, 2011,56(Suppl 1):S77-82.

[13] GE Z P, YANG P, LI G, et al. Age estimation based on pulp cavity/chamber volume of 13 types of tooth from cone beam computed tomography images[J]. Int J Legal Med, 2016,130 (4):1159-1167.

[14] 刘彦,牛忠英,王睿,等. 基于CBCT对汉族人下颌第一恒磨牙髓室高度及牙本质厚度与年龄关系的研究[J]. 实用口腔医学杂志, 2017,33(3):381-384.

[15] GE Z P, MA RH, Li G, et al. Age estimation based on pulp chamber volume of first molars from cone beam computed tomography images[J]. Forensic Sci Int, 2015,253:e131-133,e137.

[16] ABOSHI H, TAKAHASHI T, KOMURO T. Age estimation using microfocus X-ray computed tomography of lower premolars [J]. Forensic Sci Int, 2010, 200(1-3):35-40.

[17] OI T, SAKA H, IDE Y. Three-dimensional observation of pulp cavities in the maxillary first premolar tooth using micro-CT[J]. Int Endod J, 2004,37(1):46-51.