近日，厦门大学公共卫生学院陈田木教授、厦门市妇幼保健院周裕林教授等在《Epidemics》上发表了“HPV感染的流行病学特征和传播性”的一项长期回顾性研究。HPV受试者主要来自厦门市、漳州市和泉州市，它们一起被称为福建金三角，收集的数据来自因各种妇科疾病而就诊的女性患者的HPV检测结果，收集日期为2013年1月至2021年4月，共计114,824例。该研究构建了HPV传播动力学模型，模拟和比较不同HPV基因型的传播性，用基本传播系数(R0)定量描述。

 建立HPV传播动力学模型

该研究首先建立了易感-症状-宫颈癌恢复(SICR)模型，并对模型进行了校准。采用Berkeley Madonna 8.3.18软件构建模型，采用最小均方根(LRMS)方法进行拟合，再根据HPV的传播途径，优化了基于男女传播的SICR模型。该模型描述了从个人易感到HPV感染再到病毒清除或癌症的整个过程，性传播被认为是模型中HPV传播的主要途径。根据研究，绝大多数HPV感染是短暂性感染，只有持续感染高危HPV才会导致宫颈癌，因此在模型中假设HPV感染最终会带来两种结果：个体清除和宫颈癌。为了保持分析的一致性并减少混杂因素，从而更准确地估算基本传播系数，研究者作出了几项假设。首先，假定在研究期间，个体感染特定型别后获得了该型别的终身保护，对该型别不再变得易感；其次，假定不同型别感染所经历的疾病进程是一样的。模型纳入了易感、感染及清除相关的组份，对传播动力学的相关参数进行了拟合。

结果

1.HPV的流行病学特征

该研究中有11872份HPV阳性样本，总感染率为14.6% [95% CI:14.4 , 14.9 ]。感染率前五位的基因型为HPV52（3.4%）、HPV58（2.3%）、HPV16（1.9%）、HPV53（1.6%）和HPV51（1.5%）。在高危型别中，感染率排名前5位的分别是HPV52（3.4%） 、HPV58（2.3%）、HPV16（1.9%）、HPV51（1.5%）和HPV39（1.5%）。

2.阳性样本的年龄组分类

所有阳性样本均按年龄组分类，高危HPV感染率的变化趋势与总体感染率基本一致，HPV16/18型和非高危型在所有年龄段的感染率均低于高危型。卡方检验显示，χ2 = 411.475，P< 0.001，差异均有统计学意义，说明不同年龄组间HPV感染率存在差异。

所有型别感染均分布在15-83岁之间，HPV16/18有两个分布高峰，20岁以下年龄组和60-69岁年龄组。研究发现各型别的感染多分布在20~69岁年龄组，主要集中在20~29岁、30~39岁和40~49岁年龄组，20岁以下和70岁以上年龄组的感染率较低。

3.曲线拟合及传播性计算

采用SICR模型对研究区域的HPV病例进行拟合，模型整体拟合良好。观察到HPV不同型别的感染没有明显的周期性，大多数型别呈波动下降趋势，而一些型别的下降趋势并不明显，如HPV51、39、42、43等。

每种型别的传播系数(R0)由中位数测量，并显示结果。在所有HPV基因型中，13种HPV基因型的R0>1，有10种为高危型，前5名是HPV56、18、58、52和53，其中HPV56、18、58和52是高风险，而HPV53不是，R0值分别为3.35（CI：0.00–9.99）、3.20（CI：0.00–6.46）、3.19（CI：1.27–6.94）、3.19（CI：1.01–8.42）和2.99（CI：0.00–9.39）。此外，HPV52型的R0>1，持续时间约为51个月，持续时间最长。还有另外6种R0>1的高危型为HPV16、51、39、33、31和35；3种R0>1的非高危型为HPV53、6和11。因此，一般来说，高危型别的传播风险高于低危型别。

该研究基于连续9年来自医院妇女妇科检查的HPV阳性数据，建立了一个HPV性传播的传播动力学模型。通过模型拟合计算不同型别HPV的传播性，并量化各型别的传染性。这是该领域首次从HPV传播动力学的角度来讨论福建金三角地区HPV类型的流行特征和传播性。

研究显示包括非疫苗覆盖型别以内的高危型HPV有着高的传播风险。此外，非高危HPV53、11和6具有一定的传播风险，提示该领域女性发生尖锐湿疣的潜在风险不应被低估。尽管与高危HPV相比，非高危HPV感染通常会导致生殖器疣等良性病变，但它也会对女性的身心健康造成巨大危害。因此，本病例区域应密切关注HPV53、11、6的预防和筛查，以实现早期发现和治疗。