Virus transcript levels and cell growth rates after naturally occurring HPV16 integration events in basal cervical keratinocytes.

Cervical carcinogenesis is characterized by a clonal selection process in which the high-risk human papillomavirus (HRHPV) genome usually changes from the extra-chromosomal (episomal) state seen in productive infections to DNA that is integrated into host chromosomes. However, it is not clear whether all HRHPV integration events provide cells with a selective growth advantage compared with the episome-containing cells from which they originate. It is also unclear whether selection of cells containing a particular integrant from a mixed population simply reflects the highest levels of virus oncogene expression or has additional determinants. These early events in cervical carcinogenesis cannot readily be addressed by cross-sectional studies of clinical samples. We used the W12 model system to generate a panel of cervical squamous cell clones that were derived from an identical background under non-competitive conditions and differed only by the genomic site of HPV16 integration. Compared with the 'baseline' episome-containing cells from which they were isolated, only 9/17 clones (53%) showed significantly greater growth rates and only 7/17 (41%) showed significantly greater expression of the major virus oncogenes E7/E6. There were significant variations in levels of HPV16 transcription per DNA template, changes that were associated with histone modifications in the integrated virus chromatin. Cell growth rates showed only weak and non-significant associations with protein and mRNA levels for E7, E6, and the mean E7/E6 values. We conclude that HPV16 integration in basal cervical cells does not necessarily lead to increased levels of virus oncogenes, or to a competitive growth advantage, when compared with the initiating episome-containing cells.

Depletion of HPV16 early genes induces autophagy and senescence in a cervical carcinogenesis model, regardless of viral physical state.

In cervical carcinomas, high-risk human papillomavirus (HR-HPV) may be integrated into host chromosomes or remain extra-chromosomal (episomal). We used the W12 cervical keratinocyte model to investigate the effects of HPV16 early gene depletion on in vitro cervical carcinogenesis pathways, particularly effects shared by cells with episomal versus integrated HPV16 DNA. Importantly, we were able to study the specific cellular consequences of viral gene depletion by using short interfering RNAs known not to cause phenotypic or transcriptional off-target effects in keratinocytes. We found that while cervical neoplastic progression in vitro was characterized by dynamic changes in HPV16 transcript levels, viral early gene expression was required for cell survival at all stages of carcinogenesis, regardless of viral physical state, levels of early gene expression or histology in organotypic tissue culture. Moreover, HPV16 early gene depletion induced changes in host gene expression that were common to both episome-containing and integrant-containing cells. In particular, we observed up-regulation of autophagy genes, associated with enrichment of senescence and innate immune-response pathways, including the senescence-associated secretory phenotype (SASP). In keeping with these observations, HPV16 early gene depletion induced autophagy in both episome-containing and integrant-containing W12 cells, as evidenced by the appearance of autophagosomes, punctate expression of the autophagy marker LC3, conversion of LC3B-I to LC3B-II, and reduced levels of the autophagy substrate p62. Consistent with the reported association between autophagy and senescence pathways, HPV16 early gene depletion induced expression of the senescence marker beta-galactosidase and increased secretion of the SASP-related protein IGFBP3. Together, these data indicate that depleting HR-HPV early genes would be of potential therapeutic benefit in all cervical carcinogenesis pathways, regardless of viral physical state. In addition, the senescence/SASP response associated with autophagy induction may promote beneficial immune effects in bystander cells.

LIN28 Expression in malignant germ cell tumors downregulates let-7 and increases oncogene levels.

Despite their clinicopathologic heterogeneity, malignant germ cell tumors (GCT) share molecular abnormalities that are likely to be functionally important. In this study, we investigated the potential significance of downregulation of the let-7 family of tumor suppressor microRNAs in malignant GCTs. Microarray results from pediatric and adult samples (n = 45) showed that LIN28, the negative regulator of let-7 biogenesis, was abundant in malignant GCTs, regardless of patient age, tumor site, or histologic subtype. Indeed, a strong negative correlation existed between LIN28 and let-7 levels in specimens with matched datasets. Low let-7 levels were biologically significant, as the sequence complementary to the 2 to 7 nt common let-7 seed "GAGGUA" was enriched in the 3' untranslated regions of mRNAs upregulated in pediatric and adult malignant GCTs, compared with normal gonads (a mixture of germ cells and somatic cells). We identified 27 mRNA targets of let-7 that were upregulated in malignant GCT cells, confirming significant negative correlations with let-7 levels. Among 16 mRNAs examined in a largely independent set of specimens by quantitative reverse transcription PCR, we defined negative-associations with let-7e levels for six oncogenes, including MYCN, AURKB, CCNF, RRM2, MKI67, and C12orf5 (when including normal control tissues). Importantly, LIN28 depletion in malignant GCT cells restored let-7 levels and repressed all of these oncogenic let-7 mRNA targets, with LIN28 levels correlating with cell proliferation and MYCN levels. Conversely, ectopic expression of let-7e was sufficient to reduce proliferation and downregulate MYCN, AURKB, and LIN28, the latter via a double-negative feedback loop. We conclude that the LIN28/let-7 pathway has a critical pathobiologic role in malignant GCTs and therefore offers a promising target for therapeutic intervention.

Depletion of polycistronic transcripts using short interfering RNAs: cDNA synthesis method affects levels of non-targeted genes determined by quantitative PCR.

BACKGROUND: Short interfering RNAs (siRNAs) are often used to deplete viral polycistronic transcripts, such as those encoded by human papillomavirus (HPV). There are conflicting data in the literature concerning how siRNAs targeting one HPV gene can affect levels of other genes in the polycistronic transcripts. We hypothesised that the conflict might be partly explained by the method of cDNA synthesis used prior to transcript quantification. FINDINGS: We treated HPV16-positive cervical keratinocytes with siRNAs targeting the HPV16 E7 gene and used quantitative PCR to compare transcript levels of E7 with those of E6 and E2, viral genes located upstream and downstream of the target site respectively. We compared our findings from cDNA generated using oligo-dT primers alone with those from cDNA generated using a combination of random hexamer and oligo-dT primers. Our data show that when polycistronic transcripts are targeted by siRNAs, there is a period when untranslatable cleaved mRNA upstream of the siRNA binding site remains detectable by PCR, if cDNA is generated using random hexamer primers. Such false indications of mRNA abundance are avoided using oligo-dT primers. The period corresponds to the time taken for siRNA activity and degradation of the cleaved transcripts. Genes downstream of the siRNA binding site are detectable during this interval, regardless of how the cDNA is generated. CONCLUSIONS: These data emphasise the importance of the cDNA synthesis method used when measuring transcript abundance following siRNA depletion of polycistronic transcripts. They provide a partial explanation for erroneous reports suggesting that siRNAs targeting HPV E7 can have gene-specific effects.

Overexpression of the oncostatin M receptor in cervical squamous cell carcinoma cells is associated with a pro-angiogenic phenotype and increased cell motility and invasiveness.

Oncostatin M receptor (OSMR) shows frequent copy number gain and overexpression in advanced cervical squamous cell carcinoma (SCC). We used cell-based in vitro assays, RNA interference, and integrative gene expression profiling to investigate the functional significance of this observation. CaSki and SW756 were selected as representative cervical SCC cells that overexpressed OSMR, and ME180 and MS751 as cells that did not. The STAT-dependent pro-angiogenic factors VEGF-A and ID1 were rapidly induced by OSM in CaSki/SW756 but not in ME180/MS751. However, rapid induction did occur in MS751 following forced OSMR overexpression, while depleting OSMR in CaSki abrogated VEGF-A expression. Conditioned medium from both CaSki and SW756 stimulated endothelial tube formation in vitro, effects that were inhibited by depleting OSMR in the SCC cells. For both CaSki and SW756, migration in a wound healing assay and invasion through Matrigel were stimulated by OSM and consistently inhibited by OSMR depletion. The phenotype was rescued by transfection with OSMR containing a silent mutation that provided specific siRNA resistance. Overall, there was a positive correlation between OSMR levels and invasiveness. We used gene expression profiling to identify genes induced by OSM in CaSki/SW756 but not in ME180/MS751. The most prominent gene ontology category groups for the differentially expressed genes were cell motility/invasion, angiogenesis, signal transduction, and apoptosis. We also profiled 23 cervical SCC samples, identifying genes that were differentially expressed in cases with OSMR overexpression versus those without. Integration of the datasets identified 15 genes that showed consistent differential expression in association with OSMR levels in vitro and in vivo. We conclude that OSMR overexpression in cervical SCC cells provides increased sensitivity to OSM, which induces pro-malignant changes. OSMR is a potential prognostic and therapeutic target in cervical SCC. The genes that mediate OSM:OSMR effects will be valuable indicators of the effectiveness of antibody blockade in pre-clinical systems.

The two most common histological subtypes of malignant germ cell tumour are distinguished by global microRNA profiles, associated with differential transcription factor expression.

BACKGROUND: We hypothesised that differences in microRNA expression profiles contribute to the contrasting natural history and clinical outcome of the two most common types of malignant germ cell tumour (GCT), yolk sac tumours (YSTs) and germinomas. RESULTS: By direct comparison, using microarray data for paediatric GCT samples and published qRT-PCR data for adult samples, we identified microRNAs significantly up-regulated in YSTs (n = 29 paediatric, 26 adult, 11 overlapping) or germinomas (n = 37 paediatric). By Taqman qRT-PCR we confirmed differential expression of 15 of 16 selected microRNAs and further validated six of these (miR-302b, miR-375, miR-200b, miR-200c, miR-122, miR-205) in an independent sample set. Interestingly, the miR-302 cluster, which is over-expressed in all malignant GCTs, showed further over-expression in YSTs versus germinomas, representing six of the top eight microRNAs over-expressed in paediatric YSTs and seven of the top 11 in adult YSTs. To explain this observation, we used mRNA expression profiles of paediatric and adult malignant GCTs to identify 10 transcription factors (TFs) consistently over-expressed in YSTs versus germinomas, followed by linear regression to confirm associations between TF and miR-302 cluster expression levels. Using the sequence motif analysis environment iMotifs, we identified predicted binding sites for four of the 10 TFs (GATA6, GATA3, TCF7L2 and MAF) in the miR-302 cluster promoter region. Finally, we showed that miR-302 family over-expression in YST is likely to be functionally significant, as mRNAs down-regulated in YSTs were enriched for 3' untranslated region sequences complementary to the common seed of miR-302a~miR-302d. Such mRNAs included mediators of key cancer-associated processes, including tumour suppressor genes, apoptosis regulators and TFs. CONCLUSIONS: Differential microRNA expression is likely to contribute to the relatively aggressive behaviour of YSTs and may enable future improvements in clinical diagnosis and/or treatment.

In vitro progression of human papillomavirus 16 episome-associated cervical neoplasia displays fundamental similarities to integrant-associated carcinogenesis.

An important event in the development of cervical squamous cell carcinoma (SCC) is deregulated expression of high-risk human papillomavirus (HR-HPV) oncogenes, most commonly related to viral integration into host DNA. Mechanisms of development of the approximately 15% of SCCs that contain extrachromosomal (episomal) HR-HPV are poorly understood due to limited longitudinal data. We therefore used the W12 model to study mechanisms of cervical carcinogenesis associated with episomal HPV16. In vitro progression of W12 normally occurs through selection of cells containing integrated HPV16. However, in one long-term culture, keratinocytes developed a selective growth advantage and invasive phenotype while retaining HPV16 episomes at increased copy number in the absence of transcriptionally active integrants. Longitudinal investigations revealed similarities between the episome- and integrant-associated routes of neoplastic progression. Most notable were dynamic changes in viral early gene expression in episome-retaining cells, consistent with continually changing selective pressures. An early increase in viral transcription preceded elevated episome copy number and was followed by a reduction to near baseline after the development of invasiveness. Episomal transcriptional deregulation did not require selection of a specific sequence variant of the HPV16 upstream regulatory region, although increased levels of acetylated histone H4 around the late promoter implicated a role for altered chromatin structure. Interestingly, invasive episome-retaining cells showed high levels of HPV16 E2/E6 proteins (despite decreased transcript levels) and reduced expression of IFN-stimulated genes, adaptations that support viral persistence and cell survival. Our findings suggest a unified working model for events important in cervical neoplastic progression regardless of HR-HPV physical state.

Malignant germ cell tumors display common microRNA profiles resulting in global changes in expression of messenger RNA targets.

Despite their extensive clinical and pathologic heterogeneity, all malignant germ cell tumors (GCT) are thought to originate from primordial germ cells. However, no common biological abnormalities have been identified to date. We profiled 615 microRNAs (miRNA) in pediatric malignant GCTs, controls, and GCT cell lines (48 samples in total) and re-analyzed available miRNA expression data in adult gonadal malignant GCTs. We applied the bioinformatic algorithm Sylamer to identify miRNAs that are of biological importance by inducing global shifts in mRNA levels. The most significant differentially expressed miRNAs in malignant GCTs were all from the miR-371-373 and miR-302 clusters (adjusted P < 0.00005), which were overexpressed regardless of histologic subtype [yolk sac tumor (YST)/seminoma/embryonal carcinoma (EC)], site (gonadal/extragonadal), or patient age (pediatric/adult). Sylamer revealed that the hexamer GCACTT, complementary to the 2- to 7-nucleotide miRNA seed AAGUGC shared by six members of the miR-371-373 and miR-302 clusters, was the only sequence significantly enriched in the 3'-untranslated region of mRNAs downregulated in pediatric malignant GCTs (as a group), YSTs and ECs, and in adult YSTs (all versus nonmalignant tissue controls; P < 0.05). For the pediatric samples, downregulated genes containing the 3'-untranslated region GCACTT showed significant overrepresentation of Gene Ontology terms related to cancer-associated processes, whereas for downregulated genes lacking GCACTT, Gene Ontology terms generally represented metabolic processes only, with few genes per term (adjusted P < 0.05). We conclude that the miR-371-373 and miR-302 clusters are universally overexpressed in malignant GCTs and coordinately downregulate mRNAs involved in biologically significant pathways.

Characterization of naturally occurring HPV16 integration sites isolated from cervical keratinocytes under noncompetitive conditions.

As the high-risk human papillomavirus (HPV) integrants seen in anogenital carcinomas represent the end-point of a clonal selection process, we used the W12 model to study the naturally occurring integration events that exist in HPV16-infected cervical keratinocytes before integrant selection. We performed limiting dilution cloning to identify integrants present in cells that also maintain episomes. Such integrants arise in a natural context and exist in a noncompetitive environment, as they are transcriptionally repressed by episome-derived E2. We found that integration can occur at any time during episome maintenance, providing biological support for epidemiologic observations that persistent HPV infection is a major risk factor in cervical carcinogenesis. Of 24 different integration sites isolated from a single nonclonal population of W12, 12 (50%) occurred within chromosome bands containing a common fragile site (CFS), similar to observations for selected integrants in vivo. This suggests that such regions represent relatively accessible sites for insertion of foreign DNA, rather than conferring a selective advantage when disrupted. Interestingly, however, integrants and CFSs did not accurately colocalize. We further observed that local DNA rearrangements occur frequently and rapidly after the integration event. The majority of integrants were in chromosome bands containing a cancer-associated coding gene or microRNA, indicating that integration occurs commonly in these regions, regardless of selective pressure. The cancer-associated genes were generally a considerable distance from the integration site, and there was no evidence for altered expression of nine strong candidate genes. These latter observations do not support an important role for HPV16 integration in causing insertional mutagenesis.

Critical evaluation of HPV16 gene copy number quantification by SYBR green PCR.

BACKGROUND: Human papilloma virus (HPV) load and physical status are considered useful parameters for clinical evaluation of cervical squamous cell neoplasia. However, the errors implicit in HPV gene quantification by PCR are not well documented. We have undertaken the first rigorous evaluation of the errors that can be expected when using SYBR green qPCR for quantification of HPV type 16 gene copy numbers. We assessed a modified method, in which external calibration curves were generated from a single construct containing HPV16 E2, HPV16 E6 and the host gene hydroxymethylbilane synthase in a 1:1:1 ratio. RESULTS: When testing dilutions of mixed HPV/host DNA in replicate runs, we observed errors in quantifying E2 and E6 amplicons of 5-40%, with greatest error at the lowest DNA template concentration (3 ng/microl). Errors in determining viral copy numbers per diploid genome were 13-53%. Nevertheless, in cervical keratinocyte cell lines we observed reasonable agreement between viral loads determined by qPCR and Southern blotting. The mean E2/E6 ratio in episome-only cells was 1.04, but with a range of 0.76-1.32. In three integrant-only lines the mean E2/E6 ratios were 0.20, 0.72 and 2.61 (values confirmed by gene-specific Southern blotting). When E2/E6 ratios in fourteen HPV16-positive cervical carcinomas were analysed, conclusions regarding viral physical state could only be made in three cases, where the E2/E6 ratio was < or = 0.06. CONCLUSION: Run-to-run variation in SYBR green qPCR produces unavoidable inaccuracies that should be allowed for when quantifying HPV gene copy number. While E6 copy numbers can be considered to provide a useable indication of viral loads, the E2/E6 ratio is of limited value. Previous studies may have overestimated the frequency of mixed episomal/integrant HPV infections.

Pediatric malignant germ cell tumors show characteristic transcriptome profiles.

Malignant germ cell tumors (GCT) of childhood are rare and heterogeneous neoplasms thought to arise from primordial germ cells. They vary substantially in their natural history and show important clinical differences from their adult counterparts. To address the biological basis for these observations, we have undertaken a comprehensive analysis of global gene expression patterns in pediatric malignant GCTs and compared these findings with published data on adult testicular GCTs (TGCT). Our study included 27 primary tumors and assessed the principal malignant histologic types of pediatric GCT, yolk sac tumor (YST; n = 18), and seminoma (n = 9). Analysis of Affymetrix U133A GeneChip data was performed using the statistical software environment R, including gene set enrichment analysis, with cross-validation at the RNA and protein level. Unsupervised analysis showed complete separation of YSTs and seminomas by global gene expression profiles and identified a robust set of 657 discriminatory transcripts. There was no segregation of tumors of the same histology arising at different sites or at different ages within the pediatric range. In contrast, there was segregation of pediatric malignant GCTs and adult malignant TGCTs, most notably for the YSTs. The pediatric seminomas were significantly enriched for genes associated with the self-renewing pluripotent phenotype, whereas the pediatric YSTs were significantly enriched for genes associated with a differentiation and proliferation phenotype. We conclude that histologic type is the key discriminator in pediatric malignant GCTs and that the observed clinical differences between malignant GCTs of children and adults are mirrored by significant differences in global gene expression.

Interferon-beta treatment of cervical keratinocytes naturally infected with human papillomavirus 16 episomes promotes rapid reduction in episome numbers and emergence of latent integrants.

Following integration of human papillomavirus (HPV) into the host genome, overexpression of the viral oncogenes E6 and E7 requires loss of the transcriptional repressor functions of E2. A key step in HPV-related carcinogenesis is therefore clearance of residual viral episomes, which encode E2. As spontaneous loss of HPV-16 episomes in vitro is associated with increased expression of antiviral genes inducible by type I interferon (IFN), we used the W12 model to examine the effects of exogenous IFN-beta on cervical keratinocytes containing HPV-16 episomes as a result of 'natural' infection in vivo. In contrast to studies of cells transfected with HPV-31 or bovine papillomavirus, IFN-beta caused rapid reduction in numbers of HPV-16 episomes. This was associated with the emergence of cells bearing previously latent integrants, in which there was increased expression of E6 and E7. Our data indicate that integrated HPV-16 can exist in a minority of cells in a mixed population without exerting a selective advantage until episome numbers are reduced. The kinetics of cell death and changes in viral transcription and translation that we observed support a model where integrants are initially present in cells also containing episomes, with generalized episome clearance by IFN-beta resulting in integrant de-repression. We conclude that IFN-beta can hasten the transition from episomal to integrated HPV-16 in naturally infected cervical keratinocytes. Greater emphasis should be placed on episome loss in models of HPV-related carcinogenesis. We provide the strongest evidence to date that treating HPV-16 lesions by inducing an IFN response may cause clinical progression.

Selection of cervical keratinocytes containing integrated HPV16 associates with episome loss and an endogenous antiviral response.

Integration of high-risk human papillomavirus (HRHPV) into the host genome is a key event in cervical neoplastic progression. Integration is associated with deregulated expression of the viral oncogenes E6 and E7 and acquisition of a selective growth advantage for cells containing integrants. Overexpression of the viral transcriptional regulator E2 from heterologous promoters has an inhibitory effect on transcription from integrated HRHPV. Therefore, we hypothesized that loss of E2-expressing episomes from cells in which integration had previously occurred would be required for such cells to gain a growth advantage. Using the unique W12 model of cervical squamous carcinogenesis, we show that cells containing integrated HPV16 reproducibly emerged during long-term culture when there had been a rapid fall in episome numbers. During the period of emergence, it is possible to isolate single-cell clones containing an intracellular mixture of the integrant being selected and episomes at reduced load. The lower level of E2 expression seen in such cells is associated with partial inhibition of transcription from the HPV16 integrant. Full deregulation is not observed until complete loss of E2-expressing episomes occurs. Microarray analysis showed that episome loss was closely associated with endogenous activation of antiviral response genes that are also inducible by the type I IFN pathway. Taken together, our results indicate that episome loss, associated with induction of antiviral response genes, is a key event in the spontaneous selection of cervical keratinocytes containing integrated HPV16. We conclude that cervical carcinogenesis requires not only HRHPV integration, but also loss of inhibitory episomes.

Acquisition of high-level chromosomal instability is associated with integration of human papillomavirus type 16 in cervical keratinocytes.

Whereas two key steps in cervical carcinogenesis are integration of high-risk human papillomavirus (HR-HPV) and acquisition of an unstable host genome, the temporal association between these events is poorly understood. Chromosomal instability is induced when HR-HPV E7 oncoprotein is overexpressed from heterologous promoters in vitro. However, it is not known whether such events occur at the "physiologically" elevated levels of E7 produced by deregulation of the homologous HR-HPV promoter after integration. Indeed, an alternative possibility is that integration in vivo is favored in an already unstable host genome. We have addressed these issues using the unique human papillomavirus (HPV) 16-containing cervical keratinocyte cell line W12, which was derived from a low-grade squamous intraepithelial lesion and thus acquired HPV16 by "natural" infection. Whereas W12 at low passage contains HPV16 episomes only, long-term culture results in the emergence of cells containing integrated HPV16 only. We show that integration of HPV16 in W12 is associated with 3' deletion of the E2 transcriptional repressor, resulting in deregulation of the homologous promoter of the integrant and an increase in E7 protein levels. We further demonstrate that high-level chromosomal instability develops in W12 only after integration and that the forms of instability observed correlate with the physical state of HPV16 DNA and the level of E7 protein. Whereas intermediate E7 levels are associated with numerical chromosomal abnormalities, maximal levels are associated with both numerical and structural aberrations. HR-HPV integration is likely to be a critical event in cervical carcinogenesis, preceding the development of chromosomal abnormalities that drive malignant progression.