22928 results for "ega"

in 76.91 milliseconds.

• THAP11 mutations in a patient with a cblX-like phenotype implicates THAP11 in the regulation of cobalamin metabolism and early vertebrate development

  Mutations in HCFC1 are associated with cblX (MIM309541), an X-linked recessive disorder, characterized by defects in cobalamin metabolism and other developmental defects. HCFC1 is a transcriptional co-regulator, which interacts with transcription factors to regulate the expression of a myriad of genes. Notably, HCFC1 regulates cobalamin metabolism via the regulation of MMACHC expression through its interaction with THAP11, a THAP domain-containing transcription factor. In addition, the HCFC1/THAP11 complex potentially regulates genes involved in diverse cellular functions including cell cycle, proliferation, and transcription. Thus, it is likely that mutations in THAP11 may also result in biochemical and other phenotypes similar to those observed in patients with cblX. We report a patient who presented with phenotypic features that overlap cblX, but did not have any mutations in either MMACHC or HCFC1. We sequenced THAP11 by Sanger sequencing and discovered a potentially pathogenic, homozygous variant in THAP11, c.240C>G (p.Phe80Leu). Functional analysis in the developing zebrafish embryo, demonstrate that both THAP11 and HCFC1 regulate the proliferation and differentiation of neural precursors, suggesting important roles in normal brain development. Further, the loss of THAP11 in zebrafish embryos, results in craniofacial abnormalities including the complete loss of Meckel’s cartilage, the ceratohyal, and all of the ceratobranchial cartilages. These data are consistent with our previous work which demonstrated a role for HCFC1 in vertebrate craniofacial development. Furthermore, our high throughput RNA-sequencing analysis reveals several overlapping gene targets of HCFC1 and THAP11 that likely contribute to the phenotypes associated with cblX and related disorders. Thus, both HCFC1 and THAP11 play important roles in the regulation of cobalamin metabolism as well as other pathways involved in early vertebrate development.

  Study EGAS00001002201

• The application of RNA sequencing for the diagnosis and genomic classification of pediatric acute lymphoblastic leukemia

  Acute lymphoblastic leukaemia (ALL) is the most common childhood malignancy, and implementation of risk-adapted therapy has been instrumental in the dramatic improvements in clinical outcomes. A key to risk-adapted therapies includes the identification of genomic features of individual tumors, including chromosome number (for hyper- and hypodiploidy) and gene fusions, notably ETV6-RUNX1, TCF3-PBX1 and BCR-ABL1 in B-cell ALL (B-ALL). RNA-sequencing (RNA-seq) of large ALL cohorts has expanded the number of recurrent gene fusions recognized as drivers in ALL, and identification of these new entities will contribute to refining ALL risk stratification. We used RNA-seq on 126 ALL patients from our clinical service to test the utility of including RNA-seq in standard of care diagnostic pipelines to detect gene rearrangements and IKZF1 deletions. RNA-seq identified 86% of rearrangements detected by standard of care diagnostics. KMT2A-rearrangements (MLL), whilst usually identified, were the most commonly missed by RNA-seq as a result of low expression. RNA-seq identified rearrangements that were not detected by standard of care testing in nine patients. These were found in patients who were not classifiable using standard molecular assessment. We developed an approach to detect the commonest IKZF1 deletion from RNA-seq data and validated this using a RQ-PCR assay. We applied an expression classifier to identify Ph-like B-ALL patients. T-ALL proved a rich source of novel gene fusions which have clinical implications or provide insights into disease biology. Our experience shows that RNA-seq can be implemented within an individual clinical service to enhance the current molecular diagnostic risk classification of ALL.

  Study EGAS00001004212

• European BestAgeing Study on microRNA candidates for cardiovascular disease

  Background: Circulating miRNAs have emerged as promising biomarker candidates due to their stability and their role in regulating key pathological pathways in cardiovascular disease (CVD). Yet, large-scale, multicenter studies examining their diagnostic and prognostic potential are scarce. This study evaluates the potential of miRNA expression profiles to inform disease classification and risk stratification across major CVD phenotypes, including acute coronary syndrome (ACS), chronic coronary artery disease (CAD), dilated cardiomyopathy (DCM), and ischemic cardiomyopathy (ICM), in a large, multicenter European cohort. Methods: We assessed genome-wide miRNA expression profiles in a total of 1,209 cardiovascular patients and 848 controls in a uniform, standardized fashion, which renders this study one of the largest prospective miRNA studies. To focus on only the most biologically plausible miRNAs for clinical translation, we mined all original studies of miRNA candidates in CVD and performed differential miRNA expression and enrichment analysis. We then trained disease-specific binary classification models to evaluate the diagnostic potential of miRNA signatures. Finally, we evaluated prognosis and disease severity based on distinct miRNA levels. Results: 634 original abstracts were identified, detailing 166 ACS, 181 CAD, 56 DCM, and 182 ICM miRNAs. Without further optimization, the signatures of a priori miRNAs already yielded very good diagnostic performance with ROC AUC of 0.83 – 0.95. There was an improvement when considering additional miRNAs in a discovery setting. Interestingly, in ACS, CAD and DCM we observed a significantly worse prognosis in probands with higher scores in the miRNA signatures, indicating additional prognostic information. Conclusions: The European BestAgeing miRNA study reveals emerging associations of several miRNA signatures with cardiovascular disease discrimination and prognostication, providing a foundation for future external validation and potential clinical translation of this class of markers.

  Study EGAS00001008346

• Tracking the evolution of Therapy-Related Myeloid Neoplasms using chemotherapy signatures

  Patients treated with cytotoxic therapies, including autologous stem cell transplantation, are at risk for developing therapy-related myeloid neoplasms (tMN). Pre-leukemic clones (i.e., clonal hematopoiesis; CH) are detectable years before the development of these aggressive malignancies, though the genomic events leading to transformation and expansion are not well-defined. Here, leveraging distinctive chemotherapy-associated mutational signatures from whole-genome sequencing data and targeted sequencing of pre-chemotherapy samples, we reconstruct the evolutionary life-history of 39 therapy-related myeloid malignancies. A dichotomy is revealed, in which neoplasms with evidence of chemotherapy-induced mutagenesis from platinum and melphalan are hypermutated and enriched for complex structural variants (i.e., chromothripsis) while neoplasms with non-mutagenic chemotherapy exposures are genomically similar to de novo acute myeloid leukemia. Using chemotherapy-associated mutational signatures as temporal barcodes linked to a discrete clinical exposure in each patient’s life, we estimate that several complex events and genomic drivers are acquired after chemotherapy is administered. For patients with prior multiple myeloma who were treated with high-dose melphalan and autologous stem cell transplantation, we demonstrate that tMN can develop from either a reinfused CH clone that escapes chemotherapy exposure and is selected following reinfusion, or from TP53-mutant CH that survives direct myeloablative conditioning and acquires melphalan-induced DNA-damage. Overall, we reveal a novel mode of tMN progression that is not reliant on direct mutagenesis or even exposure to chemotherapy. Conversely, for tMN that evolve under the influence of chemotherapy-induced mutagenesis, distinct chemotherapies not only select pre-existing CH, but also promote the acquisition of recurrent genomic drivers.

  Study EGAS00001006903

• ctDNA monitoring using patient-specific sequencing and integration of variant reads - Lung cohort

  Circulating tumor-derived DNA (ctDNA) can be used to monitor cancer dynamics noninvasively. Detection of ctDNA can be challenging in patients with low-volume or residual disease, where plasma contains very few tumor-derived DNA fragments. We show that sensitivity for ctDNA detection in plasma can be improved by analyzing hundreds to thousands of mutations that are first identified by tumor genotyping. We describe the INtegration of VAriant Reads (INVAR) pipeline, which combines custom error-suppression methods and signal-enrichment approaches based on biological features of ctDNA. With this approach, the detection limit in each sample can be estimated independently based on the number of informative reads sequenced across multiple patient-specific loci. We applied INVAR to custom hybrid-capture sequencing data from 176 plasma samples from 105 patients with melanoma, lung, renal, glioma, and breast cancer across both early and advanced disease. By integrating signal across a median of >105 informative reads, ctDNA was routinely quantified to 1 mutant molecule per 100,000, and in some cases with high tumor mutation burden and/or plasma input material, to individual parts per million. This resulted in median Area Under the Curve (AUC) values of 0.98 in advanced cancers, and 0.80 in early stage and challenging settings for ctDNA detection. We generalized this method to whole-exome and whole-genome sequencing, showing that the INVAR may be applied without requiring personalized sequencing panels, so long as a tumor mutation list is available. As tumor sequencing becomes increasingly performed, such methods for personalized cancer monitoring may enhance the sensitivity of cancer liquid biopsies.

  Study EGAS00001004447

• Deep single-cell RNA sequencing data for 12346 T cells from tumour, adjacent normal tissue and peripheral blood of treatment-naive NSCLC patients

  Cancer immunotherapies have shown sustained clinical responses in treating non-small cell lung cancer (NSCLC), but efficacy varies between patients and is believed to depend in part on the amount and properties of tumor infiltrating lymphocytes (TILs). To comprehensively depict and dissect the baseline landscape of the composition, lineage and functional states of TILs in lung cancer, here we generated deep single-cell RNA sequencing data for 12,346 T cells from the primary tumour, adjacent normal tissues and peripheral blood of 14 treatment-naive NSCLC patients. Combined expression and TCR-based lineage tracking revealed a significant proportion of effector T cells with common origins and similar functional states across peripheral blood and tumours pointing towards a highly migratory nature of these T cells. We also observed tumour-infiltrating CD8+ T cells undergoing extensive clonal expansion and exhaustion, with two clusters of cells exhibiting states preceding exhaustion. Survival analysis on independent datasets suggested that high ratio of pre-exhausted to exhausted T cells is associated with better prognosis of lung adenocarcinoma (LUAD). In addition, we observed further heterogeneity within the tumour regulatory T cells (Tregs), characterized by the bimodal distribution of TNFRSF9, an activation marker for antigen-specific Tregs. The gene signature of this group of activated tumour Tregs, which included IL1R2, correlated with poor prognosis in LUAD. The T cell clusters revealed by our single cell analyses provide a new approach for patient stratification, and the accompanying compendium of data will help the research community to gain further insight into the functional states and dynamics of T cell responses in lung cancer.

  Study EGAS00001002430

• Biomarker Screen for Efficacy of Oncolytic Virotherapy in Patient-derived Pancreatic Cancer Cultures

  Pancreatic ductal adenocarcinoma (PDAC) is a tumour entity with unmet medical need. To assess the therapeutic potential of oncolytic virotherapy (OVT) against PDAC, different oncolytic viruses (OVs) are currently investigated in clinical trials. However, systematic comparisons of these different OVs in terms of efficacy against PDAC and biomarkers predicting therapeutic response are lacking. We screened fourteen patient-derived PDAC cultures for their sensitivity to five clinically relevant OVs, namely serotype 5 adenovirus Ad5-hTERT, herpes virus T-VEC, measles vaccine strain MV-NIS, reovirus jin-3, and protoparvovirus H 1PV using live cell analysis, quantification of viral genome/gene expression, cell viability as well as cytotoxicity assays. Sensitivity was correlated with transcriptome profiles to identify potential predictive biomarkers for response to OVT. Patient-derived PDAC cultures showed individual response patterns to OV treatment. Twelve of fourteen cultures were responsive to at least one OV, with no single OV proving superior or inferior across all cultures. Known host factors for distinct viruses were retrieved as potential biomarkers. Compared to the basal-like molecular subtype, the quasi-mesenchymal subtype of PDAC was more sensitive to H-1PV, jin-3, and T-VEC. Expression of viral entry receptors did not correlate with sensitivity to OV treatment. Rather, cellular pathways controlling immunological, metabolic, and proliferative signalling appeared to determine outcome. For instance, high baseline expression of interferon-stimulated genes (ISGs) correlated with relative resistance to oncolytic measles virus, whereas cGAS expression was associated with exceptional response. Combination treatment of MV-NIS with a cGAS inhibitor improved tumour cell killing in several PDAC cultures.

  Study EGAS00001007001

• Whole-genome Sequencing Suggests Mechanisms for 22q11.2 deletion-associated Parkinson’s disease

  Objectives: To investigate disease risk mechanisms of early-onset Parkinson’s disease (PD) associated with the recurrent 22q11.2 deletion, a genetic risk factor for early-onset PD. Methods: In a proof-of-principle study, we used whole-genome sequencing (WGS) to investigate sequence variants in nine adults with 22q11.2DS, three with neuropathologically confirmed early-onset PD and six without PD. Adopting an approach used recently to study schizophrenia in 22q11.2DS, here we tested candidate gene-sets relevant to PD. Results: No mutations common to the cases with PD were found in the intact 22q11.2 region. While all were negative for rare mutations in a gene-set comprising PD disease-causing and risk genes, another candidate gene-set of 1000 genes functionally relevant to PD presented a nominally significant (P=0.03) enrichment of rare putatively damaging missense variants in the PD cases. Polygenic score results, based on common variants associated with PD risk, were non-significantly greater in those with PD. Conclusions: The results of this first-ever pilot study of WGS in PD suggest that the cumulative burden of genome-wide sequence variants may contribute to expression of early-onset PD in the presence of threshold-lowering dosage effects of a 22q11.2 deletion. We found no evidence that expression of PD in 22q11.2DS is mediated by a recessive locus on the intact 22q11.2 chromosome or mutations in known PD genes. These findings offer initial evidence of the potential effects of multiple within-individual rare variants on the expression of PD and the utility of next generation sequencing for studying the etiology of PD.

  Study EGAS00001002275

• Multi-Layered Molecular Profiling Informs the Diagnosis and Targeted Therapy of Desmoplastic Small Round Cell Tumor

  Desmoplastic small round cell tumor (DSRCT) is an ultra-rare soft tissue sarcoma with dismal prognosis and limited therapeutic options. Here, we identified individual therapeutic targets in a cohort of 30 DSRCT patients by use of multilayered molecular profiling (genomic, transcriptomic, and (phospho-)proteomic analyses) within the DKFZ NCT DKTK MASTER program. 28 of the 30 patients received at least one personalized treatment recommendation. In line with the observation that DSRCTs are an entity with a low mutational load, most of the treatment recommendations were based on target gene or -protein expression. 13 patients received a recommended therapy, resulting in four partial responses, four stable diseases, and five progressive diseases. Eight of 13 patients received the multi tyrosine kinase inhibitor pazopanib, with mixed response. Based on high ERBB2 expression levels, two patients received the antibody drug conjugate trastuzumab deruxtecan (T-DXd), showed exceptional, long-lasting responses and are currently (8/2024) still on treatment and progression-free. We detected high ERBB2 protein expression in all DSRCT samples with proteome data (n = 9), which suggests that ERBB2 may be a particularly useful therapeutic target for DSRCT patients. For one of the two exceptional responders we did not detect a substantial ERBB2/HER2 overexpression, which may indicate that T-DXd could also be efficacious in DSRCT patients with low levels of HER2, as also described for breast cancer patients. In conclusion, multi-omic molecular profiling revealed personalized therapeutic targets in patients with DSRCT with HER2 targeting, with T-DXd as a promising therapeutic option for this ultra-rare and difficult-to-treat tumor entity.

  Study EGAS00001007934

• Whole genome sequencing of tumor samples from advanced pre-treated NSCLC patients recruited in a phase I/II single-arm trial utilising CVA21, an oncolytic coxsackie virus, in combination with pembrolizumab.

  Acquired or de novo resistance to immune checkpoint inhibitors occurs in the majority of advanced non-small cell lung cancers (NSCLC). There is an unmet need to improve outcomes for patients with this condition. Oncolytic viruses represent an attractive treatment approach due to their dual activity in inducing tumor cell lysis directly and potentially augmenting anti-tumor immunity. Here, we present the safety, efficacy and translational findings from a phase I/II single-arm trial utilising CVA21, an oncolytic coxsackie virus, in combination with pembrolizumab in patients with advanced pre-treated NSCLC. We performed paired biopsies pre- and post-treatment in 10 patients who received intravenous CVA21 and pembrolizumab, 8 of whom had prior treatment with immune checkpoint inhibitor (ICI) therapy. Whole genome sequencing and spatial proteomics were performed to comprehensively characterise the response to CVA21. Combination CVA21/pembrolizumab (anti PD-1) therapy was well tolerated with no serious treatment-related adverse events. Partial responses were seen in two patients with prior acquired anti-PD-1 resistance and disease stabilisation in 6 patients, giving a clinical benefit rate of 80%. High baseline tumor mutational burden and PD-L1 expression were observed in patients with better response to treatment. Interestingly, an increase in antigen presentation and CD8+ T cell infiltration was observed on treatment compared with baseline in patients with better progression-free survival. This study demonstrates the potential of CVA21 to modulate the immunogenicity of tumor cells and remodelling the tumor microenvironment, providing novel insights for patient selection for trials involving novel immunotherapeutic approaches.

  Study EGAS00001008258

• Sequencing of paediatric High Grade Gliomas and DIPG

  Paediatric high grade glioma (pHGG) and diffuse intrinsic pontine glioma (DIPG) are a phenotypically diverse collection of tumours which molecular profiling studies have refined into age- and location-based subgroups driven by unique genetic and epigenetic alterations. Despite this, the rarity of these tumours means individual studies are underpowered to investigate subgroup-specific changes and to identify low frequency events recurrently targeting specific pathways. We have collated genomic data from 143 unpublished cases up to 30 years of age along with those from 18 previously published datasets in an integrated analysis comprising 1067 unique patients across all anatomical compartments of the CNS, with histone mutation status available for 893 cases (n=67 H3.3G34R/V, n=316 H3.3K27M, n=68 H3.1/H3.2K27M). Distinct clinicopathological and molecular subgroups are defined by these histone mutations, with several novel co-segregating mutations identified through these aggregated analyses, including loss of FBXW7 in H3.3G34R/V, TOP3A rearrangements in H3.3K27M, and BCOR mutations in H3.1K27M. Histone wild-type subgroups can be further refined by the presence of key oncogenic events (PDGFRA or EGFR amplification/mutations) or methylation profiles which molecularly more closely resemble pleomorphic xanthoastrocytoma or other low grade gliomas. The burden of genomic aberrations decreases with age, highlighting the infant population as biologically and clinically distinct. Across the whole cohort, we identify novel previously unrecognised pathway dysregulation in small subsets of tumours (e.g. splicing, WNT, immune response) which make up the umbrella classification of paediatric diffusely infiltrating gliomas. The integrated dataset helps to further define the molecular diversity of the disease, opening up novel avenues for biological study and providing a basis for functionally defined future treatment stratification.

  Study EGAS00001002314

• Suppressors and activators of JAK-STAT signaling at diagnosis and relapse of acute lymphoblastic leukemia in Down Syndrome

  Children with Down Syndrome (DS) are prone to development of high risk B-cell precursor ALL (DS-ALL) that differs genetically from most sporadic pediatric ALLs. Increased expression of CRLF2, the receptor to thymic stromal lymphopoietin (TSLP) characterizes about half of DS-ALLs and also a subgroup of sporadic “Philadelphia-like ALLs”. To understand the pathogenesis of relapsed DS-ALL we performed integrative genomic analysis of 25 matched diagnosis remission and relapse DS-ALLs. We found that the CRLF2 rearrangements are early events during DS-ALL evolution and generally stable between diagnoses and relapse. Secondary activating signaling events in JAK-STAT/RAS pathway were ubiquitous but highly redundant between diagnosis and relapse, suggesting that signaling is essential but that no specific mutations are “relapse driving”. We further found that activated JAK2 may be naturally suppressed in 25% of CRLF2 positive DS-ALL by loss-of-function aberrations in USP9X, a deubuitinase previously shown to stabilize the activated phosphorylated JAK2. Interrogation of large ALL genomic databases extended our findings up to 25% of CRLF2 positive, “Ph like” ALLs. Pharmacological or genetic inhibition of USP9X as well as treatment with low-dose Ruxolitinib enhanced the survival of pre-B ALL cells overexpressing mutated JAK2. Thus, somehow counterintuitive, we found that suppression of JAK-STAT “hyper-signaling” may be beneficial to leukemic B cell precursors. This and the reduction of JAK mutated clones at relapse suggest that the therapeutic effect of JAK specific inhibitors may be limited. Rather, combined signaling inhibitors or direct targeting of the TSLP receptor may be a useful therapeutic strategy for DS-ALL.

  Study EGAS00001002410

• Genomic diversity of the African-descent Makranis of Pakistan

  From the eighth century onward, the Indian Ocean was the scene of extensive trade of sub-Saharan African slaves, via sea routes controlled by Muslim Arab and Swahili traders. Several populations in present-day Pakistan and India are thought to be the descendants of such slaves, yet their genetic ancestries and history of admixture and natural selection remain largely undefined. Here we studied the genome-wide diversity of the African-descent Makranis, residing on the Arabian Sea coast of Pakistan, together with that of four neighboring Pakistani populations, to investigate the genetic legacy, population dynamics and tempo of the Indian Ocean slave trade. We show that the Makranis result from an admixture event between local Baluch tribes and Bantu-speaking populations from eastern/southeastern Africa, which we dated to ~300 years ago during the Omani Empire domination. Levels of inbreeding, measured through cumulative runs of homozygosity, were found to be similar across Pakistani populations, suggesting that the admixed Makranis rapidly adopted the traditional practice of endogamous marriages. Finally, we searched for signatures of post-admixture selection at several traits evolving under positive selection, including skin color, lactase persistence and resistance to malaria. We demonstrate that the African-specific Duffy-null blood group — believed to confer resistance against Plasmodium vivax infection — was recently introduced to Pakistan through the slave trade and evolved under positive selection in this vivax malaria-endemic region. Collectively, our study reconstructs the genetic and adaptive history of a neglected episode of the African Diaspora, and illustrates the impact of recent admixture on the diffusion of adaptive traits across human populations.

  Study EGAS00001002558

• Ultra-Fast Patient-Derived Xenografts Identify Functional and Spatial Tumour Heterogeneities that Drive Therapeutic Resistance

  Treatment of patients with advanced cancers increasingly relies on agents targeting specific molecular or cellular aberrations. These expensive therapies often fail because of pre-existing drug resistance. Intra-tumoural molecular and micro-environmental heterogeneity resistance can create resistance in one or more subclones of a tumour, leading to spatio-functional heterogeneity. We created a rapid approach for quantifying tumour functional phenotypes based on tumour-implantation into the chorioallantoic membrane of chick embryos. We apply this technique to 43 separate biopsies from 6 patients with renal cell carcinoma (RCC) with 36 replicates per biopsy, achieving a 93.6% engraftment rate (1449/1548). Challenging these models with sunitinib, an anti-angiogenic therapy frequently used in RCC allowed us to quantify spatio-functional heterogeneity in individual patients. Prediction of de novo drug resistance in advanced cancer patients would facilitate precision medicine by tailoring systemic therapies and improve clinical trial design outcome measures. We describe a patient-derived xenograft platform where fresh patient tumour samples are implanted into the chorioallantoic membrane (CAM) of chick embryos at high engraftment efficiencies, permitting large scale “tumour avatar” studies to guide the selection of drugs and to predict clinical drug resistance outcomes in renal cell carcinoma (RCC) patients. Applying this “tumour avatar” model in a co-clinical trial on RCC patients, we observe intra-tumoural functional heterogeneity with sunitinib treatment, revealing the presence of drug-resistant clones prior to the initiation of therapy. Multi-regional exome Matched DNA sequencing of the primary tumour and metastases revealed DNA mutational signatures associated with sunitinib resistance at baseline and after treatment. These These studies demonstrate a rapid and affordable method of findings confirm the existence of genetic and surveying functional tumour heterogeneities heterogeneity.

  Study EGAS00001002627

• The light chain IgLV3-21 defines a new poor prognostic subgroup in Chronic Lymphocytic Leukemia: results from a multicenter study

  Purpose: Unmutated (UM) immunoglobulin heavy chain region (IgHV) status or IgHV3-21 gene usage is associated with poor prognosis in chronic lymphocytic leukemia (CLL). Interestingly, IgHV3-21 is often co-expressed with light chain IgLV3-21 which is potentially able to trigger cell-autonomous BCR-mediated signaling. However this light chain has never been characterized independently of the heavy chain IgHV3-21.Experimental Design: we performed total RNA sequencing to characterize IgLV3-21 patients (n=32) and investigated IgLV3-21 prognostic impact in terms of treatment-free (TFS) and overall (OS) survival in 3 other independent cohorts for a total of 813 patients where IgLV3-21 presence was tested by real-time PCR and confirmed by Sanger sequencing.Results: Using total RNA sequencing to characterize 32 patients with high-risk CLL, we found a high frequency (28%) of IgLV3-21 rearrangements. Gene set enrichment analysis revealed that these patients present higher expression of genes responsible for ribosome biogenesis and translation initiation (P<0.0001) and MYC target genes (P=0.0003). Patients with IgLV3-21 rearrangements displayed a significantly shorter TFS and OS (P<0.05), in particular in IgHV mutated patients. In each of the 3 independent validation cohorts (n=813), we showed that IgLV3-21 rearrangements -similar to UM IgHV status- conferred poor prognosis compared to mutated IgHV (P<0.0001). Importantly, we confirmed in multivariate analysis that this was independent of IgHV mutational status or subset #2 stereotyped receptor (P<0.0001). Conclusions: We therefore demonstrate for the first time that a light chain can impact on CLL prognosis and that IgLV3-21 light chain usage defines a new poor prognostic subgroup in CLL.

  Study EGAS00001002894

• Exploiting immune cell receptor information to quantify index switching in single cell transcriptome sequencing experiment

  By offering high sequencing speed and ultra-high-throughput at a low price, Illumina next-generation sequencing platforms has been widely adopted in recent years. However, an experiment with multiplexed library could be at risk of molecular recombination, known as “index switching”, which causes certain amount of reads to be assigned to an incorrect library. It is reported that a new advance, exclusion amplification (ExAmp) in conjunction with the patterned flow cell technology introduced on HiSeq 3000/HiSeq 4000/HiSeq X sequencing systems, potentially suffers from a higher rate of index switching than conventional bridge amplification. We took advantage of the diverse but highly cell-unique expression of immune cell receptors to quantify index switching on single cell RNA-seq data that were sequenced on HiSeq 3500 and HiSeq 4000. By utilizing the unique T-cell receptor (TCR) expression, we could quantify the spread-of-signal from many different wells (n=51 from total three batches) due to index switching. We used TraCer to reconstruct full-length TCR from all samples, and then used Kallisto to quantify TCR gene expression. We found index switching in all three batches of experiments investigated. The median percentage of incorrectly detected markers was estimated to be 4.2% (interquartile range (IQR): 2.0%-8.7%). We did not detect any consistent pattern of certain indices to be more prone for switching than others, suggesting that index switching is a stochastic process. We confirm that index switching is a problem that affects all samples run in multiplexed libraries on Illumina HiSeq 3500 and HiSeq 4000 platforms.

  Study EGAS00001002911

• A genome scan for genes underlying adult body size differences between Central African hunter-gatherers and farmers.

  The evolutionary and biological bases of the Central African pygmy phenotype, a characteristic of rainforest huntergatherers defined by reduced body size compared with neighboring farmers, remain largely unknown. Here, we perform a joint investigation in Central African hunter-gatherers and farmers of adult standing height, sitting height, leg length, and body mass index (BMI), considering 358 hunter-gatherers and 169 farmers with genotypes for 153,798 SNPs. In addition to reduced standing heights, hunter-gatherers have shorter sitting heights and leg lengths and higher sitting/standing height ratios than farmers and lower BMI for males. Standing height, sitting height, and leg length are strongly correlated with inferred levels of farmer genetic ancestry, whereas BMI is only weakly correlated, perhaps reflecting greater contributions of non-genetic factors to body weight than to height. Single- and multi-marker association tests identify one region and eight genes associated with hunter-gatherer/farmer status, and 24 genes associated with the height-related traits. Many of these genes have putative functions consistent with roles in determining their associated traits and the pygmy phenotype, and they include three associated with standing height in non-Africans (PRKG1, DSCAM, MAGI2). We find evidence that European height-associated SNPs or variants in linkage disequilibrium with them contribute to standing- and sitting-height determination in Central Africans, but not to the differential status of hunter-gatherers and farmers. These findings provide new insights into the biological basis of the pygmy phenotype, and they highlight the potential of cross-population studies for exploring the genetic basis of phenotypes that vary naturally across populations. Pemberton et al. Human Genetics 2018 https://doi.org/10.1007/s00439-018-1902-3

  Study EGAS00001002975

• Evolutionary dynamics of residual disease in human glioblastoma

  Glioblastoma is the most common and aggressive adult brain malignancy against which conventional surgery and chemoradiation provide limited benefit. Even when a good treatment response is obtained, recurrence inevitably occurs either locally (c.80%) or distally (c.20%), driven by cancer clones that are often genomically distinct from those in the primary tumour. Glioblastoma cells display a characteristic infiltrative phenotype, invading the surrounding tissue and often spreading across the whole brain. We hypothesised that cancer cells responsible for relapse reside in two compartments of residual disease that are left behind after treatment: the infiltrated normal brain parenchyma, and the sub-ventricular zone (SVZ). In this model, residual disease subclones diverge early during glioblastoma evolution, may remain dormant in the normal parenchyma and are responsible for recurrence. To test our hypothesis, we performed whole-exome sequencing of 69 multi- region samples collected using fluorescence-guided resection from 11 patients, including the infiltrating tumour margin (M) and the SVZ for each patient, as well as matched blood. We used a phylogenomic approach to dissect the spatio-temporal evolution of each tumour and unveil the relation between residual disease and the main tumour mass. We also analysed two patients with paired primary-recurrence samples with matched residual disease. Our results suggest that the infiltrative phenotype is an early evolutionary trait of glioblastoma, giving rise to the tumour mass detected at surgery. After treatment, the same infiltrative clone may seed the growth of a recurrent tumour, thus representing the ‘missing link’ between the primary tumour and recurrent disease. These results are consistent with recognised clinical phenotypic behaviour and suggest that more specific therapeutic targeting of cells in the infiltrated brain parenchyma may improve patient’s outcome.

  Study EGAS00001003043

• Genomic landscape of oral cancers (Complete Genomics WGS)

  Human papillomavirus (HPV) is a necessary but insufficient cause of a subset of oral squamous cell carcinomas (OSCC) that is increasing markedly in frequency. To identify contributory, secondary genetic alterations in these cancers, we used comprehensive genomics methods to compare 149 HPV-positive and 335 HPV-negative OSCC tumor/normal pairs. Different behavioral risk factors underlying the two OSCC types were reflected in distinctive genomic mutational signatures. In HPV-positive OSCC, the signatures of APOBEC cytosine deaminase editing, associated with anti-viral immunity, were strongly linked to overall mutational burden. By contrast, in HPV-negative OSCC, T>C substitutions in the sequence context 5’-ATN-3’ correlated with tobacco exposure. Universal expression of HPV E6*1 and E7 oncogenes was a sine qua non of HPV-positive OSCC. Significant enrichment of somatic mutations was confirmed or newly identified in PIK3CA, KMT2D, FGFR3, FBXW7, DDX3X, PTEN, TRAF3, RB1, CYLD, RIPK4, ZNF750, EP300, CASZ1, TAF5, RBL1, IFNGR1 and NFKBIA. Of these, many affect host pathways already targeted by HPV oncoproteins, including the p53 and pRB pathways, or disrupt host defenses against viral infections, including interferon and nuclear factor-B signaling. Frequent copy number changes were associated with concordant changes in gene expression. Chr. 11q (including CCND1) and 14q (including DICER1 and AKT1) were recurrently lost in HPV-positive OSCC, in contrast to their gains in HPV-negative OSCC. High-ranking variant allele fractions implicated ZNF750, PIK3CA and EP300 mutations as candidate driver events in HPV-positive cancers. We conclude that virus-host interactions cooperatively shape the unique genetic features of these cancers, distinguishing them from their HPV-negative counterparts.

  Study EGAS00001002393

• Indonesian Genome Diversity Project

  Indonesia, the world’s fourth largest country by population, comprises an archipelago of about 900 permanently inhabited islands in tropical Asia, and hosts an astonishing array of human diversity that remains largely underrepresented in modern biological surveys. The region is unique for its key role in both the early and recent evolution of Homo outside of Africa. It has one of the first traces of anatomically modern humans in Eurasia, and was the focus of extensive interaction between archaic hominins and humans – both in the likely co-existence of H. floresiensis with modern humans, and in the introgression of Denisovans into the ancestors of its present inhabitants. More recently, Indonesia was a center of the spread of Neolithic culture by Austronesian speakers. Advancing maritime technologies allowed farming populations to treat this region as a springboard to reach oceanic islands as remote as Madagascar, Hawaii, Easter Island and New Zealand. Here we present complete human genomes from healthy individuals from multiple key islands in Indonesia and Papua. The genomes were sequenced to high depth (30x). The purpose of this dataset is to address different questions of human evolution in Island Southeast Asia and beyond, including early settlement of the area by anatomically modern humans, population structure, genetic admixture and adaptation, and contacts with archaic hominins. We employed a sampling strategy to capture a broad range of diversity across the region incorporating an average of 12 samples from each of 14 islands spanning Indonesia and Papua, from Sumatra in the west to New Britain in the east. Genomes from this dataset form part of an ongoing project to describe and understand human evolutionary history in Island Southeast Asia and Eurasia in general.

  Study EGAS00001003054

• Histone chaperone CHAF1A promotes proliferation and tumorigenicity in gastric cancer and impacts prognosis via context-depedent regulation of gene expression

  The histone chaperone CHAF1A is associated with some tumors, but its mechanisms in tumor biology remain elusive. Gastric cancer (GC) has significant heterogeneity among patients, and novel prognostic markers and therapeutic targets may improve its outcomes. The role of CHAF1A in GC is unknown. In this study, we demonstrated that CHAF1A was overexpressed in GC tissues accompanied by low expression of wild-type P53 and high expression of the proliferation marker MKI67. In vitro, CHAF1A knockdown inhibited cell proliferation, induced cell cycle arrest and promoted apoptosis, while CHAF1A overexpression had contrary effects. In vivo, CHAF1A knockdown inhibited tumorigenicity of GC cells. A microarray assay with CHAF1A inhibition and a PathScan signaling antibody array with CHAF1A overexpression showed that CHAF1A inhibited the P53 pathway, activated stress response and induced glycolytic metabolism. In 665 patients, the expression level of CHAF1A protein was an independent predictor for overall survival and disease-free survival in non-cardia GC and also had survival associations in patients with small or perineural invasion-negative tumors. Tumor mRNA analyses based on next-generation sequencing in patients indicated that CHAF1A promoted the Warburg effect by depressing oxidative phosphorylation and increasing glycolysis. Furthermore, the heterogeneous clinical significance of CHAF1A among patient subgroups was associated with its differential regulation of gene expression involved with glycolysis and cell proliferation and survival, indicating that the role of CHAF1A is depended on patient characteristics. These results reveal critical roles and clinical values for CHAF1A in GC. CHAF1A may have context-dependent effects in GC, which increases understandings for GC heterogeneity.

  Study EGAS00001003064

• Exome sequencing in bipolar disorder families

  Bipolar disorder (BD) is a major psychiatric disorder affecting around 1% of the global population. BD is characterized by recurrent manic and depressive episodes, and has an estimated heritability of around 70% Research has identified the first BD susceptibility genes. However, the underlying pathways and regulatory networks remain largely unknown. Since research has shown that the cumulative impact of common alleles with small effect appears to explain only around 25-38% of the phenotypic variance for BD, rare variants of high penetrance may also contribute to BD risk. The present study investigated the role of rare, nonsynonymous, and potentially functional variants via whole exome sequencing in 15 BD cases from two large, multiply affected families from Cuba. The high prevalence of BD in these pedigrees renders them promising in terms of the identification of genetic risk variants with large effect sizes. In addition, SNP array data were used to calculate polygenic risk scores for affected and unaffected family members. After correction for multiple testing, polygenic risk scores for common, BD-associated genetic variants were not significantly higher in BD cases than in healthy relatives. Exome sequencing identified a total of 17 rare and potentially damaging variants in 17 genes. The identified variants were shared by all investigated BD cases in the respective pedigree. The most promising variant was located in the gene SERPING1 (p.L349F), which has been reported as a genome-wide significant risk gene for schizophrenia. A further promising variant was located in CAPN2, which has been implicated in synaptic plasticity. The present data suggest novel candidate genes for BD susceptibility, and may facilitate the discovery of disease-relevant pathways and regulatory networks.

  Study EGAS00001003085

• Resolving the Full Spectrum of Human Genome Variation using Linked-Reads

  Large-scale population based analyses coupled with advances in technology have demonstrated that the human genome is more diverse than originally thought. Standard short-read approaches, used primarily due to accuracy, throughput and costs, fail to give a complete picture of a genome. They struggle to identify large, balanced structural events, cannot access repetitive regions of the genome and fail to resolve the human genome into its two haplotypes. Here we describe an approach that retains long range information while harnessing the power of short reads. Starting from only 1ng of DNA, we produce barcoded short read libraries. The use of novel informatic approaches allows for the barcoded short reads to be associated with the long molecules of origin producing a novel datatype known as 'Linked-Reads'. This approach allows for simultaneous detection of small and large variants from a single Linked-Read library. We have previously demonstrated the utility of whole genome Linked-Reads (lrWGS) for performing diploid, de novo assembly of individual genomes. In this manuscript, we show the utility of reference based analysis using a single Linked-Read library for full spectrum genome analysis. We demonstrate the ability of Linked-Reads to reconstruct megabase scale haplotypes and to recover parts of the genome that are typically inaccessible to short reads, including phenotypically important genes such as STRC, SMN1 and SMN2. We demonstrate the ability of both lrWGS and Linked-Read Whole Exome Sequencing (lrWES) to identify complex structural variations, including balanced events, single exon deletions, and single exon duplications. The data presented here show that Linked-Reads provide a scalable approach for comprehensive genome analysis that is not possible using short reads alone.

  Study EGAS00001003121

• RNA seq on 19 samples of Radiation-Induced Meningiomas

  Purpose: Majority of pediatric cancers require the irradiation of the central nervous system (CNS), and as more patients survive into adulthood from improved oncological therapy the sequelae of brain radiation are increasing in prevalence. Radiation-induced meningiomas (RIMs), one such secondary effect, demonstrate a clinically more aggressive behaviour than sporadic meningiomas (SMs). We aimed to describe the genomic mutational landscape of RIMs.Methods: We analyzed a principal cohort of 18 RIMs, with 31 RIMs overall, from patients who received childhood radiation therapy and 30 SMs, as a comparator population. We performed a multiplatform integrative genomic analysis; including methylation, whole exome and RNA sequencing. Results: RIMs exhibited a five-fold increase in copy number alterations, commonly the loss of chromosome 1p (17/18 RIMs) and 22q (17/18 RIMs), which was significantly more than observed in sporadic meningiomas. Furthermore, RNA sequencing data revealed an NF2 gene fusion event in 35.3% of RIMs. In all 6 cases, there was a complete NF2 exon spliced into a complete exon of a reciprocal gene, suggesting that the breakpoints of genomic rearrangement are intronic. All tumours with the NF2 fusion also possessed monosomy of chromosome 22q, rendering the cells with homozygous disruption of NF2. Clinically, RIMs with the NF2 fusion exhibited ill-defined borders and a tendency to develop in anatomic frontal location. Also, targeted sequencing panel confirmed that RIMs had fewer nonsynonymous NF2 mutations (6.5% vs. 30% in SM) and absence of mutations in TRAF7, SMO, KLF4, PIK3CA and AKT1, genes traditionally involved in SMs.Conclusion: Our study demonstrates that RIMs have distinct genomic drivers of oncogenesis as compared to SMs, specifically NF2 inactivation through fusion event. Radiation therapy possibly triggers genomic structural rearrangements through error-prone repair of double-stranded DNA breaks.

  Study EGAS00001002318

• 1 Intratumoral genetic heterogeneity and clonal evolution following neoadjuvant chemoradiotherapy (nCRT) in locally advanced rectal tumors.

  Purpose: Neoadjuvant chemoradiotherapy (nCRT) is the standard treatment forpatients with locally advanced rectal cancer. However, little is known about the effect of nCRT in the mutational landscape and intratumoral genetic heterogeneity (ITGH) of rectal tumors. Moreover, response to nCRT vary substantially among patients and the genetic determinants of treatment response are not fully understood. Here, we sought to determine the impact of nCRT in the mutational landscape and ITGH of rectal tumors and to identify putative drivers of therapeutic resistance. Patients and Methods: We analyzed whole-exome sequencing (WES) and clinical data from 79 primary non-treated rectal cancers obtained from The Cancer Genome Atlas (TCGA). We also compared the mutational landscape of a matched set of peripheral blood cells, irradiated tumor-adjacent colonic mucosa and pre and post-treatment tumor samples to determine the iatrogenic effects of nCRT and its impact in tissue genetic heterogeneity. Finally, we performed WES of 7 matched pre and post-treatment tumor samples to examine how nCRT affects the clonal architecture of rectal tumors. Results: We show that rectal tumors exhibit remarkable ITGH, which increases with disease progression. We also show that nCRT, per se, does not introduce novel somatic mutations, nor alters the clonal architecture of the irradiated colonic mucosa. Instead, we show that it acts as a potent selective pressure, favoring the expansion of tumor cell subpopulations more prone to resist therapy and increasing ITGH in the residual tumor. Conclusion: Our results reveal the highly heterogeneous and dynamic clonal architecture of rectal tumors undergoing nCRT. Our results also uncover putative genetic determinants of response to nCRT and a possible role for ITGH in predicting treatment outcome.

  Study EGAS00001003250

• Enhanced detection of circulating tumor DNA by fragment size analysis

  Existing methods to improve detection of circulating tumor DNA (ctDNA) have focused on sensitivity for detecting genomic alterations but have rarely considered the biological properties of plasma cell-free DNA (cfDNA). We hypothesized that differences in fragment lengths of circulating DNA could be exploited to enhance sensitivity for detecting the presence of ctDNA and for non-invasive genomic analysis of cancer. We surveyed ctDNA fragment sizes in 344 plasma samples from 200 cancer patients using low-pass whole-genome sequencing (0.4×). To establish the size distribution of mutant ctDNA, tumor-guided personalized deep sequencing was performed in 19 patients. We detected enrichment of ctDNA in fragment sizes between 90–150 bp, and developed methods for in vitro and in silico size selection of these fragments. Selecting fragments between 90–150 bp improved detection of tumor DNA, with more than 2-fold median enrichment in >95% of cases, and more than 4-fold enrichment in > 10% of cases. Analysis of size-selected cfDNA identified clinically actionable mutations and copy number alterations that were otherwise not detected. Identification of plasma samples from patients with advanced cancer was improved by predictive models integrating fragment length and copy number analysis of cfDNA, with AUC> 0.99 compared to AUC < 0.80 without fragmentation features. Increased identification of cfDNA from patients with glioma, renal, and pancreatic cancer was achieved with AUC > 0.91, compared to AUC < 0.5 without fragmentation features. Fragment size analysis and selective sequencing of specific fragment sizes can boost ctDNA detection. This could reduce the required depth of cfDNA sequencing for clinical applications, earlier diagnosis and study of tumor biology.

  Study EGAS00001003258

• Genetic characterization of B-cell prolymphocytic leukemia: a hierarchical prognostic model involving MYC and TP53 abnormalities - WXS

  B-PLL is a rare disease, whose molecular pathology is largely unknown. We report here the cytogenetic and molecular findings in a large series of 34 B-PLL. Karyotype (K) was complex (≥3 abnormalities) in 73%, and highly complex (HCK≥5) in 45%. The most frequent chromosomal aberrations were: translocation targeting the MYC gene [t(MYC)] (62%), 17p deletion including TP53 gene (38%), trisomy 18/18q (30%), 13q14 deletion (29%), trisomy 3 (24%), trisomy 12 (24%) and 8p deletion (23%). Whole-Exome Sequencing performed in 16 patients revealed recurrent mutations in TP53 (6/16, 38%), MYD88 (n=4), BCOR (n=4), MYC (n=3), SF3B1 (n=3), FAT1 (n=3), SETD2 (n=2), CHD2 (n=2), CXCR4 (n=2) and BCLAF1 (n=2). The main group of patients (21/34, 62%) had a t(MYC) associated with a higher percentage of prolymphocytes (p=0.03), CD38 expression (p<0.001), lower K complexity (p=0.0004), mutations in MYC and in genes involved in RNA metabolism and chromatin remodeling. Principal component analysis of gene expression data showed that patients with t(MYC) clustered together. A second group with MYC gain (5/34, 15%), was associated with HCK≥5 (p=0,01) and trisomy 3 (p=0,008). Altogether, 26/34 patients (76%) had a MYC activation, translocation or gain, that were mutually exclusive. We identified 3 distinct cytogenetic prognostic groups (p=0.0006): lower risk: absence of MYC activation (median not reached); intermediate risk: MYC activation without del17p (125 months); high risk: MYC activation + del17p (11 months). Our results show that cytogenetic analysis is a useful diagnostic tool in B-PLL that improves prognostic stratification.

  Study EGAS00001003275

• Genomic landscape of oral cancers (Illumina RNA-Seq)

  Human papillomavirus (HPV) is a necessary but insufficient cause of a subset of oral squamous cell carcinomas (OSCC) that is increasing markedly in frequency. To identify contributory, secondary genetic alterations in these cancers, we used comprehensive genomics methods to compare 149 HPV-positive and 335 HPV-negative OSCC tumor/normal pairs. Different behavioral risk factors underlying the two OSCC types were reflected in distinctive genomic mutational signatures. In HPV-positive OSCC, the signatures of APOBEC cytosine deaminase editing, associated with anti-viral immunity, were strongly linked to overall mutational burden. By contrast, in HPV-negative OSCC, T>C substitutions in the sequence context 5’-ATN-3’ correlated with tobacco exposure. Universal expression of HPV E6*1 and E7 oncogenes was a sine qua non of HPV-positive OSCC. Significant enrichment of somatic mutations was confirmed or newly identified in PIK3CA, KMT2D, FGFR3, FBXW7, DDX3X, PTEN, TRAF3, RB1, CYLD, RIPK4, ZNF750, EP300, CASZ1, TAF5, RBL1, IFNGR1 and NFKBIA. Of these, many affect host pathways already targeted by HPV oncoproteins, including the p53 and pRB pathways, or disrupt host defenses against viral infections, including interferon and nuclear factor-B signaling. Frequent copy number changes were associated with concordant changes in gene expression. Chr. 11q (including CCND1) and 14q (including DICER1 and AKT1) were recurrently lost in HPV-positive OSCC, in contrast to their gains in HPV-negative OSCC. High-ranking variant allele fractions implicated ZNF750, PIK3CA and EP300 mutations as candidate driver events in HPV-positive cancers. We conclude that virus-host interactions cooperatively shape the unique genetic features of these cancers, distinguishing them from their HPV-negative counterparts.

  Study EGAS00001003237

• The Immune Microenvironment, Genome–Wide Copy Number Aberrations and Survival in Mesothelioma

  IntroductionResults of recent clinical studies of immune checkpoint inhibitors in malignant pleural mesothelioma (MPM) have dampened initial enthusiasm. However, the immune environment and targets of these treatments such as program cell death protein 1 (PD-1) and its ligand PD-L1 have not been well characterised in MPM. Using a large cohort of patients, we investigated PD-L1 expression, immune infiltrates and genome-wide copy number status and correlated them to clinicopathological features.Patients and methodsTissue microarrays were constructed and stained with PD-L1(clone E1L3N, CST), CD4, CD8 and FOXP3 antibodies. PD-L1 positivity was defined as ≥5% membranous staining regardless of intensity and PD-L1+hi as ≥50%. Genomic DNA from a representative subset of 113 patients was used for genome-wide copy number analysis. The percent genome alteration (PGA) was computed as a proxy for genomic instability, and statistical analyses were used to relate copy number aberrations (CNA) to other variables.ResultsAmongst 329 patients evaluated, PD-L1+ was detected in 130/311 (41.7%), but PD-L1+hi was only seen in 30/311 (9.6%). PD-L1+ correlated with non–epithelioid histology and increased infiltration with CD4+, CD8+ and FOXP3+ lymphocytes. PD-L1+hi expression correlated with worse prognosis (HR = 2.37; 95%CI= 1.57-3.56; P < 0.001) on univariate analysis, but not in multivariate analysis. Higher PGA was associated with epithelioid histology and poorer survival (HR = 1.59; 95%CI= 1.01-2.5; P = 0.04), but not PD-L1 expression.ConclusionPD-L1 expression was associated with non-epithelioid MPM, poor clinical outcome and increased immunological infiltrates. Increased genomic instability did not correlate with PD-L1 expression, but was associated with poorer survival.

  Study EGAS00001002323

• Whole exome sequencing of cell-free DNA reveals temporo-spatial heterogeneity and identifies treatment-resistant clones in neuroblastoma

  Background: Neuroblastoma (NB) displays important clinical and genetic heterogeneity, with emergence of new mutations at tumor progression. Patients and Methods: To study clonal evolution during treatment and follow-up, an innovative method based on circulating cell-free DNA (cfDNA) analysis by whole exome sequencing (WES) paired with target sequencing was realized in sequential liquid biopsy samples of 19 NB patients. Results: WES of the primary tumor and cfDNA at diagnosis showed overlap of Single Nucleotide Variants (SNVs) and Copy Number Alterations (CNAs), with 41% and 93% of all detected alterations common to the primary NB and cfDNA. CfDNA WES at a second time point indicated a mean of 22 new SNVs for patients with progressive disease. Relapse-specific alterations included genes of the MAPK pathway and targeted the protein kinase A signaling pathway. Deep coverage target sequencing of intermediate timepoints during treatment and follow-up identified distinct subclones. For 17 seemingly relapse-specific SNVs detected by sequencing detected these alterations in minor subclones, with relapse-emerging SNVs targeting genes of neuritogenesis and cell cycle. Furthermore a persisting, resistant clone with concomitant disappearance of other clones was identified by a mutation in the ubiquitin protein ligase HERC2. Conclusion: Modelisation of mutated allele fractions in cfDNA indicated distinct patterns of clonal evolution, with either a minor, treatment-resistant clone expanding to a major clone at relapse, or minor clones collaborating towards tumor progression. Identification of treatment-resistant clones will enable development of more efficient treatment strategies.cfDNA WES at relapse but not tumor or cfDNA WES at diagnosis, deep coverage target.

  Study EGAS00001002705

• Integrative genomic analysis reveals cancer-associated mutations at diagnosis of CML in patients with high risk disease

  Genomic events associated with poor outcome in chronic myeloid leukemia (CML) are poorly understood. We performed whole exome sequencing, copy number variation and/or RNA-Seq for 65 patients to discover mutations at diagnosis and blast crisis (BC). Forty-six chronic phase patients with the extremes of outcome were studied at diagnosis. Cancer gene variants were detected in 15/27 patients (56%) with subsequent BC or poor outcome and in 3/19 optimal responders (16%), P=.007. Frequently mutated genes at diagnosis were ASXL1, IKZF1 and RUNX1. The methyltransferase SETD1B was a novel recurrently mutated gene. A novel class of variant associated with the Philadelphia translocation was detected at diagnosis in 11/46 patients (24%) comprising fusions and/or rearrangement of genes on the translocated chromosomes, with evidence of fragmentation, inversion and imperfect sequence reassembly. These were more frequent at diagnosis in patients with poor outcome: 9/27 (33%) versus 2/19 optimal responders (11%), P=.07. Thirty-nine patients were tested at BC and all had cancer gene variants, including ABL1 kinase domain mutations in 58%. However, ABL1 mutations co-occurred with other mutated cancer genes in 89% of cases, and these predated ABL1 mutations in 62% of evaluable patients. Gene fusions not associated with the Philadelphia translocation occurred in 42% of patients at BC and commonly involved fusion partners that were known cancer genes (78%). Genomic analysis revealed numerous relevant variants at diagnosis in patients with poor outcome and all patients at BC. Future refined biomarker testing of specific variants will likely provide prognostic information to facilitate a risk-adapted therapeutic approach.

  Study EGAS00001003071

• Comparison of capture-based method for transcriptome profiling of formalin-fixed paraffin embedded tumor samples

  Background: The need for fresh frozen (FF) tissue limits implementing RNA sequencing (RNA-seq) in the clinic. The majority of clinical samples are stored as formalin-fixed, paraffin-embedded (FFPE) tissues. Exome capture platforms have been developed for RNA-seq from FFPE samples. However, these methods have not been systematically compared. Methods: Transcriptomic analysis of 32 FFPE tumor samples from 11 patients was performed using three exome capture-based methods: Agilent SureSelect V6, TWIST NGS Exome, and IDT XGen Exome Research Panel. We compared these methods to TruSeq RNA-seq of fresh frozen (FF-TruSeq) tumor samples from the same patients. We assessed the recovery of clinically relevant biological features. Results: The Spearman's correlation coefficients between the global expression profiles of the three capture-based methods from FFPE and matched FF-TruSeq were high (rho = 0.72-0.9, p < 0.05). A significant correlation between the expression of key immune genes between individual capture-based methods and FF-TruSeq (rho = 0.76-0.88, p < 0.05) was observed. All exome capture-based methods reliably detected outlier expression of actionable genes, including ERBB2, MET, NTRK1, and PPARG. In urothelial cancer samples, the Agilent assay was associated with the highest molecular subtype concordance with FF-TruSeq (Cohen's k = 0.7, p < 0.01). Both Agilent and IDT detected all the clinically relevant fusions that were initially identified in FF-TruSeq. Conclusion: All exome capture-based methods had comparable performance and concordance with FF-TruSeq. By enabling the interrogation of FFPE tumor samples, our findings will enable the implementation of RNA-seq in the clinic to guide precision oncology approaches.

  Study EGAS00001005255

• Are children born after medical assisted reproduction at greater risk of having an increased de novo mutation rate?

  Introduction: De novo mutations (DNMs) play a prominent role in sporadic disorders with reduced fitness such as infertility and intellectual disability. Advanced paternal age is known to increase disease risk in offspring by increasing the number of DNMs in their genome. Less is known about the effect of assisted reproduction techniques (ART) on the number of DNMs in offspring. With the on-going trend of delayed parenthood more children are now born both from older fathers and through ART. Materials and Methods: We investigated 49 trios (mother, father and child) and 2 quartets (mother, father and 2 siblings) divided into children born after spontaneous conception (n=18); born after in vitro fertilisation (IVF) (n=17) and born after intracytoplasmic sperm injection combined with testicular sperm extraction (ICSI-TESE) (n=18). Groups further divided by paternal age, young (<35) or old (>45 years of age at conception). Whole-genome sequencing was performed twice to independently detect and validate all DNMs in children. Results: A clear paternal age effect was observed, with 70 DNMs detected on average in children born to young fathers and 94 DNMs in those born to older fathers (p = 0.001). No significant differences were observed between different methods of conception (p = 1) with paternal age affecting all methods equally. Conclusions: Paternal age, not method of conception, had a major effect on the observed number of DNMs in offspring. Given the role DNMs in disease risk, this negative result is good news for IVF and ICSI-TESE born children, if replicated in larger cohorts.

  Study EGAS00001005569

• Targeted Gene Panel for 171 PTCLs

  Mature T-cell lymphomas, including peripheral T-cell lymphoma (PTCL) and extranodal NK/T-cell lymphoma (NKTL), represent a heterogeneous group of non-Hodgkin lymphomas with dismal outcomes and limited treatment options. To determine the extent of involvement of the JAK/STAT pathway in this malignancy, we performed targeted capture sequencing of 188 genes in this pathway in 171 PTCL and NKTL cases. A total of 272 nonsynonymous somatic mutations in 101 genes were identified in 73% of the samples, including 258 single-nucleotide variants and 14 insertions or deletions. Recurrent mutations were most frequently located in STAT3 and TP53 (15%), followed by JAK3 and JAK1 (6%) and SOCS1 (4%). A high prevalence of STAT3 mutation (21%) was observed specifically in NKTL. Novel STAT3 mutations (p.D427H, E616G, p.E616K, and p.E696K) were shown to increase STAT3 phosphorylation and transcriptional activity of STAT3 in the absence of cytokine, in which p.E616K induced programmed cell death-ligand 1 (PD-L1) expression by robust binding of activated STAT3 to the PD-L1 gene promoter. Consistent with these findings, PD-L1 was overexpressed in NKTL cell lines harboring hotspot STAT3 mutations, and similar findings were observed by the overexpression of p.E616K and p.E616G in the STAT3 wild-type NKTL cell line. Conversely, STAT3 silencing and inhibition decreased PD-L1 expression in STAT3 mutant NKTL cell lines. In NKTL tumors, STAT3 activation correlated significantly with PD-L1 expression. We demonstrated that STAT3 activation confers high PD-L1 expression, which may promote tumor immune evasion. The combination of PD-1/PD-L1 antibodies and STAT3 inhibitors might be a promising therapeutic approach for NKTL, and possibly PTCL.

  Study EGAS00001002740

• The Human Phenotype Project (HPP) is a large-scale deep-phenotype prospective longitudinal and ethnically diverse cohort

  The Human Phenotype Project (HPP) is a large-scale deep-phenotype prospective longitudinal and ethnically diverse cohort and biobank that we established. To date, approximately 28,000 participants have enrolled in the study, with over 13,000 having completed their initial visit. The project is aimed at identifying novel molecular signatures with diagnostic, prognostic and therapeutic value, and at developing AI-based predictive models for disease onset and progression. A unique feature of the HPP is its deep and longitudinal profiling that includes medical history, lifestyle and nutritional habits, vital signs, anthropometrics, blood tests, continuous glucose and sleep monitoring, imaging modalities and molecular profiling of the transcriptome, genetics, gut, vaginal and oral microbiome, metabolome and immune system. For several of these modalities, the HPP is the first or largest cohort of its kind. Our analyses of this data provide novel insights into the variation of phenotypes with age and ethnicity, highlighting the importance of having ethnically diverse cohorts and the need for establishing age-dependent norms. We demonstrate how the HPP can be used to unravel personalized molecular signatures of disease by comparison to patient-specific and disease-specific matched healthy controls. Leveraging the extensive dietary and lifestyle HPP data, we systematically identify associations between lifestyle factors and health outcomes. Finally, we present a multi-modal foundation AI model that is trained using self-supervised learning on diet and CGM data and outperforms existing methods in predicting future onset of disease. We extend this AI framework to integrate all data modalities of each subject as a continuous sequence of diverse medical events, creating a digital twin that can simulate interventions and predict health trajectories and outcomes.

  Study EGAS00001008040

• Integrated Genomic, Epigenetic, and Expression Analyses of Ovarian Cancer Cell Lines

  Epithelial ovarian carcinoma of all subtypes has one of the worst outcomes of any human cancer. To improve our understanding of underlying biology and therapeutic possibilities in this disease, we performed whole genome sequencing, methylation, and expression analyses of a comprehensive set of 47 ovarian cancer cell lines representing different histologic tumor subtypes. Given the challenges of genomic analyses in tumors without matched normal samples, we developed novel approaches for detection of somatic alterations and used these to identify sequence, copy number, and rearrangement changes and integrated these with genome-wide epigenetic and expression alterations. Ovarian cancer cell lines had molecular changes in genes involved in cell cycle, PI3K, RAS, BRCA, chromatin regulating, and p53 pathways, and displayed global mutation and methylation changes similar to primary ovarian carcinoma. We identified alterations not previously implicated in ovarian cancer including amplification or overexpression of ASXL1 and H3F3B, deletion or underexpression of CDC73 and members of the TGF beta receptor pathway, including TGFBR2, SMAD3 and SMAD4, in-frame rearrangements of YAP1-MAML2 and IKZF2-ERBB4, and fusions involving MET, NF1, FBXW7 and CCND1. Integrating cell line alterations with dose-response data using three targeted therapies revealed novel molecular dependencies, including increased sensitivity of tumors with PIK3CA and PPP2R1A alterations to PI3K inhibitor GNE-493, those with MYC amplifications to PARP inhibitor BMN673, and SMAD4 alterations to MEK inhibitor MEK162. Genome-wide analyses revealed intrachromosomal rearrangements as an improved measure of sensitivity to PARP inhibition compared to the homologous recombination defect (HRD) score. This study provides a comprehensive resource of molecular information for ovarian cancer cell line models and a pharmacogenomic platform for developing rational cancer therapeutic strategies.

  Study EGAS00001002998

• Whole exome sequencing of small cell neuroendocrine cancer of the cervix

  In order to improve treatment selection for high grade neuroendocrine carcinomas of the cervix (NECC), we performed a comparative genomic analysis between this rare tumor type and other cervical cancer types, as well as extra-cervical neuroendocrine small cell carcinomas of the lung and bladder. We performed whole exome sequencing on fresh-frozen tissue from 15 NECCs and matched normal tissue. We then identified mutations and copy number variants using standard analysis pipelines. Published mutation tables from cervical cancers and extra-cervical small cell carcinomas were used for comparative analysis. Descriptive statistical methods were used and a two-sided threshold of P < .05 was used for significance. In the NECC cohort, we detected a median of 1.7 somatic mutations per megabase (range 1.0-20.9). PIK3CA p.E545K mutations were the most frequency observed oncogenic mutation (4/15 tumors, 27%). Activating MAPK pathway mutations in KRAS (p.G12D) and GNAS (p.R201C) co-occurred in two tumors (13%). In total we identified PI3-kinase or MAPK pathway activating mutations in 67% of NECC. When compared to NECC, lung and bladder small cell carcinomas exhibited a statistically significant higher rate of coding mutations (P < .001 for lung; P = .001 for bladder). Mutation of TP53 was uncommon in NECC (13%) and was more frequent in both lung (103 of 110 tumors [94%], P < .001) and bladder (18 of 19 tumors [95%], P < .001) small cell carcinoma. These comparative genomics data suggest that NECC may be genetically more similar to common cervical cancer subtypes than to extra-cervical small cell neuroendocrine carcinomas of the lung and bladder. These results may have implications for the selection of cytotoxic and targeted therapy regimens for this rare disease.

  Study EGAS00001003142

• RNA-seq as a tool for evaluating human embryo competence

  The majority of embryos that are created through IVF do not implant. One possibility for this inefficiency is an incomplete understanding of the molecular factors affecting embryo competence. Currently, the process of selecting an embryo for uterine transfer utilizes an ad-hoc combination of morphological criteria, the kinetics of development, and genetic testing for aneuploidy. However, no single criterion can ensure selection of a viable embryo. In contrast, RNA-sequencing of embryos could yield highly dimensional data, which may provide additional insight and illuminate the discrepancies among current selection criteria. Indeed, recent advances enabling the production of RNA-sequencing (RNA-seq) libraries from single cells have facilitated the application of this technique to the study of some transcriptional events in early human development. However, these studies have not assessed the quality of their constituent embryos relative to commonly used embryological criteria. Here, we perform proof-of-principle advancement to clinical selection procedures by generating high quality RNA-seq libraries from a trophectoderm biopsy as well as the remaining whole embryo. We combine state-of-the-art embryological methods with low-input RNA-seq to develop the first transcriptome-wide approach for use in future predictive embryology studies. Specifically, we demonstrate the capacity of RNA-seq as a promising tool in preimplantation screening by showing that biopsies of an embryo can capture valuable information content available in the whole embryo from which they are derived. Furthermore, we show that this technique can be used to generate a RNA-based digital karyotype, and to develop a foundational dataset for identifying candidate competence-associated genes. Together, these data establish the foundation for a future RNA-based diagnostic in IVF.

  Study EGAS00001003667

• Altered oligodendrocyte heterogeneity in Multiple sclerosis revealed by single nuclei RNA sequencing

  This study is currently hosted by the European Nucleotide Archive. To access the data contained within the Study please follow the link below: https://www.ebi.ac.uk/ena/browser/view/PRJEB39323 Oligodendrocyte (OL) pathology is increasingly implicated in neurodegenerative diseases, as they are involved in metabolic support of axons and functional cross-talk with other brain cells. Rodent OLs are heterogeneous, with developmental and biological differences, but the extent of heterogeneity in the normal human brain and its contribution to any changes to disease remains unknown. Here we performed single nuclei RNA-sequencing (snRNA-seq) from white matter (WM) areas of post mortem human brain both in control (Ctr) and multiple sclerosis (MS) patients. We identified several sub-clusters of oligodendroglia in the Ctr human WM, some similar to those in mouse, and defined new markers for these cell states. Strikingly, some of these sub-clusters were under-represented in MS tissue, while others were more prevalent than in controls. We found a lack of OL precursor cells (OPCs) and an OL subcluster in an intermediate stage of differentiation in MS lesions and in normal appearing white matter (NAWM), suggesting either depletion by the disease or by a regenerative response. The differences in mature OL sub-clusters indicate different functional states of OLs in MS tissue and, as this is similar in NAWM to lesions, that MS is a more diffuse brain disease than the focal demyelinating lesions suggest. We were also able to identify new putative markers of different MS lesion subtypes. Our findings of an altered heterogeneity of oligodendroglia in MS may have an important contribution to our understanding of disease progression and may alter therapeutic approaches to MS.

  Study EGAS00001003412

• A molecular cell atlas of the human lung from single cell RNA sequencing

  Although single cell RNA sequencing studies have begun providing compendia of cell expression profiles, it has proven more difficult to systematically identify and localize all molecular cell types in individual organs to create a full molecular cell atlas. Here we describe droplet- and plate-based single cell RNA sequencing applied to ~75,000 human lung and blood cells, combined with a multi-pronged cell annotation approach, which have allowed us to define the gene expression profiles and anatomical locations of 58 cell populations in the human lung, including 41 of 45 previously known cell types or subtypes and 14 new ones. This comprehensive molecular atlas elucidates the biochemical functions of lung cell types and the cell-selective transcription factors and optimal markers for making and monitoring them; defines the cell targets of circulating hormones and predicts local signaling interactions including sources and targets of chemokines in immune cell trafficking and expression changes on lung homing; and identifies the cell types directly affected by lung disease genes and respiratory viruses. Comparison to mouse identified 17 molecular types that appear to have been gained or lost during lung evolution and others whose expression profiles have been substantially altered, revealing extensive plasticity of cell types and cell-type-specific gene expression during organ evolution including expression switches between cell types. This atlas provides the molecular foundation for investigating how lung cell identities, functions, and interactions are achieved in development and tissue engineering and altered in disease and evolution.

  Study EGAS00001004344

• We evaluate the potential for routine WGS using ONT by sequencing the well-characterised reference sample NA12878 and the genome of an individual with ataxia-pancytopenia syndrome accompanied by severe immune dysregulation.

  Whole-genome sequencing (WGS) is becoming widely used in clinical medicine in diagnostic contexts and to inform treatment choice. While current sequencing technologies have been extremely successful, their reliance on short read lengths necessarily involves some limitations: accurately assaying certain genomic regions and classes of variation can be problematic. Recent advances in throughput and cost have made WGS using the Oxford Nanopore Technologies (ONT) MinION long-read single-molecule sequencer a potential solution to these challenges. Here we evaluate the potential for routine WGS using ONT by sequencing the well-characterised reference sample NA12878 and the genome of an individual with ataxia-pancytopenia syndrome accompanied by severe immune dysregulation, to 82× and 30× respectively. For NA12878, we evaluated single-nucleotide variant (SNV) calls based on data from multiple base-calling algorithms. We demonstrate that phasing metrics from a novel, reference panel-free, long-read-based method can improve variant-calling performance from otherwise modest levels, resulting in a false discovery rate of 7.1% and false negative rate of 8.5%; remaining errors are concentrated near homopolymers and in regions of reduced sequencing coverage. In the clinical sample, we are able to identify and directly phase two non-synonymous de novo variants in SAMD9L, (OMIM #159550) inferring that they both lie on a common paternal haplotype by overlap with parental genotypes at nearby common variants. This work demonstrates that methodological innovation can substantially reduce variant-calling error rates in ONT data, and that with on-going improvements in throughput, base-calling and dedicated long-read-based SNV-calling methodology, ONT offers promise as an option for clinical WGS.

  Study EGAS00001003469

• Persistence of circulating tumor DNA in breast cancer patients during neoadjuvant treatment is a significant predictor of poor tumor response

  Purpose:Accurate assessment of response during neoadjuvant systemic treatment (NST) poses a major clinical challenge. Therefore, a minimally-invasive assessment of tumor response based on cell-free circulating tumor DNA (ctDNA) may be beneficial to guide treatment decisions.Experimental Design: We profiled 93 genes in tissue from 193 early breast cancer patients. Patient-specific assays were designed to track ctDNA during NST from 145 patients with available plasma. ctDNA presence and levels were correlated with complete pathological response (pCR) and residual cancer burden (RCB) as well as with clinicopathologic characteristics of the tumor to identify potential proxies for ctDNA release. Moreover, we tested the predictive value of driver events identified in the tumor.Results:At baseline, ctDNA could be detected in 63/145 (43.4%) patients and persisted in 25/63 (39.7%) patients at mid-therapy (MT) and 15/63 (23.8%) patients at the end of treatment. ctDNA detection at MT was significantly associated with higher RCB (OR 0.062, 95% CI 0.01-0.48, P=0.0077). Out of 31 patients with detectable ctDNA at MT, 30 patients (96.8%) were non-responders (RCB II, n=8; RCB III, n=22) and only one patient responded to the treatment (RCB I). Considering all 145 patients with baseline (BL) plasma, none of the patients with RCB 0 and only 6.7% of patients with RCB I had ctDNA detectable at MT, while 30.6% and 29.6% of patients with RCB II/III, respectively, had a positive ctDNA result. Conclusion:Overall, our results demonstrate that the detection and persistence of ctDNA at mid-therapy may have the potential to negatively predict response to neoadjuvant treatment and identify patients who will not achieve pCR or be classified with RCB II/III.

  Study EGAS00001005798

• Post_Mortem_Tissue_COVID19_RNA

  Single cell analysis of post mortem tissue samples from SARS-CoV2-infected patients We aim to directly identify the human cell types infected by SARS-CoV-2 and measure the cellular response to COVID19 infection across 20 different tissues from infected patient autopsies. We have examined the expression pattern of viral entry receptors across healthy human tissues to predict several candidate target cell types across the airway and heart. In addition, high prevalence of cardiac failure and abnormal renal function in COVID19 patients implicates heart and kidney involvement, but the pathogenesis of organ specific damage - whether via a direct cytopathic mechanism or an indirect inflammatory response - remains unknown . Currently, we lack confirmation of target cell types and cellular processes in infected tissues as autopsies are discouraged in most countries due to health and safety risks. Our collaborators Drs Michael Osborn and Brian Hanley (Imperial College) have outlined guidelines to perform post-mortem in COVID19 patients (Hanley et al., 2020) and have established a programme of autopsies for research to be performed in a high-risk facility at Westminster Public Mortuary. Here, we propose to identify infected cell types and aberrant molecular pathologies in this precious tissue resource using single cell and spatial genomics. We will prioritise three organ systems: the human airway, heart and the kidney. We will directly examine the cellular identities of SARS-CoV-2 infected cell types and identify the cellular responses to infection across these organs. This fundamental knowledge will help guide future treatment choices for COVID19. This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/

  Study EGAS00001004442

• A reference map of potential determinants for the human serum metabolome

  The collection of metabolites circulating in the human blood, termed the serum metabolome, contains a plethora of biomarkers and causative agents. Although the origin of specific compounds is known, we have a poor understanding of the key determinants of most metabolites. Here, we measured the levels of 1251 circulating metabolites in serum samples from a healthy human cohort of 491 individuals, and devised machine learning algorithms to predict their levels in held-out subjects based on a comprehensive profile consisting of host genetics, gut microbiome, clinical parameters, diet, lifestyle, anthropometric measurements and medication data. Notably, we obtained statistically significant predictions for over 76% of the profiled metabolites. Despite using the strict out-of-sample prediction metric, which is a lower bound for the explained variance, diet and microbiome each explained hundreds of metabolites, with over 50% of the variance explained in some metabolites. We further validated the robustness of the microbiome related associations by showing a high replication rate in two geographically independent cohorts that were not available to us when developing the algorithms. We also demonstrate that some of these interactions are causal, as some metabolites we predicted to be positively associated with bread increased in level following a randomized clinical trial of bread intervention. Microbiome-explained metabolites were enriched with unnamed metabolites, and we devised an algorithm that accurately predicts their biological pathway, finding that they mainly include food components, aromatic amino acids and secondary bile acid derivatives. Overall, our results unravel potential determinants of over 800 metabolites, paving the way towards mechanistic understanding of alterations in metabolites under different conditions and to designing interventions for manipulating circulating metabolite levels.

  Study EGAS00001004512

• Tumor-derived cell lines as pharmacogenomic models to predict therapeutic vulnerabilities in hepatocellular carcinoma

  Hepatocellular carcinoma (HCC) is a heterogeneous aggressive malignancy with low efficacy of current therapies at advanced stages. We integrated molecular and pharmacological profiling of a large panel of liver cancer cell lines (LCCL) to assess their clinical relevance as HCC preclinical models and identify new effective therapies and biomarkers of response. We performed multi-omic analysis including whole-exome, RNA and miRNA sequencing and quantification of 126 proteins in a series 34 LCCL coupled with a screening of 29 anti-cancer agents. Molecular profiles of LCCL and primary HCC were compared and we searched for molecular features associated with drug response. Our panel of LCCL faithfully recapitulated the most aggressive molecular "proliferation class" of HCC. Genomic alterations of the RAS-MAPK pathway correlated with trametinib sensitivity that was more potent than FGFR4 inhibitors in tumor cells depending on the FGF19/FGFR4 pathway. Moreover, we showed that FGF19 amplification was only functional in tumor cells that retained a hepatocyte differentiation program with expression of FGFR4, KLB and NR1H4 highlighting the complex interplay between the genomic and transcriptomic context in drug response. We also identified drugs/combination that could target efficiently the transcriptomic "progenitor subclass of HCC", as well as specific genetic contexts such as inactivating mutations in TSC1/TSC2 and TP53 associated with higher sensitivity to the mTOR inhibitor rapamycin and the AURKA inhibitor alisertib, respectively. LCCL represent relevant preclinical models for drug-biomarker discovery in HCC and enabled to identify molecular contexts linked with particular therapeutic vulnerabilities that may be useful to stratify patients in future clinical trials

  Study EGAS00001003536

• Survival Benefit and Genetic Profile of Pemetrexed as Initial Chemotherapy in Selected Chinese Patients with Advanced Lung Adenocarcinoma

  We retrospectively collected patients (n=1047) enrolled in the Chinese Patient Assistance Program from multiple centers who received pemetrexed alone or combined with platinum as initial chemotherapy and continued pemetrexed maintenance therapy for advanced lung adenocarcinoma from November 2014 to June 2017. The outcomes were duration of treatment (DOT) and overall survival (OS). Clinical features were analyzed for their influence on the treatment effect and prognosis. Next-generation sequencing (NGS) was performed to identify genetic biomarkers associated with the efficacy of pemetrexed. The median DOT was 9.1 months (95% CI: 8.5-9.8), and the median OS was 26.2 months (95% CI: 24.2-28.1). OS was positively correlated with DOT (r=0.403, P<0.001). Multivariable analysis showed that smoking status and Eastern Cooperative Oncology Group (ECOG) performance status (PS) were independently associated with DOT; smoking status, ECOG PS, targeted therapy, and EGFR/ALK/ROS1 status were independently associated with OS. NGS in 22 patients with available samples showed genes with high mutation rates were: TP53 (54.5%), EGFR (50.0%), MYC (18.2%), PIK3CA (13.6%). When grouped based on progression-free survival (PFS) reported in the PARAMOUNT study, the DOT >6.9 months set was associated with PIK3CA, ALK, BRINP3, CDKN2A, CSMD3, EPHA3, KRAS, and RB1 mutations, while ERBB2 mutation was observed only in the DOT ≤6.9 months set. This study shows that initial chemotherapy with pemetrexed is an effective regimen for advanced lung adenocarcinoma in selected Chinese patients. There is no specific genetic profile predicting the benefit of pemetrexed found by NGS. Biomarkers predicting the efficacy of pemetrexed need further exploration.

  Study EGAS00001004546

• A Multifactorial Tumor and Immune Cell Profile Determines Response to Immune Checkpoint blockade in Melanoma

  The complex crosstalk between tumor and immune cells during immune checkpoint blockade mandates the development of integrated models that interpret the anti-tumor immune response and predict clinical outcome. We have integrated genome-wide sequence and structural alterations with pre and on-therapy transcriptomic and T cell repertoire features for a cohort of 64 immunotherapy-naïve melanomas treated with anti-PD1 monotherapy or combined anti-PD1 and anti-CTLA4 therapy. While tumor mutation burden (TMB) was associated with improved response to therapy, expressed mutation burden was superior to TMB in predicting outcome. An increased pre-existing T cell density differentiated responding from non-responding tumors independent of therapy. Importantly, T cell repertoire reshaping determined by T cell receptor (TCR) clonotypic regressions and expansions reflected response to therapy such that a more dynamic repertoire predicted improved clinical outcome. Through whole-transcriptome analyses, we discovered differential abundance of B cell subsets in responding tumors, which highlights the importance of the interplay between pre-existing T and B cell immunity in shaping therapeutic response. Tumor temporal trajectories during therapy revealed distinct patterns of molecular response related to expressed mutation elimination or retention that accurately interpreted clinical response. High-dimensional genomic, transcriptomic and immune repertoire data were integrated by both a machine learning and a censored regression approach, resulting in a harmonized multi-modal predictor of response to immune checkpoint blockade. B cell abundance, expressed mutation load and tumor aneuploidy were combined to identify patients at high risk for recurrence, such that high risk patients had a significantly shorter progression-free survival (HR=9.18, 95% CI: 3.14-26.85, p=3.4e-06) especially in the anti-PD1/anti-CTLA4 group (HR=20.34, 95% CI: 4.08-101.32, p=1.68e-06).

  Study EGAS00001004548

• Clinical Outcomes by Tumor Mutational Burden and Inflammatory Gene Expression With Combined Nivolumab and Ipilimumab or Monotherapy in Advanced Melanoma

  Purpose Exploratory analyses of CheckMate 066 and 067 trials were conducted to investigate associations of tumor mutational burden (TMB), a 4-gene inflammatory gene expression signature, and BRAF mutation status with tumor response, progression-free survival (PFS), and overall survival (OS) in patients with advanced melanoma. Patients and Methods Patients with known programmed death ligand 1 (PD-L1) expression and BRAF mutation status received nivolumab (NIVO) or dacarbazine in CheckMate 066 and either NIVO, ipilimumab (IPI), or NIVO+IPI in CheckMate 067. Whole exome sequencing and RNA sequencing were used to determine TMB and inflammatory gene expression signature scores, respectively. These biomarkers were evaluated in terms of their association with PFS and OS. Results In the NIVO, NIVO+IPI, and IPI arms of CheckMate 067, longer survival was associated with high (> median) versus low (≤ median) TMB with hazard ratios (HRs) (95% confidence interval [CI]) for PFS of 0.45 (0.30–0.65), 0.55 (0.38–0.81), and 0.60 (0.43–0.82), and for OS of 0.46 (0.30–0.71), 0.53 (0.34–0.82), and 0.52 (0.36–0.74), respectively. For NIVO-treated patients, these results were confirmed in CheckMate 066. A survival benefit was observed with high TMB and absence of BRAF mutation. Survival was associated with high versus low inflammatory signature scores with HRs (95% CI) for PFS of 0.56 (0.34–0.94), 0.40 (0.23–0.72), and 0.43 (0.27–0.70), and for OS of 0.37 (0.20–0.66), 0.38 (0.19–0.74), and 0.46 (0.27–0.79), in the NIVO, NIVO+IPI, and IPI arms, respectively. Weak correlations were observed between PD-L1, TMB, and the inflammatory signature. Conclusions Combined assessment of TMB, inflammatory gene expression signature, and BRAF mutation status may be predictive for response to immunotherapy in advanced melanoma.

  Study EGAS00001004564

• Clinical Outcomes by Tumor Mutational Burden and Inflammatory Gene Expression With Combined Nivolumab and Ipilimumab or Monotherapy in Advanced Melanoma - CM066-WES

  Purpose Exploratory analyses of CheckMate 066 and 067 trials were conducted to investigate associations of tumor mutational burden (TMB), a 4-gene inflammatory gene expression signature, and BRAF mutation status with tumor response, progression-free survival (PFS), and overall survival (OS) in patients with advanced melanoma. Patients and Methods Patients with known programmed death ligand 1 (PD-L1) expression and BRAF mutation status received nivolumab (NIVO) or dacarbazine in CheckMate 066 and either NIVO, ipilimumab (IPI), or NIVO+IPI in CheckMate 067. Whole exome sequencing and RNA sequencing were used to determine TMB and inflammatory gene expression signature scores, respectively. These biomarkers were evaluated in terms of their association with PFS and OS. Results In the NIVO, NIVO+IPI, and IPI arms of CheckMate 067, longer survival was associated with high (> median) versus low (≤ median) TMB with hazard ratios (HRs) (95% confidence interval [CI]) for PFS of 0.45 (0.30–0.65), 0.55 (0.38–0.81), and 0.60 (0.43–0.82), and for OS of 0.46 (0.30–0.71), 0.53 (0.34–0.82), and 0.52 (0.36–0.74), respectively. For NIVO-treated patients, these results were confirmed in CheckMate 066. A survival benefit was observed with high TMB and absence of BRAF mutation. Survival was associated with high versus low inflammatory signature scores with HRs (95% CI) for PFS of 0.56 (0.34–0.94), 0.40 (0.23–0.72), and 0.43 (0.27–0.70), and for OS of 0.37 (0.20–0.66), 0.38 (0.19–0.74), and 0.46 (0.27–0.79), in the NIVO, NIVO+IPI, and IPI arms, respectively. Weak correlations were observed between PD-L1, TMB, and the inflammatory signature. Conclusions Combined assessment of TMB, inflammatory gene expression signature, and BRAF mutation status may be predictive for response to immunotherapy in advanced melanoma.

  Study EGAS00001004567

• Mutational Landscape and Tumor Burden Assessed by Cell-Free DNA in Diffuse Large B-Cell Lymphoma: a Population-based Study

  Purpose: We analyzed the utility of cfDNA in a prospective population-based cohort to determine the mutational profile, assess tumor burden, and estimate its impact in response rate and outcome in patients with DLBCL. Experimental Design: One-hundred patients were diagnosed with DLBCL during the study period. Mutational status of 112 genes was studied in cfDNA by targeted next-generation sequencing. Paired formalin-fixed paraffin-embedded samples and volumetric PET/CT were assessed when available. Results: Appropriate cfDNA to perform the analyses was obtained in 79/100 cases. At least one mutation could be detected in 69/79 cases (87%). The sensitivity of cfDNA to detect the mutations was 68% (95% CI: 56.2-78.7). The mutational landscape found in cfDNA samples was highly consistent with that shown in the tissue and allowed genetic classification in 43% of the cases. A higher amount of ctDNA significantly correlated with clinical parameters related to tumor burden (elevated LDH and β2-microglobulin serum levels, advanced stage, and high-risk IPI) and total metabolic tumor volume assessed by PET/CT. In patients treated with curative intent, high ctDNA levels (>2.5 log hGE/mL) were associated with lower complete response (65% vs. 96%, P<0.004), shorter progression-free survival (65% vs. 85%, P=0.038) and overall survival (73% vs. 100%, P=0.007) at 2 years, although it did not maintain prognostic value in multivariate analyses. Conclusions: In a population-based prospective DLBCL series, cfDNA resulted an alternative source to estimate tumor burden and to determine the tumor mutational profile and genetic classification, which have prognostic implications and may contribute to a future tailored

  Study EGAS00001004733

• Circulating Tumor DNA in Checkpoint Inhibitor treated Lung Cancer

  BackgroundImmune checkpoint inhibitors (ICI) are increasingly used to treat non-small cell lung cancer (NSCLC) and there is a critical need for novel predictive/prognostic biomarkers. We evaluated whether dynamics of circulating tumor (ctDNA) can identify NSCLC patients who benefit from single agent ICI.MethodsCirculating cell-free DNA (cfDNA) from 177 patients was analyzed prior ICI and at first treatment evaluation using a Commercial 77-gene panel. We performed de novo mutation calling and corrected for variants related clonal hematopoiesis. ctDNA changes were evaluated as predictors of tumor response estimated with PET/diagnostic CT and survival.ResultsMutations were detected in 159 of 174 (91.4%) analyzable patients. Most responders presented decreasing ctDNA levels, while non-responders had mixed molecular responses. Patients with a >50% molecular ctDNA decrease or no detectable ctDNA demonstrated improved PFS (5.9 versus 1.8 months; HR 0.49, 95%CI, 0.35-0.67, P < .001), and OS (18.9 versus 9.9 months; HR 0.41, 95%CI, 0.30-0.57, P < .001) compared to patients not achieving these endpoints. Clearance of mutations (and no detected ctDNA) had a similar effect (PFS HR 0.45, 95%CI, 0.32-0.61, P < .001; OS 0.37,95%CI, 0.27-0.52, P < .001) compared to patients having stable or additional mutations. After adjusting for clinical factors both >50% ctDNA decrease (HR 0.4, 95% CI, 0.3-0.6, P < .001) and clearance of mutations (HR 0.42, 95% CI, 0.3-0.6, P < .001) remained independent predictors for PFS and OS.ConclusionsctDNA testing reveals significant prognostic albeit not predictive information in ICI treated advanced NSCLC. This relevant information can be obtained with a commercial easy-to-implement ctDNA test that will facilitate extensive testing.

  Study EGAS00001004847

• 18 Whole Exome Sequencing for Radiation Induced-Meningiomas

  Purpose: Majority of pediatric cancers require the irradiation of the central nervous system (CNS), and as more patients survive into adulthood from improved oncological therapy the sequelae of brain radiation are increasing in prevalence. Radiation-induced meningiomas (RIMs), one such secondary effect, demonstrate a clinically more aggressive behaviour than sporadic meningiomas (SMs). We aimed to describe the genomic mutational landscape of RIMs.Methods: We analyzed a principal cohort of 18 RIMs, with 31 RIMs overall, from patients who received childhood radiation therapy and 30 SMs, as a comparator population. We performed a multiplatform integrative genomic analysis; including methylation, whole exome and RNA sequencing. Results: RIMs exhibited a five-fold increase in copy number alterations, commonly the loss of chromosome 1p (17/18 RIMs) and 22q (17/18 RIMs), which was significantly more than observed in sporadic meningiomas. Furthermore, RNA sequencing data revealed an NF2 gene fusion event in 35.3% of RIMs. In all 6 cases, there was a complete NF2 exon spliced into a complete exon of a reciprocal gene, suggesting that the breakpoints of genomic rearrangement are intronic. All tumours with the NF2 fusion also possessed monosomy of chromosome 22q, rendering the cells with homozygous disruption of NF2. Clinically, RIMs with the NF2 fusion exhibited ill-defined borders and a tendency to develop in anatomic frontal location. Also, targeted sequencing panel confirmed that RIMs had fewer nonsynonymous NF2 mutations (6.5% vs. 30% in SM) and absence of mutations in TRAF7, SMO, KLF4, PIK3CA and AKT1, genes traditionally involved in SMs.Conclusion: Our study demonstrates that RIMs have distinct genomic drivers of oncogenesis as compared to SMs, specifically NF2 inactivation through fusion event. Radiation therapy possibly triggers genomic structural rearrangements through error-prone repair of double-stranded DNA breaks.

  Study EGAS00001002317

• Breast cancer women lack normal lifelong immune response after full-term pregnancies

  Immunological tolerance during pregnancy is necessary for women’s bodies to accept the semi-allogenic fetus. Interestingly, “immune evasion”, which consists of some similar immunological changes, was recently proposed as a hallmark of cancer. Looking at these results, it seems paradoxical that increasing number of full-term pregnancies (FTP) reduces the risk of breast cancer in women. Our aim was to use a systems epidemiology approach to investigate this paradox by modelling breast cancer incidence as a function of parity in the Norwegian Women and Cancer (NOWAC) cohort, then testing the mathematical model using gene expression from blood in a subcohort of these women: the NOWAC postgenome cohort. Lastly, we undertook a gene set enrichment analysis for immunological gene sets. A linear trend nicely fitted the dataset showing an 8% decrease in risk for each FTP, independent of stratification on other risk factors and lasting for decades after a woman’s last FTP. However, when we looked at gene expression we found that hundreds of genes showed linear trends in cancer-free controls, but this was not the case in /any/ of the samples from breast cancer cases. In the gene set enrichment analysis, samples consisted of human blood or cord blood and mouse bone marrow or spleen. We found more adaptive gene sets than expected in humans and more innate gene sets than expected in mice. None of the gene sets provided information about the long-term effects of immunological tolerance during pregnancy, but the lack of relevant reductionist experiments hampers our interpretation. We suggest that, after resolution at birth, the immunological tolerance induced by pregnancy augments immune response to cancer development.

  Study EGAS00001002616

• DSRCT RNA genomic sequencing

  Purpose: To further understand the molecular pathogenesis of desmoplastic small round cell tumor (DSRCT), a rare but often fatal malignancy occurring primarily in young males, we used next-generation RNA sequencing to investigate the gene expression profiles intrinsic to this disease.Experimental Design: RNA from DSRCT tumor samples were obtained from the Cooperative Human Tissue Network biorepository. RNA sequencing was completed using the Illumina HiSeq 2000 system. Presence of the EWSR1-WT1 fusion was verified using both RNA-sequencing and PCR. Validation studies, including WT1 ChIP-seq, RNA sequencing, and EWS-WT1 knockdown experiments, were performed on an established DSRCT cell line, JN-DSRCT-1. Protein expression of clinically relevant genes determined by RNA sequencing was confirmed using immunohistochemistry. A panel of immune signature genes was also evaluated to identify possible immune therapeutic targets. Results: We analyzed the RNA from 14 tumor samples, of which two were confirmed negative for the diagnostic EWSR1-WT1 translocation. Principal component analysis of the sequencing data indicated that fusion-negative and fusion-positive tumors are transcriptionally distinct entities. Overlapping the significant ChIP-seq peaks with genes found to be highly expressed in the RNA sequencing data revealed that IGF2 and FGFR4 were both highly expressed, and targets of the EWS-WT1 fusion gene. This finding has potential clinical relevance for targeted therapeutic intervention. In addition, we identified the immune checkpoints CD200 and CD276 as potentially targetable genes whose expression is independent of the EWS-WT1 fusion gene in vitro.Conclusions: Despite the transcriptional complexity of this disease, we were able to identify IGF2, FGFR4, CD200, and CD276 as potential therapeutic targets for patients suffering from DSRCT.

  Study EGAS00001002770

• Accurate detection and classification of pediatric sarcomas based on cell-free DNA fragmentation patterns

  Sequencing of cell-free DNA in the blood of cancer patients (“liquid biopsy”) provides attractive opportunities not only for early diagnosis, but also for minimally invasive monitoring of treatment response and disease courses. To unlock liquid biopsy analysis for pediatric tumors with few genetic aberrations, we developed an integrated genetic/epigenetic analysis method and applied it to 241 deep whole-genome sequencing profiles of 95 patients with Ewing sarcoma and 31 patients with other pediatric sarcomas. We achieved sensitive detection and classification of circulating tumor DNA in peripheral blood independent of any genetic alterations. We evaluated different metrics for cell-free DNA fragmentation analysis and developed LIQUORICE, a bioinformatic tool for detecting circulating tumor DNA based on tumor-specific chromatin structure. Using machine learning methods, we combined several fragmentation-based metrics into an integrated approach for liquid biopsy analysis tailored to cancers with low mutation rates but widespread epigenetic deregulation. Clinical associations highlighted the potential value of cfDNA fragmentation patterns as prognostic biomarkers in Ewing sarcoma. Additionally, we performed low coverage whole-genome-sequencing on 43 tumor biopsy samples from patients with Ewing sarcoma, in order to compare copy number aberrations detected in cell-free DNA and biopsy samples of the same patients. For validation of the epigenetic signatures inferred from cell-free DNA, we further performed reduced representation bisulfite sequencing (RRBS) on 38 matched biopsy samples from patients with Ewing sarcoma. In summary, our study provides a comprehensive analysis of circulating tumor DNA beyond recurrent somatic mutations, and it renders the benefits of liquid biopsy more readily accessible for childhood cancers.

  Study EGAS00001005127

• Genetic diversity and continuity of the population of the UAE

  With high consanguinity rates on the Arabian Peninsula, it would not have been unexpected if the population of the UAE was shown to be relatively homogenous. However, this study of 1,000 UAE nationals, provided a contrasting perspective, one of a relatively heterogeneous population. Located at the apex of Europe, Asia and Africa, the observed diversity could be explained by a plethora of migration patterns since the first Out-of-Africa movement. A strategy to explore the extent of genetic variation of the population of the UAE is presented. The first step involved a comprehensive population stratification study that was instructive for subsequent Whole Genome Sequencing (WGS) of suitable representatives (which is described elsewhere). When this UAE data was compared to previous smaller studies from the region, the findings were consistent with a population that is a diverse and admixed group of people. However, rather than sharp and distinctive clusters, a continuous genotype distribution was observed. The analysis showed a continuous gradient of ancestral populations, suggesting that admixture on the south eastern tip of Arabian Peninsula occurred gradually. When visualized using a unique technique that combined admixture ratios and Principal Component Analysis (PCA), unappreciated diversity was revealed while mitigating projection bias of conventional PCA. The analytical strategy used here highlights the complementary nature of data from genotype array and WGS for anthropological studies. Specifically, genotype array data was instructive to select representative subjects for WGS. Furthermore, from the 2.3 million allele frequencies obtained from genotype arrays, we identified 46,481 loci with allele frequencies that were significantly different with respect to other world populations. This comparison of allele frequencies facilitates variant prioritization in common diseases. In addition, these loci bear great potential as biomarkers in anthropological and forensic studies.

  Study EGAS00001004362

• Matching of actionable mutations with therapies in cancer patients: comparison of three commercial decision support platforms

  Purpose: In oncology, precision medicine centers on tailoring the most appropriate therapy to a patient according to their tumor molecular profile. We addressed the question whether there is a consensus in the classification of actionability of somatic variants and their alignment to treatment recommendations. To this end, we evaluated three commercial clinical decision support tools to assess their variant classification and drug matching strategies.Methods: In this two centers study 48 patients with metastatic breast (n=12), colorectal (n=17) or non-small cell lung cancer (n=19) were recruited based on circulating tumor DNA levels in peripheral blood. In individual patients’ plasma samples, we established somatic copy number alterations and somatic mutations across a 77 cancer-associated gene panel and exposed the obtained molecular profiles to the decision support packages NAVIFY Mutation Profiler (NMP; Roche), Qiagen Clinical Insight (QCI) Interpret (Qiagen) and CureMatch Bionov (CureMatch).Results: In the plasma samples, altogether 492 somatic alterations were assessed. Each decision support platform varied in their format of data input, mode of variant classification and strategies for identifying druggable targets and clinical trials, which resulted in discrepancies in tier classification as well as designations of actionability. The frequency of concordant actionable events for tier I-A or tier I-B classifications, i.e. those with the strongest clinical evidence, was only 4.3%, 9.5% and 28.4% when comparing NMP with QCI, NMP with CureMatch, and CureMatch with QCI, respectively. As a consequence, the obtained treatment recommendations differed drastically, and alignment of treatment recommendations was limited to established predictive markers.Conclusions: The lack of standards for the extraction of treatment recommendations from comprehensive tumor molecular profiling data challenge the promising concepts of precision oncology.

  Study EGAS00001004383

• Combining a Universal Telomerase Based Cancer Vaccine with Ipilimumab in Patients with Metastatic Melanoma - Five-year Follow up of a Phase I/IIa Trial

  Background Ipilimumab improves survival for patients with metastatic malignant melanoma. Combining a therapeutic cancer vaccine with ipilimumab may increase efficacy by providing enhanced anti-tumor immune responses. UV1 consists of three synthetic long peptides from human telomerase reverse transcriptase (hTERT). These peptides comprise epitopes recognized by T cells from cancer patients experiencing long-term survival following treatment with a first-generation hTERT vaccine, and generate long-lasting immune responses in cancer patients when used as monotherapy. The objective of this trial was to investigate the safety and efficacy of combining UV1 with ipilimumab in metastatic melanoma. Patients and methods In this phase I/IIa, single center trial [NCT02275416], patients with metastatic melanoma received repeated UV1 vaccinations, with GM-CSF as an adjuvant, in combination with ipilimumab. Patients were evaluated for safety, efficacy and immune response. Immune responses against vaccine peptides were monitored in peripheral blood by measuring antigen-specific proliferation and IFN-γ production. Results Twelve patients were recruited. Adverse events were mainly diarrhea, injection site reaction, pruritus, rash, nausea and fatigue. Ten patients showed a Th1 immune response to UV1 peptides, occurring early and after few vaccinations. Three patients obtained a partial response and one patient a complete response. Overall survival was 50% at 5 years. Conclusion Treatment was well tolerated. The rapid expansion of UV1-specific Th1 cells in the majority of patients indicates synergy between UV1 vaccine and CTLA-4 blockade. This may have translated into clinical benefit, encouraging the combination of UV1 vaccination with standard of care treatment regimes containing ipilimumab/CTLA-4 blocking antibodies.

  Study EGAS00001005253

• Post_Mortem_Tissue_COVID19_spatial

  Single cell analysis of post mortem tissue samples from SARS-CoV2-infected patients We aim to directly identify the human cell types infected by SARS-CoV-2 and measure the cellular response to COVID19 infection across 20 different tissues from infected patient autopsies. We have examined the expression pattern of viral entry receptors across healthy human tissues to predict several candidate target cell types across the airway and heart. In addition, high prevalence of cardiac failure and abnormal renal function in COVID19 patients implicates heart and kidney involvement, but the pathogenesis of organ specific damage - whether via a direct cytopathic mechanism or an indirect inflammatory response - remains unknown . Currently, we lack confirmation of target cell types and cellular processes in infected tissues as autopsies are discouraged in most countries due to health and safety risks. Our collaborators Drs Michael Osborn and Brian Hanley (Imperial College) have outlined guidelines to perform post-mortem in COVID19 patients (Hanley et al., 2020) and have established a programme of autopsies for research to be performed in a high-risk facility at Westminster Public Mortuary. Here, we propose to identify infected cell types and aberrant molecular pathologies in this precious tissue resource using single cell and spatial genomics. We will prioritise three organ systems: the human airway, heart and the kidney. We will directly examine the cellular identities of SARS-CoV-2 infected cell types and identify the cellular responses to infection across these organs. This fundamental knowledge will help guide future treatment choices for COVID19. This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/

  Study EGAS00001004441

• Tumor-specific expression of let-7i and miR-192 is associated with resistance to cisplatin-based chemoradiotherapy in patients with larynx and hypopharynx cancer

  Background: The majority of patients with locally advanced larynx or hypopharynx squamous cell carcinoma are treated with organ preserving chemoradiotherapy (CRT). Clinical outcome following CRT varies greatly. We hypothesized tumor microRNA (miRNA) expression is predictive for outcome following CRT. Methods: Next-generation sequencing (NGS) miRNA profiling was performed on 37 formalin fixed paraffin embedded (FFPE) tumor samples. Patients with a recurrence free survival (RFS) of less than 2 years (“poor outcome”) and patients without recurrence within 2 years (“good outcome”) were compared by differential expression analysis. Tumor-specific miRNAs were selected based on normal mucosa miRNA expression data from The Cancer Genome Atlas database. A model was constructed to predict outcome using group-regularized penalized logistic ridge regression. Candidate miRNAs were validated in the initial sample set by RT-qPCR as well as in 46 additional tumor samples.Results: 1,938 known mature miRNAs were detected. Thirteen miRNAs were differentially expressed (p<0.05, FDR<0.1) according to outcome group. Initial class prediction in the NGS cohort (n=37) resulted in a model combining five miRNAs and disease stage, able to predict CRT outcome with an area under the curve (AUC) of 0.82. In the RT-qPCR cohort (n=83), 25 patients (30%) experienced early recurrence (median RFS 8 months; median follow-up 42 months). Class prediction resulted in a model combining let-7i-5p and miR-192-5p and disease stage, able to discriminate with an AUC of 0.80 patients with good versus poor clinical outcome.Conclusion: The combined miRNA expression and disease stage prediction model is superior to using either factor alone. ¬¬¬This study indicates NGS miRNA profiling using FFPE specimens is feasible, resulting in clinically relevant biomarkers.

  Study EGAS00001004001

• An analysis of humoral and cellular immune responses following COVID-19 vaccination.

  Two dose mRNA vaccination provides excellent protection against SARS-CoV-2. However, there are few data on vaccine efficacy in elderly individuals above the age of 801. Additionally, new variants of concern (VOC) with reduced sensitivity to neutralising antibodies have raised fears for vulnerable groups. Here we assessed humoral and cellular immune responses following vaccination with mRNA vaccine BNT162b22 in elderly participants prospectively recruited from the community and younger health care workers. Median age was 72 years and 51% were females amongst 140 participants. Neutralising antibody responses after the first vaccine dose diminished with increasing age, with a steep drop in participants over 80 years old. Sera from participants below and above 80 showed significantly lower neutralisation potency against B.1.1.7, B.1.351 and P.1. variants of concern as compared to wild type. Those over 80 were more likely to lack any neutralisation against VOC compared younger participants following first dose. Binding IgG and IgA antibodies were lower in the elderly, as was the frequency of SARS-CoV-2 Spike specific B- cells. We observed a trend towards lower somatic hypermutation in participants with suboptimal neutralisation, and elderly participants demonstrated clear reduction in class switched somatic hypermutation, driven by the IgA1/2 isotype. SARS-CoV-2 Spike specific T- cell IFN?? and IL-2 responses were secreted primarily by CD4 T cells, and impaired in the older age group. We conclude that the elderly are a high risk population that warrant specific measures in order to mitigate against vaccine failure, particularly where variants of concern are circulating. Specifically, the dosing interval should not be extended in this group.

  Study EGAS00001005380

• Analysis of RAD51C promoter methylation using targeted bisulfite sequencing (amplicon sequencing) in ovarian cancer pre-clinical models and patient samples.

  In high-grade serous ovarian carcinoma (HGSC), deleterious mutations in the DNA repair gene RAD51C are established drivers of defective homologous recombination and are emerging biomarkers of PARP inhibitor (PARPi) sensitivity. RAD51C promoter methylation (meRAD51C) is detected at similar frequencies to mutations, yet its effects on PARPi responses remain unresolved. In this study, three HGSC patient-derived xenograft (PDX) models with methylation at most or all examined CpG sites in the RAD51C promoter show responses to PARPi. Both complete and heterogeneous methylation patterns were associated with RAD51C gene silencing and homologous recombination deficiency (HRD). PDX models lost meRAD51C following treatment with PARPi rucaparib or niraparib, where a single unmethylated copy of RAD51C was sufficient to drive PARPi resistance. Genomic copy number profiling of one of the PDX models using SNP arrays revealed that this resistance was acquired independently in two genetically distinct lineages. In a cohort of 11 patients with RAD51C-methylated HGSC, various patterns of meRAD51C were associated with genomic 'scarring', indicative of HRD history, but exhibited no clear correlations with clinical outcome. Differences in methylation stability under treatment pressure were also observed between patients, where one HGSC was found to maintain meRAD51C after 6 lines of therapy (4 platinum-based), whilst another HGSC sample was found to have heterozygous meRAD51C and elevated RAD51C gene expression (relative to homozygous meRAD51C controls) after only neo-adjuvant chemotherapy. As meRAD51C loss in a single gene copy was sufficient to cause PARPi resistance in PDX, methylation zygosity should be carefully assessed in previously treated patients when considering PARPi therapy.

  Study EGAS00001005395

• IMCISION DNAseq

  Standard-of-care (SOC) surgery for locoregionally advanced head and neck squamous cell carcinoma (HNSCC) is intensive and results in 30‒50% five-year overall survival. Anti-PD1 immune checkpoint blockade (ICB) durably improves survival rates in recurrent/metastatic HNSCC. We report on the non-randomized phase Ib/IIa IMCISION trial (NCT03003637) of 32 HNSCC patients treated with two doses neoadjuvant nivolumab (NIVO MONO, n=6, phase Ib arm A) or two doses of nivolumab plus a single dose of ipilimumab (COMBO, n=26, 6 in phase Ib arm B and 20 in phase IIa) prior to surgery and, if indicated, adjuvant radiotherapy. Neoadjuvant ICB was feasible in all phase Ib patients, meeting the phase Ib primary feasibility endpoint. One phase IIa patient had progressive disease precluding surgery. Primary tumor pathological response (phase IIa primary endpoint), defined as the % change in viable tumor cells from baseline biopsy to on-treatment resection, was evaluable in 29 patients. . We observed a major pathological response (MPR, 90‒100% response) in 8/23 (35%) evaluable COMBO and 1/6 (17%) NIVO MONO patients. None of the patients with MPR after NIVO MONO or COMBO developed a tumor relapse after 24.0 months median follow-up. FDG-PET-based total lesion glycolysis identified MPR patients prior to surgery. A baseline AID/APOBEC-associated mutational profile and an on-treatment decrease in hypoxia signature expression were observed in MPR patients. Our data indicate that neoadjuvant NIVO MONO and COMBO are safe in HNSCC, and that particularly COMBO ICB shows encouraging effectivity. IMCISION provides rationales for future HNSCC trials aiming to pre-operatively select patients with a likely MPR after neoadjuvant ICB for de-escalation of SOC.

  Study EGAS00001005466

• Characterization of DLBCL with a PMBL gene expression signature

  Primary mediastinal large B-cell lymphoma (PMBL) is a type of aggressive B-cell lymphoma that typicallyaffects young adults, characterized by presence of a bulky anterior mediastinal mass. Lymphomas withgene expression features of PMBL have been described in non-mediastinal sites, raising questions abouthow these tumors should be classified.Here, we investigated whether these "non-mediastinal PMBLs" are indeed PMBLs or instead represent adistinct group within DLBCL. From a cohort of 325 de novo DLBCL cases, we identified tumors frompatients without evidence of anterior mediastinal involvement that expressed a PMBL expression signature(nm-PMBLsig-pos, n=16, 5%). The majority of these tumors expressed MAL and CD23 – proteins typicallyobserved in bona fide PMBL (bf-PMBL). Evaluation of clinical features of nm-PMBLsig-pos cases revealedclose associations with DLBCL, and the majority displayed a germinal center B-cell-like cell-of-origin(GCB). In contrast to bf-PMBL, nm-PMBLsig-pos patients presented at an older age, did not show pleuraldisease, and bone/bone marrow involvement was observed in three cases. However, while clinicallydistinct from bf-PMBL, nm-PMBL-sig-pos tumors resembled bf-PMBL at the molecular level with upregulationof immune response, JAK-STAT, and NF-kB signatures. Mutational analysis revealed frequent somatic genemutations in SOCS1, IL4R, ITPKB and STAT6, as well as CD83 and BIRC3, with the latter genes beingsignificantly more frequently affected than in GCB-DLBCL and bf-PMBL.Our data establish nm-PMBLsig-pos lymphomas as a group of DLBCL with distinct phenotypic and geneticfeatures, and potential implications for gene expression- and mutation-based subtyping of aggressive Bcell lymphoma and related targeted therapies.

  Study EGAS00001005057

• Whole genome sequencing of a breast cancer cohort with known functional homologous recombination status

  Optimal patient selection for PARP inhibitors and double strand break (DSB) inducing chemotherapy is of great clinical importance. In principle, the use of these therapies can be extended beyond germline BRCA1/2 mutated cancers. Although several methods for detection of BRCA-like tumors have been developed, a golden standard is lacking. We have developed the RECAP test, a functional HR assay exploiting the formation of RAD51 foci in proliferating cells after ex vivo irradiation of fresh primary breast cancer (BrC) tissue. The RECAP test seems robust and has an advantage over genetic tests, due to its functional character for exploring the HR phenotype. Another type of HRD test is the BRCAness classifier, which is based on specific genomic patterns derived from copy number data of BRCA1/2 mutated BrCs that also occur in sporadic cancers. The BRCAness classifier is a robust test that is easily applicable in the clinic and has also been validated to predict in vivo response to high dose chemotherapy in clinical trials. More recently, a new HRD test based on genomic signatures which can be extracted from whole genome sequencing has been published, HRDetect. Studies comparing these different HRD tests within the same cohort of patients do not yet exist. Here, we have the unique opportunity to compare three HRD tests within the same cohort. We aim to compare three different HRD tests (RECAP, BRCAness classifier and HRDetect) within the same cohort of primary BrCs (n=41). As a secondary aim, we will analyze additional genetic alterations in HRD tumors (both BRCA-deficient and BRCA-proficient HRD tumors).

  Study EGAS00001005572

• Genetics and therapeutic responses to TIL therapy of pancreatic cancer PDX models

  Pancreatic cancer is the 7th leading cause of cancer-related deaths worldwide. Checkpoint immunotherapy has not yet showed encouraging results in pancreatic cancer possibly owing to a poor immunogenicity and/or an immune suppressive microenvironment. The aim of this study was to develop patient-derived xenografts (PDX) models and to assess if autologous tumor-infiltrating lymphocytes (TILs) would have anti-tumoral activity in pancreatic cancer. Tumor biopsies from 29 patients were subcutaneously transplanted into NOG mice. Tumor growths were confirmed in 11 out of 29 transplantations. The PDX tumors histologically resembled their original biopsies, but since stromal cells in the PDX model tumors were from mouse, their gene expression differed from the original biopsies. Immune checkpoint ligands other than PD-L1 were expressed in pancreatic cancers but PD-L1 was rarely expressed. One of the three tumors that did express PD-L1 was an adenosquamous cancer and another had a mismatch repair deficiency. TILs were expanded from six tumors and were injected into NOG or human interleukin-2 transgenic-NOG (hIL2-NOG) mice carrying PDX tumors. Regression of tumors could be verified in hIL2-NOG mice in 3 of the 6 PDX models treated with autologous TILs, including the adenosquamous PDX model. In conclusion, PDX models of pancreatic cancer can be used to learn more about tumor characteristics and biomarkers, and to evaluate responses to ACT and combination therapies. The major benefit of the model is that modifications of T cells can be tested in an autologous humanized mouse model to gain preclinical data to support the initiation of a clinical trial.

  Study EGAS00001005596

• Proteom characterization in primary colorectal cancer and corresponding liver metastasis

  Colorectal adenocarcinomas (CRC) are one of the most commonly diagnosed tumors worldwide. Colorectal adenocarcinomas primarily metastasize into the liver and (less often) into the peritoneum. Patients suffering from CRC-liver metastasis (CRC-LM) typically present with a dismal overall survival compared to non-metastasized CRC patients. The metastasis process and metastasis promoting factors in patients with CRC are under intensive debate. However, CRC studies investigating the proteome biology are lacking. Formalin-fixed paraffin-embedded (FFPE) tissue specimens provide a valuable resource for comprehensive proteomic studies of a broad variety of clinical malignancies. The presented pilot study compares the proteome of primary CRC and patient-matched CRC-LM. The applied protocol allows a reproducible and straightforward identification and quantification of over 2,600 proteins within the dissected tumorous tissue. Subsequent unsupervised clustering reveals distinct proteome biologies of the primary CRC and the corresponding CRC-LM. Statistical analysis yields multiple differentially abundant proteins in either primary CRC or their corresponding liver metastases. A more detailed analysis of dysregulated biological processes suggests an active immune response in the liver metastases, including several proteins of the complement system. Proteins with structural roles, e.g. cytoskeleton organization or cell junction assembly appear to be less prominent in liver metastases as compared to primary CRC. Immunohistochemistry corroborates proteomic high expression levels of metabolic proteins in CRC-LM. We further assessed how the in vitro inhibition of two in CRC-LM enriched metabolic proteins affected cell proliferation and chemosensitivity. The presented proteomic investigation in a small clinical cohort promotes a more comprehensive understanding of the distinct proteome biology of primary CRC and their corresponding liver metastases.

  Study EGAS00001005641

• Genetic Risk for Subsequent Neoplasms among Long-term Survivors of Childhood Cancer in the St. Jude Lifetime Cohort

  While childhood cancer survivors are known to be at increased risk of subsequent neoplasms (SNs), the contribution of germline mutations in cancer predisposition genes to the risk of SN is largely unknown. Whole-genome sequencing (30X) was performed on 3,006 5+ year survivors of childhood cancer (median age, 35.8 [range: 7.1-69.8] years) from the St. Jude Lifetime Cohort. Survivors underwent a comprehensive clinical assessment. Germline mutations in 156 cancer predisposition genes (SJCPG156), including a subset of 60 (SJCPG60) autosomal dominant genes with moderate to high penetrance, were classified for their pathogenicity. Piecewise exponential regression, stratified by radiation exposure, was used to evaluate the relative rate (RR) and 95% confidence interval (95% CI) of SN occurrence by mutation status. Pathogenic or likely pathogenic mutations in SJCPG156 and SJCPG60 were identified in 11.7% and 5.8% of survivors, respectively. The most frequently mutated SJCPG60 genes included RB1 (n=43), NF1 (n=22), BRCA2 (n=14), BRCA1 (n=12) and TP53 (n=10). Mutations in SJCPG60 were associated with the rate of subsequent sarcoma (RR, 10.3; 95% CI, 4.2 to 25.6) and breast cancer (RR, 13.3; 95% CI, 5.8 to 30.5) among irradiated survivors, and the rate of developing any SN (RR, 4.8; 95% CI, 2.4 to 9.4) and breast cancer (RR, 7.8; 95% CI, 2.5 to 24.0) among non-irradiated survivors. The rate of developing ≥2 histologically distinct SNs was increased 17.8-fold (95% CI, 3.5 to 89.8) among SJCPG60 mutation carriers treated without radiotherapy. Similar but attenuated associations were observed for mutations in SJCPG156. Survivors who develop a SN without prior radiotherapy or who are diagnosed with a subsequent sarcoma or breast cancer within a site of prior radiotherapy should be referred to genetic counseling.

  Study EGAS00001002499

• Functional Mapping of AKT Signaling and Biomarkers of Response From the FAIRLANE Trial of Neoadjuvant Ipatasertib Plus Paclitaxel for Triple-Negative Breast Cancer

  Purpose: Despite extensive genomic and transcriptomic profiling, it remains unknown how signaling pathways are differentially activated and how tumors are differentially sensitized to certain perturbations. Here, we aim to characterize AKT signaling activity and its association with other genomic or immunohistochemistry-based PI3K/AKT pathway biomarkers as well as the clinical activity of ipatasertib (AKT inhibitor) in the FAIRLANE trial. Experimental Design: In FAIRLANE, 151 patients with early triple-negative breast cancer were randomized 1:1 to receive paclitaxel with ipatasertib or placebo for 12 weeks prior to surgery. Adding ipatasertib did not increase pathologic complete response rate and numerically improved overall response rate by magnetic resonance imaging (MRI). We used reverse-phase protein microarrays (RPPA) to examine the total level and/or phosphorylation states of over 100 proteins in various signaling or cell processes including PI3K/AKT and mTOR signaling. 125 baseline and 127 on-treatment samples were evaluable by RPPA, with 110 paired samples at both time points. Results: Tumors with genomic/protein alterations in PIK3CA/AKT1/PTEN were associated with higher levels of AKT phosphorylation. In addition, phosphorylated(p)AKT levels exhibited a significant association with enriched clinical benefit of ipatasertib, and identified patients who received benefit in the absence of PIK3CA/AKT1/PTEN alterations. Ipatasertib treatment led to a down-regulation of AKT/mTORC1 signaling, which was more pronounced among the tumors with PIK3CA/AKT1/PTEN alterations or among the responders to the treatment. Conclusions: We showed that the high baseline pAKT levels are associated with the alterations of PI3K/AKT pathway components and enriched benefit of ipatasertib in TNBC.

  Study EGAS00001005892

• Genomic landscape of oral cancers (Illumina WGS)

  Human papillomavirus (HPV) is a necessary but insufficient cause of a subset of oral squamous cell carcinomas (OSCC) that is increasing markedly in frequency. To identify contributory, secondary genetic alterations in these cancers, we used comprehensive genomics methods to compare 149 HPV-positive and 335 HPV-negative OSCC tumor/normal pairs. Different behavioral risk factors underlying the two OSCC types were reflected in distinctive genomic mutational signatures. In HPV-positive OSCC, the signatures of APOBEC cytosine deaminase editing, associated with anti-viral immunity, were strongly linked to overall mutational burden. By contrast, in HPV-negative OSCC, T>C substitutions in the sequence context 5’-ATN-3’ correlated with tobacco exposure. Universal expression of HPV E6*1 and E7 oncogenes was a sine qua non of HPV-positive OSCC. Significant enrichment of somatic mutations was confirmed or newly identified in PIK3CA, KMT2D, FGFR3, FBXW7, DDX3X, PTEN, TRAF3, RB1, CYLD, RIPK4, ZNF750, EP300, CASZ1, TAF5, RBL1, IFNGR1 and NFKBIA. Of these, many affect host pathways already targeted by HPV oncoproteins, including the p53 and pRB pathways, or disrupt host defenses against viral infections, including interferon and nuclear factor-B signaling. Frequent copy number changes were associated with concordant changes in gene expression. Chr. 11q (including CCND1) and 14q (including DICER1 and AKT1) were recurrently lost in HPV-positive OSCC, in contrast to their gains in HPV-negative OSCC. High-ranking variant allele fractions implicated ZNF750, PIK3CA and EP300 mutations as candidate driver events in HPV-positive cancers. We conclude that virus-host interactions cooperatively shape the unique genetic features of these cancers, distinguishing them from their HPV-negative counterparts.

  Study EGAS00001003228

• Monocyte and macrophage lipid accumulation results in downregulated type-I interferon responses

  Macrophages are critical components of atherosclerotic lesions and their pro- and anti-inflammatory responses influence atherogenesis. Type-I interferons (IFNs) are cytokines that play an essential role in antiviral responses and inflammatory activation and have been shown to promote atherosclerosis. Although the impact of type-I IFNs on macrophage foam cell formation is well-documented, the effect of lipid accumulation in monocytes and macrophages on type-I IFN responses remains unknown. Here we examined IFN stimulated (ISG) and non-ISG inflammatory gene expression in mouse and human macrophages that were loaded with acetylated LDL (acLDL), as a model for foam cell formation. We found that acLDL loading in mouse and human macrophages specifically suppressed expression of ISGs and IFN-β secretion, but not other pro-inflammatory genes. The downregulation of ISGs could be rescued by exogenous IFN-β supplementation. Activation of the cholesterol-sensing nuclear liver X receptor (LXR) recapitulated the cholesterol-initiated type-I IFN suppression. Additional analyses of murine in vitro and in vivo generated foam cells confirmed the suppressed IFN signalling pathways and suggest that this phenotype is mediated via downregulation of interferon regulatory factor binding at gene promoters. Finally, RNA-seq analysis of monocytes of familial hypercholesterolemia (FH) patients also showed type-I IFN suppression which was restored by lipid-lowering therapy and not present in monocytes of healthy donors. Taken together, we define type-I IFN suppression as an athero-protective characteristic of foamy macrophages. These data provide new insights into the mechanisms that control inflammatory responses in hyperlipidaemic settings and can support future therapeutic approaches focusing on reprogramming of macrophages to reduce atherosclerotic plaque progression and improve stability.

  Study EGAS00001005955

• ctDNA monitoring using patient-specific sequencing and integration of variant reads - Breast cohort

  Circulating tumor-derived DNA (ctDNA) can be used to monitor cancer dynamics noninvasively. Detection of ctDNA can be challenging in patients with low-volume or residual disease, where plasma contains very few tumor-derived DNA fragments. We show that sensitivity for ctDNA detection in plasma can be improved by analyzing hundreds to thousands of mutations that are first identified by tumor genotyping. We describe the INtegration of VAriant Reads (INVAR) pipeline, which combines custom error-suppression methods and signal-enrichment approaches based on biological features of ctDNA. With this approach, the detection limit in each sample can be estimated independently based on the number of informative reads sequenced across multiple patient-specific loci. We applied INVAR to custom hybrid-capture sequencing data from 176 plasma samples from 105 patients with melanoma, lung, renal, glioma, and breast cancer across both early and advanced disease. By integrating signal across a median of >105 informative reads, ctDNA was routinely quantified to 1 mutant molecule per 100,000, and in some cases with high tumor mutation burden and/or plasma input material, to individual parts per million. This resulted in median Area Under the Curve (AUC) values of 0.98 in advanced cancers, and 0.80 in early stage and challenging settings for ctDNA detection. We generalized this method to whole-exome and whole-genome sequencing, showing that the INVAR may be applied without requiring personalized sequencing panels, so long as a tumor mutation list is available. As tumor sequencing becomes increasingly performed, such methods for personalized cancer monitoring may enhance the sensitivity of cancer liquid biopsies.

  Study EGAS00001004446

• Profiling the unique protective properties of intracranial arterial endothelial cells

  Cardiovascular disorders, like atherosclerosis and hypertension, are increasingly known to be associated with vascular cognitive impairment (VCI). In particular, intracranial atherosclerosis is one of the main causes of VCI, although plaque development occurs later in time and is structurally different compared to atherosclerosis in extracranial arteries. Recent data suggest that endothelial cells (ECs) that line the intracranial arteries may exert anti-atherosclerotic effects due to yet unidentified pathways. To gain insights into underlying mechanisms, we isolated post-mortem endothelial cells from both the intracranial basilar artery (BA) and the extracranial common carotid artery (CCA) from the same individual (total of 15 individuals) with laser capture microdissection. RNA sequencing revealed a distinct molecular signature of the two endothelial cell populations of which the most prominent ones were validated by means of qPCR. Our data reveal for the first time that intracranial artery ECs exert an immune quiescent phenotype. Secondly, genes known to be involved in the response of ECs to damage (inflammation, differentiation, adhesion, proliferation, permeability and oxidative stress) are differentially expressed in intracranial ECs compared to extracranial ECs. Finally, Desmoplakin (DSP) and Hop Homeobox (HOPX), two genes expressed at a higher level in intracranial ECs, and Sodium Voltage-Gated Channel Beta Subunit 3 (SCN3B), a gene expressed at a lower level in intracranial ECs compared to extracranial ECs, were shown to be responsive to shear stress and/or hypoxia. With our data we present a set of intracranial-specific endothelial genes that may contribute to its protective phenotype, thereby supporting proper perfusion and consequently may preserve cognitive function. Deciphering the molecular regulation of the vascular bed in the brain may lead to the identification of novel potential intervention strategies to halt vascular associated disorders, such as atherosclerosis and vascular cognitive dysfunction.

  Study EGAS00001004479

• Clinical Outcomes by Tumor Mutational Burden and Inflammatory Gene Expression With Combined Nivolumab and Ipilimumab or Monotherapy in Advanced Melanoma - CM67 WES

  Purpose Exploratory analyses of CheckMate 066 and 067 trials were conducted to investigate associations of tumor mutational burden (TMB), a 4-gene inflammatory gene expression signature, and BRAF mutation status with tumor response, progression-free survival (PFS), and overall survival (OS) in patients with advanced melanoma. Patients and Methods Patients with known programmed death ligand 1 (PD-L1) expression and BRAF mutation status received nivolumab (NIVO) or dacarbazine in CheckMate 066 and either NIVO, ipilimumab (IPI), or NIVO+IPI in CheckMate 067. Whole exome sequencing and RNA sequencing were used to determine TMB and inflammatory gene expression signature scores, respectively. These biomarkers were evaluated in terms of their association with PFS and OS. Results In the NIVO, NIVO+IPI, and IPI arms of CheckMate 067, longer survival was associated with high (> median) versus low (≤ median) TMB with hazard ratios (HRs) (95% confidence interval [CI]) for PFS of 0.45 (0.30–0.65), 0.55 (0.38–0.81), and 0.60 (0.43–0.82), and for OS of 0.46 (0.30–0.71), 0.53 (0.34–0.82), and 0.52 (0.36–0.74), respectively. For NIVO-treated patients, these results were confirmed in CheckMate 066. A survival benefit was observed with high TMB and absence of BRAF mutation. Survival was associated with high versus low inflammatory signature scores with HRs (95% CI) for PFS of 0.56 (0.34–0.94), 0.40 (0.23–0.72), and 0.43 (0.27–0.70), and for OS of 0.37 (0.20–0.66), 0.38 (0.19–0.74), and 0.46 (0.27–0.79), in the NIVO, NIVO+IPI, and IPI arms, respectively. Weak correlations were observed between PD-L1, TMB, and the inflammatory signature. Conclusions Combined assessment of TMB, inflammatory gene expression signature, and BRAF mutation status may be predictive for response to immunotherapy in advanced melanoma.

  Study EGAS00001004555

• Germline variants in patients with rare cancers and their implications for precision cancer medicine: experiences from the Multicenter MASTER Trial by the German Cancer Consortium (HIPO_021)

  The era of precision oncology entails an increasing amount of molecular data from tumor sequencing in research and diagnostics. In this setting, germline variant analysis may add clinical benefits for patients and family members. In this study, we evaluated germline variants from genome/exome matched tumor-blood sequencing in 1.485 advanced cancer patients within a prospective observational study carried out by the German Cancer Consortium. Included patients had cancers of rare entities (79.2%) or were under the age of 51 at first diagnosis (76.6%). Ten percent of the patients harbored pathogenic germline variants in 35 tumor predisposition genes with autosomal-dominant inheritance such as BRCA1, BRCA2, TP53, RB1, CHEK2, ATM, PALB2, SDHB etc., while 4.6% of the patients were carriers of pathogenic germline variants in 18 genes with autosomal-recessive inheritance of tumor predisposition. Based on pathogenic germline variants, 14% of the patients received genetic counselling recommendations in the molecular tumor board report. Forty five percent of pathogenic germline variants supported a molecular treatment recommendation, which was implemented for 28 patients. Furthermore, we compared our approach and results to previous studies, broadening the insight into pan cancer germline findings but also revealing differences between studies on various levels, such as cohort demographics, tumor entities, gene nominations, variant detection and variant interpretation. Our data support the importance of germline variant interpretation in precision oncology studies for patients and their families, if applicable. However, guidelines for evaluation and reporting of germline variants in the context of genetic tumor profiling need to be established and further research on the relevance of germline mutations as biomarkers for treatment recommendations is needed.

  Study EGAS00001005537

• Stromal-induced epithelial-mesenchymal transition induces drug resistance in acute lymphoblastic leukemia

  Bone marrow (BM) microenvironment-induced drug resistance in acute lymphoblastic leukemia (ALL) is thought to be mediated by intercellular interaction of leukemic cells with BM niches, but the mechanistic basis of resistance is poorly understood. Here, we genomic and functional analyses, we show that ALL-BM mesenchymal stromal cells (MSC) interactions encompass direct cell-cell adhesion via integrin α4β1, uni-directional secretion of soluble factors from MSC to ALL, and bidirectional transfer of cellular vesicles through tunneling nanotubes. We show that this crosstalk drives an epithelial-mesenchymal transition (EMT)-like program in ALL cells adhering to MSC, thereby triggering stem cell-like features including drug resistance and survival advantage. Moreover, single-cell RNA sequencing of cell line and primary ALL samples identified a hybrid cell state derived from B-ALL adhered to MSC, that exhibits both B-ALL and MSC gene signatures that are orchestrated by an EMT-like program. Multiple EMT-activating pathways including WNT/β-catenin, cytokine/STAT3, TGFβ, NOTCH1, and NF-κB-dependent signaling pathways are activated in MSC-adherent ALL cells; however, only the inhibition of either WNT/β-catenin or cytokine/STAT3 pathway could attenuate the survival and drug resistance of ALL cells adherent to MSC. Notably, we identified that this hybrid cell state exhibited elevated WNT/β-catenin-mediated transcriptional activity. Finally, we show that blocking of the interaction of β-catenin and CBP interaction abrogates the survival and drug resistance of ALL cells adhered to MSC, indicating that targeting of an EMT-like program is potential therapeutic regimen in targeting cell extrinsically acquired drug resistance in ALL.

  Study EGAS00001006120

• Exploration of coding and non-coding variants in cancer using GenomePaint.

  GenomePaint (https://proteinpaint.stjude.org/genomepaint) is a dynamic visualization platform for whole-genome, whole-exome, transcriptome, and epigenomic data, featuring a novel design that captures the inter-relatedness between DNA variations and RNA expression. Regulatory non-coding variants can be inspected and discovered along with coding variants, and their functional impact further explored by examining 3D genome and/or ChIP-seq data generated from cancer cell lines. Further, GenomePaint correlates mutation and expression patterns with patient outcomes, and can display external data such as adult cancer datasets and user-provided custom tracks. We used GenomePaint to analyze multi-omics data from 3,652 pediatric cancers representing 16 histotypes, and demonstrate the visualization features through examples, including two that led to new insights into oncogenic mechanisms in pediatric cancer. The first is the discovery of a new class of pathogenic recurrent variants that cause aberrant splicing, disrupting the RING domain of CREBBP, a driver gene frequently mutated in relapsed pediatric leukemia. The second is the cis-activation of the MYC oncogene in a subset of B-lineage acute lymphoblastic leukemia (B-ALL) via duplication of the NOTCH1-MYC enhancer (N-ME), previously discovered only in T-lineage ALL. The regulatory impact of N-ME enhancer amplification was initially confirmed by allelic imbalance in published gene expression and ChIP-seq data and verified by additional Capture-C and fluorescence in situ hybridization data generated by follow-up experiments. These examples demonstrate the power of GenomePaint in enabling not only data visualization but also integrative genomic analysis that can lead to novel biological insight for follow-up experimental validation.

  Study EGAS00001004669

• Molecular and clinical diversity in primary central nervous system lymphoma: a LOC Network study

  Primary central nervous system lymphoma (PCNSL) is a distinct extranodal lymphoma presenting with limited stage disease but variable response rates to treatment despite homogenous pathological presentation. The likely underlying molecular heterogeneity and its clinical impact is poorly understood. We performed a comprehensive genome-wide analysis of 147 PCNSL from fresh-frozen tumor tissue from immunocompetent, treatment naïve PCNSL patients, employing whole-exome sequencing, assessment of somatic copy number alterations and DNA methylation, and RNA expression. These data were integrated and correlated with the clinico-radiological characteristics and outcomes of the patients. We validated our results in an independent series of 93 PCNSL formalin-fixed, paraffin-embedded (FFPE) samples. Consensus clustering of multi-omics data identified four robust, non-overlapping, prognostically significant clusters (CS) within PCNSL. The CS1 group, characterized by high proliferation and Polycomb Repressive Complex 2 (PRC2) complex activity had an intermediate outcome between CS2/CS3 and CS4. Patients who had PCNSL with an “immune-hot” (CS4) profile had the most favorable clinical outcome. In contrast, patients with the immune-cold hypermethylated CS2 and the heterogenous-immune CS3 groups had a poor prognosis. Nearly all PCNSL patients with meningeal infiltration harbored HIST1H1E mutations, enriched in the CS3 group. The integrated analysis suggests that the CS4 group may be more susceptible to immunotherapy. The integration of genome-wide data from multi-omics data revealed four molecular patterns in PCNSL with a distinctive prognostic impact that significantly improved the current clinical stratification. This molecular classification using FFPE samples facilitates routine use in clinical practice and provides potential precision-medicine strategies in PCNSL.

  Study EGAS00001006191

• Combined MEK and BCL-2/XL inhibition is effective in high-grade serous ovarian cancer patient-derived xenograft models and BIM levels are predictive of responsiveness

  Most patients with late stage high-grade serous ovarian cancer (HGSOC) initially respond to chemotherapy but inevitably relapse and develop resistance, highlighting the need for novel therapies to improve patient outcomes. The MEK/ERK pathway is activated in a large subset of HGSOC, thus making it an attractive therapeutic target. Here, we systematically evaluated the extent of MEK/ERK pathway activation and efficacy of pathway inhibition in a large panel of well-annotated HGSOC patient-derived xenograft (PDX) models. The vast majority of models were nonresponsive to the MEK inhibitor cobimetinib (GDC-0973) despite effective pathway inhibition. Proteomic analyses of adaptive responses to GDC-0973 revealed that GDC-0973 upregulated the pro-apoptotic protein BIM, thus priming the cells for apoptosis regulated by BCL2-family proteins. Indeed, combination of both MEK inhibitor and dual BCL-2/XL inhibitor (ABT-263) significantly reduced cell number, increased cell death and displayed synergy in vitro in most models. In vivo, the GDC-0973 and ABT-263 combination was well tolerated and resulted in greater tumor growth inhibition than single agents. Detailed proteomic and correlation analyses identified two subsets of responsive models – those with high BIM at baseline that was increased with MEK inhibition and those with low basal Bim and high pERK levels. Models with low BIM and low pERK were non-responsive. Our findings demonstrate that combined MEK and BCL-2/XL inhibition has therapeutic activity in HGSOC models and provide a mechanistic rationale for clinical evaluation of this drug combination as well as the assessment of the extent to which BIM and/or pERK levels predict drug combination effectiveness in chemoresistant HGSOC.

  Study EGAS00001003427

• Analysis of error profiles in deep next-generation sequencing data

  BackgroundSequencing errors are key confounding factors for detecting low-frequency genetic variants that are important for cancer molecular diagnosis, treatment, and surveillance using deep next-generation sequencing (NGS). However, there is a lack of comprehensive understanding of errors introduced at various steps of a conventional NGS workflow, such as sample handling, library preparation, PCR enrichment, and sequencing. In this study, we use current NGS technology to systematically investigate these questions.ResultsBy evaluating read-specific error distributions, we discover that the substitution error rate can be computationally suppressed to 10-5 to 10-4, which is 10- to 100-fold lower than generally considered achievable (10-3) in the current literature. We then quantify substitution errors attributable to sample handling, library preparation, enrichment PCR, and sequencing by using multiple deep sequencing datasets. We find that error rates differ by nucleotide substitution types, ranging from 10-5 for A>C/T>G, C>A/G>T, and C>G/G>C changes to 10-4 for A>G/T>C changes. Furthermore, C>T/G>A errors exhibit strong sequence context dependency, sample-specific effects dominate elevated C>A/G>T errors, and target enrichment PCR led to aprox 6-fold increase of overall error rate. We also find that more than 70% of hotspot variants can be detected at 0.1% - 0.01% frequency with the current NGS technology by applying in-silico error suppression.ConclusionsWe present the first comprehensive analysis of sequencing error sources in conventional NGS workflows. The error profiles revealed by our study highlight new directions for further improving NGS analysis accuracy both experimentally and computationally, ultimately enhancing the precision of deep sequencing.

  Study EGAS00001003444

• Identification of the dismal subtype of B-ALL with dysregulation of CDX2 and UBTF

  PURPOSE: Advance in the knowledge of genomic basis of B-cell acute lymphoblastic leukemia (B-ALL) has changed the treatment approaches and diagnostic assays and improved their outcome. Despite the recent development of next generation sequencing approaches, there are still some cases that their drivers are unknown. We aimed to define and describe the characteristics of an additional B-ALL subtype by integrating genomic, transcriptomic, and epigenomic approach. METHODS: More than 2,000 B-ALL cases of RNA-seq data from from Eastern Cooperative Oncology Group (ECOG) and the American College of Radiology Imaging Network (ACRIN) (n=764), the Children's Oncology Group (COG), Munich Leukemia Laboratory (MLL) (n=280), and other several collaborators were analyzed to identify the new subtype. Whole genome sequence was performed to detect mutations, stractural variants, and copy number alterations. The tridimentional analysis of chromatin was performed with HiChIP. RESULTS: We identified the new high-risk B-ALL subtype “CDX2/UBTF” that exhibited unique gene expression profiles, which is enriched in female (77.3%), adolescents and young adults (52.6%), and relapse cohort (3.3%). CD10 negativity and IgM positivity are hallmarks of this subtype with NTRK3 expression that can be an effective marker at diagnostic approach and a potential targeted therapy with TRK and mulkinase inbitors, larotrectinib and entrectinib. Genetically, two different alterations define CDX2/UBTF. One is UBTF-ATXN7L3 fusions caused by deletion of 17q21.31, and the other is extopic CDX2 expression through the enhancer hijacking mechanism induced by deletion of 13q12.2. Other genomic features include gain of 1q and PAX5 rearrangement (PAX5-ZCCHC7), might induce upregulation of histone cluster genes and PAX5. CONCLUSION: We described the novel B-ALL subtype “CDX2/UBTF” that has unique clinical and genomic characteristics.

  Study EGAS00001005863

• Comprehensive characterization of pre- and post-treatment samples of breast cancer reveal potential mechanisms of chemotherapy resistance

  When locally advanced breast cancer is treated with neoadjuvant chemotherapy, the recurrence risk is significantly higher if no complete pathologic response is achieved. Identification of the underlying resistance mechanisms is essential to select treatments with maximal efficacy and minimal toxicity. Here we employed gene expression profiles derived from 317 HER2-negative treatment-naïve breast cancer biopsies of patients who underwent neoadjuvant chemotherapy, deep whole exome and RNA-sequencing profiles of 22 matched pre- and post-treatment tumors, and treatment outcome data to identify biomarkers of response and resistance mechanisms. Molecular profiling of treatment-naïve breast cancer samples revealed that expression levels of proliferation, immune response and extracellular matrix (ECM) organization combined predict response to chemotherapy. Triple negative patients with high proliferation, high immune response and low ECM expression had a significantly better treatment response and survival benefit (HR 0.29, 95% CI 0.10-0.85; p=0.02), while in ER+ patients the opposite was seen (HR 4.73, 95% CI 1.51-14.8; 0=0.008). The characterization of paired pre-and post-treatment samples revealed that aberrations of known cancer genes were either only present in the pre-treatment sample (CDKN1B) or in the post-treatment sample(TP53, APC, CTNNB1). Proliferation-associated genes were frequently down-regulated in post-treatment ER+ tumors, but not in triple negative tumors. Genes involved in ECM were upregulated in the majority of post-chemotherapy samples. Genomic and transcriptomic differences between pre- and post-chemotherapy samples are common and may reveal potential mechanisms of therapy resistance. Our results show a wide range of distinct, but related mechanisms, with a prominent role for proliferation- and ECM-related genes.

  Study EGAS00001005876

• Mutational patterns and regulatory networks in epigenetic subgroups of meningioma (H033)

  DNA methylation patterns delineate clinically relevant subgroups of meningioma. We previously established the six meningioma methylation classes (MC) benign-1, 2, 3, intermediate-A, B and malignant. Here, we set out to identify subgroup-specific mutational patterns and pathway regulation. Whole-genome-sequencing was performed on 62 samples across all MCs and WHO grades from 62 patients with matched blood control, including 40 sporadic and 22 radiation-induced (Mrad) meningiomas. RNA sequencing was added for 18 cases and chromatin-immunoprecipitation for the enhancer mark H3K27ac followed by sequencing (ChIP-seq) for 16 samples. Besides the known mutations in meningioma, structural alterations were found to contribute to the spectrum of mechanisms inactivating NF2 in sporadic meningioma similar to previous reports for Mrad. Aberrations of DMD were found to be enriched in MCs with NF2 mutations and DMD was among the most differentially upregulated genes in NF2 mutant compared to NF2 wild-type cases. The mutational signature AC3 was detected in both sporadic meningioma and Mrad, but distributed across the genome in sporadic cases and enriched near genomic breakpoints in Mrad. In general, the malignant MC presented a significantly higher exposure to AC3 and higher genomic instability beyond the mutational load than the other MCs. Pathway analysis of the ChIP-seq clusters revealed that the FOXM1 is most differentially activated in high-grade cases, along with super enhancer recruitment near HOX genes and respective upregulation of expression. This data further elucidate the biological mechanisms differentiating meningiomas of different WHO grade and epigenetic subgroups and suggest leveraging the genomic instability as novel therapeutic targets.

  Study EGAS00001003481

• Genetic immune escape landscape in primary and metastatic cancer

  Immune surveillance escape is a hallmark of tumorigenesis1. Multiple studies have characterized the immune escape landscape across several untreated early-stage primary cancer types2–4. However, whether late-stage treated metastatic tumors present differences in genetic immune escape (GIE) prevalence and dynamics remains unclear. Here, we performed a pan-cancer characterization of GIE prevalence across six immune escape pathways in 6,457 uniformly processed Whole Genome Sequencing (WGS) tumor samples including 58 cancer types from 1,943 primary untreated patients and 4,514 metastatic patients. To effectively address the complexity of the Human Leukocyte Antigen (HLA-I) locus and to characterize its tumor status, we developed LILAC, an open-source integrative framework. We demonstrate that one in four tumors harbor GIE alterations, with high mechanistic and frequency variability across cancer types. GIE prevalence is highly consistent between primary and metastatic tumors for most cancer types with few exceptions such as prostate and thyroid carcinomas that have increased immune evasion frequencies in metastatic tumors. Positive selection analysis revealed that GIE alterations are frequently selected for in tumor evolution and that focal LOH of HLA-I, unlike non-focal LOH of HLA-I, tends to lose the HLA allele that presents the largest neoepitope repertoire. We also unraveled tumor genomic features contributing to immune escape incidence, including DNA repair deficiency, APOBEC activity, tobacco associated mutation load and viral DNA integration. Finally, there is a strong tendency for mid and high tumor mutation burden (TMB) tumors to preferentially select LOH of HLA-I for GIE whereas hypermutated samples favor global immune evasion strategies. Our results indicate that genetic immune escape is generally a pre-metastatic event during tumor evolution and that tumors adapt different strategies depending on their neoepitope burden.

  Study EGAS00001006123

• Combination of CDK4/6 with BET-bromodomain and PI3K/mTOR inhibitors in medulloblastoma in vitro and in vivo

  Despite improvement in the treatment of medulloblastoma (MB) over the last years, numerous patients with MYC- and MYCN-driven tumors still fail the current therapeutic regimen. The retinoblastoma (RB) pathway is often compromised in MB with RB remaining intact suggesting that CDK4/6 inhibition might represent a therapeutic strategy for which drug combination remains understudied. Using tumor lines grown in stem cell conditions referred as tumorspheres as well as tumor cells freshly dissociated from the brain of mice harboring intracranial patient-derived xenografts (PDX) of Group3 (G3) and Sonic Hedgehog (SHH) MB, we conducted high throughput drug combination screens using the CDK4/6 inhibitor (CDK4/6i) ribociclib as an anchor and 87 oncology drugs approved by FDA or currently in clinical trials. Bromodomain and extra terminal (BET) and PI3K/mTOR inhibitors potentiated ribociclib’s inhibition of proliferation in G3 and SHH MB. A reverse combination screen using BET inhibitor JQ1 as the anchor, revealed the three FDA approved CDK4/6i as the most potentiating drugs. Those synergistic combinations were tested in vivo in mice intracranially implanted with G3 and SHH MB tumor cells. Ribociclib showed single agent activity in several in vivo MB models. However, despite in vitro synergy, combinations of ribociclib with JQ1 and ribociclib with the PI3K/mTOR inhibitor paxalisib failed to improve the survival of MB bearing mice compared to ribociclib alone. Our data illustrates the difficulty in translating in vitro findings to in vivo for brain tumors partly due to insufficient free brain concentration and altered signaling pathways activation in vitro compared to in vivo.

  Study EGAS00001006286

• Long-term temporal stability of peripheral blood DNA methylation alterations in patients with inflammatory bowel disease.

  Introduction There is great current interest in the potential application of DNA methylation alterations in peripheral blood as biomarkers of susceptibility, progression, and treatment response in inflammatory bowel disease (IBD). To date, the intra-individual stability of peripheral blood leukocyte (PBL) methylation in IBD patients, a key consideration in biomarker development, has not been characterised. Here, we studied the long-term stability of all probes located on the Illumina HumanMethylation EPIC BeadChip array. Methods We followed a cohort of 46 adult IBD patients (36 Crohn’s disease (CD), 10 ulcerative colitis (UC), median age 44 (IQR: 27-56), 50% female) that received standard care without any intervention at the Amsterdam UMC. Paired PBL samples were collected at two time points with a median 7 (range: 2-9) years in between. Differential methylation and intra-class correlation (ICC) analyses were used to identify time-associated differences and temporally stable CpGs, respectively. Results Around 60% of all CpG loci on the EPIC array presented poor (ICC <0.50); 119,388 (≈14%) and 41,274 (≈5%) showed good (ICC ≥0.75) or excellent (ICC ≥0.90) intra-individual stability. Functional over-representation analyses of the stable methylated CpGs implicated genes involved in cell-cell signaling, adhesion and neurogenesis. Focusing on previously identified consistently differentially methylated positions indicated that 22 CD-, 11 UC-, and 24 IBD-associated loci demonstrated good stability (ICC ≥0.75) over time, where we observed a marked stability of CpG loci associated to the HLA genes. Conclusion Our data provide insight into the long-term stability of the PBL DNA methylome within an IBD context, facilitating the selection of biologically relevant and robust IBD-associated epigenetic biomarkers with increased potential for independent validation.

  Study EGAS00001006501

• BCR Signaling in human BM PC

  The functions of bone marrow plasma cells (BMPC) beyond antibody production are not fully elucidated and distinct subsets of BMPC suggest potential different functions. Phenotypic differences were identified for human BMPC depending on CD19 expression. Since CD19 is a co-stimulatory molecule of the B-cell-receptor (BCR), and IgA+ and IgM+ BMPC express the BCR on their surface, we here studied whether CD19 expression affects cellular responses, such as BCR signaling and the expression of checkpoint molecules. We analyzed 132 BM samples from individuals undergoing routine total hip arthroplasty. We found that both CD19+ and CD19- BMPC expressed BCR signaling molecules. Notably, the BCR-associated kinase spleen tyrosine kinase (SYK) including pSYK was higher expressed in CD19+BMPC compared to CD19- BMPC. BCR stimulation also resulted in increased kinase phosphorylation downstream of the BCR while expression of CD19 remained stable afterwards. Interestingly, the BCR response was restricted to IgA+ BMPC independently of CD19 expression. With regard to the expression of checkpoint molecules, CD19- BMPC expressed higher levels of co-inhibitory molecule programmed cell death protein-1 (PD-1) than CD19+ BMPC. IgA+ BMPC characteristically upregulated PD-1 upon BCR stimulation in contrast to other PC subsets and inhibition of the kinase SYK abrogated PD-1 upregulation. In contrast, expression of PD-1 ligand, B and T lymphocyte attenuator (BTLA) and CD28 did not change upon BCR activation of IgA+ BMPC. Here, we identify a distinct functiotype of IgA+ BMPC that is independent of the phenotypic heterogeneity of the subsets according to their CD19 expression. The data suggest that IgA+ BMPC underlie different regulatory principles and/or exert distinct regulatory functions.

  Study EGAS00001004948

• Immune microenvironment and lineage tracing help deciphering Rosette-forming GlioNeuronal Tumors: a multi-omic analysis of 9 cases

  Rosette-forming GlioNeuronal Tumor (RGNT) is a rare central nervous system neoplasm containing two components, glial and neuronal. The delineation of histological diagnosis of RGNT from similar low-grade tumors such as pilocytic astrocytoma (PA) and ganglioglioma (GG) may be challenging. We performed here a comprehensive molecular analysis of a cohort of tumors with same histology features to identify molecular characteristics of RGNT. A cohort of 9 rosette-forming tumors histologically-diagnosed were analyzed at molecular level using multimodal approaches as Whole Exome Sequencing (WES), RNAseq and methylome. In our cohort, 3 tumors were plotted within the Methylation Class-RGNT (MC-RGNT) characterized by FGFR1 mutation associated with PIK3CA or NF1 mutations. Transcriptome analysis was performed in 7 cases. RNAseq identified a “Hot” and a “Cold” transcriptomic group; the latter includes the 3 MC-RGNT and 1 MC-Pilocytic Astrocytoma and exhibited a “Cold” immune tumor microenvironment in comparison with the “Hot” group. The distinct immune cell content of both groups was confirmed by quantitative immunostainings. Get Set Enrichment Analysis showed that the “Cold” group had upregulated NOTCH pathway and mainly oligodendrocyte precursor cell and neuronal phenotypes while the “Hot” group exhibited predominantly astrocytic and neural stem cell phenotypes. In silico deconvolution identified the cerebellar granule cell lineage as a putative origin and as the cell context to understand the effects of genetic alterations and NOTCH signaling. Histological diagnosis of rosette-forming tumors encompasses heterogeneous tumor entities. Our study identified distinct tumor cell contexts and microenvironments as key features to better understand and manage these neoplasms.

  Study EGAS00001006502

• Exome sequencing from cfDNA blood samples. 159 samples at 2x101bp Illumina reads in Fastq format.

  Differential presence of exons (DPE) by next generation sequencing (NGS) is a method of interpretation of whole exome sequencing. This method has been proposed to design a predictive and diagnostic algorithm with clinical value in plasma from patients bearing colorectal cancer (CRC). The aim of the present study was to determine a common exonic signature to discriminate between different clinical pictures, such as non-metastatic, metastatic and non-disease (healthy), using a sustainable and novel technology in liquid biopsy. Through DPE analysis, we determined the differences in DNA exon levels circulating in plasma between patients bearing CRC vs. healthy, patients bearing CRC metastasis vs. non-metastatic and patients bearing CRC metastasis vs. healthy comparisons. We identified a set of 510 exons (469 up and 41 down) whose differential presence in plasma allowed us to group and classify between the three cohorts. Random forest classification (machine learning) was performed and an estimated out-of-bag (OOB) error rate of 35.9% was obtained and the predictive model had an accuracy of 75% with a confidence interval (CI) of 56.6-88.5. In conclusion, the DPE analysis allowed us to discriminate between different patho-physiological status such as metastatic, non-metastatic and healthy donors. In addition, this analysis allowed us to obtain very significant values with respect to previous published results, since we increased the number of samples in our study. These results suggest that circulating DNA in patient’s plasma may be actively released by cells and may be involved in intercellular communication and, therefore, may play a pivotal role in malignant transformation (genometastasis).

  Study EGAS00001006656

• Single_cell_measurements_to_characterise_B_cell_repopulation_in_SLE_after_rituximab_therapy__a_pilot_study

  Systemic lupus erythematosus (SLE) is a systemic autoimmune disease characterised by autoantibodies against DNA and nuclear proteins. Kidney involvement (lupus nephritis) is one of the most serious manifestations and is the leading determinant of long-term prognosis. Therapies targeting B cells (e.g. rituximab) are used in the treatment of severe lupus, including lupus nephritis. Rituximab is a monoclonal antibody that targets the CD20 surface protein on the surface of B cells. Following rituximab treatment, there is near complete depletion of B cells in peripheral blood. Peripheral blood B cells slowly return over a period of months. The speed of repopulation varies but typically B cells have returned to baseline levels between 9 and 12 months post rituximab therapy. Return of B cells may be associated with increased disease activity, but this is not always the case. For example, some patients remain in clinical remission for several years following one cycle of rituximab, despite B cell repopulation, whilst others relapse as soon as the B cell number reaches the normal range. The reasons for the different response remain unclear. The B cell subpopulations that return following B cell depleting therapy have been previously studied using flow cytometry but have never been characterised using single cell -omic technologies and thus no data on their molecular signatures are available. Previous studies have shown that naïve B cells were the predominant repopulating B-cells, indicating de novo production of B cells. Obtaining an integrated understanding of the molecular programs of these cells at the single cell level will enhance our ability to predict disease outcome, to identify patients with poorer B-cell drug responses, understand disease heterogeneity and design more effective B-cell targeted therapies approach.

  Study EGAS00001006798

• Analysis of T-cell receptor clonotypes in tumor micro-environment identifies shared cancer type-specific signatures

  Despite the conventional view that a truly random V(D)J recombination process should generate a highly diverse immune repertoire, emerging reports suggest that there is a certain bias towards the generation of shared/public immune receptor chains. These studies were performed in viral diseases where public T-cell receptors (TCR) appear to confer better protective responses. Selective pressures generating common TCR clonotypes are currently not well understood but it is believed that they confer a growth advantage. As very little is known about public TCR clonotypes in cancer, here we set out to determine the extent of shared TCR clonotypes in the intra-tumor microenvironments of virus- and non-virus- driven head and neck cancers using TCR sequencing. We report that tumor-infiltrating T-cell clonotypes were indeed shared across individuals with the same cancer type, where the majority of shared sequences were specific to the cancer type (ie viral versus non-viral). These shared clonotypes were not particularly enriched in EBV-associated nasopharynx cancer, but in both cancers, exhibited distinct characteristics, namely shorter CDR3 lengths, restricted V- and J-gene usages and also demonstrated convergent V(D)J recombination. Many of these shared TCRs were expressed in patients with a shared HLA background. Pattern recognition of CDR3 amino acid sequences revealed strong convergence to specific motifs that were unique to each cancer type, suggesting they may be enriched for specificity to common antigens found in the tumor microenvironment. The identification of shared TCRs in infiltrating tumor T-cells not only adds to our understanding of the tumor adaptive immune recognition but could also serve as disease-specific biomarkers and guide the development of future immunotherapies.

  Study EGAS00001005480

• T-cell receptor sequencing of tumor-infiltrating lymphocytes (TIL) in breast cancer

  The spatial distribution of tumor-infiltrating lymphocytes (TIL) predicts breast cancer outcome and response to systemic therapy, highlighting the importance of an intact tissue structure for characterizing patients’ tumors. Here, we present ST-FFPE, a spatial transcriptomics method for the analysis of formalin-fixed paraffin-embedded samples, which opens the possibility of interrogating widely available archival tissue. The method involves extraction, exome capture and sequencing of RNA from different tumor compartments microdissected by laser-capture, and can be used to study the cellular composition of various tumor environments. Focusing on triple-negative breast cancer (TNBC), we characterized T cells, B cells, dendritic cells, fibroblasts and endothelial cells in both stromal and intra-epithelial compartments. We found a highly variable spatial distribution of immune cell subsets among tumors. This analysis revealed that the immune repertoires of intra-epithelial T and B cells were consistently less diverse than those of stromal T and B cells. T-cell receptor (TCR) sequencing confirmed a reduced diversity and revealed higher clonality of intra-epithelial T cells relative to the corresponding stromal T cells. Analysis of the top ten dominant clonotypes in the two compartments showed a majority of shared but also some unique clonotypes both in stromal and intra-epithelial T cells. Hyperexpanded clonotypes were more abundant among intra-epithelial than stromal T cells. These findings validate the ST-FFPE method and suggest an accumulation of antigen-specific T cells within tumor core. Because ST-FFPE is readily applicable for analysis of previously collected tissue samples, it could be useful for rapid assessment of intratumoral cellular heterogeneity in multiple disease and treatment settings.

  Study EGAS00001007125

• IL7-receptor expression is frequent in T-cell acute lymphoblastic leukemia and predicts sensitivity to JAK-inhibition

  T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive hematological malignancy with a dismal prognosis related to refractory/relapsing diseases, raising the need for new targeted-therapies. Activating mutations of the IL7-receptor pathway genes (IL7Rp) play a proven leukemia-supportive role in T-ALL. JAK-inhibitors such as ruxolitinib have recently demonstrated preclinical efficacy. However, prediction markers for sensitivity to JAK-inhibitors are still lacking. Herein, we show that IL7R (CD127) expression is more frequent (~70%) than IL7Rp-mutations in T-ALL (~30%). We compared the so-called non-expressers (no IL7R-expression/IL7Rp-mutation), expressers (IL7R-expression without IL7Rp-mutation) and mutants (IL7Rp-mutations). Integrative multi-omics analysis outlined IL7R-deregulation in virtually all T-ALL subtypes, at the epigenetic-level in non-expressers, genetic-level in mutants, and post-transcriptional level in expressers. Ex-vivo data using primary-derived xenografts support that IL7Rp is functional whenever the IL7R is expressed, regardless of the IL7Rp mutational status. Consequently, ruxolitinib impaired T-ALL survival in both expressers and mutants. Interestingly, we show that expressers displayed ectopic IL7R-expression and IL7Rp-addiction conferring a deeper sensitivity to ruxolitinib. Conversely, mutants were more sensitive to venetoclax than expressers. Overall, combination of ruxolitinib and venetoclax resulted in synergistic effects in both groups. We illustrate the clinical relevance of this association by reporting achievement of complete remission in two patients with refractory/relapsed-T-ALL. This provides proof of concept for translation of this strategy into clinics as bridge to transplant. Altogether, IL7R-expression can be used as a biomarker for sensitivity to JAK-inhibition, thereby expanding the fraction of T-ALL patients eligible to ruxolitinib up to nearly ~70% of T-ALL.

  Study EGAS00001007144

• cfDNA dataset from the urine supernatant of ovarian cancer patients and healthy controls

  Background Ovarian cancer is the deadliest gynecological malignancy worldwide, due to frequent diagnosis at advanced stage. Simple and non-invasive methods for earlier and accurate detection are urgently needed. Here, the presence of tumor-derived DNA in home-collected urine, cervicovaginal self-samples and clinician-taken cervical scrapes of ovarian cancer patients was explored by DNA methylation and copy number analysis. Methods A total of 428 samples (urine, cervicovaginal self-samples, and clinician-taken cervical scrapes) of 110 healthy controls and 54 patients with benign (n=25) or malignant (n=29) ovarian masses were analyzed. Different urine fractions were examined (full void urine, urine supernatant, and urine sediment). All samples were tested for 12 methylation markers by quantitative methylation-specific PCR. Shallow whole-genome sequencing was performed to detect copy number aberrations and verify the presence of tumor-derived DNA. Results Full void urine was most discriminatory between healthy controls and ovarian cancer patients (C2CD4D, p=0.008; CDO1, p=0.022; MAL, p=0.008), followed by cervical scrapes (C2CD4D, p=0.001; CDO1, p=0.004). A significant difference between benign and malignant ovarian masses was only observed for GHSR in the urine sediment (p=0.024). Methylation levels in cervicovaginal self-samples did not discriminate between healthy controls, benign and malignant ovarian masses. Copy number changes were identified in 17% of urine supernatant samples and smaller fragment sizes were seen in urine supernatant samples with a high tumor fraction. Conclusions This study demonstrates the presence of ovarian cancer-derived DNA in home-collected urine. Additional studies are warranted to further explore the clinical applicability of this approach.

  Study EGAS00001007238

• Whole genome sequencing of 108 epileptic patients from CENet cohort

  Epileptic patients have been sequenced as part of the CENet cohort which is composed of patients with Genetic Generalized Epilepsy (GGE) or Non-Acquired Focal Epilepsy (NAFE) collected in CHUM Research Center in Montreal. Patients were diagnosed by epileptologists. The clinical epilepsy phenotype was classified according to the current classification by the International League against Epilepsy (ILAE). Libraries preparation and whole-genome sequencing: gDNA was cleaned up using ZR-96 DNA Clean & ConcentratorTM-5 Kit (Zymo) prior to being quantified using the Quant-iTTM PicoGreen dsDNA Assay Kit (Life Technologies) and its integrity assessed on agarose gels. Libraries were generated using the TruSeq DNA PCR-Free Library Preparation Kit (Illumina) according to the manufacturer’s recommendations. Libraries were quantified using the Quant-iTTM PicoGreen dsDNA Assay Kit (Life Technologies) and the Kapa Illumina GA with Revised Primers-SYBR Fast Universal kit (Kapa Biosystems). Average size fragment was determined using a LabChip GX (PerkinElmer) instrument. The libraries were denatured in 0.05N NaOH and diluted to 8pM using HT1 buffer. The clustering was done on a Illumina cBot and the flowcell was ran on a HiSeq 2500 for 2x125 cycles (paired-end mode) using v4 chemistry and following the manufacturer's instructions. A phiX library was used as a control and mixed with libraries at 0.01 level. Bioinformatics: The Illumina control software was HCS 2.2.58, the real-time analysis program was RTA v. 1.18.64. Program bcl2fastq v1.8.4 was used to demultiplex samples and generate fastq reads. The filtered reads were aligned to reference Homo_sapiens assembly b37. Each readset was aligned to create a .cram file.

  Study EGAS00001007507

• Copy numbers in resectable gastric cancer treated with surgery alone

  Purpose: Gastric cancer is a heterogeneous disease requiring stratification to recognize its underlying biology and variable clinical outcome. The Cancer Genome Atlas (TCGA) recognizes four molecular subgroups: Epstein-Barr virus positive (EBV+) and microsatellite instable (MSI) with superior survival compared to genomically stable (GS) and chromosomal instable (CIN). Since classification criteria to distinguish GS from CIN were not defined, alternative markers such as Lauren’s histopathology or CDH1/p53 immunohistochemistry are commonly applied. Here we compared survival of gastric cancer subgroups determined by four methods. Experimental Design: EBV- and microsatellite stable tumors from a total of 309 patients treated with surgery alone were assigned to subgroups by: (i) specific copy number aberrations as determined by TCGA (GS or CIN), (ii) genome instability index (GII; high or low), (iii) Lauren’s classification (diffuse or intestinal), and (iv) CDH1/p53 immunohistochemistry (IHC; aberrant CDH1, normal-, or aberrant p53 expression). Subgroups were associated with cancer-related survival (CRS). Results: Five-year CRS was 42.0% for diffuse and 49.5% for intestinal type tumors, and 57.8% for GS and 41.6% for CIN tumors. Classification by GII or CDH1/p53 IHC did not correlate with CRS. The combination of TCGA and Lauren classifications resulted in four distinct subgroups. Five-year CRS for GS-intestinal (n=24), GS-diffuse (n=57), CIN-intestinal (n=142) and CIN-diffuse (n=86) was 61.4%, 56.5%, 47.6%, and 31.5%, respectively. Conclusions: TCGA’s GS and CIN have additional prognostic value to Lauren’s diffuse and intestinal subgroups in gastric cancer. GS-intestinal, GS-diffuse, CIN-intestinal and CIN-diffuse constitute biologically distinct entities and are suggested as stratification variables for future studies.

  Study EGAS00001007394

• Evaluation of somatic mutations in urine samples as a non-invasive method for the detection and molecular classification of endometrial cancer

  Purpose: Current diagnostic methods for endometrial cancer lack specificity, leading to many women undergoing invasive procedures. The aim of this study was to evaluate somatic mutations in urine to accurately discriminate endometrial cancer patients from controls. Experimental Design: Overall, 72 samples were analyzed using next-generation sequencing with molecular identifiers targeting 47 genes. We evaluated urine supernatant samples from women with endometrial cancer (n=19) and age-matched controls (n=20). Cell pellets from urine and plasma samples from seven cases were sequenced; further, we also evaluated paired tumor samples from all cases. Finally, immunohistochemical markers for molecular profiling were evaluated in all tumor samples. Results: Overall, we were able to identify mutations in DNA from urine supernatant samples in 100% of endometrial cancers. In contrast, only one control (5%) showed variants at a variant allele frequency (VAF)≥2% in the urine supernatant samples. The molecular classification obtained by using tumor samples and urine samples showed good agreement. Analyses in paired samples revealed a higher number of mutations and VAFs in urine supernatants than in urine cell pellets and blood samples. Conclusions: Evaluation of somatic mutations using urine samples may offer a user-friendly and reliable tool for endometrial cancer detection and molecular classification. The diagnostic performance for endometrial cancer detection was very high, and cases could be molecularly classified using these noninvasive and self-collected samples. Additional multi-center evaluations using larger sample sizes are needed to validate the results and understand the potential of urine samples for the early detection and prognosis of endometrial cancer.

  Study EGAS00001007396

• Regenerative effects of MSC treatment on busulfan-induced small intestine damaged organoids.

  Allogenic hematopoietic stem cell transplantation (HSCT) is a curative treatment for leukemia and a range of non-malignant disorders. The success of the therapy is hampered by occurrence of acute graft-versus-host disease (aGvHD); an inflammatory response damaging recipient organs, with gut, liver, and skin being the most susceptible. Intestinal GvHD injury is often a life-threatening complication in patients unresponsive to steroid treatment. Second-line available therapies are immunosuppressants or mesenchymal stromal/stem cell (MSC) infusions. Data from our institution and others demonstrate rescue of approximately 40-50% of patients suffering from aGvHD with mesenchymal stromal/stem cells (MSCs) and minor side effects. Although promising, the better understanding of MSC mode of action and patient response to MSC-based therapy is essential to improve this lifesaving treatment. Here, we developed a novel 3D co-culture model of human small intestinal organoids and MSCs, which allows to study the regenerative effects of MSCs on intestinal epithelium in a more physiologically relevant setting than existing in vitro systems. Using this model we mimicked chemotherapy-mediated damage of the intestinal epithelium. The treatment with busulfan, the chemotherapeutic commonly used as conditioning regiment before the HSCT, affected pathways regulating EMT, proliferation, and apoptosis in small intestinal organoids, as shown by transcriptomic and proteomic analysis. The co-culture of busulfan treated intestinal organoids with MSCs reversed the effects of busulfan on the transcriptome and proteome of intestinal epithelium, which we also confirmed by functional evaluation of proliferation and apoptosis. Collectively, we demonstrate that our novel in vitro co-culture system is a new valuable tool to facilitate the investigation of the molecular mechanisms behind the therapeutic effects of MSCs on damaged intestinal epithelium.

  Study EGAS00001007432

• ExomeSeq Neoantigen Immunogenicity Landscapes

  Immune checkpoint blockade (ICB) therapy is a cornerstone of oncologic treatment for patients with advanced stage non-small cell lung cancer (NSCLC) and other malignancies. Neoantigen immunoediting drives ICB efficacy, yet the fundamental physiochemical characteristics of neoantigens and how neoantigen immunogenicity shapes treatment response remains poorly understood. To help address these questions, a prospective clinical trial of NSCLC patients treated with nivolumab was conducted. We assessed genomic alterations in tumors from 58 patients and performed large-scale neoepitope immunogenicity analyses, before and during treatment (CheckMate153, CA209-153). Tumors were analyzed by whole-exome and transcriptome sequencing. In responding patients, loss of mutation and neoantigen burden early during therapy associated with clinical benefit. We evaluated the immunogenicity of 1,453 candidate neoantigens and identified 502 neopeptides that bound to MHC I and 196 neopeptides that were immunogenically recognized by T cells in the setting of nivolumab treatment. These T cell reactive neoantigens were differentially present in clonal populations that underwent distinctive evolutionary trajectories across responders and nonresponders. Mapping these neoantigens to tumor clonal dynamics and clinical response revealed strong selection against immunogenic neoantigen harboring clones compared to non-immunogenic clones. Using this large collection of neoantigens, we identified position specific amino acid features related to immunogenicity, which we used to develop and validate an immunogenicity score. Changes in the genomic and neoantigen immunogenicity landscapes were accompanied by temporal changes in the tumor microenvironment. Nivolumab-induced microenvironmental evolution in NSCLC shared some similarities with that in melanoma, yet critical differences in immunologic programs were apparent from comparative network analysis between tumor types. This study provides unprecedented molecular portraits of the neoantigen landscapes underlying nivolumab’s mechanism of action.

  Study EGAS00001007508

• RNASeq Neoantigen Immunogenicity Landscape

  Immune checkpoint blockade (ICB) therapy is a cornerstone of oncologic treatment for patients with advanced stage non-small cell lung cancer (NSCLC) and other malignancies. Neoantigen immunoediting drives ICB efficacy, yet the fundamental physiochemical characteristics of neoantigens and how neoantigen immunogenicity shapes treatment response remains poorly understood. To help address these questions, a prospective clinical trial of NSCLC patients treated with nivolumab was conducted. We assessed genomic alterations in tumors from 58 patients and performed large-scale neoepitope immunogenicity analyses, before and during treatment (CheckMate153, CA209-153). Tumors were analyzed by whole-exome and transcriptome sequencing. In responding patients, loss of mutation and neoantigen burden early during therapy associated with clinical benefit. We evaluated the immunogenicity of 1,453 candidate neoantigens and identified 502 neopeptides that bound to MHC I and 196 neopeptides that were immunogenically recognized by T cells in the setting of nivolumab treatment. These T cell reactive neoantigens were differentially present in clonal populations that underwent distinctive evolutionary trajectories across responders and nonresponders. Mapping these neoantigens to tumor clonal dynamics and clinical response revealed strong selection against immunogenic neoantigen harboring clones compared to non-immunogenic clones. Using this large collection of neoantigens, we identified position specific amino acid features related to immunogenicity, which we used to develop and validate an immunogenicity score. Changes in the genomic and neoantigen immunogenicity landscapes were accompanied by temporal changes in the tumor microenvironment. Nivolumab-induced microenvironmental evolution in NSCLC shared some similarities with that in melanoma, yet critical differences in immunologic programs were apparent from comparative network analysis between tumor types. This study provides unprecedented molecular portraits of the neoantigen landscapes underlying nivolumab’s mechanism of action.

  Study EGAS00001007509

• Dynamics of circulating tumor DNA in acute myeloid leukemia (AML) patients who undergo allogeneic transplantation

  Allogeneic hematopoietic cell transplantation (allo-HCT) is the only curative option for many patients with acute myeloid leukemia (AML). With the increasing use of non-myeloablative conditioning regimens, it is widely assumed that the cytotoxic conditioning itself does not eradicate all leukemia cells. Instead, the donor-derived immune cells are believed to eliminate residual leukemia cells over time. This so-called graft-versus-leukemia (GvL) is postulated to be the primary mechanism by which allo-HCT achieves sustained remissions in patients with AML. However, the precise quantification of GvL effects and its contribution to leukemia eradication has not yet been demonstrated by sequential, highly sensitive measurements of residual disease at various time points following allo-HCT. In the current study, we employed a novel, highly sensitive and specific plasma-based assay to track mutations in the leukemia cells of 14 patients before and after allo-HCT. If the primary therapeutic benefit of allo-HCT is truly mediated through a GvL effect, we hypothesized that leukemia cells would persist following cytotoxic conditioning and would often be reduced after post-transplant immunosuppression was discontinued. In the majority of patients (11 of 14, 79%), the conditioning used for allo-HCT did not eradicate the leukemia cells. In seven of these 11 patients, the leukemia cells persisted until immunosuppression was discontinued, at which time the persistent leukemia cells were reduced or eradicated. These results demonstrate the feasibility of using circulating tumor DNA (ctDNA) as a dynamic marker of GvL as well as minimal residual disease for patients with AML undergoing alloHCT, and have substantial implications for their management.

  Study EGAS00001007969

• MethylationEPIC BeadChip samples of pre- and post-5ZA-treated head and neck cancer patients refractory to anto-PD-1 therapy

  Background Immune checkpoint blockade (ICB) is an effective treatment in a subset of patients diagnosed with head and neck squamous cell carcinoma (HNSCC); however, most patients are refractory. Methods In a non-randomized, open-label phase 1b clinical trial, participants with recurrent and/or metastatic (R/M) HNSCC were treated with low-dose 5’azacytidine (5’aza) daily for either 5- or 10-days in combination with fixed doses of durvalumab and tremelimumab after progression on ICB therapy. The primary objective was to assess the biologically effective dose of 5’aza as determined by molecular changes in paired baseline and on-treatment tumor biopsies; the secondary objective was safety. Results Thirty-eight percent (3/8) of participants with evaluable paired tissue samples had a greater than 2-fold increase from baseline in interferon-gamma (IFN-ɣ) signature and CD274 (programmed cell death protein 1 ligand, PD-L1) expression within the tumor microenvironment (TME), which was associated with increased CD8+ T cell infiltration and decreased infiltration of CD4+ T regulatory cells. The mean neutrophil-to-lymphocyte ratio (NLR) decreased by greater than 50%, from 14.2 (SD 22.6) to 6.9 (SD 5.2). Median overall survival (OS) was 16.3 months (95% CI 1.9, NA) and 2-year OS rate was 24.7% (95% CI: 4.5%, 53.2%) with 58% (7/12) of treated participants demonstrated prolonged OS of greater than 12 months. Conclusion Our findings suggest that low-dose 5’aza can reprogram systemic host immune responses and the local TME to increase IFN-ɣ and PD-L1 expression. The increased expression of these established biomarkers correlated with prolonged OS upon re-challenge with ICB.

  Study EGAS00001007998

• Influence of age on molecular changes and treatment stratification in multimodal glioblastoma therapy

  Background: The increasing life expectancy has led to a rise in high-grade gliomas, particularly in patients of older age. The impact of chronological age, associated molecular changes and contributing factors on therapy response, progression-free survival (PFS) and overall survival (OS) remain uncertain. Methods: Between January 1st, 2016, and December 31th, 2021, patients with glioblastoma IDH-wt (CNS WHO grade 4) were included. RNA- and T-cell receptor sequencing was performed in a subset of patients. Primary outcome parameters were extent of resection (EOR), progression free- and overall survival. Multivariate analysis was performed to identify independent predictors of survival. Results: A total of 286 patients were included. The mean age was 66.4 years. Chronological age did not correlate with overall survival (r=0.12). With age, tumor immune composition changed towards macrophages with M2 phenotype and decreased numbers of CD8+ and neutrophiles. Aging related genes were not enriched. Beyond 69.5 years of age, despite similar PFS (p=0.0783), OS (p<0.0001) was significantly shorter, compared to younger patients. Extent of resection and adjuvant radiochemotherapy were significant predictors of OS (p<0.0001). No difference in PFS and OS was found between younger and older patients treated with maximal safe resection, followed by concomitant RT and chemotherapy and adjuvant chemotherapy. Conclusions: This study suggests that age alone should not be the sole determinant in therapeutic decision-making for glioblastoma patients. Gross total resection and adjuvant therapy, tailored to individual clinical status, significantly impact survival outcomes. Patient-specific treatment strategies are important to minimize undertreatment and optimize therapeutic outcomes in older patients.

  Study EGAS00001008246

• Next Generation Children project

  Whole genome sequencing of sick children in neonatal and paediatric intensive care units, aligned to reference assembly GRCh38.

  Dataset EGAD00001007868

• EPICC: Evolutionary Predictions in Colorectal Cancer WGS-3

  WGS profiling bam files from colorectal carcinoma and adenoma.

  Dataset EGAD00001008642

• Clonal dynamics after allogeneic haematopoietic cell transplantation using genome-wide somatic mutations - WGS

  Allogeneic haematopoietic cell transplantation (HCT) replaces the stem cells responsible for blood production with those harvested from a donor, and is received by 40,000 patients worldwide each year. To quantify dynamics of long-term stem cell engraftment, we sequenced whole genomes of 2,824 single-cell-derived haematopoietic colonies from blood samples of 10 donor-recipient pairs taken 9-31 years after HLA-matched sibling HCT. With younger donors, 10,000-50,000 stem cells had engrafted and were still contributing to haematopoiesis at time of sampling, but estimates were 10-fold lower with older donors. Engrafted stem cells made multilineage contributions to myeloid, B-lymphoid and T-lymphoid populations, although individual clones often showed biases towards one or other mature cell type. Recipients had lower clonal diversity than matched donors, equivalent to ~10-15 years of additional ageing, arising from up to 25-fold greater expansion of stem cell clones. An HCT-related population bottleneck alone could not explain these differences: instead, phylogenetic trees evinced two distinct modes of HCT-specific selection. In 'pruning selection', cell divisions underpinning recipient-enriched clonal expansions had occurred in the donor, preceding transplant - their selective advantage derived from preferential mobilisation, harvest, survival ex vivo or initial homing. In 'growth selection', cell divisions underpinning clonal expansion occurred through proliferative advantage in the recipient's marrow after homing - clones with multiple driver mutations especially demonstrated this pattern. Uprooting stem cells from their native environment and transplanting them to foreign soil exaggerates selective pressures, distorting and accelerating the loss of clonal diversity compared to the unperturbed haematopoiesis of donors.

  Dataset EGAD00001010872

• Single cell RNA sequencing of synovial B cells in early Rheumatoid Arthritis

  This dataset includes 3432 paired single cell sequencing fastq files derived from synovial B cells of 5 early Rheumatoid Arthritis patients. Libraries were prepared using the SmartSeq2 protocol. All wells contained ERCC spike-ins. Libraries were sequenced on an Illumina NovaSeq instrument with 2x100bp paired-end reads yielding a median of 2.5M reads/well. Sample aliases ending on 368 and 384 represent empty control wells. Sample aliases starting with P11417_4 and P13157_1 belong to ACPA- patient A7, sample aliases starting with P11417_5 and P11417_6 belong to ACPA+ patient A1, sample aliases starting with P13157_2 belong to ACPA- patient A3, sample aliases starting with P13157_3 and P13157_6 belong to ACPA+ patient A2, sample aliases starting with P13157_4 and P13157_5 belong to ACPA- patient A4.

  Dataset EGAD00001009076

• A catalog of the genetic causes of Hereditary Angioedema in the Canary Islands (Spain)

  Hereditary angioedema (HAE) is a rare disease where known causes involve C1 inhibitor dysfunction or dysregulation of the kinin cascade. The updated HAE management guidelines recommend performing genetic tests to reach a precise diagnosis. Unfortunately, genetic tests are still uncommon in the diagnosis routine. Here, we characterized for the first time the genetic causes of HAE in affected families from the Canary Islands (Spain). Whole-exome sequencing data was obtained from 41 affected patients and unaffected relatives from 29 unrelated families identified in the archipelago. The Hereditary Angioedema Database Annotation (HADA) tool was used for pathogenicity classification and causal variant prioritization among the genes known to cause HAE. Manual reclassification of prioritized variants was used in those families lacking known causal variants. We detected a total of eight different variants causing HAE in this patient series, affecting essentially SERPING1 and F12 genes, one of them being a novel SERPING1 variant (c.686-12A>G) with a predicted splicing effect which was reclassified as likely pathogenic in one family. Altogether, the diagnostic yield by assessing previously reported causal genes and considering variant reclassifications according to the American College of Medical Genetics guidelines reached 66.7% (95% Confidence Interval [CI]: 30.1-91.0) in families with more than one affected member and 10.0% (95% CI: 1.8-33.1) among cases without family information for the disease. Despite the genetic causes of many patients remain to be identified, our results reinforce the need of genetic tests as first-tier diagnostic tool in this disease, as recommended by the international WAO/EAACI guidelines for the management of HAE.

  Study EGAS00001006547

• Single-cell decoding of drug induced transcriptomic reprogramming in triple negative breast cancers

  Background The encoding of cell intrinsic drug resistance states in breast cancer reflects the contributions of genomic and non-genomic variations and requires accurate estimation of clonal fitness from co-measurement of transcriptomic and genomic data. Somatic copy number (CN) variation is the dominant mutational mechanism leading to transcriptional variation and notably contributes to platinum chemotherapy resistance cell states. Here, we deploy time series measurements of triple negative breast cancer (TNBC) single-cell transcriptomes, along with co-measured single-cell CN fitness, identifying genomic and transcriptomic mechanisms in drug associated transcriptional cell states. Results We present scRNA-seq data (53,641 filtered cells) from serial passaging TNBC patient-derived xenograft (PDX) experiments spanning 2.5 years, matched with genomic single-cell CN data from the same samples. Our findings reveal distinct clonal responses within TNBC tumors exposed to platinum. Clones with high drug fitness undergo clonal sweeps and show subtle transcriptional reversion, while those with weak fitness exhibit dynamic transcription upon drug withdrawal. Pathway analysis highlights convergence on epithelial-mesenchymal transition and cytokine signaling, associated with resistance. Furthermore, pseudotime analysis demonstrates hysteresis in transcriptional reversion, indicating generation of new intermediate transcriptional states upon platinum exposure. Conclusions Within a polyclonal tumor, clones with strong genotype-associated fitness under platinum remained fixed, minimizing transcriptional reversion upon drug withdrawal. Conversely, clones with weaker fitness display non-genomic transcriptional plasticity. This suggests CN-associated and CN-independent transcriptional states could both contribute to platinum resistance. The dominance of genomic or non-genomic mechanisms within polyclonal tumors has implications for drug sensitivity, restoration and re-treatment strategies.

  Study EGAS00001007242

• H3Africa - Genomic Characterization and Surveillance of Microbial Threats in West Africa

  Viral diseases are remarkably heterogeneous in their presentation and clinical course. Although diseases such as Lassa Fever and Ebola Virus Disease are thought to be severe and rapidly fatal, they can often present with mild symptoms or remain asymptomatic. Host genetic variability is likely an important mediator of this heterogeneity. Complete characterization of the clinical, laboratory, and host genetic features of these illnesses is not only imperative to rapid diagnosis and effective clinical care, but also important for further scientific study. In this project, we propose to address this need through deep characterization of the clinical and laboratory features of viral disease and through study of genetic mediators of resistance to Lassa Fever and Ebola Virus Disease. Will carry out deep clinical characterization and laboratory assessment of viral febrile cases at two West African Collaborative Centers, and will use this data to develop prognostic models using cutting-edge tools in exploratory data visualization and machine learning. Investigate host genetic factors of Lassa Fever and Ebola Virus Disease through genome-wide association studies. We will then develop field-deployable genotyping assays for the top-associated variants, and apply them in new cohorts to replicate our findings and investigate patterns of genetic resistance in clinical subgroups. This project can shed light on the clinical manifestations of viral illness and their host genetic mediators, leading to potential new insights in disease biology and pathophysiology. We will also create a rich resource of clinical, laboratory, and genetic data for these illnesses in two countries in West Africa, which will serve as a foundation for advancing the clinical care and scientific study of viral disease in the region.

  Study EGAS00001007250

• Tubulointerstitial fibrosis is the histological hallmark of chronic kidney disease (CKD). Hypoxia and inflammation (i.e., interleukin (IL)-1β signalling) are independent mediators of tubulointerstitial fibrosis. However, the physiological response of human kidney tubular cells to IL-1β/IL-1RI signalling under the hypoxic conditions of CKD is poorly understood and remains a clinical imperative for therapeutic targeting. This study reports that hypoxia and IL-1β act in synergy to trigger cell cycle arrest/cellular senescence of ex vivo patient-derived primary proximal tubular epithelial cells (PTECs).

  Hypoxia and interleukin (IL)-1β are independent mediators of tubulointerstitial fibrosis, the histological hallmark of chronic kidney disease (CKD). Here, we examine how hypoxia and IL-1β act in synergy to augment maladaptive proximal tubular epithelial cell (PTEC) repair in human CKD. Ex vivo patient-derived PTECs were cultured under normoxic (21% O2) or hypoxic (1% O2) conditions in the absence or presence of IL-1β and examined for maladaptive repair signatures. Hypoxic PTECs incubated with IL-1β displayed a discrete transcriptomic profile distinct from PTECs cultured under hypoxia alone, IL-1β alone or under normoxia. Hypoxia+IL-1β-treated PTECs had 692 ‘unique’ differentially expressed genes (DEGs) compared to normoxic PTECs, with ‘cell cycle’ the most significantly enriched KEGG pathway based on ‘unique’ down-regulated DEGs (including CCNA2, CCNB1 and CCNB2). Hypoxia+IL-1β-treated PTECs displayed signatures of cellular senescence, with reduced proliferation, G2/M cell cycle arrest, elevated senescence-associated β-galactosidase (SA-β-gal) activity and increased production of pro-inflammatory/fibrotic senescence-associated secretory phenotype (SASP) factors compared to normoxic conditions. Treatment of Hypoxia+IL-1β-treated PTECs with either a senolytic drug combination, quercetin+dasatinib, or a type I IL-1 receptor (IL-1RI) neutralizing antibody attenuated senescent cell burden. Our data identify a mechanism whereby hypoxia in combination with IL-1β/IL-1RI signalling trigger PTEC senescence, providing novel therapeutic and diagnostic check-points for restoring tubular regeneration in human CKD.

  Study EGAS00001007904

• Identification of intronic sequences promoting natural skipping of ABCA4 exon 15 using long-read transcriptomics and midigene assays

  Stargardt disease is an inherited retinal disease that leads to progressive loss of vision, and is caused by biallelic variants in the ABCA4 gene. A significant proportion of ABCA4 variants of variable penetrance affect pre-mRNA splicing. Previously, we observed the exclusion of certain ABCA4 exons into the final transcript in the absence of any disease-causing variants. In this study, we explored this phenomenon that we coined natural exon skipping (NES) for the entire ABCA4 transcript, using targeted RT-PCR analysis as well as long-read transcriptome sequencing of control retina samples. For the most abundant NES event, which is skipping of ABCA4 exon 15, overlapping ABCA4 deletion constructs were employed to identify responsible cis-regulatory sequence motifs. Subsequently, antisense oligonucleotides (AONs) were designed and tested to reverse NES and restore fully-spliced ABCA4 mRNA. We identified intronic regions on each side of exon 15 that, when deleted, increased the inclusion of exon 15 in the final transcripts. Interestingly, these motifs appear to work in a synergistic manner. Further fine-mapping of these regions revealed several cis-regulatory motifs that may be responsible for inducing NES, on each side of the exon. AONs targeting these motifs unfortunately did not decrease the skipping of exon 15. In conclusion, we discovered an intriguing phenomenon that, in the absence of disease-causing variants, results in the skipping of one or more exons in the ABCA4 transcript. These NES events appear to be regulated by cis-regulatory motifs within the flanking introns. This work captured the extend of NES throughout ABCA4 transcript and the role of cis-regulatory motifs in pre-mRNA splicing.. The use of AONs targeting cis-regulatory motifs is an innovative strategy for splicing modulation, suitable for complementing future variant-specific splicing modulation therapies in ABCA4 retinopathies.

  Study EGAS50000000071

• Blueprint Data Access Committee.

  Access to the sequence and alignment and genotype level data produced by the Blueprint Consortium is controlled by the Blueprint Data Access Committee (DAC). Access to data will be granted to qualified investigators for appropriate use. Please email blueprint-dac@ebi.ac.uk to request the Data Access Application and Data Access Agreement.

  Dac EGAC00001000135

• Integrating molecular imaging and transcriptomic profiling in advanced HER2-positive breast cancer receiving trastuzumab emtansine (T-DM1): an analysis of the ZEPHIR clinical trial

  Purpose: The ZEPHIR clinical trial evaluated the role of [89Zr]trastuzumab-PET/CT (HER2-PET/CT) and 2-[18F]fluoro-2-deoxy-D-glucose PET/CT ([18F]FDG-PET/CT) in predicting outcomes in patients with advanced HER2-positive breast cancer treated with trastuzumab emtansine (T-DM1). Here, we combined molecular/metabolic imaging and transcriptomic data to investigate the biological processes associated with [89Zr]trastuzumab and [18F]FDG uptake, and to dissect the mechanisms involved in T-DM1 resistance. Experimental design: RNA was extracted from metastasis biopsies obtained at study inclusion in the ZEPHIR trial. HER2-PET/CT and [18F]FDG-PET/CT imaging data of biopsied lesions were integrated with RNA sequencing data. Lesions were compared based on the level of [89Zr]trastuzumab uptake (dichotomized into high/low classes) as well as on the presence/absence of metabolic response, defined comparing baseline and on-treatment [18F]FDG-PET/CT. Results: We analyzed matched transcriptomic and molecular/metabolic imaging data for 24 metastases (one lesion/patient). Pathways involved in extracellular matrix (ECM) organization, glycosylphosphatidylinositol synthesis, and collagen formation were enriched in lesions presenting low [89Zr]trastuzumab uptake. [18F]FDG uptake at baseline correlated with proliferation and immune-related processes. Hypoxia and ECM-related processes were enriched in lesions that showed no metabolic response to T-DM1. Interestingly, immune-related processes (e.g., Fc gamma receptor activation, complement cascade) were associated with high [89Zr]trastuzumab uptake and metabolic response. Conclusions: To our knowledge, this study represents the first correlative analysis between [89Zr]trastuzumab tumor uptake and gene expression profiling in humans. Our findings point toward a role of ECM in impairing [89Zr]trastuzumab tumor uptake and T-DM1 metabolic response in advanced HER2-positive breast cancer, and highlight the potential of molecular imaging to depict tumor microenvironment features.

  Study EGAS50000000470

• Genomic Profiling of Pediatric Tumors by Cell Free DNA Sequencing

  BACKGROUND Analysis of cell-free DNA (cfDNA) from peripheral blood has been studied extensively for detection and monitoring therapy response in adult cancers, but utility in pediatric tumors is unclear. METHODS We studied 233 children with hematologic, solid and brain tumors. cfDNA was extracted from 1 mL of plasma collected at diagnosis. Circulating tumor DNA (ctDNA) in cfDNA was quantified and diverse classes of somatic genomic variants were detected using PeCan-Seq. PeCan-Seq findings were compared with matched tumor whole-genome and whole-exome sequencing data. Minimal residual disease (MRD) status was evaluated by cfDNA PeCan-Seq in 14 patients with B-cell acute lymphoblastic leukemia. RESULTS Plasma cfDNA yields at diagnosis were higher for hematologic malignancies than for solid and brain tumors. Patients with hematologic malignancy were 100% ctDNA-positive, the median ctDNA fraction in cfDNA was 0.77. We detected 97% variants that were present in tumors and covered by the capture panel, and 50 subclonal variants that were present in cfDNA only. Our approach enabled detection of sequence and structural variants in B-ALL ctDNA at varying levels of MRD. Nineteen of 38 patients with solid tumors were ctDNA-positive, with high prevalence in Ewing Sarcoma, liver tumors and osteosarcoma. One of 19 brain tumor patients was ct-DNA positive, indicating lack of utility for plasma-derived cfDNA in this disease. Notably, cfDNA analysis of a single plasma sample identified genomic features in ctDNA derived from B-ALL, concurrent primary neuroblastoma, and a later neuroblastoma relapse, indicating that genomic variations from multiple tumors could be detected simultaneously. CONCLUSIONS Pediatric patients with hematologic malignancy had an abundance of plasma ctDNA at diagnosis. Plasma ctDNA levels varied between different solid tumor types, more studies are needed to determine utility. Panel-based cfDNA sequencing could serve as a non-invasive approach for pediatric cancer diagnosis and disease monitoring.

  Study EGAS50000000072

• Integrative analysis reveals a macrophage-predominant, immunosuppressive immune microenvironment and subtype-specific therapeutic vulnerabilities in advanced salivary gland cancer

  Background Salivary gland cancers (SGC) are rare and heterogeneous malignant tumors. Advanced SGC lack established treatment options and show poor response to immunotherapy. Here, an integrative multi-omics analysis in a large cohort of advanced SGC revealed insights into the tumor immune microenvironment (TIM) and distinct mechanisms of immune evasion. Methods A total of 104 patients with recurrent/metastatic SGC from the DKTK MASTER program were included in this study. Whole-exome or whole-genome sequencing and RNA-sequencing was performed on fresh frozen tumor tissue. The tumor immune microenvironment was analyzed using CIBERSORT deconvolution analysis and immune gene expression scores in bulk RNA-sequencing data. Single-nuclei sequencing and immunohistochemistry analyses were performed in selected samples. Results were validated in bulk RNA-sequencing data of a previously published independent dataset. Results Bulk transcriptome analysis revealed an immune-deserted TIM in the majority of advanced SGC samples. Immune exclusion was most prominent in adenoid cystic carcinoma (ACC) subgroup 1 exhibiting a downregulation of the antigen processing machinery. Only a small subset of advanced SGC, including few adenoid cystic carcinoma, exhibited T-cell inflammation, which was correlated with tumor mutational burden in Non-ACC samples. Subtype specific expression of immune checkpoints as well as cancer testis antigens were identified with prominent expression of VTCN1 in luminal cells within ACC. Single-cell RNA-sequencing and bulk RNA-seq deconvolution analysis validated immune cell exclusion and revealed a TIM that was dominated by M2 macrophages across SGC subtypes. Among evaluable patients treated with immune checkpoint inhibitors, a high T-cell to macrophage ratio was associated with clinical benefit. Conclusions These data support biomarker-based development of immune-checkpoint inhibition and the development of novel immune-checkpoint inhibitors and cellular therapies in SGC.

  Study EGAS50000000809

• A prospective study identifies MisMatch Repair genes as candidate predisposing genes for uveal melanoma.

  Objective: The aim of this study is to identify new predisposing genes for uveal melanoma (UM), the most common intraocular malignancy in adults. Design: Prospective cohort study from Institut Curie. Participants: 381 unselected patients diagnosed with UM between July 2021 and February 2023 at Institut Curie. Methods: All patients received genetic counseling and consented to extended genetic testing. A panel of 122 genes was analyzed by targeted sequencing on germline DNA from these patients. Main outcomes measures: UM was diagnosed by ophthalmologists at Institut Curie, and a senior pathologist confirmed the diagnosis when tumor or biopsy was available. Classification of the germline genetic variants was performed according to the ACMG guidelines. Results: We identified at least one pathogenic variant (PV) in 79 patients. Among those variants, 20 were found in clinically relevant genes for the patient’s follow-up and genetic counseling. We found a significant enrichment of germline pathogenic variants in Mismatch Repair (MMR) MLH1, MSH2 and MSH6 genes in this UM cohort, suggesting that MMR germline alterations could predispose to UM. From this cohort, one tumor was available from a patient carrying a MLH1 germline mutation. The tumor was characterized by a monosomy 3, in line with a loss of heterozygosity (LOH) of MLH1, located on chromosome 3. Loss of expression of MLH1 and PMS2 was observed by immunohistochemistry, and MMR mutational signatures SBS6, ID1 and ID2 were identified from the whole genome sequencing of the tumor. These data suggest that MLH1 is implicated in the oncogenesis of this UM. Conclusions: This study is the largest prospective constitutional study performed on UM patients, and the first demonstrating that MMR germline alterations are significantly enriched among these patients. Further analyses will be needed to estimate the risk of developing UM for patients with Lynch syndrome

  Study EGAS50000000914

• Emirati Genome Project subset of 43,608 WGS samples for population-scale variant discovery and allele frequency mapping.

  Here, we present a comprehensive genomic characterization of a cohort of 43,608 Emirati genomes sequenced as part of the Emirati Genome Program (EGP). This study identified more than 421 million single nucleotide variants and indels and more than 600 million copy number- and structural variants. Small variants had 756 million molecular effects annotated. Of the identified variants, 69.04% were observed with more than one allele. Of 7.7 million polymorphic variants having an allele frequency (AF) of more than 5% in EGP, 1,348 have a predicted deleterious effect on a protein. Characterization of global variation reveals that EGP represents a genetic continuum encompassing the range of African, Asian, and European populations. It is best described by two Arabian, an Eurasian, and an African component, with the predominant Arabian component linked to mitochondrial haplogroups J and T, which are commonly attributed to the Middle East. Various aneuploidies of sex chromosomes were detected in 93 individuals overall and of chromosome 21 in 41 individuals. The median inbreeding coefficient and cumulative runs of homozygosity (ROHs) lengths were increased due to extensive consanguinity, being largest in the groups with Arabian main ancestry components and higher than reported for Qatar. Families were identified based on genetic relatedness and classified into 264 families with unrelated parents and 247 families with third- and fourth-degree consanguineous parents. Representative consanguineous pedigrees of families in EGP were outlined. Cumulative ROHs were affected by the main ancestry component and significantly increased in offspring of consanguineous parents, with a pronounced difference between 3rd and 4th-degree relatedness. Investigation of cumulative AFs of variants causing Mendelian diseases highlighted genes related to alpha- and beta-thalassemia (HBB, HBA2). It showed a high burden of variants causing severe recessive diseases, metabolic and retinal disorders, and hearing loss. In summary, EGP represents a landmark effort in characterizing the genetic diversity of the Emirati population, leveraging the largest Middle Eastern cohort reported to date.

  Study EGAS50000001071

• Single-Gene vs. Panel Sequencing in Advanced HR+/HER2− Breast Cancer

  Background: Hormone receptor-positive (HR+)/HER2-negative (HER2−) breast cancer is the most common subtype of breast cancer, with biomarker-driven therapies improving outcomes. Liquid biopsies, which analyze circulating tumor DNA (ctDNA), are increasingly essential in clinical management, offering a minimally invasive approach to assess somatic alterations and guide therapy decisions. However, the choice of assay directly affects clinical utility and global access to validated testing platforms remains inconsistent. To address these challenges, this study aimed to compare the performance of single-gene and panel-based sequencing for detecting PIK3CA mutations in HR+/HER2− breast cancer and evaluate the added value of broader genomic profiling in guiding treatment decisions. Patients and methods: We conducted a prospective, multicenter study comparing our previously established SiMSen-Seq (SSS) assay, a single-gene method targeting PIK3CA hotspot mutations, with the AVENIO ctDNA Expanded assay, enriching for 77 clinically relevant cancer genes. Additionally, tumor fractions were estimated with mFAST-SeqS by assessing chromosomal aneuploidy. We analyzed 161 plasma samples collected from 146 patients with HR+/HER2− advanced breast cancer before starting a new line of palliative therapy. Results: High concordance (92.6%) was observed between SSS and AVENIO for PIK3CA mutation detection, with accurate variant allele frequency (VAF) quantification confirmed by Bland-Altman analysis and linear mixed-effects modeling. AVENIO identified additional PIK3CA mutations in 19% of samples as well as other actionable alterations, including ESR1 mutations (17.5%) and PI3K pathway alterations (40.6%). Panel sequencing improved the interpretation of negative liquid biopsy results by integrating tumor fraction and VAF metrics. Conclusions: Panel sequencing assays such as AVENIO provide broader genomic coverage, enhance the detection of actionable mutations and refine the interpretation of negative liquid biopsy findings. These findings support the integration of comprehensive ctDNA profiling into routine clinical practice for HR+/HER2− advanced breast cancer while highlighting the need for improved assay accessibility and expanded genomic targets.

  Study EGAS00001008200

• CDK4/6 inhibition in advanced chordoma: final results of the NCT PMO-1601

  Background. Chordoma is a rare bone tumor with no approved systemic therapy. Inactivation of the tumor suppressor p16 (encoded by CDKN2A) leads to aberrant cyclin D-CDK4/6-RB pathway activity, which can be inhibited by palbociclib in chordoma cell lines. Patients and Methods. We conducted a phase II single-arm, open-labeled trial on palbociclib in adult patients with advanced chordomas with p16 loss (by immunohistochemistry) or CDKN2A loss (by genomic analysis) and CDK4/6 and RB1 presence (by immunohistochemistry or RNA-Sequencing). Patients received 125mg palbociclib once daily for 21 days in a 28-day cycle. The study used a Simon optimal 2-stage design. The primary endpoint was disease control rate after six cycles. Response was evaluated according to Response Evaluation Criteria in Solid Tumors (RECIST) version 1.1. Secondary endpoints included progression-free survival, overall survival, and biomarker analysis. The study was considered positive if 25% of patients reached the primary endpoint. Results. Between December 2017 and April 2023, 28 patients (median age 60 years, range 31-84) were enrolled. After a median follow-up of 28 months, the disease control rate was 39%, with 11 out of 28 patients achieving stable disease. No objective responses were observed. The median progression-free survival was 5.6 months, and the median overall survival was 24.6 months. Treatment was well tolerated without new safety signals. CDKN2A deletions were the most common genomic alterations underlying p16 loss and were accompanied by MTAP co-deletion. There was no correlation between responder phenotypes by immunohistochemistry and outcome. RNA-based analysis showed heterogeneous pathways deregulation at study entry suggesting underlying primary resistance mechanisms Conclusion. This trial met its primary endpoint with a 39% disease control rate, though antitumor activity was modest. Molecularly tailored combination therapies should be considered to improve efficacy. (ClinicalTrials.gov number, NCT03110744).

  Study EGAS00001007985

• McGill Epigenomics Mapping Centre

  The McGill Epigenomics Mapping Centre (EMC) and Data Coordination Centre (EDCC) were established in 2012 at the McGill University and Genome Quebec Innovation Centre (Montreal, Canada) to support large-scale human epigenome mapping for a broad spectrum of cell types and diseases. Previous epigenetic studies have provided profiles of specific chromatin marks in relation to basic biological processes, and future studies of molecular mechanisms must build upon these proofs-of-concept by integrating sequence-based variation with multiple levels of epigenetic and transcriptional regulation across the genome of human tissues and animal disease models. This project leverages the high-throughput sequencing infrastructure and expertise in genomics and transcriptomics at the McGill Innovation Centre to carry out this research. Controlled access to the data by collaborators and the greater scientific community is gained via a portal, which takes advantage of Compute Canada high-performance computing cluster resources to manage the large volume of data associated from the generation of reference epigenome maps. McGill University is one of two Canadian mapping and data coordination centres, the other based at the Michael Smith Genome Sciences Centre in Vancouver, British Columbia. The generation of comprehensive epigenome maps at McGill University is part of a larger international effort that is coordinated by the International Human Epigenome Consortium (IHEC), whose overall long-term objective is to determine the extent to which the epigenome shapes human populations over generations and in response to the environment. This project is funded under the Canadian Epigenetics, Environment, and Health Research Consortium (CEEHRC) by the Canadian Institutes of Health Research and by Genome Quebec, with additional support from Genome Canada. The computing and networking infrastructure, and part of the software development, are provided by Compute Canada and CANARIE.

  Study EGAS00001000995