Search

According to the current guidelines, watchful waiting (WW) is a feasible option for patients with good or intermediate prognosis renal-cell carcinoma (RCC). However, some patients rapidly progress during WW, requiring the initiation of treatment. Here, we explore whether we can identify those patients using circulating cell-free DNA (cfDNA) methylation. We first defined a panel of RCC-specific circulating methylation markers by intersecting differentially methylated regions from a publicly available dataset with known RCC methylation markers from the literature. The resulting RCC-specific methylation marker panel of 22 markers was subsequently evaluated for an association with rapid progression by methylated DNA sequencing (MeD-seq) in serum from 10 HBDs and 34 RCC patients with a good or intermediate prognosis starting WW in the IMPACT-RCC study. Patients with an elevated RCC-specific methylation score compared to HBDs had a shorter progression-free survival (PFS, p = 0.018), but not a shorter WW-time (p = 0.15). Cox proportional hazards regression showed that only the International Metastatic Renal Cell Carcinoma Database Consortium (IMDC) criteria were significantly associated with WW time (HR 2.01, p = 0.01), whereas only our RCC-specific methylation score (HR 4.45, p = 0.02) was significantly associated with PFS. The results of this study suggest that cfDNA methylation is predictive of PFS but not WW.

Cellular senescence is a state of stable growth arrest and a desired outcome of tumor suppressive interventions. Treatment with many anti-cancer drugs can cause premature senescence of non-malignant cells. These therapy-induced senescent cells can have pro-tumorigenic and pro-disease functions via activation of an inflammatory secretory phenotype (SASP). Inhibitors of cyclin-dependent kinases 4/6 (CDK4/6i) have recently proven to restrain tumor growth by activating a senescence-like program in cancer cells. However, the physiological consequence of exposing the whole organism to pharmacological CDK4/6i remains poorly characterized. Here, we show that exposure to CDK4/6i induces non-malignant cells to enter a premature state of senescence dependent on p53. We observe in mice and breast cancer patients that the CDK4/6i-induced senescent program activates only a partial SASP enriched in p53 targets but lacking pro-inflammatory and NF-κB-driven components. We find that CDK4/6i-induced senescent cells do not acquire pro-tumorigenic and detrimental properties but retain the ability to promote paracrine senescence and undergo clearance. Our results demonstrate that SASP composition is exquisitely stress-dependent and a predictor for the biological functions of different senescence subsets.

Purpose: Detection of leptomeningeal metastasis is hampered by limited sensitivities of currently used techniques: MRI and cytology of cerebrospinal fluid (CSF). Detection of cell-free tumor DNA in CSF has been proposed as a tumor-specific candidate to detect leptomeningeal metastasis at an earlier stage. The aim of this study was to investigate mutation and aneuploidy status in CSF-derived cell-free DNA (cfDNA) of patients with breast cancer with a clinical suspicion of leptomeningeal metastasis. Experimental Design: cfDNA was isolated from stored remnant CSF and analyzed by targeted next-generation sequencing (NGS; n = 30) and the modified fast aneuploidy screening test-sequencing system (mFAST-SeqS; n = 121). The latter method employs selective amplification of long interspaced nuclear elements sequences that are present throughout the genome and allow for fast and cheap detection of aneuploidy. We compared these results with the gold standard to diagnose leptomeningeal metastasis: cytology. Results: Leptomeningeal metastasis was cytology proven in 13 of 121 patients. Low DNA yields resulted in insufficient molecular coverage of NGS for the majority of samples (success rate, 8/30). The mFAST-SeqS method, successful in 112 of 121 (93%) samples, detected genome-wide aneuploidy in 24 patients. Ten of these patients had cytology-proven leptomeningeal metastasis; 8 additional patients were either concurrently diagnosed with central nervous system metastases by radiological means or developed these soon after the lumbar puncture. The remaining six cases were suspected of leptomeningeal metastasis, but could not be confirmed by cytology or imaging. Aneuploidy was associated with development of leptomeningeal metastasis and significantly worse overall survival. Conclusions: Aneuploidy in CSF-derived cfDNA may provide a promising biomarker to improve timely detection of leptomeningeal metastasis.

Publisher's version (útgefin grein), Circular RNAs (circRNAs) are a class of RNAs that is under increasing scrutiny, although their functional roles are debated. We analyzed RNA-seq data of 348 primary breast cancers and developed a method to identify circRNAs that does not rely on unmapped reads or known splice junctions. We identified 95,843 circRNAs, of which 20,441 were found recurrently. Of the circRNAs that match exon boundaries of the same gene, 668 showed a poor or even negative (R < 0.2) correlation with the expression level of the linear gene. In silico analysis showed only a minority (8.5%) of circRNAs could be explained by known splicing events. Both these observations suggest that specific regulatory processes for circRNAs exist. We confirmed the presence of circRNAs of CNOT2, CREBBP, and RERE in an independent pool of primary breast cancers. We identified circRNA profiles associated with subgroups of breast cancers and with biological and clinical features, such as amount of tumor lymphocytic infiltrate and proliferation index. siRNA-mediated knockdown of circCNOT2 was shown to significantly reduce viability of the breast cancer cell lines MCF-7 and BT-474, further underlining the biological relevance of circRNAs. Furthermore, we found that circular, and not linear, CNOT2 levels are predictive for progression-free survival time to aromatase inhibitor (AI) therapy in advanced breast cancer patients, and found that circCNOT2 is detectable in cell-free RNA from plasma. We showed that circRNAs are abundantly present, show characteristics of being specifically regulated, are associated with clinical and biological properties, and thus are relevant in breast cancer., We thank the Erasmus MC Cancer Computational Biology Center for giving access to their IT infrastructure and the software that was used for the computations and data analysis in this study. We thank Sandra Albassam for her help with the first versions of the script to identify circular regions. We thank Maurice P.H.M. Jansen, Jean C. Helmijr, Inge de Kruijff, and Manouk K. Bos for their help in evaluating plasma samples that were gathered in the EU-FP7 CareMore (nr 601760) project. We thank for technical support Miriam Ragle Aure and Anita Langerød of the Oslo University Hospital, Norway; Ewan Birney of the European Bioinformatics Institute, UK; and Stefania Tommasi of the IRCCS Istituto Tumori “Giovanni Paolo II,” Bari, Italy. We thank the Oslo Breast Cancer Research Consortium (OSBREAC), Norway (https://www.ous-research.no/home/kgjebsen/home/14105) for contributing patient samples and Sabine Linn and Marleen Kok of The Netherlands Cancer Institute for contributing samples for the AI cohort. Finally, we thank all members of the ICGC Breast Cancer Working Group. This work has been funded through the ICGC Breast Cancer Working group by the Breast Cancer Somatic Genetics Study (a European research project funded by the European Community’s Seventh Framework Programme (FP7/2010-2014) under the grant agreement number 242006) and the Triple Negative project funded by the Wellcome Trust (grant reference 077012/Z/05/Z). F.G.R.-G. and S.M. were funded by BASIS. J.A.F. was funded through an ERC Advanced Grant (ERC-2012-AdG-322737) and ERC Proof-of-Concept Grant (ERC-2017-PoC-767854). K.U. was funded by the Daniel den Hoed Foundation. S.N.-Z. is a Wellcome Beit Fellow and personally funded by a Wellcome Trust Intermediate Fellowship (WT100183MA). A.L.R. is partially supported by the Dana-Farber/Harvard Cancer Center SPORE in Breast Cancer (NIH/NCI 5 P50 CA16 8504-02). A.M.S. was supported by Cancer Genomics Netherlands (CGC.nl) through a grant from the Netherlands Organization of Scientific research (NWO). M. Smid was supported by the EU-FP7-DDR response project. C.D. was supported by a grant from the Breast Cancer Research Foundation. J.E. was funded by The Icelandic Centre for Research (RANNIS).

Simple Summary Currently, the most commonly used method to monitor response to treatment in metastatic breast cancer patients is by radiological imaging. However, these imaging techniques are relatively insensitive and give little to no insight into biological tumor characteristics that might be relevant for the choice of treatment. Circulating tumor DNA (ctDNA), released by tumor cells into the blood of cancer patients, can be used to overcome these shortcomings. Besides the fact that specific alterations are known to predict response to treatment and development of resistance, the total amount of ctDNA is believed to reflect the proliferation rate of the tumor, suggesting ctDNA levels can be used as a general tool to evaluate treatment response. Different methods are available to measure ctDNA primarily based on detection of cancer-specific somatic mutations, DNA methylation, and copy number variations. In this review we have critically analyzed recently published studies using blood-derived ctDNA of metastatic breast cancer patients on multiple time points to monitor disease response in respect to analytical validity and clinical utility. Abstract Monitoring treatment response in metastatic breast cancer currently consists mainly of radiological and clinical assessments. These methods have high inter-observer variation, suboptimal sensitivity to determine response to treatment and give little insight into the biological characteristics of the tumor. Assessing circulating tumor DNA (ctDNA) over time could be employed to address these limitations. Several ways to quantify and characterize ctDNA exist, based on somatic mutations, copy number variations, methylation, and global circulating cell-free DNA (cfDNA) fragment sizes and concentrations. These methods are being explored and technically validated, but to date none of these methods are applied clinically. We systematically reviewed the literature on the use of quantitative ctDNA measurements over time to monitor response to systemic therapy in patients with metastatic breast cancer. Cochrane, Embase, PubMed and Google Scholar databases were searched to find studies focusing on the use of cfDNA to longitudinally monitor treatment response in advanced breast cancer patients until October 2020. This resulted in a total of 33 studies which met the inclusion criteria. These studies were heterogeneous in (pre-)processing procedures, applied techniques and design. An association between ctDNA and treatment response was found in most of the included studies, independent of the applied assay. To implement ctDNA-based response monitoring into daily clinical practice for metastatic breast cancer patients, sample (pre-) processing procedures need to be standardized and large prospectively collected sample cohorts with well annotated clinical follow-up are required to establish its clinical validity.

In contrast to primary colorectal cancer (CRC) little is known about the genomic landscape of metastasized CRC. Here we present whole genome sequencing data of metastases of 429 CRC patients participating in the pan-cancer CPCT-02 study (NCT01855477). Unsupervised clustering using mutational signature patterns highlights three major patient groups characterized by signatures known from primary CRC, signatures associated with received prior treatments, and metastasis-specific signatures. Compared to primary CRC, we identify additional putative (non-coding) driver genes and increased frequencies in driver gene mutations. In addition, we identify specific genes preferentially affected by microsatellite instability. CRC-specific 1kb-10Mb deletions, enriched for common fragile sites, and LINC00672 mutations are associated with response to treatment in general, whereas FBXW7 mutations predict poor response specifically to EGFR-targeted treatment. In conclusion, the genomic landscape of mCRC shows defined changes compared to primary CRC, is affected by prior treatments and contains features with potential clinical relevance., Molecular landscapes of metastatic colorectal cancers (mCRC) have often been restricted to coding regions or low numbers of patients. Here the authors present a whole-genome landscape of 429 mCRC patients, revealing the mutational impact of prior therapies and potential actionable targets.

Background There are profound individual differences in clinical outcomes between colorectal cancers (CRCs) presenting with identical stage of disease. Molecular stratification, in conjunction with the traditional TNM staging, is a promising way to predict patient outcomes. We investigated the interconnectivity between tumor stage and tumor biology reflected by the Consensus Molecular Subtypes (CMSs) in CRC, and explored the possible value of these insights in patients with stage II colon cancer. Methods We performed a retrospective analysis using clinical records and gene expression profiling in a meta-cohort of 1040 CRC patients. The interconnectivity of tumor biology and disease stage was assessed by investigating the association between CMSs and TNM classification. In order to validate the clinical applicability of our findings we employed a meta-cohort of 197 stage II colon cancers. Results CMS4 was significantly more prevalent in advanced stages of disease (stage I 9.8% versus stage IV 38.5%, p p p = 0.008). Conclusions Considerable interconnectivity between tumor biology and tumor stage in CRC exists. This implies that the TNM stage, in addition to the stage of progression, might also reflect distinct biological disease entities. These insights can potentially be utilized to optimize identification of high-risk stage II colon cancers.

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.

miRNAs represent an emerging class of promising biomarkers for cancer diagnostics. To perform reliable miRNA expression analysis using quantitative PCR, adequate data normalization is essential to remove nonbiological, technical variations. Ideal reference genes should be biologically stable and reduce technical variability of miRNA expression analysis. Herein is a new strategy for the identification and evaluation of reference genes that can be applied for miRNA-based diagnostic tests without entailing excessive additional experiments. We analyzed the expression of 11 carefully selected candidate reference genes in different types of cervical specimens [ie, tissues, scrapes, and self-collected cervicovaginal specimens (self-samples)]. To identify the biologically most stable reference genes, three commonly used algorithms (GeNorm, NormFinder, and BestKeeper) were combined. Signal-to-noise ratios and P values between control and disease groups were calculated to validate the reduction in technical variability on expression analysis of two marker miRNAs. miR-423 was identified as a suitable reference gene for all sample types, to be used in combination with RNU24 in cervical tissues, RNU43 in scrapes, and miR-30b in self-samples. These findings demonstrate that the choice of reference genes may differ between different types of specimens, even when originating from the same anatomical source. More important, it is shown that adequate normalization increases the signal-to-noise ratio, which is not observed when normalizing to commonly used reference genes.

Purpose: Epigenetic host cell changes involved in cervical cancer development following a persistent high-risk human papillomavirus (hrHPV) infection, provide promising markers for the management of hrHPV-positive women. In particular, markers based on DNA methylation of tumor suppressor gene promoters are valuable. These markers ideally identify hrHPV-positive women with precancer (CIN2/3) in need of treatment. Here, we set out to identify biologically relevant methylation markers by genome-wide methylation analysis of both hrHPV-transformed cell lines and cervical tissue specimens. Experimental Design and Results: Genome-wide discovery by next-generation sequencing (NGS) of methyl-binding domain–enriched DNA (MBD-Seq) yielded 20 candidate methylation target genes. Further verification and validation by multiplex-targeted bisulfite NGS and (quantitative) methylation-specific PCR (MSP) resulted in 3 genes (GHSR, SST, and ZIC1) that showed a significant increase in methylation with severity of disease in both tissue specimens and cervical scrapes (P < 0.005). The area under the ROC curve for CIN3 or worse varied between 0.86 and 0.89. Within the group of CIN2/3, methylation levels of all 3 genes increased with duration of lesion existence (P < 0.0005), characterized by duration of preceding hrHPV infection, and were significantly higher in the presence of a 3q gain (P < 0.05) in the corresponding tissue biopsy. Conclusions: By unbiased genome-wide DNA methylation profiling and comprehensive stepwise verification and validation studies using in vitro and patient-derived samples, we identified 3 promising methylation markers (GHSR, SST, and ZIC1) associated with a 3q gain for the detection of cervical (pre)cancer. Clin Cancer Res; 23(14); 3813–22. ©2017 AACR.

Cervical screening by high-risk HPV (hrHPV) testing requires additional risk stratification (triage), as most infections are transient and only a subset of hrHPV-positive women harbours clinically relevant disease. Molecular triage markers such as microRNAs (miRNAs) and DNA methylation markers are particularly promising, as they can be objectively tested directly on hrHPV-positive scrapes and cervicovaginal self-samples. Here, we evaluated the marker potential of 10 candidate miRNAs in 209 hrHPV-positive scrapes of women with underlying precancer (cervical intraepithelial neoplasia, grade 2–3 (CIN2-3)), cancer, or without disease (CIN0/1). A predictive miRNA classifier for CIN3 detection was built using logistic regression, which was compared to and combined with DNA methylation marker FAM19A4. Markers were correlated to histology parameters and hrHPV genotype. A miRNA classifier consisting of miR-149, miR-20a, and miR-93 achieved an area under the curve (AUC) of 0.834 for CIN3 detection, which was not significantly different to that of FAM19A4 methylation (AUC: 0.862, p = 0.591). Combining miRNA and methylation analysis demonstrated complementarity between both marker types (AUC: 0.939). While the miRNA classifier seemed more predictive for CIN2, FAM19A4 methylation was particularly high in HPV16-positive and histologically advanced CIN3, i.e. CIN3 with high lesion volume. The miRNA classifier, FAM19A4 methylation, and the miRNA/methylation combination were highest in cancer-associated scrapes. In conclusion, a panel of three miRNAs is discriminatory for CIN3 in hrHPV-positive scrapes and can complement DNA methylation analysis for the efficient detection of cervical disease. Combined analysis of the two marker types warrants further evaluation as triage strategy in hrHPV-based screening.

Summary Our aim is to improve omics based prediction and feature selection using multiple sources of auxiliary information: co-data. Adaptive group regularized ridge regression (GRridge) was proposed to achieve this by estimating additional group-based penalty parameters through an empirical Bayes method at a low computational cost. We illustrate the GRridge method and software on RNA sequencing datasets. The method boosts the performance of an ordinary ridge regression and outperforms other classifiers. Post-hoc feature selection maintains the predictive ability of the classifier with far fewer markers. Availability and Implementation GRridge is an R package that includes a vignette. It is freely available at (https://bioconductor.org/packages/GRridge/). All information and R scripts used in this study, including those on retrieval and processing of the co-data, are available from http://github.com/markvdwiel/GRridgeCodata. Supplementary information Supplementary data are available at Bioinformatics online.

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.

Cervical cancer develops from precancerous high-grade cervical intraepithelial neoplasia (CIN) harboring a transforming infection with high-risk human papillomavirus, which is characterized by p16INK4a overexpression. Once such a lesion has developed, progression toward an invasive squamous cell carcinoma (SCC) may take one or more decades, underlining the heterogeneity of these lesions in terms of duration of existence and progression risk. We performed array-based comparative genomic hybridization (array CGH) on 46 p16INK4a immunopositive CIN2/3 lesions to determine whether this heterogeneity is reflected in their chromosomal profiles. Chromosomal profiles of CIN2/3 lesions were related to those of invasive cervical SCC and promoter methylation of CADM1, a tumor suppressor gene known to be functionally involved in the tumorigenic phenotype of cervical cancer cells. Frequent alterations found in CIN2/3 lesions included gains located at chromosome 1, 3, 7, and 20 and losses located at 4, 11, 16, 17, and 19. Unsupervised hierarchical clustering identified two subsets of CIN2/3 lesions, chromosomal profiles of one of which closely resembled invasive SCCs. Gains of 1, 3q, and 20 were characteristic for CIN2/3 lesions with chromosomal signatures resembling carcinomas. In addition, dense promoter methylation of the CADM1 gene was significantly more frequent in these CIN2/3 lesions (P = 0.004). No chromosomal alterations were detected in six CIN1 lesions, five of which were completely p16INK4a immunonegative. These findings suggest that biomarkers associated with gains at chromosomes 1, 3q, and 20 are potential hallmarks of advanced p16INK4a-positive CIN2/3 lesions with a high short-term risk of progression. [Cancer Res 2009;69(2):647–55]

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.

Infection with high-risk human papillomavirus (HR-HPV) has been associated with intraepithelial neoplasia and carcinomas at various sites of the anogenital tract, including the cervix, vulva, vagina, penis and anus. Although HR-HPV is a necessary cause for cervical cancer, the majority of anal cancers and a subset of cancers at other genital sites, additional (epi)genetic events are required for malignant transformation. HPV-mediated transformation of human epithelial cells has been recognized as a multi-step process resulting from deregulated transcription of the viral oncogenes E6 and E7 in the proliferating cells. Interference of E6 and E7 with cell cycle regulators induces genetic instability, which drives the continuous selection of oncogenic alterations providing cells with a malignant phenotype. Early genetic events during cervical carcinogenesis associated with immortalization, include deletions at chromosomes 3p, 6 and 10p, whereas amongst others gain of chromosome 3q, loss of chromosome 11 and epigenetic alterations such as inactivation of the TSLC1 tumor suppressor gene represent later events associated with tumor invasion.

Little is known about the alterations in microRNA (miRNA) expression patterns during the consecutive stages of cervical cancer development and their association with chromosomal instability. In this study, miRNA expression in normal cervical squamous epithelium, high-grade precancerous lesions (cervical intraepithelial neoplasia (CIN2-3)), squamous cell carcinomas (SCCs) and adenocarcinomas (AdCAs) was integrated with previously generated chromosomal profiles of the same samples. Significantly differential expression during the consecutive stages of cervical SCC development was observed for 106 miRNAs. Of these differentially expressed miRNAs, 27 showed early transiently altered expression in CIN2-3 lesions only, 46 miRNAs showed late altered expression in SCCs only and 33 showed continuously altered expression in both CIN2-3 and SCCs. Altered expression of five significantly differentially expressed miRNAs, hsa-miR-9 (1q23.2), hsa-miR-15b (3q25.32), hsa-miR-28-5p (3q27.3), hsa-miR-100 and hsa-miR-125b (both 11q24.1), was directly linked to frequent chromosomal alterations. Functional analyses were performed for hsa-miR-9, representing a potential oncogene with increased expression linked to a chromosomal gain of 1q. Hsa-miR-9 overexpression was found to increase cell viability, anchorage-independent growth and migration in vitro. Upon organic raft culturing, hsa-miR-9 hampered differentiation and induced proliferation in all strata of the epithelial layer. These findings support a potential oncogenic function of hsa-miR-9 in cervical cancer. In summary, differential expression of 106 miRNAs, partly associated with chromosomal alterations, was observed during cervical SCC development. Altered expression of hsa-miR-9 associated with a chromosomal gain of chromosome 1q was shown to be functionally relevant, underlining the importance of deregulated miRNA expression in cervical carcinogenesis. © 2013 Macmillan Publishers Limited.

High-grade cervical intraepithelial neoplasia (CIN2/3) represents a heterogeneous disease both with respect to clinical behavior and chromosomal aberrations detected. We hypothesized that the extent of chromosomal aberrations reflects the duration of their existence. Chromosomal profiles were determined of CIN3 of women with a known 5-year history of high-risk human papillomavirus virus (hrHPV) infection, in which duration of prior hrHPV infection was considered a proxy for duration of CIN3 existence. Eleven women had a

Autofluorescence bronchoscopy (AFB) is a valid strategy for detecting premalignant endobronchial lesions. However, no biomarker can reliably predict lung cancer risk of subjects with AFB-visualized premalignant lesions.The present study set out to identify AFB-visualized squamous metaplastic (SqM) lesions with malignant potential by DNA copy number profiling.Regular AFB examinations in 474 subjects at risk of lung cancer identified six subjects with SqM lesions at baseline, and carcinoma in situ or carcinoma (carcinoma in situ or greater) at the initial SqM site at follow-up bronchoscopy. These progressive SqM lesions were compared for immunostaining pattern and array comparative genomic hybridization-based chromosomal profiles with 23 SqM lesions of subjects who remained cancer-free. Specific DNA copy number alterations (CNAs) linked to cancer risk were identified and accuracy of CNAs to predict endobronchial cancer in this series was determined.At baseline, p53, p63, and Ki-67 immunostaining were not predictive for a differential clinical outcome of SqM lesions. The mean number of CNAs in baseline SqM of cases was significantly higher compared with control subjects (P0.01). Chromosomal regions significantly more frequently altered in SqM of cases were 3p26.3-p11.1, 3q26.2-q29, 9p13.3-p13.2, and 17p13.3-p11.2 (family-wise error rate0.10). CNAs were specifically detected at the site of future cancer. In cases, baseline-detected CNAs persisted in subsequent biopsies taken from the initial site, and levels increased toward cancer progression. In this series, a model based on CNAs at 3p26.3-p11.1, 3q26.2-29, and 6p25.3-24.3 predicted cancer with 97% accuracy.The data suggest that the presence of specific CNAs in SqM lesions predict endobronchial cancer.

Background A substantial number of microRNAs (miRNAs) is subject to epigenetic silencing in cancer. Although epigenetic silencing of tumour suppressor genes is an important feature of cervical cancer, little is known about epigenetic silencing of miRNAs. Since DNA methylation-based silencing of hsa-miR-124 occurs in various human cancers, we studied the frequency and functional effects of hsa-miR-124 methylation in cervical carcinogenesis. Results Quantitative MSP analysis of all 3 loci encoding the mature hsa-miR-124 (hsa-miR-124-1/-2/-3) showed methylation in cervical cancer cell lines SiHa, CaSki and HeLa as well as in late passages of human papillomavirus (HPV) type 16 or 18 immortalised keratinocytes. Treatment of SiHa cells with a demethylating agent reduced hsa-miR-124 methylation levels and induced hsa-miR-124 expression. In HPV-immortalised keratinocytes increased methylation levels were related to reduced hsa-miR-124 expression and higher mRNA expression of IGFBP7, a potential hsa-miR-124 target gene. Ectopic hsa-miR-124 expression in SiHa and CaSki cells decreased proliferation rates and migratory capacity. Combined hsa-miR-124-1 and/or hsa-miR-124-2 methylation analysis of 139 cervical tissue specimens showed an increasing methylation frequency from 0% in normal tissues up to 93% in cervical carcinomas. Increased methylation levels of hsa-miR-124-1 and hsa-miR-124-2 were significantly correlated with reduced hsa-miR-124 expression in cervical tissue specimens. Combined hsa-miR-124-1 and/or hsa-miR-124-2 methylation analysis of 43 cervical scrapes of high-risk HPV positive women was predictive of underlying high-grade lesions. Conclusions DNA methylation-based silencing of hsa-miR-124 is functionally involved in cervical carcinogenesis and may provide a valuable marker for improved detection of cervical cancer and its high-grade precursor lesions.

Background It is well known that a persistent infection with high-risk human papillomavirus (hrHPV) is causally involved in the development of squamous cell carcinomas of the uterine cervix (CxSCCs) and a subset of SCCs of the head and neck (HNSCCs). The latter differ from hrHPV-negative HNSCCs at the clinical and molecular level. Methods To determine whether hrHPV-associated SCCs arising from different organs have specific chromosomal alterations in common, we compared genome-wide chromosomal profiles of 10 CxSCCs (all hrHPV-positive) with 12 hrHPV-positive HNSCCs and 30 hrHPV-negative HNSCCs. Potential organ-specific alterations and alterations shared by SCCs in general were investigated as well. Results Unsupervised hierarchical clustering resulted in one mainly hrHPV-positive and one mainly hrHPV-negative cluster. Interestingly, loss at 13q and gain at 20q were frequent in HPV-positive carcinomas of both origins, but uncommon in hrHPV-negative HNSCCs, indicating that these alterations are associated with hrHPV-mediated carcinogenesis. Within the group of hrHPV-positive carcinomas, HNSCCs more frequently showed gains of multiple regions at 8q whereas CxSCCs more often showed loss at 17p. Finally, gains at 3q24-29 and losses at 11q22.3-25 were frequent (>50%) in all sample groups. Conclusion In this study hrHPV-specific, organ-specific, and pan-SCC chromosomal alterations were identified. The existence of hrHPV-specific alterations in SCCs of different anatomical origin, suggests that these alterations are crucial for hrHPV-mediated carcinogenesis.

We aimed to link DNA methylation events occurring in cervical carcinomas to distinct stages of HPV-induced transformation. Methylation specific-multiplex ligation-dependent probe amplification (MS-MLPA) analysis of cervical carcinomas revealed promoter methylation of 12 out of 29 tumour suppressor genes analysed, with MGMT being most frequently methylated (92%). Subsequently, consecutive stages of HPV16/18-transfected keratinocytes (n=11), ranging from pre-immortal to anchorage-independent phenotypes, were analysed by MS-MLPA. Whereas no methylation was evident in pre-immortal cells, progression to anchorage independence was associated with an accumulation of frequent methylation events involving five genes, all of which were also methylated in cervical carcinomas. TP73 and ESR1 methylation became manifest in early immortal cells followed by RARβ and DAPK1 methylation in late immortal passages. Complementary methylation of MGMT was related to anchorage independence. Analysis of nine cervical cancer cell lines, representing the tumorigenic phenotype, revealed in addition to these five genes frequent methylation of CADM1, CDH13 and CHFR. In conclusion, eight recurrent methylation events in cervical carcinomas could be assigned to different stages of HPV-induced transformation. Hence, our in vitro model system provides a valuable tool to further functionally address the epigenetic alterations that are common in cervical carcinomas.

Summary: CGHcall achieves high calling accuracy for array CGH data by effective use of breakpoint information from segmentation and by inclusion of several biological concepts that are ignored by existing algorithms. The algorithm is validated for simulated and verified real array CGH data. By incorporating more than three classes, CGHcall improves detection of single copy gains and amplifications. Moreover, it allows effective inclusion of chromosome arm information. Availability: An R-package (GUI), a manual and an example data set are available at http://www.few.vu.nl/~mavdwiel/CGHcall.html. Contact: mark.vdwiel@vumc.nl Supplementary information: Supplementary data are available at Bioinformatics online.

Purpose: Adenoid cystic carcinoma (ACC) is a malignant salivary gland tumor with a high mortality rate due to late, distant metastases. This study aimed at unraveling common genetic abnormalities associated with ACC. Additionally, chromosomal changes were correlated with patient characteristics and survival. Experimental Design: Microarray-based comparative genomic hybridization was done to a series of 18 paraffin-embedded primary ACCs using a genome-wide scanning BAC array. Results: A total of 238 aberrations were detected, representing more gains than losses (205 versus 33, respectively). Most frequent gains (>60%) were observed at 9q33.3-q34.3, 11q13.3, 11q23.3, 19p13.3-p13.11, 19q12-q13.43, 21q22.3, and 22q13.33. These loci harbor numerous growth factor [fibroblast growth factor (FGF) and platelet-derived growth factor (PDGF)] and growth factors receptor (FGFR3 and PDGFRβ) genes. Gains at the FGF(R) regions occurred significantly more frequently in the recurred/metastasized ACCs compared with indolent ACCs. Furthermore, patients with 17 or more chromosomal aberrations had a significantly less favorable outcome than patients with fewer chromosomal aberrations (log-rank = 5.2; P = 0.02). Conclusions: Frequent DNA copy number gains at loci of growth factors and their receptors suggest their involvement in ACC initiation and progression. Additionally, the presence of FGFR3 and PDGFRβ in increased chromosomal regions suggests a possible role for autocrine stimulation in ACC tumorigenesis.

Background Prediction in high dimensional settings is difficult due to the large number of variables relative to the sample size. We demonstrate how auxiliary ‘co-data’ can be used to improve the performance of a Random Forest in such a setting. Results Co-data are incorporated in the Random Forest by replacing the uniform sampling probabilities that are used to draw candidate variables by co-data moderated sampling probabilities. Co-data here are defined as any type information that is available on the variables of the primary data, but does not use its response labels. These moderated sampling probabilities are, inspired by empirical Bayes, learned from the data at hand. We demonstrate the co-data moderated Random Forest (CoRF) with two examples. In the first example we aim to predict the presence of a lymph node metastasis with gene expression data. We demonstrate how a set of external p-values, a gene signature, and the correlation between gene expression and DNA copy number can improve the predictive performance. In the second example we demonstrate how the prediction of cervical (pre-)cancer with methylation data can be improved by including the location of the probe relative to the known CpG islands, the number of CpG sites targeted by a probe, and a set of p-values from a related study. Conclusion The proposed method is able to utilize auxiliary co-data to improve the performance of a Random Forest. Electronic supplementary material The online version of this article (doi:10.1186/s12859-017-1993-1) contains supplementary material, which is available to authorized users.

Background To determine which changes in the host cell genome are crucial for cervical carcinogenesis, a longitudinal in vitro model system of HPV-transformed keratinocytes was profiled in a genome-wide manner. Four cell lines affected with either HPV16 or HPV18 were assayed at 8 sequential time points for gene expression (mRNA) and gene copy number (DNA) using high-resolution microarrays. Available methods for temporal differential expression analysis are not designed for integrative genomic studies. Results Here, we present a method that allows for the identification of differential gene expression associated with DNA copy number changes over time. The temporal variation in gene expression is described by a generalized linear mixed model employing low-rank thin-plate splines. Model parameters are estimated with an empirical Bayes procedure, which exploits integrated nested Laplace approximation for fast computation. Iteratively, posteriors of hyperparameters and model parameters are estimated. The empirical Bayes procedure shrinks multiple dispersion-related parameters. Shrinkage leads to more stable estimates of the model parameters, better control of false positives and improvement of reproducibility. In addition, to make estimates of the DNA copy number more stable, model parameters are also estimated in a multivariate way using triplets of features, imposing a spatial prior for the copy number effect. Conclusion With the proposed method for analysis of time-course multilevel molecular data, more profound insight may be gained through the identification of temporal differential expression induced by DNA copy number abnormalities. In particular, in the analysis of an integrative oncogenomics study with a time-course set-up our method finds genes previously reported to be involved in cervical carcinogenesis. Furthermore, the proposed method yields improvements in sensitivity, specificity and reproducibility compared to existing methods. Finally, the proposed method is able to handle count (RNAseq) data from time course experiments as is shown on a real data set. Electronic supplementary material The online version of this article (doi:10.1186/1471-2105-15-327) contains supplementary material, which is available to authorized users.