HPV type-related chromosomal profiles in high-grade cervical intraepithelial neoplasia

Use of peripheral blood transcriptomic biomarkers to distinguish high-grade cervical squamous intraepithelial lesions from low-grade lesions

It is crucial to classify cervical lesions into high-grade squamous intraepithelial lesions (HSILs) and low-grade SILs (LSILs), as LSILs are conservatively treated by observation, based on an expectation of natural regression, whereas HSILs usually require electrosurgical excision. In the present study, peripheral blood gene expression profiles were analyzed to identify transcriptomic biomarkers distinguishing HSILs from LSILs. A total of 102 blood samples were collected from women with cervical SILs (66 HSIL and 36 LSIL) for microarray hybridization. Candidate gene signatures were identified using AdaBoost algorithms, and a predictive model was constructed using logistic regression to differentiate HSILs from LSILs. To correct for possible bias as a result of the limited sample size and to verify the stability of the predictive model, a two-fold cross validation and null set analysis was conducted over 1,000 iterations. The functions of the transcriptomic biomarkers were then analyzed to elucidate the pathogenesis of cervical SIL. A total of 10 transcriptomic genes (STMN3, TRPC4AP, DYRK2, AGK, KIAA0319L, GRPEL1, ZFC3H1, LYL1, ITGB1 and ARHGAP18) were identified. The predictive model based on the 10-gene panel exhibited well-discriminated power. A cross validation process using known disease status exhibited almost the same performance as that of the predictive model, whereas null-set analysis with randomly reassigned disease status exhibited much lower predictive performance for distinguishing HSILs from LSILs. These biomarkers were involved in the 'Rho GTPase cycle', 'mitochondrial protein import', 'oncogenic MAPK signaling', 'integrin cell surface interaction' and 'signaling by BRAF and RAF fusions'. In conclusion, peripheral blood gene expression analysis is a promising method for distinguishing HSILs from LSILs. The present study proposes 10 candidate genes that could be used in the future as diagnostic biomarkers and potential therapeutic targets for cervical SILs. A simple, non-invasive blood test would be clinically useful in the diagnosis and classification of patients with cervical SILs.

Expression of Pregnancy Specific β-1 Glycoprotein 1 in Cervical Cancer Cells

BACKGROUND

Cervical Cancer (CC) is a worldwide public health concern associated with genetic alterations, among these the gain of the 19q chromosome harboring the Pregnancy Specific Glycoproteins (PSG) gene family. These proteins play a critical role in pregnancy, with participation in immunotolerance, angiogenesis, and invasion processes, which are also observed in carcinogenesis. The aim of this study was to determine the molecular alterations of PSG1 and its relationship with CC.

METHODS

PSG1 Copy Number Variation (CNV) was evaluated in 31 CC and eight normal cervical tissues by qPCR. PSG1 expression was correlated with HPV detection and IL-10 and TGF-β expression in CC samples. Finally, PSG1 protein expression was evaluated by immunofluorescence in CC cell lines, by immunohistochemistry in a tissue microarray, and by immunoblotting in the sera of women with normal cervix, pre-invasive lesions, and CC.

RESULTS

PSG1 showed a gain of 25.6% in CNV and gene expression in CC. There was a lack of PSG1 expression in normal cervical epithelium and positive immunostaining in 57% of CC tissues, while all CC cell lines expressed PSG1. Finally, PSG1 was immunodetected in 90% of pre-invasive lesions and in all CC serum samples, but not in healthy women. PSG1 expression correlates with the expression of IL-10 and TGF-β in CC tissues, but not with the presence of HPV.

CONCLUSION

These data show evidence of the differential expression of PSG1 in CC that could explain its participation in tumor-biology and immunotolerance mechanisms. Further, its immunodetection could provide early detection of this cancer.

CADM1, MAL, and miR124 Promoter Methylation as Biomarkers of Transforming Cervical Intrapithelial Lesions

Background: Squamous intraepithelial lesions/cervical intraepithelial neoplasias (SIL/CIN) are high-risk human papilloma virus (hrHPV)-related lesions which are considered as high grade (HSIL/CIN2-3) or low grade (LSIL/CIN1) lesions according to their risk of progression to cervical cancer (CC). Most HSIL/CIN2-3 are considered as transforming hrHPV infections, so truly CC precursors, although some clear spontaneously. hrHPV testing has a high sensitivity for the detection of HSIL/CIN2-3 but a relatively low specificity for identifying transforming lesions. We aimed to determine whether the combination of CADM1, MAL and miR124 promoter methylation status assessed in histological samples can be used as a biomarker in the identification of transforming HSIL/CIN lesions. Design: 131 cervical biopsies, including 8 cases with no lesion and a negative hrHPV test result (control group), 19 low-grade (L)SIL/CIN1, 30 HSIL/CIN2, 60 HSIL/CIN3, and 14 CC were prospectively collected. hrHPV was detected and genotyped using the polymerase chain reaction (PCR)-based technique SPF10 HPV LIPA. A multiplex quantitative methylation-specific PCR (qMSP) was used to identify the methylation status of the CADM1, MAL, and miR124 promoter genes. Results: Significantly higher methylation levels of CADM1, MAL and miR-124 were found in HSIL/CIN2-3 and CC compared with normal and LSIL lesions. DNA methylation of at least one gene was detected in 12.5% (1/8) of normal samples, 31.5% (6/19) of LSIL/CIN1, 83.3% (25/30) of HSIL/CIN2, 81.6% (49/60) of HSIL/CIN3 and 100% (14/14) of CC (p < 0.001). The sensitivity and specificity for HSIL/CIN2-3 and CC of having at least one methylated gene were 84.6% and 74.0%, respectively. The sensitivity and specificity of the combination of at least one methylated gene and a positive hrHPV test were 80.7% and 85.1% for HSIL/CIN2-3 and CC, respectively. Conclusions: The methylation rate of CADM1, MAL and miR124 increases with the severity of the lesion. Further research is warranted to evaluate the usefulness of these biomarkers for the identification of transforming HSIL/CIN.

Molecular heterogeneity in human papillomavirus-dependent and -independent vulvar carcinogenesis

Vulvar squamous cell carcinoma (VSCC) and precancerous vulvar intraepithelial neoplasia (VIN) can develop through human papillomavirus (HPV)-dependent and -independent pathways, indicating a heterogeneous disease. Only a minority of VIN progress, but current clinicopathological classifications are insufficient to predict the cancer risk. Here we analyzed copy number alterations (CNA) to assess the molecular heterogeneity of vulvar lesions in relation to HPV and cancer risk. HPV-status and CNA by means of whole-genome next-generation shallow-sequencing were assessed in VSCC and VIN. The latter included VIN of women with associated VSCC (VINVSCC ) and women who did not develop VSCC during follow-up (VINnoVSCC ). HPV-testing resulted in 41 HPV-positive (16 VINVSCC , 14 VINnoVSCC , and 11 VSCC) and 24 HPV-negative (11 VINVSCC and 13 VSCC) lesions. HPV-positive and -negative VSCC showed a partially overlapping pattern of recurrent CNA, including frequent gains of 3q and 8q. In contrast, amplification of 11q13/cyclinD1 was exclusively found in HPV-negative lesions. HPV-negative VINVSCC had less CNA than HPV-negative VSCC (P = .009), but shared chromosome 8 alterations. HPV-positive VINnoVSCC had less CNA than VINVSCC (P = .022). Interestingly, 1pq gain was detected in 81% of HPV-positive VINVSCC and only in 21% of VINnoVSCC (P = .001). In conclusion, HPV-dependent and -independent vulvar carcinogenesis is characterized by distinct CNA patterns at the VIN stage, while more comparable patterns are present at the cancer stage. Cancer risk in VIN seems to be reflected by the extent of CNA, in particular chromosome 1 gain in HPV-positive cases.

Next generation sequencing and its applications in HPV-associated cancers

Approximately 18% of all human cancers have a viral etiology, and human papillomavirus (HPV) has been identified as one of the most prevalent viruses that plays causative role in nearly all cervical cancers and, in addition, in subset of head and neck, anal, penile and vulvar cancers. The recent introduction of next generation sequencing (NGS) and other omics approaches have resulted in comprehensive knowledge on the pathogenesis of HPV-driven tumors. Specifically, these approaches have provided detailed information on genomic HPV integration sites, disrupted genes and pathways, and common and distinct genetic and epigenetic alterations in different human HPV-associated cancers. This review focuses on HPV integration sites, its concomitantly disrupted genes and pathways and its functional consequences in both cervical and head and neck cancers. Integration of NGS data with other omics and clinical data is crucial to better understand the pathophysiology of each individual malignancy and, based on this, to select targets and to design effective personalized treatment options.

Differential in vitro immortalization capacity of eleven (probable) [corrected] high-risk human papillomavirus types

Epidemiological studies identified 12 high-risk HPV (hrHPV) types and 8 probable/possible hrHPV types that display different cancer risks. Functional studies on transforming properties of hrHPV are mainly limited to HPV16 and -18, which induce immortalization of human foreskin keratinocytes (HFKs) by successive bypass of two proliferative life span barriers, senescence and crisis. Here, we systematically compared the in vitro immortalization capacities, as well as influences on p53, pRb, hTERT, growth behavior, and differentiation capacity, of nine hrHPV types (HPV16, -18, -31, -33, -35, -45, -51, -52, and -59), and two probable hrHPV types (HPV66 and -70). By retroviral transduction, the respective E6/E7 coding sequences were expressed in HFKs from two or three independent donors. Reduced p53 levels and low-level hTERT expression in early-passage cells, as seen in HPV16-, -31-, -33-, and -35-, and to a lesser extent HPV18-transduced HFKs, was associated with continuous growth and an increased immortalization capacity. Less frequent immortalization by HPV45 and -51 and immortalization by HPV66 and -70 was preceded by an intervening period of strongly reduced growth (crisis) without prior increase in hTERT expression. Immortalization by HPV59 was also preceded by a period crisis, despite the onset of low hTERT expression at early passage. HPV52 triggered an extended life span but failed to induce immortality. Variations in p53 and pRb levels were not correlated with differences in alternative E6/E7 mRNA splicing in all hrHPV-transduced HFKs. On collagen rafts, transductants showed disturbed differentiation reminiscent of precancerous lesions. In conclusion, in vitro oncogenic capacities differ between the established hrHPV types, and both some established and probable hrHPV types display weak or moderate immortalization potential.

Chromosomal gains and losses in human papillomavirus-associated neoplasia of the lower genital tract - a systematic review and meta-analysis

BACKGROUND

Overexpression of the human papillomavirus (HPV) oncogenes E6 and E7 is necessary for the development of distinct lower genital tract cancers. However, secondary cellular genomic alterations are mandatory to promote progression of HPV-induced premalignant stages. We aimed at identifying the chromosomal regions most frequently gained and lost and the disease stage at which the latter occurs. These regions might be relevant for carcinogenesis and could serve as diagnostic markers to identify premalignant lesions with high progression risk towards invasive cancer.

METHODS

We performed a systematic literature review and meta-analysis of studies listed in PubMed that analysed chromosomal copy number alterations by comparative genomic hybridisation (CGH) in HPV-positive and -negative cancers or premalignant lesions of the anogenital tract (cervix, anus, vagina, penis and vulva).

FINDINGS

Data were extracted and analysed from 32 studies. The most common alterations in cervical squamous cell carcinoma (SCC) (12 studies, 293 samples) were gains at 3q with a rate of 0.55 (95% confidence interval (CI) 0.43-0.70), losses at 3p (0.36, 95%CI 0.27-0.48) and losses at 11q (0.33, 95%CI 0.26-0.43). Gains at 3q were particularly frequent in HPV16-positive cervical SCC (0.84, 95%CI 0.78-0.90). Also more than one quarter of high grade cervical intraepithelial neoplasia (CIN) harboured gains of 3q (0.27, 95%CI 0.20-0.36), but the rate in low grade CIN was low (0.02, 95%CI 0.00-0.09). For HPV-associated vulvar SCC (four studies, 30 samples) the same common alterations as in cervical SCC were reported. Studies on non-cervical and non-vulvar SCC and premalignant lesions of the lower genital tract are scarce.

INTERPRETATION

3q gains were most frequently found in HPV16-positive cervical SCC. The results suggest the selection of HPV-transformed cell clones harbouring 3q gains in high grade premalignant lesions, while alterations in low grade lesions are rare.

Chromosomal gains measured in cytology samples from women with abnormal cervical cancer screening results

OBJECTIVE

Chromosomal gains at 3q26, 5p15 and 20q13 have been described in cervical precancer and cancer. We evaluated a novel fluorescence in situ hybridization (FISH) assay that detects gains at these three loci simultaneously as a possible biomarker for detecting cervical precancer.

METHODS

Chromosomal copy numbers at 3q26, 5p15, 20q13 and the centromere of chromosome7 (cen7) in liquid-based cytology specimens from 168 women enrolled in the Biopsy Study were determined by FISH. The number of cells with ≥ 3 or ≥ 4 signals for a genomic locus was enumerated and diagnostic test performance measures were calculated using receiver operating characteristic (ROC) analyses. Sensitivity and specificity values were determined for the detection of CIN2+ and/or HSIL.

RESULTS

The median number of cells with ≥ 3 signals increased with the severity of cervical lesion for each genomic locus (p-trend<0.02 for each locus). ROC analysis for the number of cells with ≥ 3 signals resulted in area under the curve values of 0.70 (95% CI: 0.54-0.86), 0.67 (0.52-0.83), 0.67 (0.51-0.83) and 0.78 (0.64-0.92) for 3q26, 5p15, 20q13 and cen7, respectively, for the detection of CIN2+ and/or HSIL. Positivity for gains at multiple loci resulted in only slightly better test performance measures than those for the individual probes for four distinct combinations of probes.

CONCLUSIONS

Chromosomal gains at 3q26, 5p15, 20q13 and cen7 are associated with severity of cervical lesions. Further studies are required to quantify risk stratification of FISH assays for cervical cancer screening.

Identification of eight candidate target genes of the recurrent 3p12-p14 loss in cervical cancer by integrative genomic profiling

The pathogenetic role, including its target genes, of the recurrent 3p12-p14 loss in cervical cancer has remained unclear. To determine the onset of the event during carcinogenesis, we used microarray techniques and found that the loss was the most frequent 3p event, occurring in 61% of 92 invasive carcinomas, in only 2% of 43 high-grade intraepithelial lesions (CIN2/3), and in 33% of 6 CIN3 lesions adjacent to invasive carcinomas, suggesting a role in acquisition of invasiveness or early during the invasive phase. We performed an integrative DNA copy number and expression analysis of 77 invasive carcinomas, where all genes within the recurrent region were included. We selected eight genes, THOC7, PSMD6, SLC25A26, TMF1, RYBP, SHQ1, EBLN2, and GBE1, which were highly down-regulated in cases with loss, as confirmed at the protein level for RYBP and TMF1 by immunohistochemistry. The eight genes were subjected to network analysis based on the expression profiles, revealing interaction partners of proteins encoded by the genes that were coordinately regulated in tumours with loss. Several partners were shared among the eight genes, indicating crosstalk in their signalling. Gene ontology analysis showed enrichment of biological processes such as apoptosis, proliferation, and stress response in the network and suggested a relationship between down-regulation of the eight genes and activation of tumourigenic pathways. Survival analysis showed prognostic impact of the eight-gene signature that was confirmed in a validation cohort of 74 patients and was independent of clinical parameters. These results support the role of the eight candidate genes as targets of the 3p12-p14 loss in cervical cancer and suggest that the strong selection advantage of the loss during carcinogenesis might be caused by a synergetic effect of several tumourigenic processes controlled by these targets.

Focal aberrations indicate EYA2 and hsa-miR-375 as oncogene and tumor suppressor in cervical carcinogenesis

Cervical cancer results from persistent infection with high-risk human papillomavirus (hrHPV). Common genetic aberrations in cervical (pre)cancers encompass large genomic regions with numerous genes, hampering identification of driver genes. This study aimed to identify genes functionally involved in HPV-mediated transformation by analysis of focal aberrations (<3 Mb) in high-grade cervical intraepithelial neoplasia (hgCIN). Focal chromosomal aberrations were determined in high-resolution array comparative genomic hybridization data of 60 hgCIN. Genes located within focal aberrations were validated using 2 external gene expression datasets or qRT-PCR. Functional roles of candidate genes EYA2 (20q13) and hsa-miR-375 (2q35) were studied by siRNA-mediated knock-down and overexpression, respectively, in hrHPV-containing cell lines. We identified 74 focal aberrations encoding 305 genes. Concurrent altered expression in hgCIN and/or cervical carcinomas compared with normal cervical samples was shown for ATP13A3, HES1, OPA1, HRASLS, EYA2, ZMYND8, APOBEC2, and NCR2. Gene silencing of EYA2 significantly reduced viability, migratory capacity, and anchorage-independent growth of HPV16-transformed keratinocytes. For hsa-miR-375, a direct correlation between a (focal) loss and significantly reduced expression was found. Downregulation of hsa-miR-375 expression was confirmed in an independent series of cervical tissues. Ectopic expression of hsa-miR-375 in 2 cervical carcinoma cell lines reduced cellular viability. Our data provide a proof of concept that chromosomal aberrations are actively contributing to HPV-induced carcinogenesis and identify EYA2 and hsa-miR-375 as oncogene and tumor suppressor gene, respectively.