Genome-Wide Profiling of Human Papillomavirus DNA Integration into Human Genome and Its Influence on PD-L1 Expression in Chinese Uygur Cervical Cancer Women

Yang-chun, Feng; Sen-yu, Wang; Yuan, Zhang; Yan-chun, Huang

doi:https://doi.org/10.1155/2020/6284960

Journal of Immunology Research

On this page

Abstract Introduction Materials and Methods Discussion Conclusion Data Availability Ethical Approval Consent Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 6284960 | https://doi.org/10.1155/2020/6284960

Genome-Wide Profiling of Human Papillomavirus DNA Integration into Human Genome and Its Influence on PD-L1 Expression in Chinese Uygur Cervical Cancer Women

Feng Yang-chun,¹Wang Sen-yu,¹Zhang Yuan,²and Huang Yan-chun¹

Academic Editor: Bogdan Kolarz

Received03 Aug 2019

Revised03 Dec 2019

Accepted30 Dec 2019

Published23 Jan 2020

Abstract

Background. The Uygur is the fifth most populous ethnic group in China. Compared to other Chinese population, cervical cancer in them had high incidence, and HPV infection also was particular. Their HPV integration situation has never been reported. We aimed to investigate the integration situation of 20 subtypes of HPV gene into host cell genome in Chinese Uygur cervical cancer patients; meanwhile, we explored the influence of gene integration on PD-L1 expression. Methods. 40 frozen Chinese Uygur cervical cancer specimens with positive HPV infection were obtained from the cancer prevention and treatment institute of Tumor Hospital Affiliated to Xinjiang Medical University. The integration situation of HPV gene into host cell genome was detected by Agilent SureSelect™ Target Enrichment Chip and Next-Generation Sequencing. The related genes were analyzed by GO functional annotation and KEGG pathway enrichment. The expression levels of PD-L1 in cancer cells were tested by immunohistochemical assay (IHC). Meanwhile, the relationship between PD-L1 levels in cancer cells and gene integration were analyzed. Results. The HPV multiple infection rate by HIVID was as high as 92.5%, much higher than 35.0% by the commercial kit (). There were 13423 integration events in 40 specimens, involving 6867 human genes. These integration events were distributed on all human chromosomes, and chromosome 19 had the excessive concentration phenomenon of integration events. There were some integration hotspots in human genome such as PPP1R37, HECW2, EMBP1, ANKRD50, SPTBN4, LINC00895, LYRM4-AS1, LINC00374, RBFOX1, CSMD1, CDH13, and KLHL4. Insertion breakpoints can be found in all gene regions of the HPV genome. The actual observation of the integration times of E1 and E6 was much higher than the expected value, while the actual observation times of E5 were much lower than the expected value. The result of GO functional analysis showed that binding molecular function and cellular process biological process were the main ways to influence the cell biological behavior of HPV gene integration. The enrichment pathway analysis of KEGG showed that pathways in cancer were the most important enrichment pathways involved in the genomic integration of HPV. The positive PD-L1 rate was 62.5%. Logistic regression analysis showed that 9p24.1 existing integration sites and the number of all gene integration were risk factors for PD-L1 expression (odds ratio 17.313 and 1.012; 95% confidence interval 1.691-177.213 and 1.001-1.023). Conclusions and Relevance. Most high-frequency sites of HPV integration in Chinese Uygur cervical cancer are related to cancer progression, and the gene integration hotspots may be potential HPV carcinogenic targets. The problem of multiple HPV infection in Chinese Uygur cervical cancer patients should be paid attention. L1 and E6 genes are inapposite as the target gene of commercial HPV type detection kit, because of high-frequency breakpoints in these genes. The gene integration especially the integration existing on 9p24.1 could affect the expression level of PD-L1.

1. Introduction

Cervical cancer was one of the most common malignant tumors in China, in which its prevalence rate was 7th in all cancers and its mortality was 8th in all female cancers [1]. HPV infection is known as the main cause of cervical cancer; most of HPV infection could be cleaned by the host immune responses, but some persistent HPV infection would lead to the cytopathy of cervical cancer cells [2–4]. The previous studies found that the HPV persistent infection consisted of two phases: in the first phase of viral replication, HPV firstly infected basal cells of the cervix, bare virus DNA was the free state in the host cell nucleus. HPV DNA independently began to replicate through using the ingredients of host cells when basal cells were in the progress of differentiation. The complete virus would be assembled by E1-E7 proteins and L1-L2 proteins. The virus DNA encapsulated in a protein coat would release from the surface mature cervical epithelial cells which create the next infection of normal cervical cells. The second phase was the transformation stage, HPV DNA integrated into the host cell chromosome. In immature epithelial cells, HPV E6 and E7 gene began to have uncontrolled expression, and the excessive expression of E6 and E7 protein would, respectively, combine with the gene products of p53 and Rh, which led to the immaturity of cervical epithelial cells; meanwhile, the viruses could not complete the whole life cycle, the L1 protein expression was missing, the infection cells created immune escape, the virus could not be cleared by the immune system, the immaturity of cervical epithelial cells led to immortalization, and lastly, it caused the lesions and canceration of cervical cells. Due to the fact that the HPV gene into the human genome was the key step in cervical cancer occurrence, so to explore the integration sites of HPV genome into the human genome and the integrated position influence on human cell genome had important clinical significance [5–7].

The mainstream methods of HPV detection included two types according to the basic principle, such as the second-generation hybrid capture method (HC-2) and fluorescence PCR method, which could be used to differentiate the subtypes of HPV infection and evaluate the total amount of viral replication in the high-risk type HPV. But these tests only focused on whether HPV gene existed in host cells, which could not detect whether HPV genome was integrated into the human genome, could not provide the location of the HPV genome into the host genome, and could not find which segment of the viral genome could integrate into the host genome. Because HPV gene integrating into the human genome was the key step that led to the host cell pathological changes, therefore it needed to analyze the HPV gene integration status, in order to assess the potential effect of HPV infection on the infected cervical cells, and thus, it met the requirements of precision medical treatment.

Though HPV genome integrating into human genome research had been reported, most of the studies only focused on certain genes of HPV16 or HPV18 type, such as the genetic integration of HPV16 E6 and E7 genes [8–12]; the integration of whole HPV genome into the host cell genome was reported in few reports. The Chinese Uygur has over a population of 10 million. The HPV infection rate and HPV subtypes in them were obviously different from the other areas in China, and they also had the highest morbidity compared with other population [13]. So it was believed that the HPV gene integration sites should have their own characteristics in cervical cancer of Chinese Uygur women.

HPV infection played an important role in the process of cervical cancer progression, which was one of the important drivers of cervical cancer. The HPV gene integrated into the cervical cell genome, then changed the cervical cell genome state, eventually creating important changes in the biological behavior of cervical cells, such as changing the normal immune state of cervical cells, leading to the immune escape of cervical pathological cells. Many studies have confirmed that PD-L1 on the surface of cervical cancer cells was an important molecular marker of tumor immunity. Therefore, the relationship between HPV gene integration into the cervical cell genome and PD-L1 expression of cervical cells was explored, which could help further clarify the role of HPV in the evolution of cervical lesions and provide new ideas for the prevention and treatment of cervical cancer.

2. Materials and Methods

2.1. Sample Source

The 40 cases frozen cervical cancer pathological specimens were collected from the specimen repository of the cancer prevention research institute in the Tumor Hospital Affiliated to Xinjiang Medical University. All cases were Chinese Uygur cervical cancer, all of have existing HPV infection which were proved by flow-through hybridization and gene chip (FHGC). The ethics committee of the Tumor Hospital Affiliated to Xinjiang Medical University approved the study and the consent procedure.

2.2. Experimental Procedure

2.2.1. The Test of HPV Gene Integration

The basic technique of the study was performed with an Agilent SureSelect™ Target Enrichment System and the second-generation sequencing technology (the HiSeq 2000 sequencer of Illumina).

(1) Chip Design and Order. We designed the sequence-capture probes according to 20 types of HPV genome (6, 11, 16, 18, 31, 33, 35, 39, 42, 43, 44, 45, 51, 52, 53, 56, 58, 59, 66, and 68) sequences, and these probes were produced by Agilent Inc., USA. The enrichment chip was also ordered from Agilent Inc., USA.

(2) DNA Extraction. The DNA of 40 cases were extracted by a Genomic DNA Extraction Kit. Agarose gel electrophoresis and Qubit fluorometer were used to evaluate the DNA quality.

(3) Sample Test. DNA fragments 150–200 bp were obtained after ultrasonic fragmentation; then, these fragments were purified, end blunted, “A” tailed, and adaptor ligated. The products of ligation were amplified by PCR and constructed into a DNA sequencing library. Bioanalyzer 2100 (Agilent Technologies, USA) was used to quantify the concentration of the DNA library. Next, hybridization was performed according to the instruction of Target Enrichment Protocol. This process was carried out at 65°C for 24 h, and the untargeted fragments were then removed by washing buffer. These eluted fragments were amplified by PCR and were further processed to paired-end 100 bp read-length for sequencing by the HiSeq 2000 sequencer (Illumina Inc., USA).

(4) Data Analysis. First, clean reads were obtained through removing low-quality and duplicated reads, as well as adaptor-contaminated reads. Burrows-Wheeler Aligner (BWA) was used to align clean reads onto human (NCBI build 37, HG19) and HPV genome. The paired-end reads that can be mapped to human or HPV reference genome with both ends were removed. The remaining reads were paired-end assembled to reconstruct fragment sequences, which was used to locate the position of breakpoints more precisely. Subsequently, reference mapping was performed again to remap the paired-end assembled reads onto human and HPV genome using BWA. The position of a breakpoint was defined as the junction of human and HPV sequence in a paired-end assembled read. To minimize the impact of different sequencing data amount of each sample, the normalized support-reads number (Norm value) was introduced, which was equal to the supported reads number of each HPV breakpoint per million clean read pairs. The breakpoints with were retained. The Sanger sequencing validation rate for the selected breakpoints was 82.5%. ANNOVAR was used to do the annotation for the integrated breakpoints.

2.2.2. Immunohistochemistry of PD-L1

Monoclonal antibodies of PD-L1 were bought from Golden Bridge Biotechnology Co. Ltd (Beijing, China). H&E staining was used to assess the pathological slices. The detection method used was IHC streptavidin peroxidase and diaminobenzidine chromogenic detection. All slices were produced using the following procedures: washing in xylene, ethanol, steamed classification, boiling water to wash, dissolving endogenous peroxide physical activity by incubation with 3% H₂O₂, antigen repair, and dyeing. Phosphate-buffered saline (PBS), rather than the first antibody, was the negative control; the previous positive specimen was the positive control. For PD-L1 expression in cancer cells, we chose the nests of cervical cancer including more than 100 cancer cells and then calculated the proportion of positive cells (10% or less was negative, more than 10% was positive); the positive criterion is in the cytomembrane and cytoplasm or both at the same time appeared tan particles.

2.3. Statistical Analysis

The Genesis analysis tool of webgenstal (http://www.webgestalt.org) and DAVID (https://david-d.ncifcrf.gov/) were used to complete the Gene Ontology categories and Kyoto Encyclopedia of Genes and Genomes (KEGG) gene function pathway categories. Multivariate logistic regression analysis was used to evaluate the relationship between the number of gene integration, the integration in 9p24.1, and PD-L1 expression by SPSS 17.0 statistical software.

3. Result

3.1. The HPV Genome Detection of 40 Cases

The most frequent was HPV16 type, followed with HPV18 type, HPV33 type, HPV39 type, and HPV45 type; their detection rate was 87.5%, 45.0%, 32.5%, 32.5%, and 25.0%, respectively. The detailed HPV infection results were shown in Table 1. From Table 1, the multiple HPV infection rate was 92.5%, which was much higher than the FHGC test in clinical laboratory (35.0%). The other high-risk HPV types besides HPV16 and HPV18 also had higher detection rate than FHCG test.

3.2. HPV Integration in Different Human Genome Chromosomes

13423 gene integration events were detected in 40 cases for 20 HPV gene subtypes, which involved 6867 human genes. The HPV gene was integrated into 24 human chromosomes and the mitochondrial genome; the detailed distribution numbers of integration events were in Figure 1. The most integrated events were detected in chromosome 2, which had 1119 integration events and had 8.34% of the total examined rate. The lowest number was found in chromosome 21, only 140 integration events.

Taking the length of every human chromosome as reference, the expected value of integration events that should exist on every chromosome was obtained, then compared to the number of integration events observed on each chromosome. The results are shown in Figure 2. In Figure 2, besides the mitochondrial chromosomes, the number between observation and expectation shows obvious difference in chromosomes 3, 4, 9, 11, 13, 15, 18, 19, 21, and 22, especially in chromosome 19, which had 611 observation integrated events, it is 2.34 times than 261 expectation events.

3.3. The Hot Integration Sites of HPV Integration in Human Genome

A total of 13423 integration events were detected in 40 samples, involving a total of 6867 human genes. Among these genes, the vast majority of genes could only be detected in one sample, and only a few genes could be detected in more than two samples. The distribution of every gene in human chromosome and the distribution of sample number obtained by detection were shown in Figure 3.

According to the results of statistical analysis, the top 10 genes in the occurrence frequency of 40 samples were PPP1R37, EMBP1, HECW2, ANKRD50, SPTBN4, LINC00895, CSMD1, LINC00374 and CDH13, RBFOX1, LYRM4-AS1, and KLHL4 genes. The detection rates were 25% (10/40), 22.5% (9/40), 22.5% (9/40), 22.5% (9/40), 22.5% (9/40), 22.5% (9/40), 20% (8/40), 20% (8/40), 17.5% (7/40), 17.5% (7/40), 17.5% (7/40), and 17.5% (7/40). In this study, these top 10 genes were referred to as hot integration gene sites of HPV gene integration.

3.4. Distribution of Integration Breakpoints in the HPV16 Genome

Among the 13423 integration events in the 40 samples, 5057 events were found to be caused by the integration of HPV16 gene into the host cell genome, which is from 35 samples. According to these 5057 integration events, the gene breakpoint distribution of HPV16 genome was analyzed, which is shown in Figure 4.

In Figure 4, every location of the HPV16 genome can be integrated into the host cell genome. According to the length of each HPV16 functional gene area, there were 81 gene breakpoints in every 100 bp intervals of E6 gene, which had the most frequent gene breakpoints. But there were 46 gene breakpoints in every 100 bp intervals of E5 gene, which was the least.

Taking the length of every functional gene of HPV16 genome as reference, the expected value of gene breakpoints that should exist in every functional gene region was obtained, then compared to the actual number of observed gene breakpoints. The results were shown in Figure 5. There was a statistical difference between the expected value and the observed value of the gene breakpoints in L1, E1, E5, and E6 genes. The observed gene breakpoints of E1 and E6 genes were much higher than the expected value, which were 1.23 times and 1.26 times, respectively. However, the observed number of gene breakpoints in E5 gene was much lower than the expected value, accounting for only 72.2% of the expected value.

3.5. Functional and Pathway Enrichment Analysis of the Integration Sites

The results of Gene Ontology categories are showed in Figure 6 and Table 2; GO categories showed that binding of the molecular function was the main part in the functional level of 6867 genes, which included 142 genes. Meanwhile, cellular process of the biological process was the key component, which had 156 genes. But for the cellular component categories, there was no preponderant component.

Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was performed on 58 human genes, whose integration sites were detected more than 20 times. The top 10 enrichment pathways of gene function pathway categories are showed in Figure 7 and Table 3. The pathways in cancer was the top enrichment pathway, which included 8 genes, including RELA, BCL2, LAMA5, CBLB, FGF18, HDAC2, RXRA, and FGF14, which accounted for 14% of gene number in the pathways.

3.6. Multifactor Analysis between the Number of Gene Integration, the Integration in 9p24.1, and PD-L1 Expression

The positive and negative PD-L1 expression in cervical cancer cells was present in 25 cases (62.5%) and 15 cases (37.5%), respectively (example in Figure 8). The number of gene integration and the gene integration status in 9p24.1 were shown in Table 4. The corresponding logistic regression expression was

(a)

(b)

Table 5 shows the results of multivariate logistic regression analysis. Therefore, the number of gene integration and the gene integration in 9p24.1 were the risk factors which were significantly associated with increased risk of PD-L1 expression.

4. Discussion

The analysis methods of the Agilent SureSelect™ Target Enrichment System and the second-generation sequencing technology (the HiSeq 2000 sequencer of Illumina) included to one of the high-throughput viral integration detection (HIVID) had proved to be very useful test methods [14, 15]. In the study, the Agilent SureSelect™ Target Enrichment Chip was designed by Agilent company, which also had its patent right, which ensured the accuracy of test results. The chip included 20 HPV types; it had firstly reported so many HPV type integration detection. Because in the previous study, most HPV type integration detection only included 11 HPV types [16], the study greatly increased the detection range of HPV types.

At present, the current clinical detection methods of HPV type were FHGC, real-time fluorescence quantification, and so on, which were the principle of various commercial HPV type detection kits. For 40 specimens, the retrospective results of HPV type detection were given by commercial HPV type detection kits, which was approved by CFDA, and the detection method was FHGC. In the retrospective results, the HPV subtypes were mainly HPV16 and HPV18; the HPV multiple infection ratio was 35% (14/40), which was exactly consistent with the previous reports [17]. But by the HIVID method, for 40 specimens, HPV16 was the most frequently detected, followed by HPV33, HPV39, HPV18, and HPV45. 37 samples were HPV multiple infection, which accounted for 92.5%. So the results between the new HIVID method and FHGC were very different. The new HIVID method depended on enrichment chip and second-generation sequencing technology, which could improve the sensitivity of detection. Because the high HPV multiple infection ratio was found by HIVID, it may be possible to make a large correction for the existing opinion about HPV typing detection; we must pay more attention to the HPV multiple infection, the infection pattern may exist far beyond our imagination.

By analyzing the specific integration sites of HPV genes, it was found that the integration sites involved all chromosomes of human genome, with a total of 6867 human genes, indicating the randomness of HPV gene integration. PPP1R37, EMBP1, HECW2, ANKRD50, SPTBN4, LINC00895, CSMD1, LINC00374, CDH13, RBFOX1, LYRM4-AS1, and KLHL4 were the high-frequency integration sites; compared to the hotspots reported by Zheng et al. [16], there were significant differences in the hotspots. In their report, POU5F1B, FHIT, KLF12, HMGA2, KLF5, LRP1B, LEPREL1, DLG2, and SEMA3D were the hotspots. But the high-frequency integration sites reported by Li et al. were FHIT, CASC8, KLF5, LINC00392, RAD51B, CASC21, ERBB2, TP63, TEX41, RAP2B, and MYC [17]. These high-frequency integration sites were quite different, which indicates that there were different high-frequency integration sites in different sample sources. With the development of detection technology, new high-frequency integration sites are constantly found. Existing studies have reported that some Common Fragile Sites (CFS) widely exist in human chromosomes [18], which are often the common integration sites of various carcinogenic virus genes into the human genome including HPV. The reported CFS were 40, and two high-frequency integration sites in the study belong to CFS, such as CSMD1 and CDH13. The result verified the assertion that the CFS also were the hotspots of the integration of HPV gene.

In Figure 4, it was found that breakpoints could occur in any part of the viral genome, perhaps enabling the virus to adapt to the changing environment during carcinogenesis. The E1 gene functions had most breakpoint number in every 100 bp intervals, and E2 gene functions had the least breakpoint number in every 100 bp intervals. In the previous report, Zheng et al. also showed the same phenomenon [16]. First, E1 involves viral DNA replication, which plays a key role when the virus begins to replicate. Because E1 had the most breakpoint number, the virus must be different to replicate [19–22]. Then, L1 gene also had some breakpoints, and the detection target for most commercial HPV type detection kits was L1 gene. So the two reasons could lead to missed detection due to HPV genomic integration [23, 24].

By the functional and pathway enrichment analysis, it could be found that the HPV gene integration could mainly affect the physiological process of host cells and could also affect pathways in cancer. The GO and KEGG categories analyzed the integration gene sites, which was the first comprehensive analysis about the integration landscape. The pathways in cancer were a cell signal transduction system which played an important role in cell carcinogenesis and cancer evolution. Because HPV gene integration had the deepest level of impact about the pathways in cancer, so it could be supposed that the HPV gene integration into the cervical cells could cause the occurrence and progression of cervical cancer.

As for the relationship between HPV infection and PD-L1 expression in cervical cells, it is still a great controversy. Most reports showed that HPV infection promoted the high expression of PD-L1 [25–27]. Previous studies of our research group have also confirmed that high-risk HPV infection had higher levels of PD-L1 expression in cervical cancer patients [28]. In 2016, Keisuke et al. demonstrated that gene rearrangement on the 3-UTR of CD274 gene in tumor cells because of HPV gene integration would lead to abnormal expression of PD-L1 [29]. Also, some reports explored that the change of 9p24.1 could increase the high expression of PD-L1, because of 9q24.1 including many genes such as JAK2 and PDCD1LG2, which were closely related to PD-L1 expression [30]. Meanwhile, the latest research reported that the tumor-specific superenhancer in 9q24.1 could increase the expression of PD-L1 and PD-L2 [31]. So these genes in 9q24.1 were changed by HPV-DNA integration, which may create some kind of effect to affect PD-L1 expression. PD-L1 antibody drugs are considered as the star of tumor immunotherapy, and it has been proved that tumor mutation burden (TMB) was closely related to the efficacy of PD-L1 antibody drugs in a variety of tumor [32, 33]. The integration of viral genes into tumor cells has an effect on the genome of tumor cells, which was also a mutation burden in theory, so we called it viral gene integration burden (VIB). Existing theory proved that the potential cause of TMB to improve the efficacy of tumor immunotherapy is the production of neoantigen. But by the study, we found that the VIB similar to TMB in cervical cancer was closely related to PD-L1 expression. The number of gene integration and the presence of gene integration in 9p24.1 were both risk factors that promote the high expression of PD-L1. It was proved that this phenomenon of VIB not only promoted the production of neoantigen but also improved the expression of PD-L1 in tumor cells, thereby improving the efficacy of PD-L1 antibody immunotherapy.

5. Conclusion

The chip enrichment combined with the second-generation sequencing technology could significantly improve the detection sensitivity of HPV integration sites; multiple infectious problem was worthy of in-depth study, the HPV gene integrated into the human genome was universal and random, but also, there were some hotspots in human genome, and the hotspots in Chinese Uygur cervical cancer women were different from the other Chinese population. The breakpoints of HPV genome could lead to missed detection of the commercial HPV type detection kits. The integration of HPV genes into the human genome could affect pathways in cancer, which was a key step in the development of cervical cancer. The gene integration especially the integration existing on 9p24.1 could affect the expression level of PD-L1. The VIB should regulate the high expression of PD-L1 in some cancer cells because of oncogenic virus such as HBV and EB.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Ethical Approval

All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. The ethics committee of the Tumor Hospital Affiliated to Xinjiang Medical University approved the study and the consent procedure.

Informed consent was obtained from all individual participants included in the study.

Conflicts of Interest

The authors declare that they have no conflict of interest.

Acknowledgments

The study was funded by the Clinical Research Project of Tumor Hospital Affiliated to Xinjiang Medical University (no. 2019CL01002). The authors are grateful to the Department of Gynaecology Tumour III for the specimen source.

References

W. Chen, R. Zheng, P. D. Baade et al., “Cancer statistics in China, 2015,” CA: A Cancer Journal for Clinicians, vol. 66, no. 2, pp. 115–132, 2016.
View at: Publisher Site | Google Scholar
P. Guan, R. Howell-Jones, N. Li et al., “Human papillomavirus types in 115,789 HPV-positive women: a meta-analysis from cervical infection to cancer,” International Journal of Cancer, vol. 131, no. 10, pp. 2349–2359, 2012.
View at: Publisher Site | Google Scholar
F. Yang-chun, Z. Yuan, L. Cheng-ming, H. Yan-chun, and M. Xiu-min, “Increased HPV L1 gene methylation and multiple infection status lead to the difference of cervical epithelial cell lesion in different ethnic women of Xinjiang, China,” Medicine, vol. 96, no. 12, article e6409, 2017.
View at: Publisher Site | Google Scholar
F. Yang-chun, Y. Jia, and H. Yan-chun, “The application of CYFRA21-1 in cervical lesions screening in high-risk human papillomavirus infected women,” Ginekologia Polska, vol. 87, no. 9, pp. 617–620, 2016.
View at: Publisher Site | Google Scholar
R. Chetty, “70 years of the JCP-highly cited papers: the causal relation between human papillomavirus and cervical cancer,” Journal of Clinical Pathology, vol. 70, no. 12, p. 997, 2017.
View at: Publisher Site | Google Scholar
W. J. Shi, H. Liu, D. Wu, Z. H. Tang, Y. C. Shen, and L. Guo, “E6/E7 proteins are potential markers for the screening and diagnosis of cervical pre-cancerous lesions and cervical cancer in a Chinese population,” Oncology Letters, vol. 14, no. 5, pp. 6251–6258, 2017.
View at: Publisher Site | Google Scholar
G. Hudelist, M. Manavi, K. I. Pischinger et al., “Physical state and expression of HPV DNA in benign and dysplastic cervical tissue: different levels of viral integration are correlated with lesion grade,” Gynecologic Oncology, vol. 92, no. 3, pp. 873–880, 2004.
View at: Publisher Site | Google Scholar
B. Xu, S. Chotewutmontri, S. Wolf et al., “Multiplex identification of human Papillomavirus 16 DNA integration sites in cervical carcinomas,” PLoS One, vol. 8, no. 6, article e66693, 2013.
View at: Publisher Site | Google Scholar
K. L. Dall, C. G. Scarpini, I. Roberts et al., “Characterization of naturally occurring HPV16 integration sites isolated from cervical keratinocytes under noncompetitive conditions,” Cancer Research, vol. 68, no. 20, pp. 8249–8259, 2008.
View at: Publisher Site | Google Scholar
M. Schmitz, C. Driesch, K. Beer-Grondke, L. Jansen, I. B. Runnebaum, and M. Dürst, “Loss of gene function as a consequence of human papillomavirus DNA integration,” International Journal of Cancer, vol. 131, no. 5, pp. E593–E602, 2012.
View at: Publisher Site | Google Scholar
A. Gradíssimo Oliveira, C. Delgado, N. Verdasca, and A. Pista, “Prognostic value of human papillomavirus types 16 and 18 DNA physical status in cervical intraepithelial neoplasia,” Clinical Microbiology and Infection, vol. 19, no. 10, pp. E447–E450, 2013.
View at: Publisher Site | Google Scholar
Y. Shi, L. Li, Z. Hu et al., “A genome-wide association study identifies two new cervical cancer susceptibility loci at 4q12 and 17q12,” Nature Genetics, vol. 45, no. 8, pp. 918–922, 2013.
View at: Publisher Site | Google Scholar
A. Abudukadeer, S. Azam, A. Z. Mutailipu, L. Qun, G. Guilin, and S. Mijiti, “Knowledge and attitude of Uyghur women in Xinjiang province of China related to the prevention and early detection of cervical cancer,” World Journal of Surgical Oncology, vol. 13, no. 1, p. 110, 2015.
View at: Publisher Site | Google Scholar
A. Gnirke, A. Melnikov, J. Maguire et al., “Solution hybrid selection with ultra-long oligonucleotides for massively parallel targeted sequencing,” Nature Biotechnology, vol. 27, no. 2, pp. 182–189, 2009.
View at: Publisher Site | Google Scholar
L.-H. Zhao, X. Liu, H.-X. Yan et al., “Genomic and oncogenic preference of HBV integration in hepatocellular carcinoma,” Nature Communications, vol. 7, no. 1, article 12992, 2016.
View at: Publisher Site | Google Scholar
Z. Hu, D. Zhu, W. Wang et al., “Genome-wide profiling of HPV integration in cervical cancer identifies clustered genomic hot spots and a potential microhomology-mediated integration mechanism,” Nature Genetics, vol. 47, no. 2, pp. 158–163, 2015.
View at: Publisher Site | Google Scholar
W. Li, Y. Qi, X. Cui et al., “Characteristic of HPV integration in the genome and transcriptome of cervical cancer tissues,” BioMed Research International, vol. 2018, Article ID 6242173, 7 pages, 2018.
View at: Publisher Site | Google Scholar
G. Gao and D. I. Smith, “Role of the common fragile sites in cancers with a human papillomavirus etiology,” Cytogenetic and Genome Research, vol. 150, no. 3-4, pp. 217–226, 2016.
View at: Publisher Site | Google Scholar
Y.-C. Feng, Y.-c. Huang, J. Yang, C.-m. Liu, and Z.-z. Chen, “DNA ploidy of cervical epithelial cells should be a cure criterion of high-risk HPV infection in Xinjiang Uygur women,” OncoTargets and Therapy, vol. 8, pp. 827–833, 2015.
View at: Publisher Site | Google Scholar
D. Das, N. W. Smith, X. Wang, S. L. Richardson, M. C. T. Hartman, and I. M. Morgan, “Calcein represses human papillomavirus 16 E1-E2 mediated DNA replication via blocking their binding to the viral origin of replication,” Virology, vol. 508, pp. 180–187, 2017.
View at: Publisher Site | Google Scholar
L. E. King, J. C. Fisk, E. S. Dornan, M. M. Donaldson, T. Melendy, and I. M. Morgan, “Human papillomavirus E1 and E2 mediated DNA replication is not arrested by DNA damage signalling,” Virology, vol. 406, no. 1, pp. 95–102, 2010.
View at: Publisher Site | Google Scholar
W. A. A. Tjalma and C. E. Depuydt, “Cervical cancer screening: which HPV test should be used—L1 or E6/E7?” European Journal of Obstetrics & Gynecology and Reproductive Biology, vol. 170, no. 1, pp. 45-46, 2013.
View at: Publisher Site | Google Scholar
G. H. O. Cavalcante, J. M. G. de Araújo, J. V. Fernandes, and D. C. F. Lanza, “A seminested PCR assay for detection and typing of human papillomavirus based on E1 gene sequences,” Diagnostic Microbiology and Infectious Disease, vol. 91, no. 1, pp. 20–26, 2018.
View at: Publisher Site | Google Scholar
M. Satra, D. N. Vamvakopoulou, D. O. Sioutopoulou et al., “Sequence-based genotyping HPV L1 DNA and RNA transcripts in clinical specimens,” Pathology - Research and Practice, vol. 205, no. 12, pp. 863–869, 2009.
View at: Publisher Site | Google Scholar
A. M. Hong, P. Ferguson, T. Dodds et al., “Significant association of PD-L1 expression with human papillomavirus positivity and its prognostic impact in oropharyngeal cancer,” Oral Oncology, vol. 92, no. 5, pp. 33–39, 2019.
View at: Publisher Site | Google Scholar
Y. Meng, H. Liang, J. Hu et al., “PD-L1 expression correlates with tumor infiltrating lymphocytes and response to neoadjuvant chemotherapy in cervical cancer,” Journal of Cancer, vol. 9, no. 16, pp. 2938–2945, 2018.
View at: Publisher Site | Google Scholar
M. Choschzick, A. Gut, and D. Fink, “PD-L1 receptor expression in vulvar carcinomas is HPV-independent,” Virchows Archiv, vol. 473, no. 4, pp. 513–516, 2018.
View at: Publisher Site | Google Scholar
F. Yang-Chun, C. Zhen-Zhen, H. Yan-Chun, and M. Xiu-Min, “Association between PD-L1 and HPV status and the prognostic value for HPV treatment in premalignant cervical lesion patients,” Medicine, vol. 96, no. 25, article e7270, 2017.
View at: Publisher Site | Google Scholar
K. Kataoka, Y. Shiraishi, Y. Takeda et al., “Aberrant _PD-L1_ expression through 3-UTR disruption in multiple cancers,” Nature, vol. 534, no. 7607, pp. 402–406, 2016.
View at: Publisher Site | Google Scholar
M. R. Green, S. Monti, S. J. Rodig et al., “Integrative analysis reveals selective 9p24.1 amplification, increased PD-1 ligand expression, and further induction via JAK2 in nodular sclerosing Hodgkin lymphoma and primary mediastinal large B-cell lymphoma,” Blood, vol. 116, no. 17, pp. 3268–3277, 2010.
View at: Publisher Site | Google Scholar
Y. Xu, Y. Wu, S. Zhang et al., “A tumor-specific super-enhancer drives immune evasion by guiding synchronous expression of PD-L1 and PD-L2,” Cell Reports, vol. 29, no. 11, pp. 3435–3447.e4, 2019.
View at: Publisher Site | Google Scholar
T. Jiang, J. Shi, Z. Dong et al., “Genomic landscape and its correlations with tumor mutational burden, PD-L1 expression, and immune cells infiltration in Chinese lung squamous cell carcinoma,” Journal of Hematology & Oncology, vol. 12, no. 1, p. 75, 2019.
View at: Publisher Site | Google Scholar
S. Lu, J. E. Stein, D. L. Rimm et al., “Comparison of biomarker modalities for predicting response to PD-1/PD-L1 checkpoint Blockade,” JAMA Oncology, vol. 5, no. 8, p. 1195, 2019.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Feng Yang-chun et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies