Human Papilloma Virus (HPV) Integration Signature in Cervical Cancer: Identification of MACROD2 Gene As HPV Hot Spot Integration Site
Total Page:16
File Type:pdf, Size:1020Kb
www.nature.com/bjc ARTICLE Molecular Diagnostics Human papilloma virus (HPV) integration signature in Cervical Cancer: identification of MACROD2 gene as HPV hot spot integration site Maud Kamal1,2, Sonia Lameiras3, Marc Deloger4, Adeline Morel5, Sophie Vacher5, Charlotte Lecerf1,2, Célia Dupain1,2, Emmanuelle Jeannot 5,6, Elodie Girard4, Sylvain Baulande3, Coraline Dubot1,2, Gemma Kenter7, Ekaterina S. Jordanova 7,8, Els M. J. J. Berns9, Guillaume Bataillon6, Marina Popovic10, Roman Rouzier11,12, Wulfran Cacheux13, Christophe Le Tourneau 1,2,4,12, Alain Nicolas14, Nicolas Servant4, Suzy M. Scholl1,2, Ivan Bièche5,15 and RAIDs Consortium BACKGROUND: Cervical cancer (CC) remains a leading cause of gynaecological cancer-related mortality with infection by human papilloma virus (HPV) being the most important risk factor. We analysed the association between different viral integration signatures, clinical parameters and outcome in pre-treated CCs. METHODS: Different integration signatures were identified using HPV double capture followed by next-generation sequencing (NGS) in 272 CC patients from the BioRAIDs study [NCT02428842]. Correlations between HPV integration signatures and clinical, biological and molecular features were assessed. RESULTS: Episomal HPV was much less frequent in CC as compared to anal carcinoma (p < 0.0001). We identified >300 different HPV-chromosomal junctions (inter- or intra-genic). The most frequent integration site in CC was in MACROD2 gene followed by MIPOL1/TTC6 and TP63. HPV integration signatures were not associated with histological subtype, FIGO staging, treatment or PFS. HPVs were more frequently episomal in PIK3CA mutated tumours (p = 0.023). Viral integration type was dependent on HPV genotype (p < 0.0001); HPV18 and HPV45 being always integrated. High HPV copy number was associated with longer PFS (p = 0.011). CONCLUSIONS: This is to our knowledge the first study assessing the prognostic value of HPV integration in a prospectively annotated CC cohort, which detects a hotspot of HPV integration at MACROD2; involved in impaired PARP1 activity and chromosome instability. British Journal of Cancer (2021) 124:777–785; https://doi.org/10.1038/s41416-020-01153-4 BACKGROUND genome. Both patterns may be present jointly (episomal/ Cervical cancer (CC) remains a leading cause of gynaecological integrated).5 It is thought that the longer half-life of integrated cancer-related mortality worldwide and constitutes the second viral transcripts compared to half-life of episomal transcripts most common malignancy in women.1 Although patients with CC favours cellular immortalisation and transformation into cancer exhibit differences in clinical behaviour, infection by high-risk cells while also providing a selective growth advantage.6 Most human papilloma virus (HPV) remains an important initiating often, the integration of HPV DNA leads to a breakpoint in the E2 event in CC tumorigenesis,2 and one of the most important risk gene, resulting in de-repression of the E6 and E7 viral oncogenes. factors for developing CC.3 Most HPV infections are cleared When the virus remains episomal, expression of E6 and E7 spontaneously by the immune system, yet in some cases, it proteins may result from leaky expression or epigenetics persists leading to cancer.4 Following infection, the virus can dysregulation. E6 and E7 proteins impact the function of p53 remain in its episomal form, or become integrated into the host and pRb proteins, allowing squamous cell tumorigenesis.6 1Department of Drug Development and Innovation, Institut Curie, PSL Research University, 75005 Paris & 92210, Saint-Cloud, France; 2Department of Drug Development and Innovation, Institut Curie, PSL Research University, 92210 Saint-Cloud, France; 3Institut Curie, Genomics of Excellence (ICGex) Platform, PSL Research University, 75005 Paris, France; 4Bioinformatics and Computational Systems Biology of Cancer, PSL Research University, Mines Paris Tech, INSERM U900, 75005 Paris, France; 5Department of Genetics, Institut Curie, PSL Research University, 75005 Paris, France; 6Department of Pathology, Institut Curie, PSL Research University, 75005 Paris, France; 7Center for Gynaecologic Oncology Amsterdam, Amsterdam UMC and The Netherlands Cancer Institute - Antoni van Leeuwenhoek Hospital, Amsterdam, The Netherlands; 8Department of Pathology, Leiden University Medical Center, Leiden, The Netherlands; 9Department of Medical Oncology, Erasmus MC, 3000 CA Rotterdam, The Netherlands; 10Oncology Institute of Vojvodina, Put doktora Goldmana, 421204 Sremska Kamenica, Serbia; 11Department of Surgery, Institut Curie, PSL Research University, 92210 Saint-Cloud, France; 12Paris-Saclay University, Paris, France; 13Hopital Privé Pays de Savoie, Service d’oncologie médicale, 19 avenue Pierre Mendès France, 74100 Annemasse, France; 14Institut Curie, PSL Research University, CNRS UMR3244, 75248 Paris, France and 15Faculty of Pharmaceutical and Biological Sciences, INSERM U1016, Paris Descartes University, 75005 Paris, France Correspondence: Maud Kamal ([email protected]) Members of the RAIDs Consortium are listed end of the paper. These authors contributed equally: Suzy M. Scholl, Ivan Bièche. Received: 18 March 2020 Revised: 2 October 2020 Accepted: 15 October 2020 Published online: 16 November 2020 © The Author(s) 2020 Published by Springer Nature on behalf of Cancer Research UK Human papilloma virus (HPV) integration signature in Cervical Cancer:. M Kamal et al. 778 Several mechanisms of integration have been reported in the Table 1. Clinical and biological characteristics of 272 patients with literature; the “looping” model of HPV integration following DNA HPV-positive cervical cancer, in relation to progression-free survival. replication and recombination (resulting in DNA concatemers)7 is the most widely accepted but not experimentally reconstituted. Patients (%) Events (%)a PFS (p HPV DNA integration into the human genome triggers various value)b genetic alterations, such as oncogenes amplification, tumour suppressor gene inactivation, inter- or intra- chromosomal Total 272 (100.0) 84 (30.9) 6,8 rearrangements as well as genetic instability. Genes localised Age near the integration sites of viral genomes can experience ≤50 140 (51.5) 41 (29.3) 0.40 (NS) changes in RNA and protein expression levels, leading to over- >50 132 (48.5) 43 (32.6) or under-expression. In 2015, whole-genome sequencing and high-throughput viral integration methods identified as many as Histologic subtype 3667 HPV integration breakpoints in cervical neoplastic lesions. Squamous cell carcinoma 230 (84.6) 71 (30.9) 0.56 (NS) Frequent integration sites have been reported in genes relevant to Adenocarcinoma 27 (9.9) 8 (29.6) 9 the neoplastic process, such as the MYC oncogene. Loss of Adenosquamous 10 (3.7) 4 (40.0) function (LOF) in the RAD51B tumour suppressor gene following carcinoma HPV DNA insertion was reported to affect the DNA repair pathway 10 Mixed form or 5 (1.8) 1 (20.0) and genomic instability in CC. undifferentiated HPV DNA integration occurs as a single copy or in multiple HPV status repeats (in tandem or dispersed).11 In 2016, Holmes et al. developed a Capture HPV method to identify five different HPV Genotype 16 155 (57.0) 50 (32.3) 0.13 (NS) signatures in 72 CC. The first two signatures contain two hybrid Genotype 18 36 (13.2) 14 (38.9) chromosomal–HPV junctions which are co-linear (2 Junctions Genotype 45 27 (9.9) 10 (37.0) “ ” “ ” Colinear 2J-COL ) or non-linear (2 Junctions Non-Linear 2J-NL ) Genotype 31 9 (3.3) 0 (0) depending on their relative orientations. It reflects two modes of Genotype 33 11 (4.0) 0 (0) viral integration, associated with chromosomal deletion or c amplification events, respectively. The third and fourth signatures Other genotypes 34 (12.5) 10 (29.4) 1234567890();,: exhibit several hybrid junctions either clustered in one chromo- FIGO stage somal region (Multiple Junctions Clustered “MJ-CL”) or scattered at I/II 205 (75.4) 50 (24.4) <0.0001 “ ” fi distinct loci (Multiple Junctions Scattered MJ-SC ) while the fth III/IV 67 (24.6) 34 (50.7) signature consists of episomal forms of HPV (EPI).12 Nodal involvement On the assumption that HPV integration types/signatures/ pattern might predict clinical outcomes, we analysed the Yes 167 (61.4) 62 (37.1) 0.0028 association between the different viral integration signatures, No 105 (38.6) 22 (21.0) clinical and pathological parameters and outcome in the large Pelvic lymph nodes cohort of 272 HPV-positive CC patients enrolled in the prospective Yes 165 (60.7) 62 (37.6) 0.0015 BioRAIDs study [NCT02428842]. No 107 (39.3) 22 (20.6) Para-aortic lymph nodes METHODS Yes 43 (15.8) 22 (51.2) 0.0001 Patients and samples No 229 (84.2) 62 (27.1) Patients included in this study were enrolled in the EU-funded Initial therapy RAIDs Network (Rational Molecular Assessment and Innovative Surgery 54 (19.9) 10 (18.5) 0.0008 Drug Selection, www.raids-fp7.eu) prospective CC BioRAIDs study [NCT02428842]. The clinical protocol together with tumour Radiotherapy 176 (64.7) 52 (29.5) sampling procedures, quality control of samples and treatment Neoadjuvant 42 (15.4) 22 (52.4) in 18 European centres (seven European countries) as well as chemotherapy – study results have been previously published.13 15 PIK3CA mutational statusd WT 182 (67.7) 60 (33.0) 0.23 (NS) HPV typing Mutated 87 (32.3) 23 (26.4) All samples included in this study were analysed for HPV type, using the SPF10 primer set and INNO-LiPA HPV genotyping extra Significant results are displayed