Molecular Epidemiological Investigation of HPV 56 in Precancerous Cervical Cancer Patients in Yunnan Province
-
摘要:
目的 研究单一感染HPV56型的癌前病变患者中HPV56 E6、E7和L1基因变异特征,获取云南地区宫颈癌癌前病变患者中HPV56型分子流行病学数据。 方法 随机选取300例宫颈癌癌前病变患者,通过HPV分型试剂盒对患者感染HPV的型别进行检测,并对HPV56型感染患者中HPV E6 E7和L1基因进行测序,最后完成对研究对象所感染的HPV56的谱系分布及HPV56 E6、E7和L1基因的遗传变异进行分析。 结果 研究对象所感染的HPV56病毒主要为A谱系,90%以上为A谱系,约10%为B谱系,HPV56 E6、E7、L1基因中共存在13个变异位点,其中E6基因中6个、E7基因中4个、L1基因中3个。HPV56 E6基因的6个变异位点中有2个为同义突变(T158C(S19S)、G437A(P112P)),4个为非同义突变(A141C(S14R)、A241G(N47S)、A263C(K54N)、G279A(D60N))。HPV56 E7基因的4个变异位点中均为非同义突变(C601A(D10E)、G605A(V12I)、G770C(E67Q)、G802C(Q77H))。HPV56 L1基因的3个变异位点中有2个为同义突变(G6454T(P321P)、A6868C(I459I)),1个为非同义突变(A5867G(N126D))。 结论 初步获得了云南地区宫颈癌癌前病变患者中HPV56 E6、E7和L1基因变异情况,获得了HPV56 E6、E7和L1基因的分子流行学特征,为进一步分析云南地区宫颈癌癌前病变患者中HPV56 E6、E7和L1基因的分子流行病学特征提供了必要的数据。 Abstract:Objective To study the genetic variation characteristics of HPV56 E6, E7 and L1 in patients with precancerous lesions infected with HPV56, and to obtain the molecular epidemiological data of HPV56 in patients with precancerous lesions of cervical cancer in Yunnan. Methods A total of 300 cervical intraepithelial neoplasia patients were enrolled and HPV genotyping kit was used to detect the types of HPV infection in the patients. The HPV E6, E7, and L1 genes in patients infected with HPV 56 were sequenced. Finally, the lineage distribution of HPV 56 and the genetic variations of HPV 56 E6, E7, and L1 genes in the study subjects were analyzed. Results In this study, the HPV56 virus infected by the study subjects mainly belonged to the A lineage, with more than 90% belonged to the A lineage, and about 10% belonged to the B lineage. There were 13 mutation loci in the E6, E7 and L1 genes of HPV56, including 6 in the E6 gene, 4 in the E7 gene and 3 in the L1 gene. Among the 6 mutation sites of HPV56 E6 gene, 2 were synonymous mutations (T158C (S19S), G437A (P112P)), and 4 were non-synonymous mutations (A141C (S14R), A241G (N47S), A263C (K54N), G279A (D60N)). The four mutation sites of HPV56 E7 gene were all non-synonymous mutations (C601A (D10E), G605A (V12I), G770C (E67Q), G802C (Q77H)). Among the three mutation sites of HPV56 L1 gene, two were synonymous mutations (G6454T (P321P), A6868C (I459I)), and one was non-synonymous mutation (A5867G (N126D)). Conclusion This study obtained preliminary data on the genetic variations of HPV 56 E6, E7, and L1 genes in CIN patients in Yunnan area, and obtained the molecular epidemiological characteristics of HPV 56 E6, E7, and L1 genes. It provides necessary data for further analysis of the molecular epidemiological characteristics of HPV 56 E6, E7, and L1 genes in CIN patients in Yunnan area. -
Key words:
- Cervical intraepithelial neoplasia /
- HPV56 /
- Variation /
- Lineage
-
宫颈癌是全球女性中第四大常见肿瘤,主要致病因素是高危型人乳头瘤病毒的持续性感染[1]。80%以上的女性在一生中经历过至少1次人乳头瘤病毒(human pappilomavirus,HPV)感染,但90%以上的HPV感染可在2 a内自然清除,仅不足1%的患者可能发展成为癌前病变(cervical intraepithelial neoplasia,CIN),甚至引发宫颈癌(cervical cancer,CC)[2]。癌前病变一般依据宫颈上皮内瘤变的严重程度划分为CIN1、CIN2、CIN3。目前,与宫颈癌发生相关的高危型HPV型别主要包括HPV16、HPV18、HPV31、HPV33、HPV35、HPV39、HPV45、HPV51、HPV52、HPV56、HPV58、HPV59、HPV68等[3]。
1989年,Lörincz等[4]在华盛顿特区1名妇女的CIN组织样本中首次发现了HPV56基因型。尽管HPV56感染导致宫颈癌发生的比例在全球范围来看相对较小,但在亚洲地区,特别是中国,HPV56具有相对较高的检出率。2014年,Yang等[5]在中国西南地区74例宫颈癌组织标本中发现HPV56阳性率为1.64%(1/74)。2019年,Zhao等[6]对华北地区HPV56的感染率研究表明,在宫颈正常组中单一感染HPV56型阳性率为0.66%(61/9256);在癌前病变组中单一感染HPV56型在CIN1期的阳性率为6.80%(7/103),在CIN2期为0%(0/97)。以上研究数据表明,HPV56在癌前病变组和宫颈癌组中的感染率远高于宫颈正常组, HPV56可能在宫颈癌的发生发展中发挥重要作用。
HPV基因组分为3个功能区,分别是长控制区,早期基因区和晚期基因区[7]。长控制区主要调控病毒的转录和复制,控制病毒蛋白和感染颗粒的产生;早期基因区编码产生E1,E2,E4,E5,E6和E7蛋白,其中E6和E7蛋白为病毒癌蛋白,在细胞癌变过程中发挥极其重要的作用;晚期基因编码产生在病毒颗粒自我组装以及侵入目标细胞时起协调和识别作用的L1和L2病毒衣壳蛋白[7]。研究表明[8-11]HPV基因变异可导致致癌风险发生改变,特别是E6和E7。E6和E7基因变异可导致编码氨基酸改变而使致癌能力发生改变。另外,L1基因变异导致的L1氨基酸序列变化可能会对病毒的抗原性和感染能力造成影响[12]。研究表明[13- 14],根据L1基因的序列变异特征,HPV56分为A谱系和B谱系。
本研究以云南省宫颈癌癌前病变患者为研究对象,将患者组织样本进行HPV分型,选取单一感染HPV56的样本提取基因组,扩增得到HPV56 E6,E7和L1基因,并对获得基因序列进行测序,以此分析云南地区HPV56谱系分布特征及E6,E7和L1基因的遗传变异数据,为进一步HPV56疫苗研发及相应的诊断和治疗提供参考数据。
1. 资料与方法
1.1 样本采集
选取2022年2月至2023年2月昆明医科大学第三附属医院妇瘤科宫颈癌癌前病变患者作为研究对象,采集宫颈组织样本,将组织样本保存于组织保存液中,-80℃保存。诊断标准依据《临床诊疗指南-妇产科学分册》[15],研究对象纳入标准为:(1)所有患者有性生活史;(2)单一宫颈、无其它宫颈病变。研究对象排除标准为:(1)未接受过化疗、放疗、手术等抗肿瘤治疗;(2)未合并或继发其他恶性肿瘤。所有研究对象均为彼此之间无血缘关系并且长期居住于云南地区的汉族个体,研究方案经昆明医科大学第三附属医院伦理审查委员会批准(KYCS2021193),所有样本采集均征得本人知情同意。
1.2 样本DNA提取
取-80℃冷冻保存的宫颈组织样本30~60 mg,于室温融化。加入钢珠及600 mL Buffer RLT,置于超声破碎仪,50 HZ破碎2 min,取出钢珠,8 000 r/min离心3 min。取上清置于新的离心柱,8000 r/min离心1 min,而后加入500 μL AW1于8 000 r/min离心30 s,弃上清。每个样本中再加入500 μL AW2后14000 r/min离心2 min。加入100 μL EB,室温孵育10 min,然后8 000 r/min离心1 min。得到的溶液即含释放的DNA,取出备用。
1.3 HPV检测分型
使用基于核酸分子导流杂交系统的HPV分型检测试剂盒(凯普生物科技有限公司)对样本DNA进行HPV分型。
1.4 HPV56 E6、E7、L1基因扩增测序
根据分型结果,选取检测出单一HPV56阳性的患者作为研究对象,扩增HPV56 E6、E7、L1基因。根据参考文献[13]设计针对HPV56 E6、E7、L1基因的扩增和测序引物,交由生工生物工程(上海)有限公司合成。用于设计扩增引物和测序引物的HPV56序列为X74483.1,各基因的扩增引物和测序引物,见表1。PCR反应体积为50 μL,反应条件为98℃预变性2 min,98℃变性10 s,Tm℃退火5 s,72℃延伸45 s,72℃延伸5 min,共35个循环。PCR产物经2%琼脂糖凝胶电泳成像验证后,把PCR原液送交生工生物工程(上海)有限公司测序。
表 1 HPV56 E6、E7、L1基因的扩增引物和测序引物Table 1. The amplification primers and sequencing primers of E6 E7 and L1 genes of HPV56基因 方向 序列5′-3′ 退火温度℃ E6 Forward ATTGGGAGTGACCGAAAAGG 63 Reverse ACAACACGCAGGTCCTCTTT E7 Forward GCTACAGATGTCAAAGTCCG 61 Reverse GCCTCTACTTCAAACCATCC L1-1 Forward GCCCCTTTAGGTAATGTGTGG 63 Reverse CCACATAGAATCACCATAGGCAT L1-1 Forward ATGATAGACACAGGATTTGGCG 63 Reverse ATACAACACACAAACACAGTTACAC 1.5 HPV56 E6、E7、L1基因变异和系统发生树分析
将测序所得的HPV56 E6、E7、L1基因序列分别用Molecular Evolutionary Genetics Analysis(MEGA)v 7.0软件中的Clustal W多重比对工具进行比对,而后分析基因序列的变异情况。用于分析基因变异情况的参考序列下载于GenBank,登录号为X74483.1。选用Molecular Evolutionary Genetics Analysis(MEGA)v 7.0软件中的ML法来构建系统发育树,通过bootstrap值为1000对树的可信度进行评估,选用模型为Kimura 2-parameter(Kimura-2参数)。构树参考序列从GenBank序列数据库下载,参考序列的登录号分别为:EF177176.1、EF177177.1、EF177178.1、EF177179.1、EF177180.1、EF177181.1、KU298915.1、KU298916.1、KU298917.1,KU298918.1和KU298919.1。
1.6 统计学处理
本研究HPV56型E6、E7和L1基因变异发生数的录入及基因变异发生频率的计算等处理均采用Microsoft Excel软件完成。
2. 结果
2.1 HPV56 E6、E7、L1基因扩增情况
本研究通过对300例癌前病变患者组织样本进行检测,发现单一感染HPV56的样本有32例,其中CIN1期16例,CIN2期16例。从32个样本中分别扩增出完整HPV56 E6序列25条(CIN1期14条,CIN2期11条)、HPV56 E7序列23条(CIN1期11条,CIN2期12条)、HPV56 L1序列28条(CIN1期16条,CIN2期12条)。
2.2 HPV56 E6、E7、L1基因的系统发生树分析
本研究用从样本中扩增得到的HPV56 E6,E7和L1序列来构建系统发生树,结果表明,扩增得到的HPV56 E6序列中A谱系23例(CIN1期14例,CIN2期9例),B谱系2例(CIN1期0例,CIN2期2例),见图1;HPV56 E7序列中A谱系21例(CIN1期11例,CIN2期10例),B谱系2例(CIN1期0例,CIN2期2例),见图2;28例HPV56 L1序列中A谱系26例(CIN1期16例,CIN2期10例),B谱系2例(CIN1期0例,CIN2期2例),见图3。
2.3 HPV56 E6、E7、L1基因的变异分析
经比对扩增得到的HPV56 E6,E7和L1基因系列,发现共存在13个遗传变异位点,其中E6基因中6个、E7基因中4个、L1基因中3个,见表2。E6基因中发现6个核苷酸变异,分别是A141C(S14R)、T158C(S19S)、A241G(N47S)、A263C(K54N),G279A(D60N)和G437A(P112P),其中A141C(S14R)、A241G(N47S)、A263C(K54N)变异占扩增序列的8.0%,T158C(S19S)、G279A(D60N)、G437A(P112P)变异占扩增序列的92.0%,见表2。E7基因中发现4个核苷酸变异,分别为C601A(D10E)、G605A(V12I),G770C(E67Q)和G802C(Q77H),其中C601A(D10E)和G770C(E67Q)变异占扩增序列的91.3%,G605A(V12I)变异占扩增序列的8.7%,G802C(Q77H)变异占扩增序列的13.0%,见表2。L1基因中发现3个核苷酸变异,分别为A5867G(N126D),G6454T(P321P)和A6868C(I459I),其中G6454T(P321P)发生100%变异,A5867G(N126D)和A6868C(I459I)变异占扩增序列的96.4%,见表2。
表 2 HPV56 E6、E7、L1基因变异情况Table 2. E6,E7 and L1 gene variation of HPV56核苷酸位置 E6 E7 L1 141 158 241 263 279 437 601 605 770 802 5867 6454 6868 X74483.1-E6 A T A A G G C G G G A G A cin1-2 − C − − A A A − C C G T C cin1-4 − C − − A A A − C − G T C cin1-5 − C − − A A A − C − 0 0 0 cin1-6 − C − − A A A − C − 0 0 0 cin1-7 − C − − A A A − C − G T C cin1-8 − C − − A A A − C − G T C cin1-9 − C − − A A A − C − G T C cin1-10 − C − − A A A − C − G T C cin1-12 − C − − A A A − C − G T C cin1-14 − C − − A A A − C − G T C cin1-16 − C − − A A A − C − G T C cin1-17 − C − − A A 0 0 0 0 G T C cin1-19 − C − − A A 0 0 0 0 0 0 0 cin1-20 − C − − A A 0 0 0 0 G T C cin1-21 0 0 0 0 0 0 0 0 0 0 G T C cin1-22 0 0 0 0 0 0 0 0 0 0 G T C cin2-2 0 0 0 0 0 0 A − C − 0 0 0 cin2-5 0 0 0 0 0 0 A − C − 0 0 0 cin2-12 0 0 0 0 0 0 A − C − 0 0 0 cin2-16 − C − − A A A − C − G T C cin2-18 − C − − A A A − C − 0 0 0 cin2-19 − C − − A A A − C − G T C cin2-22 0 0 0 0 0 0 A − C − 0 0 0 cin2-28 − C − − A A 0 0 0 0 G T C cin2-29 − C − − A A 0 0 0 0 0 0 0 cin2-32 − C − − A A 0 0 0 0 G T C cin2-33 − C − − A A A − C − 0 0 0 cin2-34 − C − − A A 0 0 0 0 0 0 0 cin2-35 − C − − A A A − C − G T C cin2-37 0 0 0 0 0 0 A − C − G T C cin2-42 C − G C − − − A − C 0 0 0 cin2-46 C − G C − − − A − C 0 0 0 氨基酸变化 S14R − N47S K54N D60N − D10E V12I E67Q Q77H N126D − − 标记为“0”的部分表示未扩增/测出基因序列; 标记为“-”的部分表示未检测到变异。 3. 讨论
HPV56是导致宫颈癌发生的高危型HPV之一[16],其感染率在1.1%~10.9%之间,尤其在亚洲地区具有相对较高检出率[6, 17-19]。因此,本研究对云南地区CIN患者中HPV56型E6、E7和L1基因进行了测序,并基于以上3个基因的序列数据对云南地区HPV56谱系特征及其变异特征进行了分析。
在本研究中,HPV56 L1基因系统发生树结果显示,92.9%的患者感染的HPV56病毒属于A谱系(26/28),而仅有7.1%的患者属于B谱系(2/28)。2018年,Jing等[13]对四川地区HPV56分布特征进行研究,结果显示,50.0%的患者感染属于A谱系(8/16),50.0%的患者属于B谱系(8/16)。该结果表明,在四川地区和云南地区流行的HPV56型的L1基因存在一定的差异,所以造成了基于L1获得的谱系分析结果存在差异。2015年,Wu等[20]对四川地区宫颈癌患者中HPV16进行分子流行病学调查后发现,四川地区HPV16类型主要属于欧洲谱系(病例组:86.7%),但是在云南地区HPV16以亚洲谱系为主导地位(病例组:72.7%)[21]。该研究结果之间的差异也充分说明四川和云南人群所感染的HPV谱系是存在一定差异的。
E6蛋白是HPV主要的2个致癌蛋白之一,在宫颈癌的发生发展中发挥着极其重要的作用。HR-HPV E6蛋白可通过直接或间接调控PI3K/AKT/mTOR、p53等细胞凋亡信号传导中的关键蛋白、从而促进肿瘤的发生发展[22-23]。HPV E6蛋白含有2个CXXC-X29 -CXXC锌指结构,分别位于E6蛋白的33-69和106-142氨基酸序列区域。本研究发现A241G(N47S),A263C(K54N)和G279A(D60N)非同义突变均位于HPV56 E6蛋白的第1个锌指结构区域内,这些位点的变异可能会导致HPV56 E6蛋白的锌指结构区域折叠结构的改变或是与E6靶标的结合而影响E6蛋白的功能[23-24]。同时,本研究发现A141C(S14R)和A263C(K54N)变异,但2008年,Hiller等[25]研究发现携带A141C(S14R)和A263C(K54N)突变的HPV56变异株的E6蛋白与HPV56原型的E6蛋白在促使p53 降解的功能上并无差异。因此,尽管A141C(S14R)和A263C(K54N)在云南地区CIN患者中存在变异,但这2个变异可能与CIN和宫颈癌的进展无关。
E7蛋白是HPV另一个重要的致癌蛋白,主要是通过与视网膜母细胞瘤蛋白质(Rb)家族结合促使Rb家族蛋白水解而刺激细胞周期失控,进而促进宫颈癌的发生发展。HPV E7蛋白N端一般含有2个保守区域(CR1和CR2),与Rb家族蛋白结合的LXCXE基序和酪蛋白激酶II磷酸化位点位于CR2区域内,另外,在E7蛋白的C端含有1个和E6蛋白相似的CXXC-X29 -CXXC锌指结构[26]。本研究发现4个非同义突变中的C601A(D10E)和G605A(V12I)位于E7蛋白的CR1区域内,G770C(E67Q)和G802C(Q77H)位于E7锌指结构区域,CR2区域的LXCXE基序未发现突变,提示CR2区域可能在E7基因中较CR1及锌指结构保守。
HPV L1蛋白是HPV病毒的衣壳蛋白,是HPV病毒的主要结构蛋白,在病毒DNA导入宿主细胞并整合进入宿主细胞染色体过程中起着重要作用。研究表明L1蛋白可以作为人体细胞免疫反应靶蛋白被用于HPV疫苗的研究中[26]。有研究表明,HPV L1基因变异而引起的氨基酸改变会影响L1蛋白的自我组装和免疫细胞对VLPS的识别从而导致免疫原性产生差异以及造成免疫逃逸发生[27- 28]。Jing等[13]发现HPV56 L1基因中引起K116M 和N126D氨基酸变化的变异可能具有破坏性,而本研究所发现的HPV56 L1基因的3个变异位点中有仅A5867G(N126D)为非同义突变,而未发现存在导致K116M氨基酸改变的变异。
在HPV疫苗研究中,预防性疫苗主要关注点在于HPV L1的特征研究,而治疗性疫苗则主要关注于HPV E6和E7的特征研究,为此,本研究为了解HPV56的变异特征和疫苗研发提供了参考数据。本研究的局限性主要在于由于HPV56在人群中感染比例较小,故纳入本研究的样本数量有限,从而导致在分析基因变异及HPV56谱系分布特征代表性可能不足。
-
表 1 HPV56 E6、E7、L1基因的扩增引物和测序引物
Table 1. The amplification primers and sequencing primers of E6 E7 and L1 genes of HPV56
基因 方向 序列5′-3′ 退火温度℃ E6 Forward ATTGGGAGTGACCGAAAAGG 63 Reverse ACAACACGCAGGTCCTCTTT E7 Forward GCTACAGATGTCAAAGTCCG 61 Reverse GCCTCTACTTCAAACCATCC L1-1 Forward GCCCCTTTAGGTAATGTGTGG 63 Reverse CCACATAGAATCACCATAGGCAT L1-1 Forward ATGATAGACACAGGATTTGGCG 63 Reverse ATACAACACACAAACACAGTTACAC 表 2 HPV56 E6、E7、L1基因变异情况
Table 2. E6,E7 and L1 gene variation of HPV56
核苷酸位置 E6 E7 L1 141 158 241 263 279 437 601 605 770 802 5867 6454 6868 X74483.1-E6 A T A A G G C G G G A G A cin1-2 − C − − A A A − C C G T C cin1-4 − C − − A A A − C − G T C cin1-5 − C − − A A A − C − 0 0 0 cin1-6 − C − − A A A − C − 0 0 0 cin1-7 − C − − A A A − C − G T C cin1-8 − C − − A A A − C − G T C cin1-9 − C − − A A A − C − G T C cin1-10 − C − − A A A − C − G T C cin1-12 − C − − A A A − C − G T C cin1-14 − C − − A A A − C − G T C cin1-16 − C − − A A A − C − G T C cin1-17 − C − − A A 0 0 0 0 G T C cin1-19 − C − − A A 0 0 0 0 0 0 0 cin1-20 − C − − A A 0 0 0 0 G T C cin1-21 0 0 0 0 0 0 0 0 0 0 G T C cin1-22 0 0 0 0 0 0 0 0 0 0 G T C cin2-2 0 0 0 0 0 0 A − C − 0 0 0 cin2-5 0 0 0 0 0 0 A − C − 0 0 0 cin2-12 0 0 0 0 0 0 A − C − 0 0 0 cin2-16 − C − − A A A − C − G T C cin2-18 − C − − A A A − C − 0 0 0 cin2-19 − C − − A A A − C − G T C cin2-22 0 0 0 0 0 0 A − C − 0 0 0 cin2-28 − C − − A A 0 0 0 0 G T C cin2-29 − C − − A A 0 0 0 0 0 0 0 cin2-32 − C − − A A 0 0 0 0 G T C cin2-33 − C − − A A A − C − 0 0 0 cin2-34 − C − − A A 0 0 0 0 0 0 0 cin2-35 − C − − A A A − C − G T C cin2-37 0 0 0 0 0 0 A − C − G T C cin2-42 C − G C − − − A − C 0 0 0 cin2-46 C − G C − − − A − C 0 0 0 氨基酸变化 S14R − N47S K54N D60N − D10E V12I E67Q Q77H N126D − − 标记为“0”的部分表示未扩增/测出基因序列; 标记为“-”的部分表示未检测到变异。 -
[1] Arbyn M,Weiderpass E,Bruni L,et al. Estimates of incidence and mortality of cervical cancer in 2018: A worldwide analysis[J]. Lancet Glob Health,2020,8(2):e191-e203. doi: 10.1016/S2214-109X(19)30482-6 [2] 李双,李明珠,丛青,等. 人乳头瘤病毒疫苗临床应用中国专家共识[J]. 中国妇产科临床杂志,2021,22(2):225-234. doi: 10.13390/j.issn.1672-1861.2021.02.045 [3] Burd E M. Human papillomavirus and cervical cancer[J]. Clin Microbiol Rev.,2003,16(1):1-17. doi: 10.1128/CMR.16.1.1-17.2003 [4] Lorincz A T, Quinn A P, Goldsborough M D, et al. Human papillomavirus type 56: A new virus detected in cervical cancers [J]. J Gen Virol, 1989, 70 ( Pt 11): 3099-3104. [5] Yang L,Yang H,Wu K,et al. Prevalence of HPV and variation of HPV 16/HPV 18 E6/E7 genes in cervical cancer in women in South West China[J]. J Med Virol,2014,86(11):1926-1936. doi: 10.1002/jmv.24043 [6] Zhao S,Zhao X,Hu S,et al. Distribution of high-risk human papillomavirus genotype prevalence and attribution to cervical precancerous lesions in rural North China[J]. Chin J Cancer Res,2019,31(4):663-672. doi: 10.21147/j.issn.1000-9604.2019.04.10 [7] Harari A,Chen Z,Burk R D. Human papillomavirus genomics: Past,present and future[J]. Curr Probl Dermatol,2014,45:1-18. [8] Zhou Z,Yang H,Yang L,et al. Human papillomavirus type 16 E6 and E7 gene variations associated with cervical cancer in a Han Chinese population[J]. Infect Genet Evol,2019,73:13-20. doi: 10.1016/j.meegid.2019.04.008 [9] Wu Y, Chen Y, Li L, et al. Analysis of mutations in the E6/E7 oncogenes and L1 gene of human papillomavirus 16 cervical cancer isolates from China [J]. J Gen Virol, 2006, 87(Pt 5): 1181-1188. [10] Lee C W,Bae J H,Lee S J,et al. Distribution of human papillomavirus type 16 E6 and E7 gene variants in the progression of cervical dysplasia in Korean women[J]. J Obstet Gynaecol Res,2011,37(10):1320-1326. doi: 10.1111/j.1447-0756.2010.01520.x [11] Sun Z,Ren G,Cui X,et al. Genetic diversity of HPV-16 E6,E7,and L1 genes in women with cervical lesions in Liaoning Province,China[J]. Int J Gynecol Cancer,2011,21(3):551-558. doi: 10.1097/IGC.0b013e3182112023 [12] Buck C B,Day P M,Trus B L. The papillomavirus major capsid protein L1[J]. Virology,2013,445(1-2):169-174. doi: 10.1016/j.virol.2013.05.038 [13] Jing Y,Wang T,Chen Z,et al. Phylogeny and polymorphism in the long control regions E6,E7,and L1 of HPV Type 56 in women from southwest China[J]. Mol Med Rep,2018,17(5):7131-7141. [14] Burk R D,Harari A,Chen Z. Human papillomavirus genome variants[J]. Virology,2013,445(1-2):232-243. doi: 10.1016/j.virol.2013.07.018 [15] 中华医学会. 临床诊疗指南. 妇产科学分册[M]. 北京: 人民卫生出版社, 2007: 69-73. [16] Munoz N,Bosch FX,de Sanjose S,et al. Epidemiologic classification of human papillomavirus types associated with cervical cancer[J]. N Engl J Med,2003,348(6):518-527. doi: 10.1056/NEJMoa021641 [17] Franceschi S,Rajkumar R,Snijders P J,et al. Papillomavirus infection in rural women in southern India[J]. Br J Cancer,2005,92(3):601-606. doi: 10.1038/sj.bjc.6602348 [18] Nielsen A,Kjaer S K,Munk C,et al. Type-specific HPV infection and multiple HPV types: Prevalence and risk factor profile in nearly 12,000 younger and older Danish women[J]. Sex Transm Dis,2008,35(3):276-282. doi: 10.1097/OLQ.0b013e31815ac5c7 [19] Liu J,Ma S,Qin C,et al. Prevalence and genotype distribution of human papillomavirus in Zhengzhou,China,in 2016[J]. Arch Virol,2020,165(3):731-736. doi: 10.1007/s00705-019-04515-3 [20] Wu H,Wu E,Ma L,et al. Lineage distribution and E2 sequence variation of high-risk human papillomavirus types isolated from patients with cervical cancer in Sichuan Province,China[J]. Arch Virol,2015,160(11):2845-2855. doi: 10.1007/s00705-015-2559-z [21] Dai S,Yao Y,Yan Z,et al. The association of human papillomavirus type 16 E2 variations with cervical cancer in a Han Chinese population[J]. Infect Genet Evol,2018,64:241-248. doi: 10.1016/j.meegid.2018.07.002 [22] Morgan E L,Scarth J A,Patterson M R,et al. E6-mediated activation of JNK drives EGFR signalling to promote proliferation and viral oncoprotein expression in cervical cancer[J]. Cell Death Differ,2020,28(5):1669-1687. [23] Scarth J A,Patterson M R,Morgan E L,et al. The human papillomavirus oncoproteins: A review of the host pathways targeted on the road to transformation[J]. J Gen Virol,2021,102(3):001540. [24] Howie H L,Katzenellenbogen R A,Galloway D A. Papillomavirus E6 proteins[J]. Virology,2009,384(2):324-334. doi: 10.1016/j.virol.2008.11.017 [25] Hiller T,Stubenrauch F,Iftner T. Isolation and functional analysis of five HPVE6 variants with respect to p53 degradation[J]. J Med Virol,2008,80(3):478-483. doi: 10.1002/jmv.21093 [26] Roman A,Munger K. The papillomavirus E7 proteins[J]. Virology,2013,445(1-2):138-168. doi: 10.1016/j.virol.2013.04.013 [27] Kirnbauer R,Taub J,Greenstone H,et al. Efficient self-assembly of human papillomavirus type 16 L1 and L1-L2 into virus-like particles[J]. J Virol,1993,67(12):6929-6936. doi: 10.1128/jvi.67.12.6929-6936.1993 [28] Tan W,Felber B K,Zolotukhin A S,et al. Efficient expression of the human papillomavirus type 16 L1 protein in epithelial cells by using Rev and the Rev-responsive element of human immunodeficiency virus or the cis-acting transactivation element of simian retrovirus type 1[J]. J Virol,1995,69(9):5607-5620. doi: 10.1128/jvi.69.9.5607-5620.1995 -