A MYCN-driven de-differentiation profile identifies a subgroup of aggressive retinoblastoma.

Retinoblastoma are childhood eye tumors arising from retinal precursor cells. Two distinct retinoblastoma subtypes with different clinical behavior have been described based on gene expression and methylation profiling. Using consensus clustering of DNA methylation analysis from 61 retinoblastomas, we identify a MYCN-driven cluster of subtype 2 retinoblastomas characterized by DNA hypomethylation and high expression of genes involved in protein synthesis. Subtype 2 retinoblastomas outside the MYCN-driven cluster are characterized by high expression of genes from mesodermal development, including NKX2-5. Knockdown of MYCN expression in retinoblastoma cell models causes growth arrest and reactivates a subtype 1-specific photoreceptor signature. These molecular changes suggest that removing the driving force of MYCN oncogenic activity rescues molecular circuitry driving subtype 1 biology. The MYCN-RB gene signature generated from the cell models better identifies MYCN-driven retinoblastoma than MYCN amplification and can identify cases that may benefit from MYCN-targeted therapy. MYCN drives tumor progression in a molecularly defined retinoblastoma subgroup, and inhibiting MYCN activity could restore a more differentiated and less aggressive tumor biology.

Pancreatic cancer acquires resistance to MAPK pathway inhibition by clonal expansion and adaptive DNA hypermethylation.

BACKGROUND: Pancreatic ductal adenocarcinoma (PDAC) is an aggressive cancer with poor prognosis. It is marked by extraordinary resistance to conventional therapies including chemotherapy and radiation, as well as to essentially all targeted therapies evaluated so far. More than 90% of PDAC cases harbor an activating KRAS mutation. As the most common KRAS variants in PDAC remain undruggable so far, it seemed promising to inhibit a downstream target in the MAPK pathway such as MEK1/2, but up to now preclinical and clinical evaluation of MEK inhibitors (MEKi) failed due to inherent and acquired resistance mechanisms. To gain insights into molecular changes during the formation of resistance to oncogenic MAPK pathway inhibition, we utilized short-term passaged primary tumor cells from ten PDACs of genetically engineered mice. We followed gain and loss of resistance upon MEKi exposure and withdrawal by longitudinal integrative analysis of whole genome sequencing, whole genome bisulfite sequencing, RNA-sequencing and mass spectrometry data. RESULTS: We found that resistant cell populations under increasing MEKi treatment evolved by the expansion of a single clone but were not a direct consequence of known resistance-conferring mutations. Rather, resistant cells showed adaptive DNA hypermethylation of 209 and hypomethylation of 8 genomic sites, most of which overlap with regulatory elements known to be active in murine PDAC cells. Both DNA methylation changes and MEKi resistance were transient and reversible upon drug withdrawal. Furthermore, MEKi resistance could be reversed by DNA methyltransferase inhibition with remarkable sensitivity exclusively in the resistant cells. CONCLUSION: Overall, the concept of acquired therapy resistance as a result of the expansion of a single cell clone with epigenetic plasticity sheds light on genetic, epigenetic and phenotypic patterns during evolvement of treatment resistance in a tumor with high adaptive capabilities and provides potential for reversion through epigenetic targeting.

Non-responsiveness to cardioprotection by ischaemic preconditioning in Ossabaw minipigs with genetic predisposition to, but without the phenotype of the metabolic syndrome.

The translation of successful preclinical and clinical proof-of-concept studies on cardioprotection to the benefit of patients with reperfused acute myocardial infarction has been difficult so far. This difficulty has been attributed to confounders which patients with myocardial infarction typically have but experimental animals usually not have. The metabolic syndrome is a typical confounder. We hypothesised that there may also be a genuine non-responsiveness to cardioprotection and used Ossabaw minipigs which have the genetic predisposition to develop a diet-induced metabolic syndrome, but before they had developed the diseased phenotype. Using a prospective study design, a reperfused acute myocardial infarction was induced in 62 lean Ossabaw minipigs by 60 min coronary occlusion and 180 min reperfusion. Ischaemic preconditioning by 3 cycles of 5 min coronary occlusion and 10 min reperfusion was used as cardioprotective intervention. Ossabaw minipigs were stratified for their single nucleotide polymorphism as homozygous for valine (V/V) or isoleucine (I/I)) in the γ-subunit of adenosine monophosphate-activated protein kinase. Endpoints were infarct size and area of no-reflow. Infarct size (V/V: 54 ± 8, I/I: 54 ± 13% of area at risk, respectively) was not reduced by ischaemic preconditioning (V/V: 55 ± 11, I/I: 46 ± 11%) nor was the area of no-reflow (V/V: 57 ± 18, I/I: 49 ± 21 vs. V/V: 57 ± 21, I/I: 47 ± 21% of infarct size). Bioinformatic comparison of the Ossabaw genome to that of Sus scrofa and Göttingen minipigs identified differences in clusters of genes encoding mitochondrial and inflammatory proteins, including the janus kinase (JAK)-signal transducer and activator of transcription (STAT) pathway. The phosphorylation of STAT3 at early reperfusion was not increased by ischaemic preconditioning, different from the established STAT3 activation by cardioprotective interventions in other pig strains. Ossabaw pigs have not only the genetic predisposition to develop a metabolic syndrome but also are not amenable to cardioprotection by ischaemic preconditioning.

Detection and Validation of Circular DNA Fragments Using Nanopore Sequencing.

Occurrence of extra-chromosomal circular DNA is a phenomenon frequently observed in tumor cells, and the presence of such DNA has been recognized as a marker of adverse outcome across cancer types. We here describe a computational workflow for identification of DNA circles from long-read sequencing data. The workflow is implemented based on the Snakemake workflow management system. Its key step uses a graph-theoretic approach to identify putative circular fragments validated on simulated reads. We then demonstrate robustness of our approach using nanopore sequencing of selectively enriched circular DNA by highly sensitive and specific recovery of plasmids and the mitochondrial genome, which is the only circular DNA in normal human cells. Finally, we show that the workflow facilitates detection of larger circular DNA fragments containing extrachromosomal copies of the MYCN oncogene and the respective breakpoints, which is a potentially useful application in disease monitoring of several cancer types.

Fast lightweight accurate xenograft sorting.

MOTIVATION: With an increasing number of patient-derived xenograft (PDX) models being created and subsequently sequenced to study tumor heterogeneity and to guide therapy decisions, there is a similarly increasing need for methods to separate reads originating from the graft (human) tumor and reads originating from the host species' (mouse) surrounding tissue. Two kinds of methods are in use: On the one hand, alignment-based tools require that reads are mapped and aligned (by an external mapper/aligner) to the host and graft genomes separately first; the tool itself then processes the resulting alignments and quality metrics (typically BAM files) to assign each read or read pair. On the other hand, alignment-free tools work directly on the raw read data (typically FASTQ files). Recent studies compare different approaches and tools, with varying results. RESULTS: We show that alignment-free methods for xenograft sorting are superior concerning CPU time usage and equivalent in accuracy. We improve upon the state of the art sorting by presenting a fast lightweight approach based on three-way bucketed quotiented Cuckoo hashing. Our hash table requires memory comparable to an FM index typically used for read alignment and less than other alignment-free approaches. It allows extremely fast lookups and uses less CPU time than other alignment-free methods and alignment-based methods at similar accuracy. Several engineering steps (e.g., shortcuts for unsuccessful lookups, software prefetching) improve the performance even further. AVAILABILITY: Our software xengsort is available under the MIT license at http://gitlab.com/genomeinformatics/xengsort . It is written in numba-compiled Python and comes with sample Snakemake workflows for hash table construction and dataset processing.

Proteomic and bioinformatic profiling of neutrophils in CLL reveals functional defects that predispose to bacterial infections.

Patients with chronic lymphocytic leukemia (CLL) typically suffer from frequent and severe bacterial infections. Although it is well known that neutrophils are critical innate immune cells facilitating the early defense, the underlying phenotypical and functional changes in neutrophils during CLL remain largely elusive. Using a murine adoptive transfer model of CLL, we demonstrate aggravated bacterial burden in CLL-bearing mice upon a urinary tract infection with uropathogenic Escherichia coli. Bioinformatic analyses of the neutrophil proteome revealed increased expression of proteins associated with interferon signaling and decreased protein expression associated with granule composition and neutrophil migration. Functional experiments validated these findings by showing reduced levels of myeloperoxidase and acidification of neutrophil granules after ex vivo phagocytosis of bacteria. Pathway enrichment analysis indicated decreased expression of molecules critical for neutrophil recruitment, and migration of neutrophils into the infected urinary bladder was significantly reduced. These altered migratory properties of neutrophils were also associated with reduced expression of CD62L and CXCR4 and correlated with an increased incidence of infections in patients with CLL. In conclusion, this study describes a molecular signature of neutrophils through proteomic, bioinformatic, and functional analyses that are linked to a reduced migratory ability, potentially leading to increased bacterial infections in patients with CLL.

Epitope similarity cannot explain the pre-formed T cell immunity towards structural SARS-CoV-2 proteins.

The current pandemic is caused by the SARS-CoV-2 virus and large progress in understanding the pathology of the virus has been made since its emergence in late 2019. Several reports indicate short lasting immunity against endemic coronaviruses, which contrasts studies showing that biobanked venous blood contains T cells reactive to SARS-CoV-2 S-protein even before the outbreak in Wuhan. This suggests a preformed T cell memory towards structural proteins in individuals not exposed to SARS-CoV-2. Given the similarity of SARS-CoV-2 to other members of the Coronaviridae family, the endemic coronaviruses appear likely candidates to generate this T cell memory. However, given the apparent poor immunological memory created by the endemic coronaviruses, immunity against other common pathogens might offer an alternative explanation. Here, we utilize a combination of epitope prediction and similarity to common human pathogens to identify potential sources of the SARS-CoV-2 T cell memory. Although beta-coronaviruses are the most likely candidates to explain the pre-existing SARS-CoV-2 reactive T cells in uninfected individuals, the SARS-CoV-2 epitopes with the highest similarity to those from beta-coronaviruses are confined to replication associated proteins-not the host interacting S-protein. Thus, our study suggests that the observed SARS-CoV-2 pre-formed immunity to structural proteins is not driven by near-identical epitopes.

Machine learning reveals a PD-L1-independent prediction of response to immunotherapy of non-small cell lung cancer by gene expression context.

OBJECTIVE: Current predictive biomarkers for PD-1 (programmed cell death protein 1)/PD-L1 (programmed death-ligand 1)-directed immunotherapy in non-small cell lung cancer (NSCLC) mostly focus on features of tumour cells. However, the tumour microenvironment and immune context are expected to play major roles in governing therapy response. Against this background, we set out to apply context-sensitive feature selection and machine learning approaches on expression profiles of immune-related genes in diagnostic biopsies of patients with stage IV NSCLC. METHODS: RNA expression levels were determined using the NanoString nCounter platform in formalin-fixed paraffin-embedded tumour biopsies obtained during the diagnostic workup of stage IV NSCLC from two thoracic oncology centres. A 770-gene panel covering immune-related genes and control genes was used. We applied supervised machine learning methods for feature selection and generation of predictive models. RESULTS: Feature selection and model creation were based on a training cohort of 55 patients with recurrent NSCLC treated with PD-1/PD-L1 antibody therapy. Resulting models identified patients with superior outcomes to immunotherapy, as validated in two subsequently recruited, separate patient cohorts (n = 67, hazard ratio = 0.46, p = 0.035). The predictive information obtained from these models was orthogonal to PD-L1 expression as per immunohistochemistry: Selecting by PD-L1 positivity at immunohistochemistry plus model prediction identified patients with highly favourable outcomes. Independence of PD-L1 positivity and model predictions were confirmed in multivariate analysis. Visualisation of the models revealed the predictive superiority of the entire 7-gene context over any single gene. CONCLUSION: Using context-sensitive assays and bioinformatics capturing the tumour immune context allows precise prediction of response to PD-1/PD-L1-directed immunotherapy in NSCLC.

Improved risk-stratification for posterior fossa ependymoma of childhood considering clinical, histological and genetic features - a retrospective analysis of the HIT ependymoma trial cohort.

INTRODUCTION: Risk stratification of children with ependymomas of the posterior fossa in current therapeutic protocols is mainly based on clinical criteria. We aimed to identify independent outcome predictors for this disease entity by a systematic integrated analysis of clinical, histological and genetic information in a defined cohort of patients treated according to the German HIT protocols. METHODS: Tumor samples of 134 patients aged 0.2-15.9 years treated between 1999 and 2010 according to HIT protocols were analyzed for histological features including mitotic activity, necrosis and vascular proliferation and genomic alterations by SNP and molecular inversion probe analysis. Survival analysis was performed by Kaplan-Meier method with log rank test and multivariate Cox regression analysis. RESULTS: Residual tumor after surgery, chromosome 1q gain and structural genomic alterations were identified as predictors of significantly shorter event-free (EFS) and overall survival (OS). Furthermore, specific histological features including vascular proliferation, necrosis and high mitotic activity were predictive for shorter OS. Multivariate Cox regression revealed residual tumor, chromosome 1q gain and mitotic activity as independent predictors of both EFS and OS. Using these independent predictors of outcome, we were able to build a 3-tiered risk stratification model that separates patients with standard, intermediate and high risk, and which outperforms current stratification procedures. CONCLUSION: The integration of defined clinical, histological and genetic parameters led to an improved risk-stratification model for posterior fossa ependymoma of childhood. After validation in independent cohorts this model may provide the basis for risk-adapted treatment of children with ependymomas of the posterior fossa.

Impact of RAS mutation subtype on clinical outcome-a cross-entity comparison of patients with advanced non-small cell lung cancer and colorectal cancer.

Mutated RAS onco-proteins are key drivers across many cancers. The distribution of somatic RAS mutations varies between cancer entities. Retrospective analyses have associated some RAS mutations with distinct clinical outcomes. However, the clinical impact of the full spectrum of RAS mutations in their disease contextuality remains to be defined. To improve upon this situation, we studied genomically and clinically annotated, prospectively recruited cohorts of patients with RAS-mutated metastatic lung cancer and colorectal cancer. Mutational spectra were compared with predictions derived from analyzing the mutagenic impact at the genome level for each entity. Interestingly, we found concordance of predicted signatures with those actually observed in our patients. Thus, composition of the functionally active RAS mutational subtypes is primarily determined by the mutagenic context. Most RAS mutations seemed dominant oncogenic drivers with entity-dependent clinical outcomes. RAS comutations were enriched in tumors harboring class 2/3 BRAF mutations, highlighting the functional dependency of some mutated BRAF isoforms on RAS. With our dataset, we established a probabilistic model for cross-entity comparison of the prognostic impact of specific RAS mutational subtypes. The resulting prognostic clusters showed largely consistent clinical categorizations in both entities. This suggests mutant subtype-specific functional properties leading to similar clinical effects. A notable exception is KRAS G12C, which imparted an adverse prognosis only in colorectal cancer. Our findings provide a framework for risk stratification of specific RAS mutations across several cancer entities, which is required to guide the analysis of clinical findings in patients treated with direct RAS inhibitors or agents targeting downstream pathways.

A mechanistic classification of clinical phenotypes in neuroblastoma.

Neuroblastoma is a pediatric tumor of the sympathetic nervous system. Its clinical course ranges from spontaneous tumor regression to fatal progression. To investigate the molecular features of the divergent tumor subtypes, we performed genome sequencing on 416 pretreatment neuroblastomas and assessed telomere maintenance mechanisms in 208 of these tumors. We found that patients whose tumors lacked telomere maintenance mechanisms had an excellent prognosis, whereas the prognosis of patients whose tumors harbored telomere maintenance mechanisms was substantially worse. Survival rates were lowest for neuroblastoma patients whose tumors harbored telomere maintenance mechanisms in combination with RAS and/or p53 pathway mutations. Spontaneous tumor regression occurred both in the presence and absence of these mutations in patients with telomere maintenance-negative tumors. On the basis of these data, we propose a mechanistic classification of neuroblastoma that may benefit the clinical management of patients.

Circulating microRNA biomarkers for metastatic disease in neuroblastoma patients.

In this study, the circulating miRNome from diagnostic neuroblastoma serum was assessed for identification of noninvasive biomarkers with potential in monitoring metastatic disease. After determining the circulating neuroblastoma miRNome, 743 miRNAs were screened in 2 independent cohorts of 131 and 54 patients. Evaluation of serum miRNA variance in a model testing for tumor stage, MYCN status, age at diagnosis, and overall survival revealed tumor stage as the most significant factor impacting miRNA abundance in neuroblastoma serum. Differential abundance analysis between patients with metastatic and localized disease revealed 9 miRNAs strongly associated with metastatic stage 4 disease in both patient cohorts. Increasing levels of these miRNAs were also observed in serum from xenografted mice bearing human neuroblastoma tumors. Moreover, murine serum miRNA levels were strongly associated with tumor volume. These findings were validated in longitudinal serum samples from metastatic neuroblastoma patients, where the 9 miRNAs were associated with disease burden and treatment response.

Cancer evolution, mutations, and clonal selection in relapse neuroblastoma.

The notion of cancer as a complex evolutionary system has been validated by in-depth molecular analyses of tumor progression over the last years. While a complex interplay of cell-autonomous programs and cell-cell interactions determines proliferation and differentiation during normal development, intrinsic and acquired plasticity of cancer cells allow for evasion of growth factor limitations, apoptotic signals, or attacks from the immune system. Treatment-induced molecular selection processes have been described by a number of studies already, but understanding of those events facilitating metastatic spread, organ-specific homing, and resistance to anoikis is still in its early days. In principle, somatic events giving rise to cancer progression should be easier to follow in childhood tumors bearing fewer mutations and genomic aberrations than their counterparts in adulthood. We have previously reported on the genetic events accompanying relapsing neuroblastoma, a solid tumor of early childhood. Our results indicated significantly higher single nucleotide variants in relapse tumors, gave hints for branched tumor evolution upon treatment and clonal selection as deduced from shifts in allelic frequencies between primary and relapsing neuroblastoma. Here, we will review these findings and give an outlook on dealing with intratumoral heterogeneity and sub-clonal diversity in neuroblastoma for future targeted treatments.

Recurrent alterations of TNFAIP3 (A20) in T-cell large granular lymphocytic leukemia.

The pathogenesis of T-cell large granular lymphocytic leukemia (T-LGL) is poorly understood, as STAT3 mutations are the only known frequent genetic lesions. Here, we identified non-synonymous alterations in the TNFAIP3 tumor suppressor gene in 3 of 39 T-LGL. In two cases these were somatic mutations, in one case the somatic origin was likely. A further case harbored a SNP that is a known risk allele for autoimmune diseases and B cell lymphomas. Thus, TNFAIP3 mutations represent recurrent genetic lesions in T-LGL that affect about 8% of cases, likely contributing to deregulated NF-κB activity in this leukemia.

Human TLR8 senses UR/URR motifs in bacterial and mitochondrial RNA.

Toll-like receptor (TLR) 13 and TLR2 are the major sensors of Gram-positive bacteria in mice. TLR13 recognizes Sa19, a specific 23S ribosomal (r) RNA-derived fragment and bacterial modification of Sa19 ablates binding to TLR13, and to antibiotics such as erythromycin. Similarly, RNase A-treated Staphylococcus aureus activate human peripheral blood mononuclear cells (PBMCs) only via TLR2, implying single-stranded (ss) RNA as major stimulant. Here, we identify human TLR8 as functional TLR13 equivalent that promiscuously senses ssRNA. Accordingly, Sa19 and mitochondrial (mt) 16S rRNA sequence-derived oligoribonucleotides (ORNs) stimulate PBMCs in a MyD88-dependent manner. These ORNs, as well as S. aureus-, Escherichia coli-, and mt-RNA, also activate differentiated human monocytoid THP-1 cells, provided they express TLR8. Moreover, Unc93b1(-/-)- and Tlr8(-/-)-THP-1 cells are refractory, while endogenous and ectopically expressed TLR8 confers responsiveness in a UR/URR RNA ligand consensus motif-dependent manner. If TLR8 function is inhibited by suppression of lysosomal function, antibiotic treatment efficiently blocks bacteria-driven inflammatory responses in infected human whole blood cultures. Sepsis therapy might thus benefit from interfering with TLR8 function.

Mutational dynamics between primary and relapse neuroblastomas.

Neuroblastoma is a malignancy of the developing sympathetic nervous system that is often lethal when relapse occurs. We here used whole-exome sequencing, mRNA expression profiling, array CGH and DNA methylation analysis to characterize 16 paired samples at diagnosis and relapse from individuals with neuroblastoma. The mutational burden significantly increased in relapsing tumors, accompanied by altered mutational signatures and reduced subclonal heterogeneity. Global allele frequencies at relapse indicated clonal mutation selection during disease progression. Promoter methylation patterns were consistent over disease course and were patient specific. Recurrent alterations at relapse included mutations in the putative CHD5 neuroblastoma tumor suppressor, chromosome 9p losses, DOCK8 mutations, inactivating mutations in PTPN14 and a relapse-specific activity pattern for the PTPN14 target YAP. Recurrent new mutations in HRAS, KRAS and genes mediating cell-cell interaction in 13 of 16 relapse tumors indicate disturbances in signaling pathways mediating mesenchymal transition. Our data shed light on genetic alteration frequency, identity and evolution in neuroblastoma.

N6-adenosine methylation in MiRNAs.

Methylation of N6-adenosine (m6A) has been observed in many different classes of RNA, but its prevalence in microRNAs (miRNAs) has not yet been studied. Here we show that a knockdown of the m6A demethylase FTO affects the steady-state levels of several miRNAs. Moreover, RNA immunoprecipitation with an anti-m6A-antibody followed by RNA-seq revealed that a significant fraction of miRNAs contains m6A. By motif searches we have discovered consensus sequences discriminating between methylated and unmethylated miRNAs. The epigenetic modification of an epigenetic modifier as described here adds a new layer to the complexity of the posttranscriptional regulation of gene expression.