TGF-ß neutralization attenuates tumor residency of activated T cells to enhance systemic immunity in mice.

A tissue resident-like phenotype in tumor infiltrating T cells can limit systemic anti-tumor immunity. Enhanced systemic anti-tumor immunity is observed in head and neck cancer patients after neoadjuvant PD-L1 immune checkpoint blockade (ICB) and transforming growth factor ß (TGF-ß) neutralization. Using T cell receptor (TCR) sequencing and functional immunity assays in a syngeneic model of oral cancer, we dissect the relative contribution of these treatments to enhanced systemic immunity. The addition of TGF-ß neutralization to ICB resulted in the egress of expanded and exhausted CD8+ tumor infiltrating lymphocytes (TILs) into circulation and greater systemic anti-tumor immunity. This enhanced egress associated with reduced expression of Itgae (CD103) and its upstream regulator Znf683. Circulating CD8+ T cells expressed higher Cxcr3 after treatment, an observation also made in samples from patients treated with dual TGF-ß neutralization and ICB. These findings provide the scientific rationale for the use of PD-L1 ICB and TGF-ß neutralization in newly diagnosed patients with carcinomas prior to definitive treatment of locoregional disease.

The tumor microenvironment state associates with response to HPV therapeutic vaccination in patients with respiratory papillomatosis.

Recurrent respiratory papillomatosis (RRP) is a rare, debilitating neoplastic disorder caused by chronic infection with human papillomavirus (HPV) type 6 or 11 and characterized by growth of papillomas in the upper aerodigestive tract. There is no approved medical therapy, and patients require repeated debulking procedures to maintain voice and airway function. PRGN-2012 is a gorilla adenovirus immune-therapeutic capable of enhancing HPV 6/11-specific T cell immunity. This first-in-human, phase 1 study (NCT04724980) of adjuvant PRGN-2012 treatment in adult patients with severe, aggressive RRP demonstrates the overall safety and clinically meaningful benefit observed with PRGN-2012, with a 50% complete response rate in patients treated at the highest dose. Responders demonstrate greater expansion of peripheral HPV-specific T cells compared with nonresponders. Additional correlative studies identify an association between reduced baseline papilloma HPV gene expression, greater interferon responses and expression of CXCL9 and CXCL10, and greater papilloma T cell infiltration in responders. Conversely, nonresponders were characterized by greater HPV and CXCL8 gene expression, increased neutrophilic cell infiltration, and reduced T cell papilloma infiltration. These results suggest that papilloma HPV gene expression may regulate interferon signaling and chemokine expression profiles within the tumor microenvironment that cooperate to govern clinical response to therapeutic HPV vaccination in patients with respiratory papillomatosis.

Oncogenic Ras and ΔNp63α cooperate to recruit immunosuppressive polymorphonuclear myeloid-derived suppressor cells in a mouse model of squamous cancer pathogenesis.

Introduction: Amplification of human chromosome 3q26-29, which encodes oncoprotein ΔNp63 among other isoforms of the p63 family, is a feature common to squamous cell carcinomas (SCCs) of multiple tissue origins. Along with overexpression of ΔNp63, activation of the protooncogene, RAS, whether by overexpression or oncogenic mutation, is frequently observed in many cancers. In this study, analysis of transcriptome data from The Cancer Genome Atlas (TCGA) demonstrated that expression of TP63 mRNA, particularly ΔNp63 isoforms, and HRAS are significantly elevated in advanced squamous cell carcinomas of the head and neck (HNSCCs), suggesting pathological significance. However, how co-overexpressed ΔNp63 and HRAS affect the immunosuppressive tumor microenvironment (TME) is incompletely understood. Methods: Here, we established and characterized an immune competent mouse model using primary keratinocytes with retroviral-mediated overexpression of ΔNp63α and constitutively activated HRAS (v-rasHa G12R) to evaluate the role of these oncogenes in the immune TME. Results: In this model, orthotopic grafting of wildtype syngeneic keratinocytes expressing both v-rasHa and elevated levels of ΔNp63α consistently yield carcinomas in syngeneic hosts, while cells expressing v-rasHa alone yield predominantly papillomas. We found that polymorphonuclear (PMN) myeloid cells, experimentally validated to be immunosuppressive and thus representing myeloid-derived suppressor cells (PMN-MDSCs), were significantly recruited into the TME of carcinomas arising early following orthotopic grafting of ΔNp63α/v-rasHa-expressing keratinocytes. ΔNp63α/v-rasHa-driven carcinomas expressed higher levels of chemokines implicated in recruitment of MDSCs compared to v-rasHa-initiated tumors, providing a heretofore undescribed link between ΔNp63α/HRAS-driven carcinomas and the development of an immunosuppressive TME. Conclusion: These results support the utilization of a genetic carcinogenesis model harboring specific genomic drivers of malignancy to study mechanisms underlying the development of local immunosuppression.

SMYD3 represses tumor-intrinsic interferon response in HPV-negative squamous cell carcinoma of the head and neck.

Cancers often display immune escape, but the mechanisms are incompletely understood. Herein, we identify SMYD3 as a mediator of immune escape in human papilloma virus (HPV)-negative head and neck squamous cell carcinoma (HNSCC), an aggressive disease with poor response to immunotherapy with pembrolizumab. SMYD3 depletion induces upregulation of multiple type I interferon (IFN) response and antigen presentation machinery genes in HNSCC cells. Mechanistically, SMYD3 binds to and regulates the transcription of UHRF1, encoding for a reader of H3K9me3, which binds to H3K9me3-enriched promoters of key immune-related genes, recruits DNMT1, and silences their expression. SMYD3 further maintains the repression of immune-related genes through intragenic deposition of H4K20me3. In vivo, Smyd3 depletion induces influx of CD8+ T cells and increases sensitivity to anti-programmed death 1 (PD-1) therapy. SMYD3 overexpression is associated with decreased CD8 T cell infiltration and poor response to neoadjuvant pembrolizumab. These data support combining SMYD3 depletion strategies with checkpoint blockade to overcome anti-PD-1 resistance in HPV-negative HNSCC.

Phenotypic plasticity and reduced tissue retention of exhausted tumor-infiltrating T cells following neoadjuvant immunotherapy in head and neck cancer.

Neoadjuvant immunotherapies (NITs) have led to clinical benefits in several cancers. Characterization of the molecular mechanisms underlying responses to NIT may lead to improved treatment strategies. Here we show that exhausted, tumor-infiltrating CD8+ T (Tex) cells display local and systemic responses to concurrent neoadjuvant TGF-ß and PD-L1 blockade. NIT induces a significant and selective increase in circulating Tex cells associated with reduced intratumoral expression of the tissue-retention marker CD103. TGF-ß-driven CD103 expression on CD8+ T cells is reversed following TGF-ß neutralization in vitro, implicating TGF-ß in T cell tissue retention and impaired systemic immunity. Transcriptional changes implicate T cell receptor signaling and glutamine metabolism as important determinants of enhanced or reduced Tex treatment response, respectively. Our analysis illustrates physiological and metabolic changes underlying T cell responses to NIT, highlighting the interplay between immunosuppression, tissue retention, and systemic anti-tumor immunity and suggest antagonism of T cell tissue retention as a promising neoadjuvant treatment strategy.

Tumor cell HLA class I expression and pathologic response following neoadjuvant immunotherapy for newly diagnosed head and neck cancer.

OBJECTIVES: Biomarkers are needed to identify patients likely to respond to neoadjuvant immunotherapy (NIT) prior to receiving definitive treatment. MATERIALS AND METHODS: We hypothesized that expression of tumor cell HLA class I would correlate with pathologic response (PR) following NIT for primary untreated head and neck cancer. Multispectral immunofluorescence of pre- and post-treatment biopsy specimens from a neoadjuvant study of bintrafusp alfa, a dual TGF-ß and PD-L1 inhibitor, was performed. RESULTS: Discordant expression of tumor cell HLA class I and PD-L1 measured by multispectral immunofluorescence was observed with most positive tumor cells expressing HLA class I or PD-L1 but not both. Spatial analysis revealed colocalization between tumor parenchyma T cells and HLA class I positive tumors cells, but no clear colocalization between T cells and PD-L1 positive tumor cells. Greater pre-treatment tumor cell HLA class I expression, but not PD-L1 expression or tumor T cell infiltration, correlated with the development of a PR. Additionally, increased tumor cell HLA class I expression after NIT compared to before NIT correlated with development of a PR, whereas inconsistent changes in PD-L1 and T cell infiltration were observed after treatment in all patients. CONCLUSIONS: These data provide the rationale for the study of tumor cell HLA class I expression in larger prospective studies powered to determine the performance of biomarkers of PR in newly diagnosed HNSCC patients receiving NIT.

Durable response in a patient with recurrent respiratory papillomatosis treated with immune checkpoint blockade.

BACKGROUND: Immune checkpoint blockade can provide clinical benefit for patients with advanced cancer. Here, we report durable disease control over many years following PD-L1 blockade through induction of a viral antigen-specific T cell response in an adult patient with recurrent respiratory papillomatosis. METHODS: Antigen-specific T cell response assays, single cell RNA-sequencing, and RNA-scope was used to study clinical tissues. RESULTS: An HPV6 E2-specific T cell clone restricted to HLA-B*55, present at low frequency in the pre-treatment papilloma, significantly expanded after six doses of PD-L1 blockade and remained present and functional at the site of initial response in the larynx as a tissue resident memory T cell for 4 years. An associated reduction in E2 target gene was observed following treatment. CONCLUSIONS: Although demonstrated in a single exceptional responder, these results highlight that immune checkpoint blockade may induce durable, viral antigen-specific immunity of sufficient magnitude to control disease in patients with nonmalignant disorders.

Enhanced neoepitope-specific immunity following neoadjuvant PD-L1 and TGF-ß blockade in HPV-unrelated head and neck cancer.

BACKGROUNDHead and neck squamous cell carcinoma not associated with HPV (HPV-unrelated HNSCC) is associated with a high rate of recurrence and poor survival.METHODSWe conducted a clinical trial in 14 patients with newly diagnosed HPV-unrelated HNSCC to evaluate the safety and efficacy of neoadjuvant bintrafusp alfa, a bifunctional fusion protein that blocks programmed death ligand 1 (PD-L1) and neutralizes TGF-ß.RESULTSBintrafusp alfa was well tolerated, and no treatment-associated surgical delays or complications occurred. Objective pathologic responses (PRs) were observed, and 12 of the 14 (86%) patients were alive and disease free at 1 year. Alterations in Treg infiltration and spatial distribution relative to proliferating CD8+ T cells indicated a reversal of Treg immunosuppression in the primary tumor. Detection of neoepitope-specific tumor T cell responses, but not virus-specific responses, correlated with the development of a PR. Detection of neoepitope-specific responses and PRs in tumors was not correlated with genomic features or tumor antigenicity but was associated with reduced pretreatment myeloid cell tumor infiltration. These results indicate that dual PD-L1 and TGF-ß blockade can safely enhance tumor antigen-specific immunity and highlight the feasibility of multimechanism neoadjuvant immunotherapy for patients with HPV-unrelated HNSCC.CONCLUSIONOur studies provide insight into the ability of neoadjuvant immunotherapy to induce polyclonal neoadjuvant-specific T cell responses in tumors and suggest that features of the tumor microenvironment, such as myeloid cell infiltration, may be a major determinant of enhanced antitumor immunity following such treatment.TRIAL REGISTRATIONClinicalTrials.gov NCT04247282.FUNDINGThis work was funded by the Center for Cancer Research, the NCI, and the Intramural Research Program of the NIDCD, NIH. Bintrafusp alfa was provided by the health care business of Merck KGaA (Darmstadt, Germany), through a Cooperative Research and Development Agreement with the NCI. Additional funding was provided by ImmunityBio through a Cooperative Research and Development Agreement with the NIDCD.

Cure of syngeneic carcinomas with targeted IL-12 through obligate reprogramming of lymphoid and myeloid immunity.

Therapeutic IL-12 has demonstrated the ability to reduce local immune suppression in preclinical models, but clinical development has been limited by severe inflammation-related adverse events with systemic administration. Here, we show that potent immunologic tumor control of established syngeneic carcinomas can be achieved by i.t. administration of a tumor-targeted IL-12 antibody fusion protein (NHS-rmIL-12) using sufficiently low doses to avoid systemic toxicity. Single-cell transcriptomic analysis and ex vivo functional assays of NHS-rmIL-12-treated tumors revealed reinvigoration and enhanced proliferation of exhausted CD8+ T lymphocytes, induction of Th1 immunity, and a decrease in Treg number and suppressive capacity. Similarly, myeloid cells transitioned toward inflammatory phenotypes and displayed reduced suppressive capacity. Cell type-specific IL-12 receptor-KO BM chimera studies revealed that therapeutic modulation of both lymphoid and myeloid cells is required for maximum treatment effect and tumor cure. Study of single-cell data sets from human head and neck carcinomas revealed IL-12 receptor expression patterns similar to those observed in murine tumors. These results describing the diverse mechanisms underlying tumor-directed IL-12-induced antitumor immunity provide the preclinical rationale for the clinical study of i.t. NHS-IL-12.

Comprehensive multiomic characterization of human papillomavirus-driven recurrent respiratory papillomatosis reveals distinct molecular subtypes.

Recurrent respiratory papillomatosis (RRP) is a debilitating neoplastic disorder of the upper aerodigestive tract caused by chronic infection with low-risk human papillomavirus types 6 or 11. Patients with severe RRP can require hundreds of lifetime surgeries to control their disease and pulmonary papillomatosis can be fatal. Here we report the comprehensive genomic and transcriptomic characterization of respiratory papillomas. We discovered and characterized distinct subtypes with transcriptional resemblance to either a basal or differentiated cell state that associate with disease aggressiveness and differ in key molecular, immune and APOBEC mutagenesis profiles. Through integrated comparison with high-risk HPV-associated head and neck squamous cell carcinoma, our analysis revealed divergent molecular and immune papilloma subtypes that form independent of underlying genomic alterations. Cumulatively our results support the development of dysregulated cellular proliferation and suppressed anti-viral immunity through distinct programs of squamous cell differentiation and associated expression of low-risk HPV genes. These analyses provide insight into the pathogenesis of respiratory papillomas and provide a foundation for the development of therapeutic strategies.

Dual PD-L1 and TGF-b blockade in patients with recurrent respiratory papillomatosis.

BACKGROUND: Recurrent respiratory papillomatosis (RRP) is a human papillomavirus (HPV) driven neoplastic disorder of the upper aerodigestive tract that causes significant morbidity and can lead to fatal airway obstruction. Prior clinical study demonstrated clinical benefit with the programmed death-ligand 1 (PD-L1) monoclonal antibody avelumab. Bintrafusp alpha is a bifunctional inhibitor of PD-L1 and transforming growth factor-beta (TGF-b) that has shown clinical activity in several cancer types. METHODS: We conducted a phase II clinical trial evaluating bintrafusp alpha in adults with RRP. Papilloma samples before and after treatment with bintrafusp alpha were assessed for correlates of response with multiplex immunofluorescence as well as immunological and genomic analyses. Post hoc analyses of papilloma samples before and after treatment with avelumab were assessed for comparison. RESULTS: Dual PD-L1/TGF-b inhibition failed to abrogate papilloma growth in most subjects and increased the frequency of clinically indicated interventions after treatment in four of eight subjects based on each subject's own historical control. TGF-b neutralization consistently decreased pSMAD3 and p21 and increased Ki67 expression within the basal layers of papillomas, indicating that TGF-b restrained proliferation. These alterations were not observed in papillomas treated with PD-L1 blockade alone. Dual PD-L1/TGF-b inhibition did not enhance anti-HPV immunity within papillomas beyond that observed with PD-L1 blockade. Genomic alterations in TGF-b superfamily genes were infrequent in papillomas and normal mucosa but present in a significant fraction of head and neck carcinomas. CONCLUSIONS: Intact TGF-b signaling restrains proliferation within papillomas, and the use of clinical agents that abrogate this pathway should be avoided in patients with RRP. TRIAL REGISTRATION NUMBERS: NCT03707587 and NCT02859454.

Chimeric antigen receptor engineered NK cellular immunotherapy overcomes the selection of T-cell escape variant cancer cells.

BACKGROUND: As heterogeneous tumors develop in the face of intact immunity, tumor cells harboring genomic or expression defects that favor evasion from T-cell detection or elimination are selected. For patients with such tumors, T cell-based immunotherapy alone infrequently results in durable tumor control. METHODS: Here, we developed experimental models to study mechanisms of T-cell escape and demonstrated that resistance to T-cell killing can be overcome by the addition of natural killer (NK) cells engineered to express a chimeric antigen receptor (CAR) targeting programmed death ligand-1 (PD-L1). RESULTS: In engineered models of tumor heterogeneity, PD-L1 CAR-engineered NK cells (PD-L1 t-haNKs) prevented the clonal selection of T cell-resistant tumor cells observed with T-cell treatment alone in multiple models. Treatment of heterogenous cancer cell populations with T cells resulted in interferon gamma (IFN-γ) release and subsequent upregulation of PD-L1 on tumor cells that escaped T-cell killing through defects in antigen processing and presentation, priming escape cell populations for PD-L1 dependent killing by PD-L1 t-haNKs in vitro and in vivo. CONCLUSIONS: These results describe the underlying mechanisms governing synergistic antitumor activity between T cell-based immunotherapy that results in IFN-γ production, upregulation of PD-L1 on T-cell escape cells, and the use of PD-L1 CAR-engineered NK cells to target and eliminate resistant tumor cell populations.

Antigen processing and presentation in cancer immunotherapy.

BACKGROUND: Knowledge about and identification of T cell tumor antigens may inform the development of T cell receptor-engineered adoptive cell transfer or personalized cancer vaccine immunotherapy. Here, we review antigen processing and presentation and discuss limitations in tumor antigen prediction approaches. METHODS: Original articles covering antigen processing and presentation, epitope discovery, and in silico T cell epitope prediction were reviewed. RESULTS: Natural processing and presentation of antigens is a complex process that involves proteasomal proteolysis of parental proteins, transportation of digested peptides into the endoplasmic reticulum, loading of peptides onto major histocompatibility complex (MHC) class I molecules, and shuttling of peptide:MHC complexes to the cell surface. A number of T cell tumor antigens have been experimentally validated in patients with cancer. Assessment of predicted MHC class I binding and total score for these validated T cell antigens demonstrated a wide range of values, with nearly one-third of validated antigens carrying an IC50 of greater than 500 nM. CONCLUSIONS: Antigen processing and presentation is a complex, multistep process. In silico epitope prediction techniques can be a useful tool, but comprehensive experimental testing and validation on a patient-by-patient basis may be required to reliably identify T cell tumor antigens.

The role of protein methyltransferases as potential novel therapeutic targets in squamous cell carcinoma of the head and neck.

Squamous cell carcinoma of the head and neck is a lethal disease with suboptimal survival outcomes and standard therapies with significant comorbidities. Whole exome sequencing data recently revealed an abundance of genetic and expression alterations in a family of enzymes known as protein methyltransferases in a variety of cancer types, including squamous cell carcinoma of the head and neck. These enzymes are mostly known for their chromatin-modifying functions through methylation of various histone substrates, though evidence supports their function also through methylation of non-histone substrates. This review summarizes the current knowledge on the function of protein methyltransferases in squamous cell carcinoma of the head and neck and highlights their promising potential as the next generation of therapeutic targets in this disease.

Feasibility, Safety, and Prediction of Complications for Minimally Invasive Myomectomy in Women With Large and Numerous Myomata.

OBJECTIVE: To assess perioperative outcomes and identify predictors of complications for minimally invasive surgery (MIS) myomectomy in a cohort of women with large and numerous myomata. DESIGN: Case-control study (Canadian Task Force classification II-2). SETTING: Academic tertiary care medical center. PATIENTS: Women undergoing MIS myomectomy performed by 3 high-volume surgeons between April 2011 and December 2014. INTERVENTIONS: Characteristics were compared between women who experienced complications and those who did not. Factors predictive of complications were then identified. MEASUREMENTS AND MAIN RESULTS: A total of 221 women underwent an MIS myomectomy, 47.5% via a laparoscopic approach and 52.5% via robotic surgery. The mean ± SD specimen weight was 408.1 ± 384.9 g, uterine volume was 586.1 ± 534.1 cm3, dominant myoma diameter was 9.6 ± 5.1 cm, and number of myomata removed was 4.5 ± 4.1. The most common complications were hemorrhage >1000 mL (8.6%) and blood transfusion (4.1%). The conversion rate was 1.8%. A dominant myoma diameter of ≥12 cm and a uterine volume of ≥750 cm3 increased the odds of complications (odds ratio [OR], 7.44; 95% confidence interval [CI], 2.03-31.84; p = .004 and OR, 6.15; 95% CI, 1.55-30.02; p = .014 respectively). A receiver operating characteristic curve considering dominant myoma diameter and uterine volume had an area under the curve of 0.81. A combination of dominant myoma diameter of ≥10 cm and uterine volume of 600 cm3 predicted complications with 79% sensitivity and 79% specificity. CONCLUSION: Our cohort had large and numerous myomata with high specimen weights, but complications were comparable to those reported in previous studies of MIS myomectomy with less complex pathology. Hemorrhage and transfusion accounted for the majority of complications, and a combination of dominant myoma diameter and uterine volume was predictive of complications. Both factors can be easily defined before surgery and may be used to guide patient counseling, referrals, and implementation of preventative measures for hemorrhage and transfusion.

Adaptive Chromatin Remodeling Drives Glioblastoma Stem Cell Plasticity and Drug Tolerance.

Glioblastoma, the most common and aggressive malignant brain tumor, is propagated by stem-like cancer cells refractory to existing therapies. Understanding the molecular mechanisms that control glioblastoma stem cell (GSC) proliferation and drug resistance may reveal opportunities for therapeutic interventions. Here we show that GSCs can reversibly transition to a slow-cycling, persistent state in response to targeted kinase inhibitors. In this state, GSCs upregulate primitive developmental programs and are dependent upon Notch signaling. This transition is accompanied by widespread redistribution of repressive histone methylation. Accordingly, persister GSCs upregulate, and are dependent on, the histone demethylases KDM6A/B. Slow-cycling cells with high Notch activity and histone demethylase expression are present in primary glioblastomas before treatment, potentially contributing to relapse. Our findings illustrate how cancer cells may hijack aspects of native developmental programs for deranged proliferation, adaptation, and tolerance. They also suggest strategies for eliminating refractory tumor cells by targeting epigenetic and developmental pathways.

A Multiplexed System for Quantitative Comparisons of Chromatin Landscapes.

Genome-wide profiling of histone modifications can provide systematic insight into the regulatory elements and programs engaged in a given cell type. However, conventional chromatin immunoprecipitation and sequencing (ChIP-seq) does not capture quantitative information on histone modification levels, requires large amounts of starting material, and involves tedious processing of each individual sample. Here, we address these limitations with a technology that leverages DNA barcoding to profile chromatin quantitatively and in multiplexed format. We concurrently map relative levels of multiple histone modifications across multiple samples, each comprising as few as a thousand cells. We demonstrate the technology by monitoring dynamic changes following inhibition of p300, EZH2, or KDM5, by linking altered epigenetic landscapes to chromatin regulator mutations, and by mapping active and repressive marks in purified human hematopoietic stem cells. Hence, this technology enables quantitative studies of chromatin state dynamics across rare cell types, genotypes, environmental conditions, and drug treatments.

Sensitive and highly resolved identification of RNA-protein interaction sites in PAR-CLIP data.

BACKGROUND: PAR-CLIP is a recently developed Next Generation Sequencing-based method enabling transcriptome-wide identification of interaction sites between RNA and RNA-binding proteins. The PAR-CLIP procedure induces specific base transitions that originate from sites of RNA-protein interactions and can therefore guide the identification of binding sites. However, additional sources of transitions, such as cell type-specific SNPs and sequencing errors, challenge the inference of binding sites and suitable statistical approaches are crucial to control false discovery rates. In addition, a highly resolved delineation of binding sites followed by an extensive downstream analysis is necessary for a comprehensive characterization of the protein binding preferences and the subsequent design of validation experiments. RESULTS: We present a statistical and computational framework for PAR-CLIP data analysis. We developed a sensitive transition-centered algorithm specifically designed to resolve protein binding sites at high resolution in PAR-CLIP data. Our method employes a Bayesian network approach to associate posterior log-odds with the observed transitions, providing an overall quantification of the confidence in RNA-protein interaction. We use published PAR-CLIP data to demonstrate the advantages of our approach, which compares favorably with alternative algorithms. Lastly, by integrating RNA-Seq data we compute conservative experimentally-based false discovery rates of our method and demonstrate the high precision of our strategy. CONCLUSIONS: Our method is implemented in the R package wavClusteR 2.0. The package is distributed under the GPL-2 license and is available from BioConductor at http://www.bioconductor.org/packages/devel/bioc/html/wavClusteR.html .

A deterministic analysis of genome integrity during neoplastic growth in Drosophila.

The development of cancer has been associated with the gradual acquisition of genetic alterations leading to a progressive increase in malignancy. In various cancer types this process is enabled and accelerated by genome instability. While genome sequencing-based analysis of tumor genomes becomes increasingly a standard procedure in human cancer research, the potential necessity of genome instability for tumorigenesis in Drosophila melanogaster has, to our knowledge, never been determined at DNA sequence level. Therefore, we induced formation of tumors by depletion of the Drosophila tumor suppressor Polyhomeotic and subjected them to genome sequencing. To achieve a highly resolved delineation of the genome structure we developed the Deterministic Structural Variation Detection (DSVD) algorithm, which identifies structural variations (SVs) with high accuracy and at single base resolution. The employment of long overlapping paired-end reads enables DSVD to perform a deterministic, i.e. fragment size distribution independent, identification of a large size spectrum of SVs. Application of DSVD and other algorithms to our sequencing data reveals substantial genetic variation with respect to the reference genome reflecting temporal separation of the reference and laboratory strains. The majority of SVs, constituted by small insertions/deletions, is potentially caused by erroneous replication or transposition of mobile elements. Nevertheless, the tumor did not depict a loss of genome integrity compared to the control. Altogether, our results demonstrate that genome stability is not affected inevitably during sustained tumor growth in Drosophila implying that tumorigenesis, in this model organism, can occur irrespective of genome instability and the accumulation of specific genetic alterations.