William P. J. Leenders

Persistent infections with high-risk human papillomavirus (hrHPV) can cause cervical squamous intraepithelial lesions (SIL) that may progress to cancer. The cervicovaginal microbiome (CVM) correlates with SIL, but the temporal composition of the CVM after hrHPV infections has not been fully clarified. To determine the association between the CVM composition and infection outcome, we applied high-resolution microbiome profiling using the circular probes-based RNA sequencing technology on a longitudinal cohort of cervical smears obtained from 141 hrHPV DNA-positive women with normal cytology at first visit, of whom 51 were diagnosed by cytology with SIL six months later. Here we show that women with a microbial community characterized by low diversity and high Lactobacillus crispatus abundance exhibit low risk to SIL development at both visits, while women with a microbial community characterized by high diversity and Lactobacillus depletion at first visit have a higher risk of developing SIL. At the level of individual species we observed that an increased abundance for Gardnerella vaginalis and Atopobium vaginae associate with SIL outcomes at both visits. These species together with Dialister micraerophilus showed a moderate discriminatory power for hrHPV infection progression. Our results suggest that the CVM can potentially be used as a biomarker for cervical disease and SIL development after hrHPV infection diagnosis with implications on cervical cancer prevention strategies and treatment of SIL.

Background: Because most cervical cancers are caused by high-risk human papillomaviruses (hrHPVs), cervical cancer prevention programs increasingly employ hrHPV testing as a primary test. The high sensitivity of HPV tests is accompanied by low specificity, resulting in high rates of overdiagnosis and overtreatment. Targeted circular probe-based RNA next generation sequencing (ciRNAseq) allows for the quantitative detection of RNAs of interest with high sequencing depth. Here, we examined the potential of ciRNAseq-testing on cervical scrapes to identify hrHPV-positive women at risk of having or developing high-grade cervical intraepithelial neoplasia (CIN). Methods: We performed ciRNAseq on 610 cervical scrapes from the Dutch cervical cancer screening program to detect gene expression from 15 hrHPV genotypes and from 429 human genes. Differentially expressed hrHPV- and host genes in scrapes from women with outcome “no CIN” or “CIN2+” were identified and a model was built to distinguish these groups. Results: Apart from increasing percentages of hrHPV oncogene expression from “no CIN” to high-grade cytology/histology, we identified genes involved in cell cycle regulation, tyrosine kinase signaling pathways, immune suppression, and DNA repair being expressed at significantly higher levels in scrapes with high-grade cytology and histology. Machine learning using random forest on all the expression data resulted in a model that detected ‘no CIN’ versus CIN2+ in an independent data set with sensitivity and specificity of respectively 85 ± 8% and 72 ± 13%. Conclusions: CiRNAseq on exfoliated cells in cervical scrapes measures hrHPV-(onco)gene expression and host gene expression in one single assay and in the process identifies HPV genotype. By combining these data and applying machine learning protocols, the risk of CIN can be calculated. Because ciRNAseq can be performed in high-throughput, making it cost-effective, it can be a promising screening technology to stratify women at risk of CIN2+. Further increasing specificity by model improvement in larger cohorts is warranted.

The R132H mutation in the metabolic enzyme isocitrate dehydrogenase 1 (IDH1) is the most important prognostic factor for the survival of glioma patients. Subsequent studies led to the discovery of a panel of enzymes mainly involved in glutamate anaplerosis and aerobic glycolysis that change in abundance as a result of the IDH1 mutation. To further study these changes, appropriate glioma models are required that accurately mimic in vivo metabolism. To investigate how metabolism is affected by in vitro cell culture, we here compared surgically obtained snap-frozen glioma tissues with their corresponding primary glioma cell culture models with a previously developed targeted mass spectrometry proteomic assay. We determined the relative abundance of a panel of metabolic enzymes. Results confirmed increased glutamate use and decreased aerobic glycolysis in resected IDH1 R132H glioma tissue samples. However, these metabolic profiles were not reflected in the paired glioma primary cell cultures. We suggest that culture conditions and tumor microenvironment play a crucial role in maintaining the in vivo metabolic situation in cell culture models. For this reason, new models that more closely resemble the in vivo microenvironment, such as three-dimensional cell co-cultures or organotypic multicellular spheroid models, need to be developed and investigated.

Background: The cervicovaginal microbiome (CVM) plays a significant role in women’s cervical health and disease. Microbial alterations at the species level and characteristic community state types (CST) have been associated with acquisition and persistence of high-risk human papillomavirus (hrHPV) infections that may result in progression of cervical lesions to malignancy. Current sequencing methods, especially most commonly used multiplex 16S rRNA gene sequencing, struggle to fully clarify these changes because they generally fail to provide sufficient taxonomic resolution to adequately perform species-level associative studies. To improve CVM species designation, we designed a novel sequencing tool targeting microbes at the species taxonomic rank and examined its potential for profiling the CVM. Results: We introduce an accessible and practical circular probe-based RNA sequencing (CiRNAseq) technology with the potential to profile and quantify the CVM. In vitro and in silico validations demonstrate that CiRNAseq can distinctively detect species in a mock mixed microbial environment, with the output data reflecting its ability to estimate microbes’ abundance. Moreover, compared to 16S rRNA gene sequencing, CiRNAseq provides equivalent results but with improved sequencing sensitivity. Analyses of a cohort of cervical smears from hrHPV-negative women versus hrHPV-positive women with high-grade cervical intraepithelial neoplasia confirmed known differences in CST occurring in the CVM of women with hrHPV-induced lesions. The technique also revealed variations in microbial diversity and abundance in the CVM of hrHPV-positive women when compared to hrHPV-negative women. Conclusions: CiRNAseq is a promising tool for studying the interplay between the CVM and hrHPV in cervical carcinogenesis. This technology could provide a better understanding of cervicovaginal CST and microbial species during health and disease, prompting the discovery of biomarkers, additional to hrHPV, that can help detect high-grade cervical lesions.

To study head and neck squamous cell carcinomas (HNSCC) in vitro, a large variety of HNSCC cell lines have been developed. Here, we characterize a panel of 22 HNSCC cell lines, thereby providing a tool for research into tumor-specific treatment options in HNSCC. Both human papillomavirus (HPV) positive and HPV negative tumor cell lines were collected from commercial and collaborative sources. Short tandem repeat profiling was used to confirm or characterize the identity of the cell lines. Targeted sequencing was performed using a standard pathology single molecule Molecular Inversion Probe panel to detect mutations for 23 tumor suppressors and oncogenes. HPV status, p16 status, radiosensitivity data, and hypoxia data are summarized from all cell lines. We detected HPV transcripts in five cell lines, all of which overexpressed p16. One HPV negative cell line was also p16 positive. We detected mutations in KIT (SCCNij185), PIK3CA (SCCNij185), and CDKN2A (UT-SCC-5 and UT-SCC-38). TP53 mutations were the most frequent, occurring in 16/22 cell lines. HPV infection and TP53 mutations were almost mutually exclusive, with the exception of 93-VU-147T. The cell lines exhibited a wide range of sensitivities towards hypoxia and irradiation. Here, we provide a description of a set of frequently used HNSCC cell lines with diverse characteristics as found in HNSCC patients.

Ovarian cancer is the most lethal gynecological malignancy due to late detection associated with dissemination throughout the abdominal cavity. Targeted photodynamic therapy (tPDT) aimed at epithelial cell adhesion molecule (EpCAM), overexpressed in over 90% of ovarian cancer metastatic lesions, is a promising novel therapeutic modality. Here, we tested the specificity and activity of conjugates of EpCAM-directed designed ankyrin repeat proteins (DARPins) with the photosensitizer IRDye 700DX in in vitro and in vivo ovarian cancer models. EpCAM-binding DARPins (Ec1: K_d = 68 pM; Ac2: K_d = 130 nM) and a control DARPin were site-specifically functionalized with fluorophores or IRDye 700DX. Conjugation of anti-EpCAM DARPins with fluorophores maintained EpCAM-specific binding in cell lines and patient-derived ovarian cancer explants. Penetration of DARPin Ec1 into tumor spheroids was slower than that of Ac2, indicative of a binding site barrier effect for Ec1. DARPin-IRDye 700DX conjugates killed EpCAM-expressing cells in a highly specific and illumination-dependent fashion in 2D and 3D cultures. Furthermore, they effectively homed to EpCAM-expressing subcutaneous OV90 xenografts in mice. In conclusion, the high activity and specificity observed in preclinical ovarian cancer models, combined with a high specificity in patient material, warrant a further investigation of EpCAM-targeted PDT for ovarian cancer.

IDH1^R132H (isocitrate dehydrogenase 1) mutations play a key role in the development of low-grade gliomas. IDH1^wt converts isocitrate to α-ketoglutarate while reducing nicotinamide adenine dinucleotide phosphate (NADP⁺), whereas IDH1^R132H uses α-ketoglutarate and NADPH to generate the oncometabolite 2-hydroxyglutarate (2-HG). While the effects of 2-HG have been the subject of intense research, the 2-HG independent effects of IDH1^R132H are still ambiguous. The present study demonstrates that IDH1^R132H expression but not 2-HG alone leads to significantly decreased tricarboxylic acid (TCA) cycle metabolites, reduced proliferation, and enhanced sensitivity to irradiation in both glioblastoma cells and astrocytes in vitro. Glioblastoma cells, but not astrocytes, showed decreased NADPH and NAD⁺ levels upon IDH1^R132H transduction. However, in astrocytes IDH1^R132H led to elevated expression of the NAD-synthesizing enzyme nicotinamide phosphoribosyltransferase (NAMPT). These effects were not 2-HG mediated. This suggests that IDH1^R132H cells utilize NAD⁺ to restore NADP pools, which only astrocytes could compensate via induction of NAMPT. We found that the expression of NAMPT is lower in patient-derived IDH1-mutant glioma cells and xenografts compared to IDH1-wildtype models. The Cancer Genome Atlas (TCGA) data analysis confirmed lower NAMPT expression in IDH1-mutant versus IDH1-wildtype gliomas. We show that the IDH1 mutation directly affects the energy homeostasis and redox state in a cell-type dependent manner. Targeting the impairments in metabolism and redox state might open up new avenues for treating IDH1-mutant gliomas.

Radiotherapy is an important treatment modality of head and neck squamous cell carcinomas (HNSCC). Multiple links have been described between the metabolic activity of tumors and their clinical outcome. Here we test the hypothesis that metabolic features determine radiosensitivity, explaining the relationship between metabolism and clinical outcome. Radiosensitivity of 14 human HNSCC cell lines was determined using colony forming assays and the expression profile of approximately 200 metabolic and cancer-related genes was generated using targeted RNA sequencing by single molecule molecular inversion probes. Results: Correlation between radiosensitivity data and expression profiles yielded 18 genes associated with radiosensitivity or radioresistance, of which adenosine triphosphate (ATP) citrate lyase (ACLY) was of particular interest. Pharmacological inhibition of ACLY caused an impairment of DNA damage repair, specifically homologous recombination, and lead to radiosensitization in HNSCC cell lines. Examination of a The Cancer Genome Atlas (TCGA) cohort of HNSCC patients revealed that high expression of ACLY was predictive for radiotherapy failure, as it was only associated with poor overall survival in patients who received radiotherapy (hazard ratio of 2.00, 95% CI: 1.12–3.55; p = 0.0184). These data were further validated in an independent cohort of HNSCC patients treated with chemoradiation. Furthermore, patients with poor locoregional control after radiotherapy have significantly higher nuclear ACLY protein levels. Together, we here show that ACLY affects DNA damage repair, and is a predictive factor for radiotherapy outcome in HNSCC.

Many biology-based precision drugs are available that neutralize aberrant molecular pathways in cancer. Molecular heterogeneity and the lack of reliable companion diagnostic biomarkers for many drugs makes targeted treatment of cancer inaccurate for many individuals. Identifying actionable hyperactive biological pathways in individual cancers may improve this situation.To achieve this we applied a novel targeted RNA next generation sequencing (t/RNA-NGS) technique to surgically obtained glioma tissues. The test combines mutation detection with analysis of biological pathway activities that are involved in tumour behavior in many cancer types (e.g. tyrosine kinase signaling, angiogenesis signaling, immune response, metabolism), via quantitative measurement of transcript levels and splice variants of hundreds of genes. We here present proof of concept that the technique, which uses molecular inversion probes, generates a histology-independent molecular diagnosis and identifies classifiers that are strongly associated with conventional histopathology diagnoses and even with patient prognosis. The test not only confirmed known glioma-associated molecular aberrations but also identified aberrant expression levels of actionable genes and mutations that have so far been considered not to be associated with glioma, opening up the possibility of drug repurposing for individual patients. Its cost-effectiveness makes t/RNA-NGS to an attractive instrument to aid oncologists in therapy decision making.

Glioblastoma (GB) is a fatal disease in which most targeted therapies have clinically failed. However, pharmacological reactivation of tumor suppressors has not been thoroughly studied as yet as a GB therapeutic strategy. Tumor suppressor Protein Phosphatase 2A (PP2A), is inhibited by non-genetic mechanisms in GB, and thus it would be potentially amendable for therapeutic reactivation. Here we demonstrate, that s mall m olecule a ctivators of P P2A (SMAPs), NZ-8-061 and DBK-1154, effectively cross the in vitro model of blood-brain barrier (BBB), and in vivo partition to mouse brain tissue after oral dosing. In vitro , SMAPs exhibit robust cell killing activity against five established GB cell lines, and nine patient-derived primary glioma cell lines. Collectively these cell lines have heterogenous genetic background, kinase inhibitor resistance profile, and stemness properties; and they represent different clinical GB subtypes. Oral dosing of either of the SMAPs significantly reduced growth of infiltrative intracranial GB tumors. DBK-1154, with both higher degree of brain/blood distribution, and more potent in vitro activity against all tested GB cell lines, also significantly increased survival of mice bearing orthotopic GB xenografts. In summary, this report presents a proof-of-principle data for BBB-permeable tumor suppressor reactivation therapy for glioblastoma cells of heterogenous molecular background.