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.

Gene dosage, expression, and ontology analysis identifies driver genes in the carcinogenesis and chemoradioresistance of cervical cancer.

Integrative analysis of gene dosage, expression, and ontology (GO) data was performed to discover driver genes in the carcinogenesis and chemoradioresistance of cervical cancers. Gene dosage and expression profiles of 102 locally advanced cervical cancers were generated by microarray techniques. Fifty-two of these patients were also analyzed with the Illumina expression method to confirm the gene expression results. An independent cohort of 41 patients was used for validation of gene expressions associated with clinical outcome. Statistical analysis identified 29 recurrent gains and losses and 3 losses (on 3p, 13q, 21q) associated with poor outcome after chemoradiotherapy. The intratumor heterogeneity, assessed from the gene dosage profiles, was low for these alterations, showing that they had emerged prior to many other alterations and probably were early events in carcinogenesis. Integration of the alterations with gene expression and GO data identified genes that were regulated by the alterations and revealed five biological processes that were significantly overrepresented among the affected genes: apoptosis, metabolism, macromolecule localization, translation, and transcription. Four genes on 3p (RYBP, GBE1) and 13q (FAM48A, MED4) correlated with outcome at both the gene dosage and expression level and were satisfactorily validated in the independent cohort. These integrated analyses yielded 57 candidate drivers of 24 genetic events, including novel loci responsible for chemoradioresistance. Further mapping of the connections among genetic events, drivers, and biological processes suggested that each individual event stimulates specific processes in carcinogenesis through the coordinated control of multiple genes. The present results may provide novel therapeutic opportunities of both early and advanced stage cervical cancers.

Integrative Analysis of DCE-MRI and Gene Expression Profiles in Construction of a Gene Classifier for Assessment of Hypoxia-Related Risk of Chemoradiotherapy Failure in Cervical Cancer.

PURPOSE: A 31-gene expression signature reflected in dynamic contrast enhanced (DCE)-MR images and correlated with hypoxia-related aggressiveness in cervical cancer was identified in previous work. We here aimed to construct a dichotomous classifier with key signature genes and a predefined classification threshold that separated cervical cancer patients into a more and less hypoxic group with different outcome to chemoradiotherapy. EXPERIMENTAL DESIGN: A training cohort of 42 patients and two independent cohorts of 108 and 131 patients were included. Gene expression data were generated from tumor biopsies by two Illumina array generations (WG-6, HT-12). Technical transfer of the classifier to a reverse transcription quantitative PCR (RT-qPCR) platform was performed for 74 patients. The amplitude ABrix in the Brix pharmacokinetic model was extracted from DCE-MR images of 64 patients and used as an indicator of hypoxia. RESULTS: Classifier candidates were constructed by integrative analysis of ABrix and gene expression profiles in the training cohort and evaluated by a leave-one-out cross-validation approach. On the basis of their ability to separate patients correctly according to hypoxia status, a 6-gene classifier was identified. The classifier separated the patients into two groups with different progression-free survival probability. The robustness of the classifier was demonstrated by successful validation of hypoxia association and prognostic value across cohorts, array generations, and assay platforms. The prognostic value was independent of existing clinical markers, regardless of clinical endpoints. CONCLUSIONS: A robust DCE-MRI-associated gene classifier has been constructed that may be used to achieve an early indication of patients' risk of hypoxia-related chemoradiotherapy failure. Clin Cancer Res; 22(16); 4067-76. ©2016 AACR.

Hypoxia-induced alterations of G2 checkpoint regulators.

Hypoxia promotes an aggressive tumor phenotype with increased genomic instability, partially due to downregulation of DNA repair pathways. However, genome stability is also surveilled by cell cycle checkpoints. An important issue is therefore whether hypoxia also can influence the DNA damage-induced cell cycle checkpoints. Here, we show that hypoxia (24 h 0.2% O2) alters the expression of several G2 checkpoint regulators, as examined by microarray gene expression analysis and immunoblotting of U2OS cells. While some of the changes reflected hypoxia-induced inhibition of cell cycle progression, the levels of several G2 checkpoint regulators, in particular Cyclin B, were reduced in G2 phase cells after hypoxic exposure, as shown by flow cytometric barcoding analysis of individual cells. These effects were accompanied by decreased phosphorylation of a Cyclin dependent kinase (CDK) target in G2 phase cells after hypoxia, suggesting decreased CDK activity. Furthermore, cells pre-exposed to hypoxia showed increased G2 checkpoint arrest upon treatment with ionizing radiation. Similar results were found following other hypoxic conditions (∼0.03% O2 20 h and 0.2% O2 72 h). These results demonstrate that the DNA damage-induced G2 checkpoint can be altered as a consequence of hypoxia, and we propose that such alterations may influence the genome stability of hypoxic tumors.

Aberrant methylation-mediated silencing of microRNAs contributes to HPV-induced anchorage independence.

Cervical cancer and a subset of anogenital and head-and-neck carcinomas are caused by high-risk types of the human papillomavirus (hrHPV). During hrHPV-induced malignant transformation keratinocytes become able to grow anchorage independently, a tumorigenic trait at least partly associated with inactivation of tumor suppressor genes. We used hrHPV-containing keratinocytes to investigate the role of DNA methylation-mediated silencing of microRNAs (miRNAs) in the acquisition of anchorage independence.Anchorage dependent (n=11) and independent passages (n=19) of 4 hrHPV-immortalized keratinocyte cell lines were treated with 2'-deoxy-5-azacytidine (DAC). Genome-wide miRNA expression profiles before and after treatment were compared to identify miRNAs silenced by methylation. Bisulfite sequencing and methylation-specific PCR showed increased methylation of hsa-mir-129-2/-137/-935/-3663/-3665 and -4281 in anchorage independent HPV-transformed keratinocytes and cervical cancer cell lines. Mature miRNAs derived from hsa-mir-129-2/-137/-3663 and -3665 showed functional relevance as they decreased anchorage independence in cervical cancer cell lines. Cervical (pre)cancerous lesions demonstrated increased methylation of hsa-mir-129-2/-935/-3663/-3665 and -4281, underlining the clinical relevance of our findings.In conclusion, methylation-mediated silencing of tumor suppressive miRNAs contributes to acquisition of an anchorage independent phenotype. This study further substantiates the importance of miRNAs during early stages of carcinogenesis and underlines their potential as both disease markers and therapeutic targets.

Interplay between promoter methylation and chromosomal loss in gene silencing at 3p11-p14 in cervical cancer.

Loss of 3p11-p14 is a frequent event in epithelial cancer and a candidate prognostic biomarker in cervical cancer. In addition to loss, promoter methylation can participate in gene silencing and promote tumor aggressiveness. We have performed a complete mapping of promoter methylation at 3p11-p14 in two independent cohorts of cervical cancer patients (n = 149, n = 121), using Illumina 450K methylation arrays. The aim was to investigate whether hyperm-ethylation was frequent and could contribute to gene silencing and disease aggressiveness either alone or combined with loss. By comparing the methylation level of individual CpG sites with corresponding data of normal cervical tissue, 26 out of 41 genes were found to be hypermethylated in both cohorts. The frequency of patients with hypermethylation of these genes was found to be higher at tumor stages of 3 and 4 than in stage 1 tumors. Seventeen of the 26 genes were transcriptionally downregulated in cancer compared to normal tissue, whereof 6 genes showed a significant correlation between methylation and expression. Integrated analysis of methylation, gene dosage, and expression of the 26 hypermethylated genes identified 3 regulation patterns encompassing 8 hypermethylated genes; a methylation driven pattern (C3orf14, GPR27, ZNF717), a gene dosage driven pattern (THOC7, PSMD6), and a combined methylation and gene dosage driven pattern (FHIT, ADAMTS9, LRIG1). In survival analysis, patients with both hypermethylation and loss of LRIG1 had a worse outcome compared to those harboring only hypermethylation or none of the events. C3orf14 emerged as a novel methylation regulated suppressor gene, for which knockdown was found to promote invasive growth in human papilloma virus (HPV)-transformed keratinocytes. In conclusion, hypermethylation at 3p11-p14 is common in cervical cancer and may exert a selection pressure during carcinogenesis alone or combined with loss. Information on both events could lead to improved prognostic markers.

A new and improved host-independent plasmid system for RK2-based conjugal transfer.

Bacterial conjugation is a process that is mediated either by a direct cell-to-cell junction or by formation of a bridge between the cells. It is often used to transfer DNA constructs designed in Escherichia coli to recipient bacteria, yeast, plants and mammalian cells. Plasmids bearing the RK2/RP4 origin of transfer (oriT) are mostly mobilized using the E. coli S17-1/SM10 donor strains, in which transfer helper functions are provided from a chromosomally integrated RP4::Mu. We have observed that large plasmids were occasionally modified after conjugal transfer when using E. coli S17-1 as a donor. All modified plasmids had increased in size, which most probably was a result of co-transfer of DNA from the chromosomally located oriT. It has earlier also been demonstrated that the bacteriophage Mu is silently transferred to recipient cells by these donor strains, and both occurrences are very likely to lead to mutations within the recipient DNA. Here we report the construction of a new biological system addressing both the above mentioned problems in which the transfer helper functions are provided by a plasmid lacking a functional oriT. This system is compatible with all other replicons commonly used in conjugation experiments and further enables the use of diverse bacterial strains as donors. Plasmids containing large inserts were successfully conjugated and the plasmid modifications observed when E. coli S17-1 was used as donor were eliminated by the use of the new host-independent vector system.

Gene network reconstruction reveals cell cycle and antiviral genes as major drivers of cervical cancer.

Although human papillomavirus was identified as an aetiological factor in cervical cancer, the key human gene drivers of this disease remain unknown. Here we apply an unbiased approach integrating gene expression and chromosomal aberration data. In an independent group of patients, we reconstruct and validate a gene regulatory meta-network, and identify cell cycle and antiviral genes that constitute two major subnetworks upregulated in tumour samples. These genes are located within the same regions as chromosomal amplifications, most frequently on 3q. We propose a model in which selected chromosomal gains drive activation of antiviral genes contributing to episomal virus elimination, which synergizes with cell cycle dysregulation. These findings may help to explain the paradox of episomal human papillomavirus decline in women with invasive cancer who were previously unable to clear the virus.

Hypoxia-induced gene expression in chemoradioresistant cervical cancer revealed by dynamic contrast-enhanced MRI.

Knowledge of the molecular background of functional magnetic resonance (MR) images is required to fully exploit their potential in cancer management. We explored the prognostic impact of dynamic contrast-enhanced MR imaging (DCE-MRI) parameters in cervical cancer combined with global gene expression data to reveal their underlying molecular phenotype and construct a representative gene signature for the relevant parameter. On the basis of 78 patients with cervical cancer subjected to curative chemoradiotherapy, we identified the prognostic DCE-MRI parameter A(Brix) by pharmacokinetic analysis of pretreatment images based on the Brix model, in which tumors with low A(Brix) appeared to be most aggressive. Gene set analysis of 46 tumors with pairwise DCE-MRI and gene expression data showed a significant correlation between A(Brix) and the hypoxia gene sets, whereas gene sets related to other tumor phenotypes were not significant. Hypoxia gene sets specific for cervical cancer created in cell culture experiments, including both targets of the hypoxia inducible factor (HIF1α) and the unfolded protein response, were the most significant. In the remaining 32 tumors, low A(Brix) was associated with upregulation of HIF1α protein expression, as assessed by immunohistochemistry, consistent with increased hypoxia. On the basis of the hypoxia gene sets, a signature of 31 genes that were upregulated in tumors with low A(Brix) was constructed. This DCE-MRI hypoxia gene signature showed prognostic impact in an independent validation cohort of 109 patients. Our findings reveal the molecular basis of an aggressive hypoxic phenotype and suggest the use of DCE-MRI to noninvasively identify patients with hypoxia-related chemoradioresistance.

Phosphorylation of EGFR measured with in situ proximity ligation assay: relationship to EGFR protein level and gene dosage in cervical cancer.

PURPOSE: We have applied the sensitive and specific in situ proximity ligation assay (PLA) to characterize Tyr1068 phosphorylation of the epidermal growth factor receptor (EGFR) in cervical cancer in relation to the protein level and gene dosage. MATERIALS AND METHODS: Pretreatment tumor biopsies from 178 patients were analyzed. EGFR protein level was determined by immunohistochemistry, and Tyr1068 phosphorylation was detected with PLA in 97 EGFR positive tumors. EGFR gene dosage was derived from array comparative genomic hybridization of 86 cases. RESULTS: EGFR was expressed in most tumors, whereas phosphorylation was seen in about half of the EGFR positive ones. A correlation was found between the expression of EGFR and phosphorylated EGFR (p=0.016, membrane; p=0.012, cytoplasm). However, tumor regions with high protein level without phosphorylation were occasionally seen and the percentage of EGFR positive cells was higher than the phosphorylated percentage (p<0.001). Moreover, an increase in the phosphorylation in both the membrane (p=0.014) and cytoplasm (p=0.002) was seen in 11 tumors with gain of EGFR. The protein level was not correlated with gene dosage. CONCLUSION: In contrast to gain of the EGFR chromosomal region, high EGFR protein level may not necessarily indicate Tyr1068 phosphorylation and thereby receptor activation in cervical cancer.

Membranous expression of ectodomain isoforms of the epidermal growth factor receptor predicts outcome after chemoradiotherapy of lymph node-negative cervical cancer.

PURPOSE: We compared the prognostic significance of ectodomain isoforms of the epidermal growth factor receptor (EGFR), which lack the tyrosine kinase (TK) domain, with that of the full-length receptor and its autophosphorylation status in cervical cancers treated with conventional chemoradiotherapy. EXPERIMENTAL DESIGN: Expression of EGFR isoforms was assessed by immunohistochemistry in a prospectively collected cohort of 178 patients with squamous cell cervical carcinoma, and their detection was confirmed with Western blotting and reverse transcriptase PCR. A proximity ligation immunohistochemistry assay was used to assess EGFR-specific autophosphorylation. Pathways associated with the expression of ectodomain isoforms were studied by gene expression analysis with Illumina beadarrays in 110 patients and validated in an independent cohort of 41 patients. RESULTS: Membranous expression of ectodomain isoforms alone, without the coexpression of the full-length receptor, showed correlations to poor clinical outcome that were highly significant for lymph node-negative patients (locoregional control, P = 0.0002; progression-free survival, P < 0.0001; disease-specific survival, P = 0.005 in the log-rank test) and independent of clinical variables. The ectodomain isoforms were primarily 60-kD products of alternative EGFR transcripts. Their membranous expression correlated with transcriptional regulation of oncogenic pathways including activation of MYC and MAX, which was significantly associated with poor outcome. This aggressive phenotype of ectodomain EGFR expressing tumors was confirmed in the independent cohort. Neither total nor full-length EGFR protein level, or autophosphorylation status, showed prognostic significance. CONCLUSION: Membranous expression of ectodomain EGFR isoforms, and not TK activation, predicts poor outcome after chemoradiotherapy for patients with lymph node-negative cervical cancer.

GeneCount: genome-wide calculation of absolute tumor DNA copy numbers from array comparative genomic hybridization data.

Absolute tumor DNA copy numbers can currently be achieved only on a single gene basis by using fluorescence in situ hybridization (FISH). We present GeneCount, a method for genome-wide calculation of absolute copy numbers from clinical array comparative genomic hybridization data. The tumor cell fraction is reliably estimated in the model. Data consistent with FISH results are achieved. We demonstrate significant improvements over existing methods for exploring gene dosages and intratumor copy number heterogeneity in cancers.