23014 results for "ega"

in 51.37 milliseconds.

• 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

• 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

• 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

• DNA methylation-based prognostic subtypes of chordoma tumors in tissue

  Chordomas are rare malignant bone cancers of the skull-base and spine. Patient survival is variable and not reliably predicted using clinical factors or molecular features. This study identifies prognostic epigenetic chordoma subtypes that are detected noninvasively using plasma methylomes. Methylation profiles of 68 chordoma surgical samples were obtained between 1996 and 2018 across three international centers along with matched plasma methylomes where available. Consensus clustering identified two stable tissue clusters with a disease-specific survival difference that was independent of clinical factors in a multivariate Cox analysis (HR = 14.2, 95%CI: 2.1–94.8, P = 0.0063). Immune-related pathways with genes hypomethylated at promoters and increased immune cell abundance were observed in the poorperforming “Immune-infiltrated” subtype. Cell-to-cell interaction plus extracellular matrix pathway hypomethylation and higher tumor purity were observed in the better-performing “Cellular” subtype. The findings were validated in additional DNA methylation and RNA sequencing datasets as well as with immunohistochemical staining. Plasma methylomes distinguished chordomas from other clinical differential diagnoses by applying fifty chordoma-versusother binomial generalized linear models in random 20% testing sets (mean AUROC = 0.84, 95%CI: 0.52–1.00). Tissue-based and plasma-based methylation signals were highly correlated in both prognostic clusters. Additionally, leave-one-out models accurately classified all tumors into their correct cluster based on plasma methylation data. Here, we show the first identification of prognostic epigenetic chordoma subtypes and first use of plasma methylome based biomarkers to noninvasively diagnose and subtype chordomas. These results may transform patient management by allowing treatment aggressiveness to be balanced with patient risk according to prognosis.

  Study EGAS00001006406

• 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

• Integrative analysis of exome-seq, RNA-seq, ATAC-seq (bulk and single-cell), and Hi-C data generated from 3-D spatially mapped samples acquired during surgical resection from 10 patients diagnosed with IDH-WT glioblastoma

  Glioblastoma (GBM), the most common primary brain cancer in adults, remains incurable with no targeted therapies approved despite decades of investigation into its molecular landscape. Treatment failure is attributed to intratumoral heterogeneity in the GBM genome and epigenome, which foster tumor evolution and selection of resistant clones. However, tumor evolution and intratumoral heterogeneity remain poorly understood on the level of the whole tumor as most studies are based on single samples per patient and lack spatial context. Here, we have used 3-D neuronavigation during surgical resection for 10 primary IDH-WT GBM patients to collect 102 samples representing maximal tumor diversity, each mapped by 3-D spatial coordinates. We have applied a strategic set of genomic and epigenomic assays spanning multiple levels of resolution to discover, orthogonally validate, and functionally assess drivers of tumor evolution and intratumoral heterogeneity. These include extrachromosomal DNA amplifications, chromothripsis events, inversions, and translocations that disrupt both the GBM genome and epigenome while generating fusion transcripts and opportunities for therapeutic intervention. We define epigenomic programs that contribute to GBM evolution and intratumoral heterogeneity, revealing their 3-D spatial patterning within whole tumors and their cell type(s) of origin in single-cell data from the same tumor samples. Notably, we distinguish neuronal, glial, and immune programs aberrantly active in tumor cells from their counterparts in normal cells and discover NEUROD1, JUN/FOS, and NF1 transcription factors as key drivers of GBM evolution and growth. Collectively, these data provide unprecedented insight into GBM evolution and intratumoral heterogeneity from single-cell to whole-tumor resolution, redefining current understanding and providing a rich resource of targets for therapeutic investigation.

  Study EGAS00001006785

• 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

• 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

• 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

• WGS of cell-free DNA derived from plasma of patients with pediatric sarcoma and healthy controls, and lcWGS/RRBS of matched tumor tissue

  This dataset contains: i) 241 deep (median 12x) whole-genome sequencing profiles of 95 patients with Ewing sarcoma, 31 patients with other pediatric sarcomas, and 22 additional profiles from healthy controls. Sequencing was performed on a NovaSeq 6000 instrument using the NovaSeq S4 2x100 bp configuration. In addition, pilot experiments for 18 cfDNA samples were performed using Illumina HiSeq 2000/2500 machines (2x75 bp configuration). Data is provided as .fastq.gz files (2 files, .R1 and .R2, per sample). i) Low coverage whole-genome-sequencing on 43 tumor biopsy samples from patients with Ewing sarcoma with matched cfDNA samples. The samples were sequenced using a NovaSeq 6000 instrument with the NovaSeq S4 1x100 bp configuration. Data are provided as unmapped (raw) .bam files. iii) Reduced-representation bisulfite sequencing data for 38 tumor biopsy samples from patients with Ewing sarcoma with matched cfDNA samples, and 2 control samples. RRBS libraries were sequenced on Illumina HiSeq 2000/2500 machines. Data are provided as unmapped (raw) .bam files.

  Dataset EGAD00001007080

• Next Generation Children project

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

  Dataset EGAD00001007868

• panALL exome sequencing, data set2, 700 samples

  panALL exome sequencing, data set2, 700 samples

  Dataset EGAD00001008268

• EPICC: Evolutionary Predictions in Colorectal Cancer WGS-3

  WGS profiling bam files from colorectal carcinoma and adenoma.

  Dataset EGAD00001008642

• A study of the genetic basis of evation by Acute Myeloid Leukaemia of Graft vs Leukaemia effects after allogeneic bone marrow transplantation

  Acute myeloid leukaemia (AML) is an aggressive and molecularly diverse disease with a poor overall survival of 20-25%. With an annual incidence of 2.9 per 100,000, AML is currently the commonest myeloid malignancy in Europe, yet the two main therapeutic options for this disease, anthracyclines and purine analogues, have remained unchanged for over 20 years. Currently patients are stratified at diagnosis according to a series of clinicopathological parameters (e.g. age, white cell count and presence/absence of previous clonal haematological disease) and molecular markers (e.g. chromosomal translocations/deletions, aneuploidy and mutations in genes such as FLT3 and NPM1). Patients with adverse prognostic features, whose prognosis is particularly poor (e.g. <15% long-term survival) are offered treatment with allogeneic bone marrow transplantation (allo-BMT) if a sibling or unrelated donor is available. This can significantly improve survival (e.g. up to 40% long-term survival in some contexts), albeit at the expense of significant toxicity and transplant-related mortality (TRM). Allo-BMT is thought to work in part by allowing the delivery of large doses of chemotherapy followed by haemopoietic "rescue" with donor haemopoietic stem cells (haemopoietic failure would otherwise ensue). However, potentially the most potent effect of allo-BMT is the cytotoxic effect of donor lymphocytes against AML blasts, a phenomenon known as graft-vs-leukaemia (GVL) effect. Increasingly, transplants using reduced chemotherapy intensity (mini-allografts) are being used that partially circumvent the toxicity from chemotherapy and rely on GVL to effect cure. Nevertheless, AML relapse after allo-BMT still occurs at a significant rate of up to 80% depending on the type of transplant. There is accumulating evidence that genetic events in residual leukaemic cells enable them to evade immunodetection and therefore survive the GVL effect and expand to cause relapse. The most striking example of this is the loss of HLA antigens after transplants in which donor and recipient are not fully HLA-matched. In these cases, the leukaemia "deletes" the genomic region containing the disparate HLA antigen which was preferentially targeted as "foreign" by the GVL effect. However, the genetic basis of immune evasion in the majority of transplants, which are fully HLA matched, is not known. One possibility is that loss of genes coding for antigens outside the HLA locus but which are also targets of GVL may operate, alternatively genetic events that affect processes downstream of immunological cytotoxicity may be responsible. The identification of genetic events that mediate immune evasion would not only facilitate the understanding of this process but can help plan therapeutic interventions that improve the outcomes of allogeneic transplantation for AML and other disorders. We intend to study this by conducting exome sequencing on 6 cases of AMLs from patients that attend my clinic at Addenbrooke's hospital and have relapsed after allogeneic transplantation. Samples from AML diagnosis, remission/normal and AML relapse (total n=18) will be studied to identify somatic mutations in the primary AML and those acquired by the relapsed clone. The 18 samples will also be studied by array CGH to detect regions of genomic amplification or deletion.

  Dataset EGAD00001000404

• 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

• Targeted sequencing of genes recurrently mutated in AML

  Background Massively parallel sequencing technology has transformed cancer genomics. It is now feasible, in a clinically relevant time-frame, for a clinically manageable cost, to screen DNA from patient tumours for mutations essentially genome-wide. The challenge for personalised medicine will be to increase the sample size to thousands or tens of thousands of well-characterised cases in order to attain sufficient statistical power to stratify patients accurately across the complexity and genomic heterogeneity expected for most of the common tumour types. Currently, whole genome sequencing on this scale is not feasible, and targeted sequencing of relevant portions of the genome will be required. Pilot data We have developed protocols for large-scale, multiplexed sequencing of 100-200 genes in thousands of samples. Essentially, using robotic technology, genomic DNA from the cancer specimen is processed into sequencing libraries with unique DNA barcodes, thereby allowing sequencing reads to be attributed to the sample they derive from. Currently, these sequencing libraries can be generated in a 96-well format using fully automated protocols, and we are exploring methods to expand this to a 384-well format. The sequencing libraries are pooled and hybridized to custom sets of RNA baits representing the genomic regions of interest. Sequencing of the pulled-down libraries is done in pools of 48-96 samples per lane of an Illumina Hi-Seq. This protocol is already implemented at the Sanger Institute. We have published proof that somatic mutations in novel cancer genes can be identified from exome-wide sequencing. In unpublished pilot data, we have established the feasibility of robotic library production, custom pull-down, and multiplexed sequencing of barcoded libraries for 100 known myeloid cancer genes across 760 myelodysplasia samples. Highlights of the data thus far analysed reveal that the coverage is remarkably even between samples; when 96 samples are run, average coverage per lane of sequencing is ~250, with 90-95% of targeted exons covered by >25 reads; known mutations can be discovered in the data set; and the protocol is amenable to whole genome amplified DNA. The bioinformatic algorithms for identification of substitutions and indels in pull-down data are well-established; we have pilot data proving that copy number changes, LOH and genomic rearrangements in specific regions of interest can also be identified by tiling of baits across the relevant loci. Proposal We propose to apply this methodology to 10000 samples from patients with AML enrolled in clinical trials over the last 10-20 years. Oncogenic point mutations and potentially genomic rearrangements will be identified, and linked to clinical outcome data, with a view to undertaking the following sorts of analyses: ? Identification of co-occurrence, mutual exclusivity and clusters of driver mutations. ? Correlation of prognosis with driver mutations and potentially gene-gene interactions ? Exploration of genomic markers of drug response Ultimately, we would like to be in a position to release the mutation data together with matched clinical outcome data to genuine medical researchers via a controlled access approach, possibly within the COSMIC framework (www.sanger.ac.uk/genetics/CGP/cosmic/). The vision here is to generate a portal whereby a clinician faced with an AML patient and his / her mutational profile can obtain a ?personalised? prediction of outcome, together with a fair assessment of the uncertainty of the estimate. With a sufficient sample size, there would also be the potential to develop decision support algorithms for therapeutic choices based on such data.

  Dataset EGAD00001000606

• 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

• 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

• Limited Association between HRR Gene Alterations and HRD in Molecular Tumor Board Cancer Samples: Who should be tested for HRD?

  Alterations in Homologous Recombination Repair (HRR) Pathway genes have been found to be associated with HR-Deficiency (HRD), which is an approved biomarker for PARP Inhibitor (PARPi) treatment. The aim of a Molecular Tumor Board (MTB) is to identify molecular alterations in cancer patients with advanced tumors that may suggest off-label treatment options. So far, few studies have analyzed the presence of HRR gene mutations and their association with HRD outside of clinical studies. Currently, no data on HRD testing in the setting of a MTB have been published. For the present study, a cohort of 237 patients encompassing 24 different tumor entities was collected from the MTB of the Comprehensive Cancer Center Erlangen-EMN. We show that an elevated Genomic Instability Score (GIS ≥ 42) can occur in samples with and without mutations in HRR-related genes. Overall, 38.1% of cancer samples with BRCA1/2 mutations, 10.9% of tumors with alterations in HRR genes other than BRCA1/2, and 4.3% of cancer samples without HRR gene mutations harbored an elevated GIS. Notably, our data show that various inactivating BRCA1/2 mutations are not associated with an elevated GIS. Taken together, panCancer assessment of HRD in addition to BRCA1/2 and other HRR gene mutational analysis is recommended to guide decisions regarding PARPi treatment. Further studies are needed to establish thresholds for GIS in non-ovarian cancer entities. Finally, HRD can be observed in 4.3% of BRCA1/2 and other HRR gene wildtype cancer samples, and may emerge as an independent biomarker for PARPi in the future.

  Study EGAS00001008063

• 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

• Molecular subtypes of pancreatic neuroendocrine tumors mutated in MEN1/DAXX/ ATRX explain biological variability

  Pancreatic Neuroendocrine Tumors (PanNETs) most probably arise from the islets of Langerhans, although direct evidence is still lacking. About 40% of non-functioning (NF) PanNETs harbour mutations in MEN1, often in conjunction with alterations in DAXX/ATRX. Despite the homogeneous genetic background, the ADM (mutation in DAXX/ATRX and MEN1) group is heterogeneous regarding molecular features and response to treatment. To explore the molecular and clinical heterogeneity within this group, we integrated transcriptomic and DNA methylation data from 36 and 93 samples, respectively. Using a multi-omic approach we described a two-level hierarchy separating PanNET subtypes mutated in MEN1, DAXX and ATRX. First, DNA methylation discriminates ADM-PanNET from PanNET mutated in MEN1 only (α-like). Compared to α-like, all ADM tumors exhibit DNA methylation changes at enhancer and lowly methylated regions as well as at peri-centromeric and telomeric regions, associated with alternative lengthening of telomeres, increased chromosomal instability, and enhanced proliferative capacity. At a second level, transcriptomic analysis revealed three distinct ADM subtypes: ADM hypoxic, ADM NST (No Special Type), and ADM immunosuppressive. The ADM hypoxic subtype is characterized by strong hypoxia signature, likely regulated epigenetically. The ADM NST subtype appears to be primarily driven by epigenetic changes that promote proliferation. The ADM immunosuppressive is highly enriched in immune component and strong metastasis-like signature. While all the three subtypes are highly associated with liver metastasis formation and invasion, the ADM immunosuppressive are significantly smaller (<2.5 cm p.value=0.023) despite displaying aggressive molecular features. By defining these three novel ADM subtypes, our study provides a refined framework for understanding PanNET heterogeneity. This classification underscores potential diagnostic markers and highlights distinct biological vulnerabilities that may inform the development of subtype-tailored therapeutic strategies.

  Study EGAS00001008272

• 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

• 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

• Genetic heritage of the Baphuthi highlights an over-ethnicised notion of 'Bushman' in the Maloti-Drakensberg, southern Africa

  Using contemporary people as proxies for ancient communities is a contentious but necessary practice in anthropology. In southern Africa, the distinction between the Cape KhoeSan and eastern KhoeSan remains unclear as ethnicity labels have been changed through time and most communities were decimated if not extirpated. The eastern KhoeSan may have had genetic distinctions from neighbouring communities who speak Bantu languages and KhoeSan further away, alternatively the identity may not have been tied to any notion of biology, instead denoting communities with a nomadic 'lifeway' distinct from African agro-pastoralism. The Baphuthi of the 1800s in the Maloti-Drakensberg, southern Africa had a substantial KhoeSan constituency and a lifeway of nomadism, cattle raiding, and horticulture. Baphuthi heritage could provide insights into the history of the eastern KhoeSan. We examine genetic affinities of 23 Baphuthi to discern whether the narrative of KhoeSan descent reflects distinct genetic ancestry. Genome-wide SNP data (Illumina GSA) were merged with 52 global populations, for 160K SNPs. Genetic analyses show no support for a unique eastern KhoeSan ancestry distinct from other KhoeSan or southern Bantu-speakers. The Baphuthi have strong affinities with early-arriving southern Bantu-speaking (Nguni) communities, as the later-arriving non-Nguni show strong evidence of recent African admixture possibly related to Late-Iron Age migrations. The references to communities as 'San',and 'Bushman' in historic literature has often been misconstrued as notions of ethnic/biological distinctions. The terms may have reflected ambiguous references to non-sedentary polities instead, as seems to be the case for the eastern ‘Bushman’ heritage of the Baphuthi.

  Study EGAS00001007080

• 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

• A transcriptomic approach to understand patient susceptibility to pneumonia after abdominal surgery

  Rationale Nosocomial infections are a major healthcare challenge, developing in over 20% of patients aged 45 or over undergoing major-abdominal surgery, with postoperative pneumonia associated with an almost five-fold increase in 30-day mortality. Objectives To describe immune-pathways and gene-networks altered following major-abdominal surgery and identify transcriptomic patterns associated with postoperative pneumonia. Methods and Measurements From a prospective consecutive cohort (n=150) undergoing major-abdominal surgery whole-blood RNA was collected preoperatively and at three time-points postoperatively (2-6, 24 and 48hrs). Twelve patients diagnosed with postoperative pneumonia and 27 matched patients remaining infection-free were identified for analysis with RNA-sequencing. Main Results Compared to preoperative sampling, 3,639 genes were upregulated and 5,043 downregulated at 2-6hrs. Pathway-analysis demonstrated innate-immune activation with neutrophil-degranulation and Toll-like-receptor signalling upregulation alongside adaptive-immune suppression. Cell-type deconvolution of preoperative RNA-sequencing revealed elevated S100A8/9-high neutrophils alongside reduced naïve CD4 T-cells in those later developing pneumonia. Preoperatively, a gene-signature characteristic of neutrophil-degranulation was associated with postoperative pneumonia acquisition (P=0.00092). A previously reported Sepsis Response Signature (SRSq) score, reflecting neutrophil-dysfunction and a more dysregulated host response, at 48hrs postoperatively, differed between patients subsequently developing pneumonia and those remaining infection-free (P=0.045). Analysis of the novel neutrophil gene-signature and SRSq scores in independent major-abdominal surgery and polytrauma cohorts indicated good predictive performance in identifying patients suffering later infection. Conclusions Major-abdominal surgery acutely upregulates innate-immune pathways while simultaneously suppressing adaptive-immune pathways. This is more prominent in patients developing postoperative pneumonia. Preoperative transcriptomic signatures characteristic of neutrophil-degranulation and postoperative SRSq scores may be useful predictors of subsequent pneumonia risk.

  Study EGAS00001007229

• 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

• Multi-omic analysis of the tumor microenvironment shows clinical correlations in Ph1 study of atezolizumab +/- SoC in MM

  Multiple myeloma (MM) remains incurable, and treatment of relapsed/refractory (R/R) disease is challenging. There is an unmet need for more targeted therapies in this setting; deep cellular and molecular phenotyping of the tumor and microenvironment in MM could help guide such therapies. This phase 1b study (NCT02431208) evaluated the safety and efficacy of the anti-programmed death-ligand 1 monoclonal antibody atezolizumab (Atezo) alone or in combination with the standard of care (SoC) treatments lenalidomide (Len) or pomalidomide (Pom) and/or daratumumab (Dara) in patients with R/R MM. Study endpoints included incidence of adverse events (AEs) and overall response rate (ORR). A novel unsupervised integrative multi-omic analysis was performed using RNA sequencing, mass cytometry immunophenotyping, and proteomic profiling of baseline and on-treatment bone marrow samples from patients receiving Atezo monotherapy or Atezo+Dara. A similarity network fusion (SNF) algorithm was applied to preprocessed data. Eighty-five patients were enrolled. Treatment-emergent deaths occurred in 2 patients; both deaths were considered unrelated to study treatment. ORRs ranged from 11.1% (Atezo+Len cohorts, n=18) to 83.3% (Atezo+Dara+Pom cohort, n=6). High-dimensional multi-omic profiling of the tumor microenvironment and integrative SNF analysis revealed novel correlations between cellular and molecular features of the tumor and immune microenvironment, patient selection criteria, and clinical outcome. Atezo monotherapy and SoC combinations were safe in this patient population and demonstrated some evidence of clinical efficacy. Integrative analysis of high dimensional genomics and immune data identified novel clinical correlations that may inform patient selection criteria and outcome assessment in future immunotherapy studies for myeloma.

  Study EGAS00001007286

• 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

• 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

• Ultra-fast deep-learned pediatric CNS tumor classification.

  The primary treatment of CNS tumors starts with a neurosurgical resection in order to obtain tumor tissue for diagnosis and to reduce tumor load and mass effect. The neurosurgeon has to decide between radical resection versus a more conservative strategy to prevent surgical morbidity. The prognostic impact of a radical resection varies between tumor types. However due to a lack of pre-operative tissue-based diagnostics, limited knowledge of the precise tumor type is available at the time of surgery. Current standard practice includes preoperative imaging and intraoperative histological analysis, but these are not always conclusive. After surgery, histopathological and molecular tests are performed to diagnose the precise tumor type. The results may indicate that an additional surgery is needed or that the initial surgery could have been less radical. Using rapid Nanopore sequencing, a sparse methylation profile can be directly obtained during surgery, making it ideally suited to enable intraoperative diagnostics. We developed a state-of-the-art neural-network approach called Sturgeon, to deliver trained models that are lightweight and universally applicable across patients and sequencing depths. We demonstrate our method to be accurate and fast enough to provide a correct diagnosis with as little as 20 to 40 minutes of sequencing data in 45 out of 49 pediatric samples, and inconclusive results in the other four. In four intraoperative cases we achieved a turnaround time of 60-90 minutes from sample biopsy to result; well in time to impact surgical decision making. We conclude that machine-learned diagnosis based on intraoperative sequencing can assist neurosurgical decision making, allowing neurological comorbidity to be avoided or preventing additional surgeries.

  Study EGAS00001007475

• 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