王巧祎 赵希芹 郭雨辰 张璐 李晓欣 王云

[摘要] 目的 探討人乳头瘤病毒52型（HPV52）E6和E7基因多态性与宫颈上皮内瘤变（CIN）的相关性。

方法 收集来自青岛地区的177例HPV52阳性且进行宫颈活检病人的宫颈脱落细胞样本，采用PCR直接测序法检测E6和E7基因序列，在GenBank中剪切HPV52标准株和变异株序列，应用Lasergene和MEGA软件对序列进行比对和系统进化树分析。比较正常、低级别上皮内瘤变（CIN1）和高级别上皮内瘤变（CIN2/3）样本中E6和E7基因变异型的分布。

结果 138例样本成功获得E6和E7基因序列。与标准株比较，132例（95.7%）样本在E6区同时出现1处错义突变（K93R）和1处同义突变（g350t），经系统进化树分析，确定为B变异型；其余样本有5例（3.6%）属于A变异型，1例（0.7%）属于C变异型。135例（97.8%）样本在E7区同时出现2处同义突变（c751t和a801g），为B变异型；其余样本有2例（1.5%）属于A变异型，1例（0.7%）属于C变异型。正常、CIN1和CIN2/3样本之间E6和E7变异型的分布差异均无显著性（P＞0.05）。

结论 来自青岛地区正常及CIN样本中HPV52 E6和E7基因均以亚洲株B变异型为主，E6和E7基因多态性可能与CIN发生发展无明显相关性。

[关键词] 人乳头瘤病毒52；乳头瘤病毒E6蛋白质类；乳头瘤病毒E7蛋白质类；宫颈上皮内瘤样病变；多态性，单核苷酸

[中图分类号] R373.9；R711.74

[文献标志码] A

[文章编号] 2096-5532（2022）05-0667-05doi：10.11712/jms.2096-5532.2022.58.136

[开放科学（资源服务）标识码（OSID）]

[网络出版] https：//kns.cnki.net/kcms/detail/37.1517.r.20220713.1040.006.html;2022-07-14 13：37：34

RELATIONSHIP BETWEEN E6 AND E7 GENE POLYMORPHISMS OF HUMAN PAPILLOMAVIRUS TYPE 52 AND CERVICAL INTRAEPITHELIAL NEOPLASIA

WANG Qiaoyi， ZHAO Xiqin， GUO Yuchen， ZHANG Lu， LI Xiaoxin， WANG Yun

（Department of Pathogenic Biology， School of Basic Medicine， Qingdao University， Qingdao 266071， China）;

[ABSTRACT] Objective To explore the relationship between the E6 and E7 gene polymorphisms of human papillomavirus type 52 （HPV52） and cervical intraepithelial neoplasia （CIN）.

Methods We collected the cervical exfoliated cell samples from 177 patients who were positive for HPV52 and underwent cervical biopsy in Qingdao， China. The E6 and E7 gene sequences of the samples were determined by PCR-based direct sequencing. With the sequences of standard and variant HPV52 strains obtained from GenBank， sequence alignment and phylogenetic analysis were carried out using Lasergene and MEGA software. The distributions of E6 and E7 variants were compared between normal， low-grade intraepithelial neoplasia （CIN1）， and high-grade intraepithelial neoplasia （CIN2/3） samples.

Results The E6 and E7 sequences of 138 samples were obtained. Compared with the standard strain， 132（95.7%） samples had a missense mutation （K93R） and a synonymous mutation （g350t） both in the E6 region， belonging to lineage B variant as demonstrated by the phylogenetic tree; in addition， 5（3.6%） and 1（0.7%） samples belonged to lineage A and C variants， respectively. One hundred and thirty-five （97.8%） samples had two synonymous mutations （c751t and a801g） both in the E7 region， belonging to lineage B variant; 2（1.5%） and 1（0.7%） samples belonged to lineage A and C variants， respectively. There were no significant differences in the distributions of E6 and E7 variants between normal， CIN1， and CIN2/3 samples （P>0.05）.

Conclusion The HPV52 E6 and E7 variants in normal and CIN samples from Qingdao predominantly belong to Asian lineage B. E6 and E7 gene polymorphisms may not be related to the development of CIN.

[KEY WORDS] human papillomavirus 52; papillomavirus E6 proteins; papillomavirus E7 proteins; cervical intraepithelial neoplasia; polymorphism， single nucleotide

宫颈癌是严重威胁女性生命的恶性肿瘤之一，2020年研究数据显示，宫颈癌的发病率和病死率在全球女性恶性肿瘤中位居第4位。宫颈上皮内瘤变（CIN）是反映宫颈癌持续发展的重要指标，而持续人乳头瘤病毒（HPV）感染是导致CIN和宫颈癌发生发展的关键因素。HPV有200余种基因型，不同地区及不同宫颈疾病人群中HPV的感染率和基因型分布均存在差异。以往的个例报道和Meta分析数据均表明，在我国HPV 52型（HPV52）是与宫颈病变发生有关的一种重要的HPV基因型。HPV编碼基因存在多态性，E6和E7是HPV致癌的关键蛋白，其多态性与宫颈病变发生的相关性备受关注。但目前对HPV52 E6和E7基因多态性的研究相对较少，其与宫颈病变的关系尚无明确结论。本研究分析比较青岛地区正常和不同级别CIN样本中HPV52 E6和E7基因多态性，为丰富HPV52分子流行病学资料和明确HPV52基因变异在宫颈病变中的作用提供依据。

1 资料与方法

1.1 研究对象

选取2019年6月—2020年12月于青岛大学附属医院中心实验室行HPV基因型检测并进行阴道镜下宫颈活检的门诊、住院病人及体检女性为研究对象。首先采集宫颈脱落细胞样本，采用HPV分型检测试剂盒（亚能生物技术深圳有限公司）进行HPV检测，样本采集和检测步骤均严格按照试剂盒说明书的要求操作，并留取HPV52单独阳性标本-20 ℃冻存。然后利用病人门诊号（唯一号）通过医院医渡云科研平台检索留取样本病例的病理检测结果，选取其中首次进行宫颈活检的病例，并进一步检索或咨询其年龄、户籍和民族等相关信息。最后纳入177例样本进行HPV52 E6和E7基因序列检测和分析，其中病理结果显示正常者68例，CIN1者57例，CIN2/3者52例；病人年龄21～77岁，平均（43.02±11.40）岁；病人均为青岛地区户籍，汉族，有性生活史，无子宫切除、宫颈手术及HPV感染治疗史。177例样本经序列检测后共138例样本成功获得E6和E7基因序列，其中病理结果正常者47例，CIN1者46例，CIN2/3者45例。正常组年龄为22～77岁，平均（44.45±12.10）岁；CIN1组年龄21～67岁，平均（43.44±12.36）岁；CIN2/3组年龄21～65岁，平均（40.91±10.53）岁，3组之间年龄差异无统计学意义（F=1.100，P=0.336）。

1.2 E6和E7基因序列检测

采用Primer-BLAST（https：//www.ncbi.nlm.nih.gov/tools/primer-blast/）在线设计引物，扩增区域包括E6编码区（nt 102～548）和E7编码区（nt 553～852）。引物名称、序列、扩增片段大小及其在基因组中的位置见表1。先以HPV52-F1和HPV52-R1为引物进行PCR扩增，出现扩增条带的直接测序；如果没有出现扩增条带，则以该扩增产物为模板进行第2轮扩增，扩增引物为HPV52-F2和HPV52-R2，出现扩增条带者进行测序。PCR扩增反应体系内含有：Taq Green PCR Master Mix（美国赛默飞世尔科技公司）15.0 μL，上下游引物各1.2 μL，DNA模板2.0 μL，用水补足至30.0 μL。循环条件：94 ℃预变性3 min；94 ℃变性30 s，55 ℃退火30 s，72 ℃延伸1 min，共35个循环；72 ℃延伸10 min。取3 μL的PCR产物进行20 g/L琼脂糖凝胶电泳，观察扩增条带。每次PCR反应均设阴性对照和阳性对照，阴性对照为HPV阴性的宫颈癌细胞系C-33A提取的DNA，阳性对照为已知HPV52阳性标本DNA。

取上述PCR产物25 μL送北京华大基因有限公司测序，测序所用引物为HPV52-F2、HPV52-R2和HPV52-R3，采用末端终止法进行双向测序。以HPV52标准株X74481（GenBank收录号X74481）序列作参比，采用Lasergene软件（version 7.0）进行序列剪接和对排。

1.3 系统进化树分析

从GenBank中获取HPV52变异型或亚型A2、B1、B2、C1、C2和D代表株E6、E7序列，这些代表株GenBank收录号分别为HQ537739、HQ537740、HQ537743、HQ537744、HQ537746和HQ537748，标准株X74481为A1变异亚型。用MEGA 7.0软件，将本研究所获得的E6和E7序列及上述代表株序列，用Clastal W方法进行比对和同源性分析，临接法（Neighbor-joining）绘制系统进化树。

1.4 统计学分析

应用SAS 9.4统计软件进行数据分析。采用Fisher确切概率法比较不同组间HPV52 E6和E7基因变异型的分布。以P＜0.05为差异有显著性。

2 结 果

2.1 E6基因序列多态性

与标准株比较，138例标本中除2例在E6基因编码区未出现变异，其余标本均出现单核苷酸突变，未见插入或缺失。共出现6处错义突变和10处同义突变。错义突变中，K93R见于134例样本（4例样本核苷酸a378c和a379g突变导致氨基酸K93R突变，130例a379g突变导致K93R突变），其中1例样本还出现K93G突变（核苷酸a378g和a379g突变）；N122K（c467a）突变见于2例样本；Q38H（a215t）、F45R（t235a）和L67V（c300g）同时出现于1例样本；W132R（t495c）见于1例样本。同义突变中，g350t见于133例样本，g356a见于15例样本，t191c和g356a突变各见于2例样本，a180c、a278t、t309c、g371a、a422g和a530g各见于1例样本。

本文138例样本共出现15条不同序列，将这15条不同基因序列与GenBank中A1、A2、B1、B2、C1、C2和D代表序列共同绘制系统进化树（图1），结果显示，2例样本序列属于A1变异亚型，与代表株X74481-A1序列相同；3例样本序列属于A2变异亚型，与代表株HQ537739-A2比较，均出现错义突变K93R，但是未出现同义突变a251g；132例样本序列属于B变异型，与代表株HQ537740-B1和HQ537743-B2的序列相同，均同时出现错义突变K93R（1例为K93G）和同义突变g350t；1例样本序列属于C2变异亚型，与代表株HQ537746-C2比较，出现同义突变a530g，但在氨基酸83仅出现核苷酸g350t突变，而代表株同时出现核苷酸g350t和c348g突变导致氨基酸L83V突变。

正常、CIN1和CIN2/3组E6基因各变异型分布差异无统计学意义（P=0.741），3组均以B变异型为主。见表2。

2.2 E7基因序列多态性

与标准株比较，138例标本中除2例標本在E7基因编码区未出现变异，其余标本均出现单核苷酸突变，未见插入或缺失。共出现12处错义突变和7处同义突变。在错义突变中，M1T（t554c）和Y11C（a584g）同时出现于1例样本；T37I（c662t）、S52D（a706g、g707a）、Y59D（a728g）、H61Y（c733t）、D64N（g742a）和L99R（t848g）同时见于1例样本，其中D64N还见于另外2例样本；M84I（g804t）见于1例样本。在同义突变中，a801g出现于136例样本，c751t出现于133例样本，t583g、a681g、t765c、a793t和c837g各见于1例样本。

本文138例样本共出现12条不同序列，将这12条不同基因序列与GenBank中A1、A2、B1、B2、C1、C2和D代表序列共同绘制系统进化树（图2），结果显示，2例样本序列与代表株X74481-A1和HQ537739-A2序列均相同，在系统树中属于A变异型，这2条序列与E6基因属于A1变异亚型的2条序列来源于相同样本；135例样本序列与代表株HQ537740-B1和HQ537743-B2比较，均同时出现同义突变c751t和a801g，属于B变异型；1例样本序列属于C变异型，其与E6基因属于C变异型的序列来源于同一样本，该序列与HQ537746-C2的序列相同，出现6处错义突变（T37I、S52D、Y59D、H61Y、D64N和L99R）和2处同义突变（t583g和a801g）。除3例在E6区属于A2变异亚型的序列在E7区属于B变异型外，其余序列在E6和E7区的变异型分类一致。

正常、CIN1和CIN2/3组中E7基因各变异型分布差异无统计学意义（P=0.472），3组均以B变异型为主。见表3。

3 讨 论

HPV52编码基因具有多态性，根据其编码的E6、E7、主要衣壳蛋白L1和长控制区（LCR）基因多态性，将HPV52分为A、B、C、D等4种基因变异型，进一步分为A1、A2、B1、B2、C1、C2和D等7种变异亚型。以往研究显示，HPV52变异型的分布具有地域差异，亚洲以B变异型为主，欧洲、美洲和非洲均以A变异型为主，C和D变异型在各地域检出率均较低。本研究首次对来自青岛地区正常和CIN样本中HPV52 E6和E7基因进行检测，结果显示本地区样本E6和E7基因均以B变异型为主，分别占95.7%和97.8%。这与日本、我国四川、上海和北京等地的检测结果相似，2篇来自日本的研究中B变异型检出率分别为96.2%（50/52）和92.3%（36/39），四川为97.6%（41/42），上海为96.3%（52/54），北京为90.0%（45/50）。但本研究样本中HPV52 E6和E7变异型的分布与台湾和湖北荆州的检测数据存在明显差异，台湾280例样本中A、B和C变异型的检出率分别为0.7%、88.2%和11.1%，湖北荆州168例样本中B、C和D变异型的检出率分别为85%、5%和10%，这两个地区B变异型的检出率相对较低，C或者D变异型的检出率相对较高。上述各地区HPV52 B变异型的检出率均高于越南（73.7%，28/38）和菲律宾（62.5%，20/32）。可见亚洲不同地域或国内不同地域HPV52 E6和E7基因变异型的分布存在差异，对不同地域样本的检测可为全面了解HPV52的分子流行病学特征提供有价值的信息。

HPV基因变异与宫颈病变的关系是HPV致癌机制研究的热点问题之一。有研究比较同一地区不同程度宫颈病变或细胞学异常人群与正常人群中HPV52 E6和（或）E7基因变异情况，结果显示不同宫颈病变人群及正常人群之间基因变异型或基因突变分布无明显差异，从而认为HPV52 E6和E7基因变异与宫颈病变发病风险或病变严重程度不相关。本研究结果显示，HPV52 E6和E7基因变异型分布在正常、CIN1和CIN2/3组中差异无显著性，且多数样本序列相同，提示HPV52 E6和E7基因多态性与CIN发生发展无明显相关性。但也有研究结果显示，HPV52 E6和E7基因某变异型或特定突变增加了CIN或宫颈癌患病风险，如韩国的一项研究中，HPV52 E6基因以B变异型和C变异型为主，而且B变异型E6基因的K93R（A379G）突变与宫颈恶性病变相关；加拿大的一项研究亦表明，LCR的t7744c和E6的K93R突变增加了CIN2/3（与CIN1+正常比较）的发生风险；台湾一项研究显示，E6和E7的C变异型增加了高级别鳞状细胞内病变及宫颈癌的发生风险；美国最新的一项病例对照研究采用全基因组测序检测50例CIN3和39例正常女性宫颈样本，发现C2变异亚型与CIN3高风险有关；ZHANG等的研究显示，与A型比较，B型和C型增加了CIN3及宫颈癌的发生风险。与HPV52 E6和E7变异增加患病风险不同，江苏地区的一项研究和广州地区的一项研究分别显示，宫颈病变（CIN1/2/3和宫颈癌合并）样本HPV52 E6基因的突变频率显著低于正常样本，从而认为E6基因突变可能阻碍宫颈病变和宫颈癌的发生发展或降低HPV的致癌作用。由此可见，HPV52 E6和E7基因多态性与宫颈病变的相关性目前尚不能得出明确结论。导致不同研究结论不一致的原因，可能与研究样本例数偏少、纳入研究病种不一致、疾病分组不一致和地域局限等有关。

总之，本研究表明，来自青岛地区宫颈细胞样本中HPV52 E6和E7基因主要为B变异型，且E6和E7基因多态性可能与CIN发生发展无明显相关性。本研究结果为了解该地区HPV52分子流行病学特征以及进一步明确HPV52基因变异与宫颈病变的关系提供了数据支持。

[参考文献]

[1]SUNG H， FERLAY J， SIEGEL R L， et al. Global cancer statistics 2020： GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries[J]. CA： A Cancer Journal for Clinicians， 2021，71（3）：209-249.

[2]BRUSSELAERS N， SHRESTHA S， VAN DE WIJGERT J， et al. Vaginal dysbiosis and the risk of human papillomavirus and cervical cancer： systematic review and meta-analysis[J]. American Journal of Obstetrics and Gynecology， 2019，221（1）：9-18.e8.

[3]GUO C L， DU H， BELINSON J L， et al. The prevalence and distribution of human papillomavirus among 10，867 Chinese Han women[J]. Infectious Agents and Cancer， 2021，16（1）：21.

[4]BAO H L， JIN C， WANG S， et al. Prevalence of cervicovaginal human papillomavirus infection and genotypes in the pre-vaccine era in China： a nationwide population-based study[J]. The Journal of Infection， 2021，82（4）：75-83.

[5]韋梦娜，余艳琴，徐慧芳，等. 中国大陆地区不同宫颈病变人群中人乳头瘤病毒感染率及型别分布的系统研究[J]. 中国肿瘤临床与康复， 2020，27（2）：133-137.

[6]赵希芹，宋海洋，孙有湘，等. 青岛市部分地区女性HPV感染现状分析[J]. 精准医学杂志， 2021，36（5）：428-432.

[7]CHEN Z G， SCHIFFMAN M， HERRERO R， et al. Evolution and taxonomic classification of human papillomavirus 16 （HPV16）-related variant genomes： HPV31， HPV33， HPV35， HPV52， HPV58 and HPV67[J]. PLoS One， 2011，6（5）：e20183.

[8]ZHANG C Q， PARK J S， GRCE M， et al. Geographical distribution and risk association of human papillomavirus genotype 52-variant lineages[J]. The Journal of Infectious Diseases， 2014，210（10）：1600-1604.

[9]ISHIZAKI A， MATSUSHITA K， HOANG H T， et al. E6 and E7 variants of human papillomavirus-16 and-52 in Japan， the Philippines， and Vietnam[J]. Journal of Medical Virology， 2013，85（6）：1069-1076.

[10]CHOI Y J， KI E Y， ZHANG C Q， et al. Analysis of sequence variation and risk association of human papillomavirus 52 va-riants circulating in Korea[J]. PLoS One， 2016，11（12）：e0168178.

[11]王本敬，刘敏娟，李文静，等. 人乳头瘤病毒52型E6和E7基因突变与宫颈癌及癌前病变关系的研究[J]. 中国计划生育和妇产科， 2018，10（6）：24-30.

[12]秦燕萍，夏小璐，杨丽芳. HPV52型E6和E7基因突变位点分布及与宫颈癌和癌前病变的相关性[J]. 中国病案， 2020，21（5）：83-86.

[13]KUMAR S， STECHER G， TAMURA K. MEGA7： molecular evolutionary genetics analysis version 7.0 for bigger datasets[J]. Molecular Biology and Evolution， 2016，33（7）：1870-1874.

[14]TENJIMBAYASHI Y， ONUKI M， HIROSE Y， et al. Whole-genome analysis of human papillomavirus genotypes 52 and 58 isolated from Japanese women with cervical intraepithelial neoplasia and invasive cervical cancer[J]. Infectious Agents and Cancer， 2017，12：44.

[15]SONG Z L， CUI Y R， LI Q F， et al. The genetic variability， phylogeny and functional significance of E6， E7 and LCR in human papillomavirus type 52 isolates in Sichuan， China[J]. Virology Journal， 2021，18（1）：94.

[16]ZHAO J W， ZHAN Q， GUO J H， et al. Phylogeny and polymorphism in the E6 and E7 of human papillomavirus： alpha-9 （HPV16，31，33，52，58）， alpha-5 （HPV51）， alpha-6 （HPV53，66）， alpha-7 （HPV18，39，59，68） and alpha-10 （HPV6，44） in women from Shanghai[J]. Infectious Agents and Cancer， 2019，14：38.

[17]李穎，周玲，廖秦平，等. 北京地区收集的人乳头瘤病毒52型基因突变和型下分类的研究[J]. 中华实验和临床病毒学杂志， 2016，30（6）：558-561.

[18]CHANG Y J， CHEN H C， LEE B H， et al. Unique variants of human papillomavirus genotypes 52 and 58 and risk of cervical neoplasia[J]. International Journal of Cancer， 2011，129（4）：965-973.

[19]LI S， YE M X， CHEN Y L， et al. Genetic variation of E6 and E7 genes of human papillomavirus 52 from Central China[J]. Journal of Medical Virology， 2021，93（6）：3849-3856.

[20]FRATI E R， BIANCHI S， AMENDOLA A， et al. Genetic characterization of variants of HPV-16， HPV-18 and HPV-52 circulating in Italy among general and high-risk populations[J]. Molecular Medicine Reports， 2020，21（2）：894-902.

[21]BEE K J， GRADISSIMO A， CHEN Z G， et al. Genetic and epigenetic variations of HPV52 in cervical precancer[J]. International Journal of Molecular Sciences， 2021，22（12）：6463.

[22]DING T， WANG X Y， YE F， et al. Distribution of human papillomavirus 58 and 52 E6/E7 variants in cervical neoplasia in Chinese women[J]. Gynecologic Oncology， 2010，119（3）：436-443.

[23]SUN Z R， LU Z T， LIU J H， et al. Genomic polymorphism of human papillomavirus type 52 in women from Northeast China[J]. International Journal of Molecular Sciences， 2012，13（11）：14962-14972.

[24]FORMENTIN A， ARCHAMBAULT J， KOUSHIK A， et al. Human papillomavirus type 52 polymorphism and high-grade lesions of the uterine cervix[J]. International Journal of Can-cer， 2013，132（8）：1821-1830.

（本文编辑 马伟平）