22928 results for "ega"

in 78.12 milliseconds.

• A novel breast cancer cell line derived from a metastatic bone lesion of a breast cancer patient

  Purpose: We sought to generate and characterize a novel cell line from a breast cancer bone metastasis in order to better study the progression of the disease.Methods: The cell line, P7731, was derived from a metastatic bone lesion of a breast cancer patient and assessed for marker expression. P7731 was analyzed for DNA copy number variation, somatic mutations and gene expression and compared to the primary tumor. Results: P7731 cells are negative for estrogen receptor alpha (ERα), progesterone receptor (PR) and HER2 (triple-negative), strongly express vimentin (100% of cells positive) and also express cytokeratins 8/18 and 19 but at lower frequencies. Flow cytometry indicates P7731 cells are predominantly CD44+/CD49f+/EpCAM-, consistent with a primitive, mesenchymal-like phenotype. The cell line is tumorigenic in immunocompromised mice. Exome sequencing identified a total of 45 and 76 somatic mutations in the primary tumor and cell line respectively, of which 32 were identified in both samples and included mutations in known driver genes PIK3CA, TP53 and ARID1A. P7731 retains the DNA copy number alterations present in the matching primary tumor. Homozygous deletions detected in the cell line and in the primary tumor were found in regions containing three known (CDKN2A, CDKN2B and CDKN1B) and 23 putative tumor suppressor genes. Cell line specific gene amplification coupled with mRNA expression analysis revealed genes and pathways with potential pro-metastatic functions.Conclusion: This novel human breast cancer-bone metastasis cell line will be a useful model to study aspects of breast cancer biology, particularly metastasis-related changes from breast to bone.

  Study EGAS00001002840

• Colorectal adenomas, carcinomas and adjacent normal NKI-AvL TGO series NGS-ProToCol

  Cancer is caused by changes in the genome of cells, resulting in aberrant transcripts and proteins. In order to understand why normal cells become cancer cells and how we can identify and destroy them, the CTMM-NGS-ProToCol project aimed to understand prostate and colorectal cancer through DNA and RNA molecular profiling, to further improve diagnosis, prognosis and treatment.The state-of-the-art within the field of genome sequencing makes it possible to identify all the modifications in our genome that encode for RNA and of the transcripts themselves. The latest development in RNA sequencing is that if one does not select for polyA RNA, but first removes rRNA from the total RNA and then performs ‘random’ RNA sequencing, 30-60% of the reads are intronic. Since about 40% of the genome consists of genes (and only 1.6% is exonic), deep RNA sequencing will cover not only the exons, but also a large portion of the introns and therefore much of the genome. The power of the combined sequencing of genome and transcriptome is the full view of copy number aberrations (CNA), exon mutations and the consequence and relevance of the mutations at the RNA level in the form of differential gene expression, alternative splicing, alternative promoter usage, novel transcripts, fusion transcripts, read-through transcripts and mutations. For the colorectal (tumor) tissue dataset, we applied molecular profiling to 30 colorectal adenomas, 30 colorectal carcinomas and 18 normal adjacent colon tissues.

  Study EGAS00001002854

• Clonal fitness inferred from timeseries modeling of single cell cancer genomes

  Progress in defining genomic fitness landscapes in cancer by copy number alterations (CNA) has been impeded by lack of single cell and timeseries sampling. We generated 42,000 single cell whole genomes (scWGS) from breast epithelium and primary triple negative breast cancer (TNBC) patient-derived xenografts (PDX) collected during multi-year time series. Using a Wright-Fisher population genetics model, we inferred reproducible CNA-defined clonal fitness dynamics induced by TP53 mutation and cisplatin chemotherapy, with accurate forecasting of experimentally enforced clonal competition dynamics. Drug treatment in three long-term serially passaged TNBC PDX resulted in cisplatin resistant clones that had shown low fitness in the untreated setting. By contrast high fitness clones from treatment naive controls were eradicated reflecting an inversion of the fitness landscape. Upon drug release selective pressure dynamics were reversed indicating a fitness cost of treatment resistance. Taken together, our findings reveal clonal fitness dynamics linked to CNA and therapeutic resistance in polyclonal tumours.

  Study EGAS00001004448

• Versatile workflow for cell-type resolved transcriptional and epigenetic profiles from cryopreserved human lung

  The complexity of the lung microenvironment and changes in cellular composition during disease progression make it exceptionally hard to understand molecular mechanisms leading to the development of chronic lung diseases. Although recent advances in cell-type resolved and single-cell sequencing approaches hold great promise for studying complex diseases, their implementation relies on local access to fresh tissue, as traditional tissue storage methods do not allow viable cell isolation. To overcome these hurdles, we developed a novel, versatile workflow that allows long-term storage of human lung tissue with high cell viability, permits thorough sample quality check before cell isolation, and is compatible with sequencing-based profiling, including single-cell approaches. We demonstrate that cryopreservation is suitable for the isolation of multiple cell types from human lung and is applicable to both healthy and diseased tissue, including COPD and tumour samples. Basal cells isolated from cryopreserved airways retain the ability to differentiate, indicating that cellular identity is not altered by cryopreservation. Importantly, using RNA sequencing and EPIC Array, we show that genome-wide gene expression and DNA methylation signatures are preserved upon cryopreservation, emphasizing the suitability of our workflow for -omics profiling of human lung cells. In addition, we obtained high-quality single-cell RNA sequencing data of cells isolated from the cryopreserved human lung, demonstrating that cryopreservation empowers single-cell approaches. Overall, thanks to its simplicity, our cryopreservation workflow is well-suited for prospective tissue collection by academic collaborators and biobanks, opening worldwide access to human tissue.

  Study EGAS00001004477

• Identification of Therapeutic Targets in Rhabdomyosarcoma Through Integrated Genomic, Epigenomic, and Proteomic Analyses

  Personalized cancer therapy based on the somatic mutations identified in patient tumors is becoming an increasingly emphasized approach to improve outcomes of patients with cancer. There are also examples of therapeutic vulnerabilities in cancer that result from changes in gene expression that are a direct or indirect result of tumor specific epigenetic perturbations. These genomic and epigenomic changes are ultimately manifest in the tumor proteome and phosphoproteome. In this study, we integrated genomic, epigenomic and proteomic data for rhabdomyosarcoma (RMS) to identify therapeutic vulnerabilities. RMS was selected for this analysis because RAS pathway mutations in rhabdomyosarcoma (RMS) are the most common potentially actionable lesions across pediatric solid tumors. The epigenomic data was useful for identifying deregulated developmental pathways in RMS including the WNT, HH, BMP, adenyl cyclase, p38/MAPK and PI3K pathways. Perturbations in those 6 myogenic signal transduction pathways were also evident in the proteome and phosphoteome data. In addition, the proteomic/phosphoproteomic data revealed that the cell cycle checkpoint, unfolded protein response and RB/E2F pathways were deregulated in RMS relative to normal muscle. Recent success targeting CDK4/6 and MEK in adult cancers with RAS mutations led us to test the value of these targets in RMS tumors. We also targeted the unfolded protein response pathway and cell cycle checkpoint pathway using molecular targeted therapeutics in orthotopic patient derived xenografts. Taken together, these data demonstrate the value of integrating epigenomic and proteomic data to identify tumor vulnerabilities that extend beyond somatic mutations identified in the genome.

  Study EGAS00001002967

• Detection and characterization of lung cancer using cell-free DNA fragmentomes

  Lung cancer remains the leading cause of cancer death world-wide, largely due to its late diagnosis. Non-invasive approaches for assessment of cell-free DNA (cfDNA) provide an opportunity for detection and intervention that may have broader accessibility than current imaging approaches. Using a machine learning model for detecting tumor-derived cfDNA through genome-wide analyses of cfDNA fragmentation, we examined a prospective study of 365 individuals at risk for lung cancer (Lung Cancer Diagnostic Study, LUCAS), including 129 individuals ultimately diagnosed with lung cancer and 236 individuals determined to not have lung cancer. We externally validated the cancer detection model using an independent cohort of 385 non-cancer individuals and 46 predominantly early stage lung cancer patients. Combining fragmentation features with clinical risk factors and CEA levels followed by CT imaging detected 94% of patients with cancer across stages and subtypes, including 91% of stage I/II and 96% of stage III/IV, at 80% specificity. Genome-wide fragmentation profiles across ~13,000 ASCL1 transcription factor binding sites could be used to distinguish individuals with small cell lung cancer (SCLC) from those with non-small cell lung cancer (NSCLC) with high accuracy (AUC=0.98). Among individuals with lung cancer, a higher cfDNA fragmentation score was associated with tumor size and invasion, and represented an independent prognostic indicator of survival. These studies provide a facile approach for non-invasive detection of lung cancer and clinical management of this disease.

  Study EGAS00001005340

• Open_Targets_020_Epigenomes_of_Cell_Lines

  There is currently a drive to establish cell based assay systems of greater human biological and disease relevance through the use of well characterised transformed cell lines, primary cells and complex cellular models (e.g. co-culture, 3D models). However, although the field is gaining valuable experience in running more non-standard & complex cell assays for target validation and compound pharmacology studies, there is the lack of a systematic approach to determine if this expansion in cell assay models is reflected in increased human biological and disease relevance. The increasing wealth of publically available transcriptomic, and epigenome (ENCODE and Epigenome Roadmap) data represents an ideal reference mechanism for determining the relationship between cell types used for target & compound studies to primary human cells and tissues from both healthy volunteers & patients. The CTTV020 epigenomes of cell line project aims to generate epigenetic and transcriptomic profiles of cell lines and compare these with existing and newly generated reference data sets from human tissue and cell types. The aim is to identify assay systems which will provide greater confidence in translating target biology and compound pharmacology to patients. Multiple cell types commonly used within research have been grouped according to biology. Examples include erythroid, lung epithelial, hepatocyte cell types and immortalised models of monocyte / macrophage biology. 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 EGAS00001003137

• Open_Targets_020_Epigenomes_of_Cell_Lines

  There is currently a drive to establish cell based assay systems of greater human biological and disease relevance through the use of well characterised transformed cell lines, primary cells and complex cellular models (e.g. co-culture, 3D models). However, although the field is gaining valuable experience in running more non-standard & complex cell assays for target validation and compound pharmacology studies, there is the lack of a systematic approach to determine if this expansion in cell assay models is reflected in increased human biological and disease relevance. The increasing wealth of publically available transcriptomic, and epigenome (ENCODE and Epigenome Roadmap) data represents an ideal reference mechanism for determining the relationship between cell types used for target & compound studies to primary human cells and tissues from both healthy volunteers & patients. The CTTV020 epigenomes of cell line project aims to generate epigenetic and transcriptomic profiles of cell lines and compare these with existing and newly generated reference data sets from human tissue and cell types. The aim is to identify assay systems which will provide greater confidence in translating target biology and compound pharmacology to patients. Multiple cell types commonly used within research have been grouped according to biology. Examples include erythroid, lung epithelial, hepatocyte cell types and immortalised models of monocyte / macrophage biology. 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 EGAS00001003139

• Healthy_ageing_thymus

  Single cell RNA sequencing of immune and non-immune cells from healthy ageing thymus Thymus is a primary lymphoid organ that creates an environment for the T lymphocyte precursors differentiation into the naive T cells, which involves TCR gene rearrangement, negative and positive selection. Despite thymus critical function for the recognition of the pathogens it starts to atrophy very early in life. Thymic Epithelial Space (TES), where thymocyte education occurs, starts to decline rapidly just after birth, while with the onset of puberty thymus is also progressively replaced with adipose tissue. This leads to a reduced naive T cell output,which decreases an organism's ability to recognize pathogens. In our lab, we have already collected and characterized cell type composition and gene regulatory networks of the fetal tand adult human thymus (Park et al, Science 2020). In the current project we plan to expand the number of paediatric and adult samples in the study with the aim to understand how the process of thymic involution happens. We will use single-cell transcriptomics to zoom on changes that occur in immune and non-immune cell types in different phases of thymic involution with the hope to understand which gene expression changes are causal for the thymic involution process and how naïve T cell production can be increased in elderly. 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 EGAS00001004311

• Novel optineurin frameshift insertion causing familial frontotemporal dementia and parkinsonism without amyotrophic lateral sclerosis

  Frontotemporal Dementia (FTD) is a common cause of Young Onset Dementia and has diverse clinical manifestations involving behaviour, executive function, language and motor function, including parkinsonism. Up to 50% of FTD patients report a positive family history, supporting a strong genetic basis, particularly in cases with both FTD and amyotrophic lateral sclerosis (FTD-ALS). Mutations in three genes are associated with the majority of familial FTD (fFTD) cases - microtubule associated protein tau gene (MAPT), progranulin gene (PGRN), and hexanucleotide repeat expansions in chromosome 9 open reading frame 72 (C9orf72), whilst mutations in other genes such as optineurin (OPTN) have rarely been reported. Mutations in OPTN have been reported mostly in familial and sporadic cases of ALS, or in rare cases of FTD-ALS, but not in association with pure or predominant FTD and/or parkinsonian phenotype. Here, we report for the first time, a family from the Philippines with 4 affected members harbouring the same OPTN frameshift insertion, presenting with FTD-related phenotypes, including one sibling with predominant parkinsonism resembling corticobasal syndrome. Notably, none of the affected members showed any evidence of motor neuron disease or ALS at the time of writing, both clinically and on electrophysiological testing, expanding the phenotypic spectrum of OPTN mutations. Close follow-up of mutation carriers for the development of new clinical features and wider investigation of additional family members with further genetic analyses will be conducted to investigate the possibility of other genetic modifiers in this family which could explain phenotypic heterogeneity.

  Study EGAS00001005220

• Exome-wide mutation analysis of cell-free DNA to simultaneously monitor the full spectrum of cancer treatment outcomes

  Purpose: Cell-free DNA (cfDNA) offers a non-invasive approach to monitor cancer. Here we develop a method using whole-exome sequencing (WES) of cfDNA for simultaneously monitoring the full spectrum of cancer treatment outcomes, including MRD, recurrence, evolution, and second primary cancers. Experimental Design: Three simulation datasets were generated from 26 cancer patients to benchmark the detection performance of MRD/recurrence and second primary cancers. For further validation, cfDNA samples (n=76) from cancer patients (n=35) with six different cancer types were used for performance validation during various treatments.Results: We present a cfDNA-based cancer monitoring method, named cfTrack. Taking advantage of the broad genome coverage of WES data, cfTrack can sensitively detect MRD and cancer recurrence by integrating signals across known clonal tumor mutations of a patient. In addition, cfTrack detects tumor evolution and second primary cancers by de novo identifying emerging tumor mutations. A series of machine learning and statistical denoising techniques are applied to enhance the detection power. On the simulation data, cfTrack achieved an average AUC of 99% on the validation dataset and 100% on the independent dataset in detecting recurrence in samples with tumor fraction ≥0.05%. In addition, cfTrack yielded an average AUC of 88% in detecting second primary cancers in samples with tumor fraction ≥0.2%. On real data, cfTrack accurately monitors tumor evolution during treatment, which cannot be accomplished by previous methods.Conclusion: Our results demonstrated that cfTrack can sensitively and specifically monitor the full spectrum of cancer treatment outcomes using exome-wide mutation analysis of cfDNA.

  Study EGAS00001005906

• A higher ctDNA fraction decreases survival in regorafenib-treated metastatic colorectal cancer patients

  The predictive effect of circulating tumor DNA (ctDNA) in colorectal cancer (CRC) treatment is still highly discussed. The primary objective of this study was to investigate a possible prognostic/predictive value of ctDNA under regorafenib treatment. This prospective multicenter translational biomarker phase II pilot study enrolled 30 metastatic CRC patients (67% men, 33% women) treated with regorafenib. ctDNA was assessed in plasma before treatment start and at defined time points during administration. Measurement of tumor fraction as well as mutation and copy number analysis of CRC driver genes were performed by next generation sequencing approaches. Multivariate analyses for survival and treatment efficacy were adjusted to age, gender and ECOG. Disease control rate was 30%. Median tumor fraction at baseline was 18.5% (0-49.9). Mutations in CRC driver genes or genes involved in angiogenesis were identified in 25 patients (83.3%). KRAS mutations were detected in 13 out of 14 KRAS positive tumors, in three patients without KRAS mutation in the respective tumor acquired mutations as a consequence of prior anti-EGFR treatment were detected. In a subset of patients, novel occurring mutations or focal amplifications were detected. A tumor fraction of 5% and higher at baseline was significantly associated with a decreased OS (p=0.022; HHR 3.110 (95%CI: 1.2-8.2). ctDNA is detectable in a high proportion of mCRC patients. Higher ctDNA levels are associated with survival among regorafenib treatment. Moreover, our data highlight the benefit of a combined evaluation of mutations and somatic copy number alterations (SCNA) in advanced cancer patients.

  Study EGAS00001004491

• Butyrate producers as potential next-generation probiotics: safety assessment of the administration of Butyricicoccus pullicaecorum to healthy volunteers

  Advances in gut microbiota research have triggered interest in developing colon butyrate producers as niche-specific next-generation probiotics, targeted at increasing colon butyrate production and countering disease-associated microbiota alterations. Crucial steps in the development of next-generation probiotics are the design of formulations with a reasonable shelf life as well as the safety demonstration of an intervention in healthy volunteers. One such potential next-generation butyrate-producing probiotic is Butyricicoccus pullicaecorum 25-3T, with demonstrated safety in in vitro as well as in animal models. Here, we examined the strain's safety, tolerability, and impact on microbiota composition and metabolic activity in healthy volunteers in a randomized, double blind, placebo-controlled, crossover study in 30 healthy volunteers. Study design consisted of two four-week intervention periods (108 CFU B. pullicaecorum [treatment] or maltodextrin [placebo] per day) with a three-week washout in between. We assessed adverse events, blood parameters (primary endpoints), and fecal microbiota composition and metabolite profiles (secondary endpoints).The number of reported adverse events during the B. pullicaecorum treatment was similar to placebo intervention, as were observed changes in blood chemistry parameters, bowel habits, and fecal calprotectin concentrations. Administration of the strain did not induce any disruptive effect in microbiota composition or metabolic activity. In this first human intervention trial with a butyrate-producing Clostridium cluster IV isolate, we demonstrated B. pullicaecorum 25-3T administration to be both safe and well tolerated by healthy participants. This safety study paves the way for the further development of the strain as a next-generation probiotic.

  Study EGAS00001003276

• Transcriptomic intra-tumor heterogeneity of colorectal cancer varies depending on tumor location within the colorectum

  Background; Intra-tumor heterogeneity (ITH) of colorectal cancer (CRC) complicates molecular tumor classification, such as transcriptional subtyping. Differences in cellular states, biopsy cell composition, and tumor microenvironment may all lead to ITH. Here we analyze ITH at the transcriptomic and proteomic levels to ascertain whether subtype discordance between multiregional biopsies reflects relevant biological ITH or lack of classifier robustness. Further, we study the impact of tumor location on ITH.MethodsMultiregional biopsies from stage II and III CRC tumors were analyzed by RNA sequencing (41 biopsies, 14 tumors) and multiplex immune protein analysis (89 biopsies, 29 tumors). CRC subtyping were performed using consensus molecular subtypes (CMS), CRC intrinsic subtypes (CRIS), and TUMOR types. ITH-scores and network maps were defined to determine the origin of heterogeneity. A validation cohort was used with one biopsy per tumor (162 tumors).ResultsOverall, inter-tumor transcriptional variation exceeded ITH, and subtyping calls were frequently concordant between multiregional biopsies. Still, some tumors had high ITH and were classified discordantly. Proximal tumors showed ITH related to differences in the microenvironment, while ITH within distal tumors were cancer-cell related. This subtype discordancy reflected actual molecular ITH within the tumors. The relevance of the subtypes was reflected at protein level where several inflammation markers were significantly increased in immune related transcriptional subtypes, which was verified in an independent cohort (Wilcoxon rank sum test; p<0.05).ConclusionOur transcriptomic and proteomic analyses show that the tumor location along the colorectum influence the ITH of CRC, which again influence the concordance of subtyping.

  Study EGAS00001004668

• Exome sequencing

  Human telomere biology disorders (TBD), a heterogeneous group of disorders characterized by telomere attrition and premature aging phenotypes, are caused by inherited loss-of-function mutations in telomere-associated genes. Here, we identify three germline heterozygous missense variants in RPA1 gene in four unrelated probands presenting with short telomeres and varying clinical features of TBD including bone marrow failure, myelodysplastic syndrome, T- and B-cell lymphopenia, pulmonary fibrosis, and skin manifestations. All variants cluster to DNA binding domain A of RPA1 protein. RPA1 is an essential single-strand DNA-binding protein essential for DNA replication and repair. RPA1 has been implicated in telomere maintenance but the mechanism remains elusive. Thus far, no genetic disorder had been linked to RPA1. Using FRET assays, we showed that RPA1E240K and RPA1V227A exhibit increased binding to single-strand and telomeric DNA, implying a gain in DNA-binding function while RPA1T270A has binding properties similar to wild type. To study the mutational effect, RPA1E240K, with the highest binding capacity, was knocked-in iPSCs using CRISPR/Cas9, which showed severe telomere shortening and impaired hematopoietic differentiation. Furthermore, in the RPA1E240K patient, stable blood counts over two decades due coincided with clonal rescue hematopoiesis with an RNA-degrading truncating RPA1 mutation and a uniparental isodisomy 17p with loss of RPA1E240K allele. Using single cell DNA sequencing, the two somatic genetic rescue events were proven to be independently acquired in hematopoietic stem cells. In summary, we describe the first human disease caused by germline RPA1 variants in individuals with TBD.

  Study EGAS00001005761

• Mutation-specific non-canonical pathway of PTEN as a distinct therapeutic target for glioblastoma

  PTEN is one of most frequently altered tumor suppressor genes in malignant tumors. The dominant negative effect of PTEN alteration suggests that aberrant function of PTEN mutation might be more disastrous than deletion, the most frequent genomic event in glioblastoma (GBM). This study aimed to understand the functional properties of various PTEN missense mutations and to investigate their clinical relevance. The genomic landscape of PTEN alteration was analyzed using the Samsung Medical Center GBM cohort and validated via The Cancer Genome Atlas dataset. Several hot-spot mutations were identified, and their subcellular distributions and phenotypes were evaluated. We established a library of cancer cell lines that over-express these mutant proteins using the U87MG and patient-derived cell models lacking functional PTEN. PTEN mutations were categorized into two major subsets: missense mutations in the phosphatase domain and truncal mutations in the C2 domain. We determined the subcellular compartmentalization of four mutant proteins (H93Y, C124S, R130Q, and R173C) from the former group and found that they had distinct localizations; those associated with invasive phenotypes (‘edge mutations’) localized to the cell periphery, while the R173C mutant localized to the nucleus. Invasive phenotypes derived from edge substitutions were unaffected by an anti-PI3K/Akt agent but were disrupted by microtubule inhibitors. PTEN mutations exhibit distinct functional properties with regard to their subcellular localization. Further, some missense mutations (‘edge mutations’) in the phosphatase domain caused enhanced invasiveness associated with dysfunctional cytoskeletal assembly, thus suggesting it to be a potent therapeutic target.

  Study EGAS00001004753

• RaScALL: Rapid screening of RNA-seq in acute lymphoblastic leukaemia

  RNA-sequencing (RNA-seq) efforts in acute lymphoblastic leukaemia (ALL) have identified numerous prognostically significant genomic alterations which can guide diagnostic risk stratification and treatment choices when detected early. However, a full RNA-seq Bioinformatics workflow is time-consuming and costly in a clinical setting where rapid detection and accurate reporting of clinically relevant alterations are essential. To accelerate the identification of ALL-specific alterations (including gene fusions, single nucleotide variants and focal gene deletions), we developed the rapid screening tool RaScALL, capable of identifying more than 100 prognostically significant lesions directly from raw sequencing reads. RaScALL uses the k-mer based targeted detection tool km and known ALL variant information to achieve a high degree of accuracy for reporting subtype defining genomic alterations compared to standard alignment-based pipelines. Gene fusions, including difficult to detect fusions involving EPOR and DUX4, were accurately identified in 98% (164 samples) of reported cases in a 180-patient Australian study cohort and 95% (n=63) of samples in a North American validation cohort. Pathogenic sequence variants were correctly identified in 75% of tested samples, including all cases involving subtype defining variants PAX5 p.P80R (n=12) and IKZF1 p.N159Y (n=4). Accurate detection of intragenic IKZF1 deletions resulting in aberrant transcript isoforms was also detectable with 98% accuracy. Importantly, the median analysis time for detection of all targeted alterations averaged 22 minutes per sample, significantly shorter than standard alignment-based approaches, ensuring accerelated risk-stratification and therapeutic triage.

  Study EGAS00001006460

• Integrative genomic analysis identifies key pathogenic mechanisms in primary mediastinal large B-cell lymphoma

  Primary mediastinal large B-cell lymphoma (PMBL) represents a clinically and pathologically distinct subtype of large B-cell lymphomas. Furthermore, molecular studies, including global gene expression profiling, have provided evidence that PMBL is more closely related to classical Hodgkin lymphoma (cHL). Although targeted sequencing studies have revealed a number of mutations involved in PMBL pathogenesis, a comprehensive description of disease-associated genetic alterations and perturbed pathways is still lacking. Here, we performed whole-exome sequencing of 95 PMBL tumors to inform on oncogenic driver genes and recurrent copy number alterations. The integration of somatic gene mutations with gene expression signatures provides further insights into genotype-phenotype interrelation in PMBL. We identified highly recurrent oncogenic mutations in the JAK-STAT and NF-kB pathways, and provide additional evidence of the importance of immune evasion in PMBL (CIITA, CD58, B2M, CD274, PDCD1LG2). Our analyses highlight the IRF-pathway as a putative novel hallmark with frequent alterations in multiple pathway members (IRF2BP2, IRF4, IRF8). In addition, our integrative analysis illustrates the importance of JAK1, RELB and EP300 mutations driving oncogenic signaling. The identified driver genes were significantly more frequently mutated in PMBL as compared to diffuse large B-cell lymphoma, whereas only a limited number of genes were significantly different between PMBL and cHL, emphasizing the close relationship between these entities. Our study, performed on a large cohort of PMBL, highlights the importance of distinctive genetic alterations for disease taxonomy with relevance for diagnostic work-up and therapeutic decision-making.

  Study EGAS00001003746

• Colorectal cancer study

  Matrix metalloproteinase-11 (MMP11) is an enzyme with proteolytic activity against matrix and non-matrix proteins. Although most MMPs are secreted as inactive proenzymes and they are later activated extracellularly, MMP11 is activated intracellularly by furin within the constitutive secretory pathway. It is a key factor in physiological tissue remodeling and its alteration may play an important role in the progression of epithelial malignancies and other diseases. TCGA colon and colorectal adenocarcinoma data showed that upregulation of MMP11 expression correlates with tumorigenesis and malignancy. Here, we provide evidence that a germline variant in the MMP11 gene (NM_005940: c.232C>T; p.(Pro78Ser)), identified by whole exome sequencing, can increase the tumorigenic properties of colorectal cancer (CRC) cells. P78S is located in the prodomain region, which is responsible for blocking MMP11’s protease activity. This variant was detected in the proband and all the cancer-affected family members analyzed, while it was not detected in healthy relatives. In silico analyses predict that P78S could have an impact on the activation of the enzyme. Furthermore, our in vitro analyses show that the expression of P78S in HCT116 cells increases tumor cell invasion and proliferation. In summary, our results show that this variant could modify the structure of the MMP11 prodomain, producing a premature or uncontrolled activation of the enzyme that may contribute to an early CRC onset in these patients. The study of this gene in other CRC cases will provide further information about its role in CRC development, which might improve patient treatment in the future.

  Study EGAS00001006489

• Robust methylation-based classification of brain tumours using nanopore sequencing

  Background: DNA methylation-based classification of cancer provides a comprehensive molecular approach to diagnose tumours. In fact, DNA methylation profiling of human brain tumours already profoundly impacts clinical neuro-oncology. However, current implementations using hybridization microarrays are time-consuming and costly. We recently reported on shallow nanopore whole-genome sequencing for rapid and cost-effective generation of genome-wide 5-methylcytosine profiles as input to supervised classification. Here, we demonstrate that this approach allows to discriminate a wide spectrum of primary brain tumours. Results: Using public reference data of 82 distinct tumour entities, we performed nanopore genome sequencing on N=382 tissue samples covering 46 brain tumour (sub)types. Using bootstrap sampling in a cohort of N = 56 cases, we find that a minimum set of 1,000 random CpG features is sufficient for high-confidence classification by ad hoc random forests. We implemented score recalibration as confidence measure for interpretation in a clinical context and empirically determined a platform-specific threshold in a randomly sampled discovery cohort (N = 185). Applying this cut-off to an independent validation series (N = 184) yielded 148 classifiable cases (sensitivity 80.4%) and demonstrated 100 % specificity. Cross-lab validation demonstrated robustness with concordant results across four laboratories in 10/11 (90.9%) cases. In a prospective benchmarking (N = 15), median time to results was 21.1 hours. Conclusions: In conclusion, nanopore sequencing allows robust and rapid methylation-based classification across the full spectrum of brain tumours. Platform-specific confidence scores facilitate clinical implementation for which prospective evaluation is warranted and ongoing.

  Study EGAS00001006540

• PD-L1 blockade in combination with carboplatin as immune induction in metastatic lobular breast cancer: the GELATO-trial

  Invasive lobular breast cancer (ILC) is the second most common histological breast cancer subtype but published data on trials specific for ILC are so far lacking. Translational research revealed that a subset of ILCs may be immune-related and more sensitive to DNA-damaging agents such as platinum. In murine ILC models, synergy between immune checkpoint blockade and platinum has been observed. Here, we tested this concept in the phase II, GELATO-trial (NCT03147040), in which patients with metastatic ILC were treated with weekly carboplatin (AUC 1.5) as immune induction treatment for 12 weeks and atezolizumab (PD-L1 blockade; every three weeks) from the third week onwards until disease progression. Four out of 23 evaluable patients had a partial response (17%, 95%CI 5-39%) and two patients had stable disease for at least 24 weeks, resulting in a clinical benefit rate of 26% (95%CI 10-48%). Out of these six patients, four patients had triple-negative ILC (TN-ILC). In serial biopsies of metastatic lesions, we observed higher CD8 T-cell infiltration, expression of immune checkpoints, and exhausted T cells upon carboplatin/PD-L1 blockade. This is the first report of a clinical trial specifically for ILC and we demonstrate promising anti-tumor activity of atezolizumab with carboplatin as immune induction, in particular for TN-ILC. While activity of carboplatin/PD-L1 blockade in classical ER+ ILC was limited, our translational data yield important insights for the design of highly needed clinical trials in ILC.

  Study EGAS00001006902

• Chronic lymphocytic leukemia IGHV somatic hypermutation detection by targeted capture next-generation sequencing

  Somatic hypermutation (SHM) status of the immunoglobulin heavy variable (IGHV) gene locus plays a crucial role in determining the prognosis and treatment of patients with chronic lymphocytic leukemia (CLL). A common approach for determining SHM status is multiplex PCR of the IGH gene followed by Sanger sequencing; however, this technique is low throughput, vulnerable to failure, and does not allow multiplexing with other molecular diagnostic assays. Here we designed and validated a DNA targeted capture approach that allows the assessment of IGHV SHM status as a sub-module of a larger next-generation sequencing (NGS) assay. We developed an end-to-end assay that takes as input raw FASTQ files and outputs a report containing IGHV SHM status and V-allele usage following the guidelines of the European Research Initiative on CLL (ERIC). We validated the approach on 35 CLL patient samples, 34 of which were previously characterized using Sanger sequencing. The targeted capture NGS panel identified the SHM status of 34 of 35 CLL patients. We showed 100% sensitivity and specificity among the 33 CLL samples with both NGS and Sanger calls. Furthermore, we demonstrated that this panel can be combined with additional targeted capture panels to detect common and prognostically important CLL single nucleotide variants (SNV) and insertion/deletion (indel) variants as well as copy number variants (del17p). Thus, a targeted capture approach to IGHV SHM detection can be easily integrated into broader sequencing panels, allowing broad CLL prognostication in a single molecular assay.

  Study EGAS00001006887

• A novel transcriptional signature identifies T-cell infiltration in high-risk paediatric cancer

  Background Molecular profiling of the tumour immune microenvironment (TIME) has enabled the rational choice of immunotherapies in some adult cancers. In contrast, the TIME of paediatric cancers is relatively unexplored. We speculated that a more refined appreciation of the TIME in childhood cancers, rather than a reliance on commonly used biomarkers such as tumour mutation burden (TMB), neoantigen load and PD-L1 expression, is an essential prerequisite for improved immunotherapies in childhood solid cancers. Methods We combined immunohistochemistry (IHC) with RNA-sequencing and whole genome sequencing across a diverse spectrum of high-risk paediatric cancers to develop an alternative, expression-based signature associated with CD8+ T-cell infiltration of the TIME. Further, we explored transcriptional features of immune archetypes, T-cell receptor sequencing diversity, assessed the relationship between CD8+ and CD4+ abundance by IHC and deconvolution predictions, and assessed common adult biomarkers such as neoantigen load and TMB. Results A novel 15-gene immune signature, Immune Paediatric Signature Score (IPASS), was identified. Using this signature, we estimate up to 31% of high-risk cancers harbour infiltrating T-cells. In addition, we showed that PD-L1 protein expression is poorly correlated with PD-L1 RNA expression and TMB and neoantigen load are not predictive of T-cell infiltration in paediatrics. Furthermore, deconvolution algorithms are only weakly correlated with IHC measurements of T-cells. Conclusions Our data provides new insights into the variable immune-suppressive mechanisms dampening responses in paediatric solid cancers. Effective immune-based interventions in high-risk paediatric cancer will require individualised analysis of the TIME.

  Study EGAS00001007029

• Treatment stratification and biomarker validation using patient-derived head and neck cancer organoids

  Background Organoids are three-dimensional structures that can be grown from patient-derived material including carcinomas. These in vitro models can be cryopreserved to generate “living” biobanks that reflect patient heterogeneity. Head and neck cancer (HNC) is a collective term used for tumors arising in the head and neck area, including squamous cell carcinomas (HNSCC) and adenocarcinoma of the salivary gland. Methods to culture patient-derived HNC organoids have previously been described. Methods Organoids were derived from resection or biopsy material of patients included between 2016 and 2022. Organoids were characterised with immunohistochemistry and whole exome sequencing (WES) and exposed to chemo-, radio- and targeted-therapies. CRISPR/Cas9-based gene editing was applied to validate proposed biomarkers. Results We have extended our previously published organoid biobank to a total of 110 organoid models, some of which were characterised by histology and WES. Organoids recapitulated both single nucleotide variants and copy number alterations as found in patient HNSCC samples. Organoid and patient response to radiotherapy (RT; primary (n=6) and adjuvant (n=15)) could be compared in 21 HNC cases. In the adjuvant RT group, organoids derived from patients that relapsed were more resistant to RT in vitro, showing potential for guiding treatment options in this subset of patients. Conclusion Taken together, these results indicate illustrate the potential of organoids as a diagnostic tool in personalised medicine for HNC and indicate the potential of this technology for biomarker discovery and validation within the context of relevant patient heterogeneity.

  Study EGAS00001007076

• A Machine Learning Tool Integrating Circulating Cell-Free DNA Methylation with Clinical Variables for Non-Invasive Diagnosis of Liver Graft Pathology

  Liver Transplant Recipients (LTRs) with elevated liver enzymes can develop graft injury due to various etiologies, such as T-cell mediated rejection (TCMR) and NASH (NASH-LT). While liver biopsy is the gold standard method for diagnosis of graft pathology, it is invasive and is associated with risks. The goal of our study was to develop a Machine Learning (ML) tool integrating clinical variables with methylation patterns on circulating DNA in plasma as a non-invasive diagnostic tool of graft injury. We generated methylation profiles of circulating DNA in a pilot study of 43 LTRs (11 NASH-LT, 19 TCMR, and 13 Control-LT), and developed an L2 multinomial logistic regression ML approach across 101 bootstrapped models to distinguish between the graft conditions. NASH was associated with distinctive methylation patterns on genes involved in fatty acid metabolism, while methylation of platelet-derived growth factors was identified in patients with TCMR. Our ML model achieved a mean multi-classification accuracy of 0.91, with mean specificity and sensitivity of 0.94 and 0.91, respectively. The model was found to be particularly adept at detecting TCMR and Control-LT, with true positive rates (TPRs) of 95% and 90%, and areas under the curve (AUROC) of 0.992 and 0.985, respectively. For NASH-LT, the models achieved a performance with a TPR of 82% and an AUROC of 0.991. These preliminary results suggest that our newly developed ML tool, leveraging cell-free DNA methylation and clinical variables, is a promising non-invasive classifier of graft pathology.

  Study EGAS00001007171

• scRNAseq for patients with immunodeficiency and HCs

  The NFAT family of transcription factors plays central roles in adaptive immunity in murine models, however, their contribution to human immune homeostasis remains poorly defined. In a multigenerational pedigree, we identified three patients carrying germline biallelic missense variants in NFATC1, presenting with recurrent infections, hypogammaglobulinemia and decreased antibody responses. The compound heterozygous NFATC1 variants identified in the patients caused decreased stability and reduced binding of DNA and interacting proteins. We observed defects in early activation and proliferation of T and B cells from these patients, amenable to reconstitution upon genetic rescue. Following stimulation, T-cell activation and proliferation were impaired, reaching that of healthy controls with delay indicative of an adaptive capacity of the cells. Assessment of the metabolic capacity of patient T cells, revealed that NFATc1-dysfunction rendered T cells unable to engage in glycolysis following stimulation, although oxidative metabolic processes were intact. We hypothesized that NFATc1-mutant T cells could compensate for the energy deficit due to defective glycolysis by enhanced lipid metabolism as an adaptation, leading to a delayed, but not lost activation responses. Indeed, we observed increased 13C-labelled palmitate incorporation into citrate indicating higher fatty acid oxidation and we demonstrated that metformin and rosiglitazone improved patient T-cell effector functions. Collectively, enabled by our molecular dissection of NFATC1 mutations and extending the role of NFATc1 in human immunity beyond receptor signaling, and reveal evidence of metabolic plasticity in the context of impaired glycolysis observed in patient T cells to remedy delayed effector responses.

  Study EGAS00001007271

• Analysis of the Elements Involved in the Enrichment of a Panel of Genomic Regions by Nanopore Sequencing Using Adaptive Sampling

  DNA sequencing using nanopore technologies with the affordable MinION device is useful for the identification and characterization of structural variants, long haplotypes, sequencing of repetitive regions and identification of epigenetic modifications. The main limitation of this approach is the low coverage obtained, which might be avoided by adaptive sampling, a computationally controlled method of enrichment for targeted genomic regions. This study dissects the factors involved in the enrichment by adaptive sampling of a panel of 18 human genome regions containing 20 genes implicated in breast cancer sequenced in 16 patients with familial breast cancer negative for NGS screening. An average coverage of 2.0x was obtained for the whole genome and 5.1x for selected regions. Sequencing time was the main factor improving coverage. The selection of long reads (>1 Kb) did not improve the enrichment. The length of the selected region, which in our study ranged from 126 to 565 Kb, did not play a significant role in enrichment. However, the region containing PMS2 showed significantly lower coverage, which could be explained by the high number of PMS2 pseudogenes (N = 14), which were also enriched. Our study shows new evidence of enrichment obtained by adaptive sampling in a panel of genomic regions and shows parameters, the relevance of sequencing time and the role of pseudogenes, that improve the enrichment yield with no library reloading or GPU use, data useful for a more efficient application of this procedure in future studies.

  Study EGAS00001007375

• Efficacy of nivolumab in pediatric cancers with high mutation burden and mismatch-repair deficiency

  Purpose: Checkpoint-inhibitors have limited efficacy for children with unselected solid and brain tumors. We report the first prospective pediatric trial (NCT02992964) using nivolumab exclusively for refractory non-hematological cancers harboring tumor mutation burden (TMB) ≥5 mutations/megabase (mut/Mb) and/or mismatch-repair deficiency (MMRD). Patients and methods: Twenty patients were screened, and ten ultimately included in the response cohort of whom nine had TMB >10mut/Mb (three initially eligible based on MMRD) and one patient had TMB between 5-10 mut/Mb. Results: Delayed immune responses contributed to best overall response of 50%, improving on initial objective responses (20%) and leading to 2-year overall survival (OS) of 50% (95% CI; 27, 93). Four children, including three with refractory malignant gliomas are in complete remission at a median follow-up of 37-months (range: 32.4-60), culminating in 2-year OS of 43% (95% CI; 18.2, 100). Biomarker analyses confirmed benefit in children with germline MMRD, microsatellite instability, higher activated and lower regulatory circulating T-cells. Stochastic mutation accumulation driven by underlying germline MMRD impacted the tumor microenvironment, contributing to delayed responses. No benefit was observed in the single patient with a MMR-proficient tumour and TMB 7.4 mut/Mb. Conclusions: Nivolumab resulted in durable responses and prolonged survival for the first time in a pediatric trial of refractory hypermutated cancers including malignant gliomas. Novel biomarkers identified here need to be translated rapidly to clinical care to identify children who can benefit from checkpoint-inhibitors, including for upfront management of their cancers.

  Study EGAS00001007393

• shallow WGS and deep targeted panel on ctDNA in rhabdomyosarcomas

  The genomic spectrum of rhabdomyosarcoma (RMS) progression from primary to relapse is not fully understood. In this study we investigate 35 patients with relapsed RMS from two contributing institutions, 18 fusion-positive (FP-RMS) and 17 fusion-negative RMS (FN-RMS). Targeted DNA or whole exome sequencing (WES) was used to detect alterations in paired primary/relapsed samples. In 10 cases, circulating tumor DNA (ctDNA) from multiple timepoints through clinical care and progression was analyzed for feasibility of liquid biopsy in monitoring treatment response/relapse. ctDNA alterations were evaluated using a targeted custom RMS panel (36 genes) at high coverage for single nucleotide variation and fusion detection, and a shallow whole genome sequencing for copy number variation. FP-RMS had a stable genome with relapse, with the most common secondary alterations : CDKN2A/B, MYCN and CDK4 alterations, being already present at diagnosis and impacting overall survival. FP-RMS lacking major secondary events at baseline acquired recurrent MYCN and AKT1 alterations. FN-RMS acquired a higher number of new alterations, most commonly SMARCA2 missense mutations. ctDNA analyses detected pathognomonic variants in all RMS patients at diagnosis, regardless of FP/FN or type of alterations, while at relapse selected alterations were confirmed in 86% of FP-RMS and 100% FN-RMS. Moreover, a higher number of fusion reads was detected with increased disease burden and at relapse in patients following a fatal outcome. These results underscore patterns of tumor progression within a relatively stable genomic landscape and provide rationale for using liquid biopsy to monitor treatment response.

  Study EGAS00001007399

• Clinical and genetic factors associated with tumor response to neoadjuvant (chemo)radiotherapy, survival and recurrence risk in rectal cancer

  Rectal cancer poses challenges in preoperative treatment response, with up to 30% achieving a complete response (CR). Personalized treatment relies on accurate identification of responders at diagnosis. This study aimed to unravel CR determinants, overall survival (OS), and time to recurrence (TTR) using clinical and targeted sequencing data. Analyzing 402 patients undergoing preoperative treatment, tumor stage, size, and treatment emerged as robust response predictors. CR rates were higher in smaller, early-stage, and intensively treated tumors. Targeted sequencing analyzed 216 cases, while 120 patients provided hotspot mutation data. KRAS mutation dramatically reduced CR odds by over 50% (odds ratio [OR]=0.3 in the targeted sequencing and OR=0.4 hotspot cohorts, respectively). In contrast, SMAD4 and SYNE1 mutations were associated with higher CR rates (OR=6.0 and 6.8, respectively). Favorable OS was linked to younger age, CR, and low baseline carcinoembryonic antigen levels. Notably, CR and an APC mutation increased TTR, while a BRAF mutation negatively affected TTR. Beyond tumor burden, SMAD4 and SYNE1 mutations significantly influenced CR. KRAS mutations independently correlated with radiotherapy resistance, and BRAF mutations heightened recurrence risk. Intriguingly, nonresponding tumors with initially small sizes carried a higher risk of recurrence. These findings offer insights into rectal cancer treatment response, guiding personalized therapeutic strategies. By uncovering factors impacting CR, OS, and TTR, this study underscores the importance of tailored approaches for rectal cancer patients. These findings, based on extensive analysis and mutation data, pave the way for personalized interventions, optimizing outcomes in the challenges of rectal cancer preoperative treatment.

  Study EGAS00001007501

• Peripheral blood DNA methylome in adalimumab-treated patients with rheumatoid arthritis

  Background and aims Rheumatoid arthritis (RA) patients are currently treated with biological agents mostly aimed at cytokine blockade such as anti-TNFα therapy. Currently, there are no biomarkers to predict therapy response to these agents. Here, we aimed to predict response to adalimumab (ADA) treatment in RA patients using DNA methylation in peripheral blood Methods DNA methylation profiling on whole peripheral blood from 92 RA patients before the start of adalimumab (ADA) was determined using Illumina HumanMethylationEPIC BeadChip Array. After 6 months, treatment response was assessed according to the Alliance of Associations for Rheumatology (EULAR) criteria for disease response. Patients were classified as responders (DAS28<3.2 or decrease of 1.2 points) or as non-responders (DAS28>5.1 or decrease of less than 0.6 points). Machine learning models were built with gradient boosting and stability selection models to predict response prior to ADA treatment with predictor DNA methylation markers. Results We demonstrated a 27-feature panel of DNA methylation markers or CpG classifiers that predict response in RA patients treated with ADA. Forty-nine patients were assigned as responder and 43 patients assigned as non-responders. We differentiate responders from non-responders with a high sensitivity (AUC 0.76). The predictor CpGs annotated to genes involved in immunological- and pathophysiological pathways related to RA such as T-cell signaling, B-cell pathology and angiogenesis. Conclusion Our findings indicate that DNA methylation signatures discriminate responders and non-responders to ADA treatment and can serve as a tool for therapy prediction.

  Study EGAS00001007578

• Pharmacogenomic Analysis Reveals New Therapeutic Options for Pleural Mesothelioma

  Background: Pleural mesothelioma (PM) urgently requires novel effective treatments. Drug repurposing is a relevant strategy to discover new therapeutic options among well-known drugs used in clinical practice. This study aimed to identify new therapies in the context of molecular heterogeneity of PM. Methods: We performed a multi-omics study of the largest cohort of patient-derived primary PM cell lines (n=58) to date, using whole exome sequencing, RNA sequencing and methylome analysis. We also conducted a large-scale drug screening of over 1300 compounds on a subset of cell lines (n=11), as well as an extension screen (n=48) for selected drugs of interest. The best candidates were validated using 3D culture and mouse models. Results: We demonstrated that the molecular landscape of the cell lines cohort recapitulates the main findings in tumours and revealed new specific features of PM. The large-scale drug screening and subsequent validations highlighted notably strong sensitivity to entinostat and fluvastatin. Importantly, both compounds displayed same efficacy in 2D and 3D culture models and single treatment prolonged survival in vivo similarly or even better than standard chemotherapy. Strikingly, combining entinostat treatment with anti-PD1 led to cures in some mice and immunized them against tumour cells. Conclusions: Among several candidates, we identified entinostat, notably in combination with immune checkpoint inhibitor, and to a lesser extent fluvastatin, as promising therapeutic options for patients. We believe our study provides a valuable resource that could contribute to advances in the field and warrants further clinical investigation.

  Study EGAS00001007866

• Capturing disease severity in LIS1-lissencephaly reveals proteostasis dysregulation in patient-derived forebrain organoids

  Lissencephaly is a developmental cortical malformation characterized by reduced to absent gyri and a disorganized cortex, often leading to severe impairments in affected individuals and a reduced life expectancy. Heterozygous mutations in the LIS1 gene, encoding a regulator of the microtubule motor dynein, cause lissencephaly with different clinical severities. While the clinical disease spectrum correlates with the degree of lissencephaly, location and type of mutation may not. We leveraged forebrain-type organoids from LIS1-lissencephaly patients, diagnosed with mild, moderate or severe lissencephaly to investigate, in a cytoarchitecture and multi-omics approach, how the severity degree in patients might relate to specific mutations in the LIS1 gene. We questioned which processes during cortical development might be differentially affected by severity grade, and whether they could be pharmacologically targeted. We found alterations in neurodevelopment often with a severity-dependent gradient. Specifically, we identified alterations of the cytoarchitecture, progenitor cell homeostasis and neurogenesis. Particularly important disease-linked molecular mechanisms were microtubule destabilization, WNT-signaling, and perturbed cadherin- and unfolded protein-binding. Some mechanisms exhibited a severity-dependent gradient, or were specific to a severe grade. We present strategies to reverse phenotypic changes in LIS1-patient organoids, and an in silico approach with therapeutic potential. Thus, we show that different LIS1-severity grades can be recapitulated in vitro, that there is a direct link between the phenotype and genotype, that organoid-based disease modeling can identify molecular underpinnings of malformations of cortical development and that organoids provide a valid platform to develop and test therapeutic strategies.

  Study EGAS00001008227

• Immunodeficiency syndrome caused by LCP1 mutations

  Severe congenital neutropenia (SCN) is characterized by neutropenia, recurrent infections and an increased leukemia risk. Multiple genetic defects underlying SCN have been identified, but a genetic diagnosis is still lacking in a significant proportion of patients. Here, we report three independent pedigrees with autosomal dominant variants in LCP1. Variants c.740-1G>T and c.740-20_744del resulted in the same alternatively spliced RNA product, causing an in-frame deletion (p.A247_E254del). Variant c.509C>T in the third pedigree resulted in p.S170L. Variant carriers suffer from neutropenia, poor or complete lack of response to G-CSF treatment, and variable degrees of lymphopenia. Patients with A247_E254del presented with tetraploid cells in the bone marrow, indicative of disturbed cytokinesis. In one of these kindreds, two individuals developed acute leukemia. G-CSF non-responsiveness and defective cell cycling were repaired upon correction of the LCP1 A247_E254del variant in patient-derived induced pluripotent stem cells, supporting the monogenic origin of the disease. Indicative of their gain-of-function, both A247_E254del and S170L variants increased F-actin bundling and formation of abnormal protrusions. Single-cell transcriptome analysis of A247_E254del bone marrow-derived hematopoietic stem and progenitor cells (HSPCs) showed deregulation of signaling pathways controlling mitosis in multi-lineage and lymphoid-primed HSPC subsets. We conclude that activating LCP1 variants cause a new autosomal dominant hematopoietic disorder. Depending on the consequences of the LCP1 variants for protein structure, patients may suffer from G-CSF refractory severe neutropenia, lymphopenia, defective cytokinesis and predisposition to leukemia.

  Study EGAS00001008293

• Gut metagenome from Estonian Microbiome Cohort

  Stool samples were collected from 2,509 Estonian Biobank participants. The shotgun metagenomic paired-end sequencing was performed by Novogene Bioinformatics Technology Co., Ltd. using the Illumina NovaSeq6000 platform, resulting in 4.62 ± 0.44 Gb of data per sample (insert size, 350 bp; read length, 2 × 250 bp). A total of 2,513 samples belonging to 2,509 individuals were sequenced, including 4 biological replicates from one individual. First, the reads were trimmed for quality and adapter sequences. The host reads that aligned to the human genome were removed using SOAP2.21 (parameters: -s 135 -l 30 -v 7 -m 200 -x 400).

  Dataset EGAD00001008448

• Myeloma Genome Project Targeted Panel Validation dataset

  372 samples consisting of 185 patient paired CD138+ tumor and non-involved DNA pairs, plus 5 Horizon Diagnostic known mutation standards (HD). Samples were processed using the KAPA HyperCap protocol and hybridized onto a targeted panel for multiple myeloma and associated diseases. Reference Sudha et al Clinical Cancer Research, 2022.

  Dataset EGAD00001008689

• Clonal selection after gene therapy in sickle cell disease

  While gene therapy (GT) provides a potentially curative treatment option for patients with sickle cell disease (SCD), the occurrence of myeloid malignancies in clinical trials has prompted concern. To interrogate potential mechanisms underlying increased cancer risk, we used hematopoietic stem cell (HSC) clonal tracking by whole genome sequencing (WGS) to map the somatic mutation and clonal landscape of 2,592 gene modified as well as unmodified single stem and progenitor cells from six SCD patients undergoing gene therapy (7-26 years old, average 12.7× depth). Pre-GT phylogenetic trees in SCD were highly polyclonal and mutation burdens per cell were elevated in some, but not all, patients. Post-GT, no clonal expansions were identified. However, an increased frequency of driver mutations associated with myeloid neoplasms or clonal hematopoiesis (DNMT3A- and EZH2-mutated clones in particular) were seen in both genetically modified and unmodified cells suggested positive selection of mutant clones during gene therapy. This work sheds light on the mutation landscape and HSC clonal dynamics in gene therapy for SCD and highlights enhanced fitness of some HSCs harboring pre-existing driver mutations following gene therapy. Future studies should define the long-term fate of mutant clones including any contribution to expansions associated with myeloid neoplasms.

  Dataset EGAD00001010913

• RNAseq data for Failure of Differentiation of the Rhombic Lip Constitutes Medulloblastoma

  We used single-cell transcriptomics to study cells from the developing human cerebellum, and show that different molecular subgroups of medulloblastoma resemble distinct glutamatergic progenitors.

  Dataset EGAD00001008458

• KiCS WGS data for academic and for-profit use

  Data from samples that are marked for both academic and for-profit use.

  Dataset EGAD00001009699

• Irish Covid19 WGS Alignment Data

  Dataset contains 483 Irish origin Individuals with Covid19. For WGS, Alignment has been done using BWA-mem (sention v 201808.03 ) and BAM file generated.

  Dataset EGAD00001010098

• Clinical relevance of somatic copy-number alterations in plasma circulating tumor DNA from advanced EGFR-mutated NSCLC patients treated with osimertinib

  Background: Somatic copy-number alterations (SCNAs) are associated with drug resistance in patients with advanced epidermal growth factor receptor (EGFR)-mutated non-small cell lung cancer (NSCLC) treated with EGFR-tyrosine kinase inhibitors (TKIs). The purpose of our study was to assess the clinical relevance and development of genome-wide SCNAs in plasma circulating tumor DNA (ctDNA) under osimertinib therapy.Methods: We included 43 patients with advanced EGFR T790M-positive lung adenocarcinoma who were treated with osimertinib after progression under previous EGFR-TKI therapy. We performed genomic profiling of ctDNA in plasma samples from each patient obtained pre-osimertinib and after patients developed resistance to osimertinib. SCNAs were detected by shallow whole-genome plasma sequencing and EGFR mutations were assessed by droplet digital PCR.Results: SCNAs in resistance-related genes (rrSCNAs) were detected in 10 out of 39 (25.6%) evaluable patients before start of osimertinib. The presence of rrSCNAs in plasma before initiation of osimertinib therapy was associated with a lower response rate to osimertinib (50% versus 86%, p=0.02) and was an independent predictor for shorter progression-free survival (adjusted HR 3.26, 95% CI 1.37-7.75, p=0.008) and overall survival (adjusted HR 2.41, 95% CI 1.09-5.34, p=0.03).Conclusions: Genomic profiling of plasma ctDNA is clinically relevant and affects the efficacy and clinical outcome of osimertinib. Our approach allows comprehensive assessment of SCNAs in plasma samples of lung adenocarcinoma patients and may help to guide genotype-specific therapeutic strategies in the future.

  Study EGAS00001004539

• Combination of ribociclib and gemcitabine for the treatment of medulloblastoma

  Group 3 (G3) medulloblastoma (MB) is one of the deadliest forms of the disease for which novel treatment is desperately needed. Here we evaluate ribociclib, a highly selective CDK4/6 inhibitor, with gemcitabine in mouse and human G3 MBs. Ribociclib CNS penetration was assessed by in vivo microdialysis and by immunohistochemistry and gene expression studies. Survival studies to determine the efficacy of ribociclib and gemcitabine combination were performed on mice orthotopically implanted with luciferase labelled mouse and human G3 MB. Pharmacokinetic-pharmacodynamic outcomes and univariable survival models were analyzed to estimate survival. Gene activity inference using NetBID and tumor differentiation analysis investigated the effects of the combination after short and long-term treatments. Tumors from mice treated with oral ribociclib displayed inhibited RB phosphorylation, downregulated E2F target genes, and decreased proliferation. Treatment of mice with the combination of ribociclib and gemcitabine was well tolerated, slowed tumor progression and metastatic spread, and increased survival. Molecular analysis of treated versus untreated tumors showed a significant decrease in the activity and expression of genes involved in cell cycle progression and DNA damage response, and an increase in activity and expression of genes implicated in neuronal identity and neuronal differentiation. Ribociclib is CNS-penetrant. When administered/combined with gemcitabine in orthotopic G3 MB models resulted in improved survival. Our findings, with both mouse and human patient-derived-orthotopic xenograft models, suggest that this combination therapy has promise for children with G3 MB and may represent an effective treatment strategy for other CNS malignancies.

  Study EGAS00001006001

• Implementation of pediatric precision oncology into clinical practice: The individualized Therapies for Children with cancer program “iTHER”

  iTHER is a prospective national precision oncology program aiming to define tumor molecular profiles in children and adolescents with primary very high-risk, relapsed, or refractory pediatric tumors in order to identify relevant aberrations to inform treatment. From April 2017 to April 2021, 302 samples of 253 patients were included. Comprehensive molecular profiling utilizing lcWGS, WES, RNAseq, Affymetrix, and/or 850k methylation profiling was performed for 226 samples. Germline pathogenic variants were identified in 16% of patients. At least one somatic alteration was detected in 204 (90.3%), and 185 (81.9%) were considered druggable, with priority levels very high (6.1%), high (21.3%), moderate (26.0%), intermediate (36.1%), and borderline (10.5%). Diagnosis was revised or refined in 8 patients (3.5%). Spatial and temporal heterogeneity was observed in paired samples of 15 patients, indicating the value of sequential analysis. Of 137 patients with follow-up beyond 12 months, 21 molecularly matched treatments were applied in 19 patients (13.9%) with modest response. Most relevant barriers to not applying targeted therapies included poor performance status, as well as limited access to drugs within clinical trial. iTHER demonstrates feasibility of comprehensive molecular profiling across all ages, tumor types and stages in pediatric cancers, informing of diagnostic, prognostic, and targetable alterations as well as reportable germline variants within a clinically relevant timeframe. Therefore, WES and RNAseq is implemented into standard clinical care at the Princess Máxima Center for all children. Improved access to targeted treatments within biology-driven combination trials is required to eventually improve survival.

  Study EGAS00001007099

• Whole_genome_sequencing_of_100_DDD_trios_with_suspected_noncoding_causal_mutations

  We would like to whole-genome sequence 100 trios (300 individuals) from the DDD project in which we suspect the proband has a pathogenic noncoding mutation, to 30X coverage. These have been carefully selected to maximise the chances of finding pathogenic noncoding mutations. Specifically, we have selected probands who have no likely causal mutation within the coding region but who have a heterozygous loss-of-function mutation inherited from one parent in a known recessive developmental disorder gene, and whose phenotype is a good match for the typical phenotype caused by recessive variants in that gene. We suspect that these probands may have a heterozygous noncoding mutation on the other allele which reduces expression of the gene so that, together, these two mutations are effectively equivalent to a homozygous loss-of-function mutation. About a quarter of the probands have affected siblings and, if those siblings have been exome sequenced, we have checked that they are identical-by-descent in the region of the gene. We propose to search in the promoter region, intronic and annotated regulatory regions around these genes of interest for rare variants inherited from the other parent (i.e. not the one who transmitted the loss-of-function variant), or de novo variants. If we find such variants, and demonstrate an enrichment of this compound heterozygous LoF noncoding mechanism above that expected by chance, we intend to follow up specific candidates by RNAseq and potentially CRISPR/Cas9 experiments. 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 please see http://www.sanger.ac.uk/datasharing/

  Study EGAS00001002452

• A proof-of-concept study of sequential treatment with the HDAC inhibitor vorinostat following BRAF and MEK inhibitors in BRAFV600 mutated melanoma

  "Purpose: Development of resistance limits the clinical benefit of BRAF and MEK inhibitors (BRAFi/MEKi) in BRAFV600 mutated melanoma. It has been shown that short-term treatment (14 days) with vorinostat was able to initiate apoptosis of the resistant tumor cells. We aimed to assess the anti-tumor activity of sequential treatment with vorinostat following BRAFi/MEKi in patients with BRAFV600 melanoma who progressed after initial response to BRAFi/MEKi. Patients and Methods: Patients with BRAFi/MEKi resistant BRAFV600 melanoma were treated with vorinostat 360 mg QD for 14 days followed by BRAFi/MEKi. The primary endpoint was an objective response rate of progressive lesions of at least 30% according to RECIST 1.1. Secondary endpoints included progression-free survival (PFS), overall survival (OS), safety, pharmacokinetics of vorinostat and translational molecular analyses using ctDNA and tumor biopsies. Results: Twenty-six patients with progressive BRAFi/MEKi resistant BRAFV600 mutated melanoma received treatment with vorinostat. Twenty-two patients were evaluable for response. The ORR was 9% (one complete response for 31.2 months and one partial response for 14.9 months. Median PFS and OS were 1.4 and 5.4 months, respectively. Common adverse events were fatigue (23%) and nausea (19%). ctDNA analysis showed emerging secondary mutations in NRAS and MEK in eight patients at time of BRAFi/MEKi resistance. Elimination of these mutations by vorinostat treatment was observed in three patients. Conclusions: Intermittent treatment with vorinostat in patients with resistant BRAFV600 mutated melanoma is well tolerated. Although the primary endpoint of this study was not met, durable anti-tumor responses were observed in a minority of patients (9%). "

  Study EGAS00001007709

• Toxigenic Clostridium perfringens isolated from at-risk pediatric inflammatory bowel disease patients

  With the goal of studying mucosal-associated bacteria and products detrimental to inflammatory bowel disease (IBD) patients, we directly cultured bacterial communities from mucosal biopsies and examined for detrimental traits. Focusing on red blood cell lysis, a virulence trait, communities from 5 of 9 patients exhibited hemolytic activity. Clostridium perfringens was isolated from each hemolytic community, and whole-genome sequencing confirmed presence of hemolytic toxins – alpha toxin in all 5 strains, perfringolysin O (PFO) in 4 strains, alveolysin in 2 strains – and absence of enterotoxins. In vitro cytotoxicity experiments demonstrated strain supernatants were toxic to various cell types located beneath the intestinal epithelial barrier, including endothelial, neuroblasts, and neutrophils, while impact on epithelial cells was less pronounced. These findings suggest C. perfringens may be damaging particularly when intestinal barrier integrity is compromised, such as in IBD patients. Further characterization using purified toxins and genetic insertion mutants confirmed PFO toxin was sufficient for cell toxicity. C. perfringens supernatants also induced activation of neuroblast and dorsal root ganglion (DRG) neuron cells in vitro, suggesting C. perfringens in inflamed mucosal tissue may directly contribute to abdominal pain, a frequent IBD symptom. Strikingly, C. perfringens was detected in biopsies of 8 of 9 pediatric patients while not in all corresponding stool samples, and in other datasets, C. perfringens was in stool samples of IBD adults (18.7-27.1%) much more frequently versus healthy (5.1%). Our findings showcase C. perfringens may have an important pathogenic impact on IBD patients. Similar to Clostridioides difficile, routine monitoring of C. perfringens and PFO toxins in patients may be beneficial.

  Study EGAS50000000334

• METABRIC

  Solid tumors are complex tissues composed of a mixture of cancer and normal cells, which complicates the interpretation of their molecular profiles. Normal cell contamination can dilute cancer cell information and tissue architecture is generally not reflected in molecular assays. To address these challenges, we developed a computational approach based on standard Haematoxylin and Eosin-stained sections and demonstrated its power in a discovery cohort of 323 breast tumors and an independent validation cohort of 241 tumors. First, to deconvolute cellular heterogeneity and detect subtle genomic aberrations, we introduced an algorithm based on tumor cellularity to increase the comparability of copy-number profiles between samples. Second, we demonstrated that a predictor for survival integrating image-based and molecular features significantly outperforms classifiers based on single data types. Third, we described and validated a novel, independent prognostic factor based on quantitative analysis of spatial patterns between stromal cells, which are not detectable by molecular assays. Our quantitative methods refine and complement molecular assays of tumor samples and could benefit all large-scale cancer studies.First, to deconvolute cellular heterogeneity and detect subtle genomic aberrations, we introduced an algorithm based on tumor cellularity to increase the comparability of copy-number profiles between samples. Second, we demonstrated that a predictor for survival integrating image-based and molecular features significantly outperforms classifiers based on single data types. Third, we described and validated a novel, independent prognostic factor based on quantitative analysis of spatial patterns between stromal cells, which are not detectable by molecular assays. Our quantitative methods refine and complement molecular assays of tumor samples and could benefit all large-scale cancer studies.

  Study EGAS00000000098

• Cell-Free DNA Genomic and Fragmentomic Features for Early Outcome Prediction in Diffuse Large B-Cell Lymphoma

  Purpose: Diffuse large B-cell lymphoma (DLBCL) patients need an accurate and early risk stratification strategy as prompt therapy escalation may improve outcome. Patients and methods: We evaluated cell free DNA (cfDNA) genomic and fragmentomic features in 190 patients with large B-cell lymphoma from more than 40 hospitals in the Benelux using Whole Genome Sequencing (WGS). We defined a new metric called ACT score (Aberrations, Contribution of short fragments, Terminal motif analyses) by combining four cfDNA signals and evaluated its performance to predict end-of-treatment (EOT) response and survival. Results: Individual cfDNA features and ACT score were altered per EOT response after one cycle of treatment (T1). The ACT score at T1 outperformed all individual cfDNA features in predicting EOT response (AUC = 0.74). Patients with a positive ACT score had an inferior outcome compared to ACT score negative patients [progression-free survival (PFS): HR 5.0 (95% CI 3.1-8.2), log-rank test, p < 0.0001) and overall survival (OS): HR 6.5 (95% CI 3.5-12.2), log-rank test, p < 0.0001]. The 2-year PFS in ACT score positive and negative patients were 29.8% and 78.1%, respectively. The prognostic value of the ACT score is independent of the International Prognostic Index and interim imaging. Conclusions: ACT score computed from a single plasma sample collected after one cycle of treatment can predict clinical outcome. This low-cost and easy-to-interpret test does not require tissue biopsies or a priori knowledge of mutations and has the potential to aid patient selection in interventional clinical trials and risk-adapted treatment strategies.

  Study EGAS50000000412

• Phase 1/2 Study of the Indoleamine 2,3-Dioxygenase 1 Inhibitor Linrodostat Mesylate Combined With Nivolumab or Nivolumab and Ipilimumab in Advanced Solid Tumors or Hematologic Malignancies

  Purpose: To evaluate the safety, efficacy, pharmacokinetics, pharmacodynamics, and biomarker data of linrodostat mesylate, a selective, oral indoleamine 2,3-dioxygenase 1 (IDO1) inhibitor, combined with nivolumab ± ipilimumab in advanced solid tumors and hematologic malignancies. Patients and Methods: In this phase 1/2 study, patients received once-daily (QD) linrodostat (part 1 [escalation], 25-400 mg; part 2 [expansion], 100 or 200 mg) plus nivolumab (480 mg every [Q] 4 weeks [W] or 240 mg Q2W) or triplet therapy (part 3, linrodostat 20-100 mg QD; nivolumab 360 mg Q3W or 480 mg Q4W; ipilimumab 1 mg/kg Q6W or Q8W). Endpoints included safety and efficacy (co-primary; parts 2, 3), pharmacokinetics, pharmacodynamics, biomarkers, and efficacy (part 1). Results: Fifty-five, 494, and 41 patients were enrolled in parts 1, 2, and 3, respectively. Linrodostat exposures exceeded predicted therapeutic target concentrations starting at 50 mg. Rates of grade 3/4 adverse events were 50.1%-63.4%. The maximum tolerated linrodostat dose was 200 mg; dose-limiting toxicities were primarily immune related. Responses were observed across different cohorts, study parts, and tumor types, particularly in immunotherapy-naive patients. Kynurenine decreased with linrodostat + nivolumab, regardless of response. In contrast, interferon gamma (IFN-γ) gene expression signature was associated with response; in nonmelanoma patients, a composite of low tryptophan 2,3-dioxygenase (TDO2) gene expression plus high IFN-γ signature was associated with response. Conclusions: Linrodostat + nivolumab ± ipilimumab demonstrated a manageable safety profile. Kynurenine changes supported IDO1 pathway inhibition but did not correlate with response. A composite biomarker of low TDO2 expression plus high IFN-γ gene expression may predict response to linrodostat + nivolumab.

  Study EGAS50000000710

• Luminal breast epithelial cells from wildtype and BRCA mutation carriers harbor copy number alterations commonly associated with breast cancer

  Cancer-associated mutations have been documented in normal tissues, but the prevalence and nature of somatic copy number alterations and their role in tumor initiation and evolution is not well understood. Here, using single cell DNA sequencing, we describe the landscape of CNA variants in >48,000 breast epithelial cells from women with normal or high risk of developing breast cancer. Individual cells with one or two of a specific subset of CNAs (e.g. 1q gain and 16q, 22q, 7q, and 10q loss) are detectable in almost all breast tissues and, in those from BRCA1 or BRCA2 mutations carriers, occurs prior to loss of heterozygosity (LOH) of the wildtype alleles. These CNAs, which are among the most common associated with ductal carcinoma in situ (DCIS) and malignant breast tumors, are enriched almost exclusively in luminal cells and are generally absent from basal myoepithelial cells. Allele-specific analysis of the enriched CNAs reveals that each allele was independently altered, demonstrating convergent evolution of these CNAs in an individual breast. Tissues from BRCA1 or BRCA2 mutation carriers contained a small percentage of cells with extreme aneuploidy, featuring loss of TP53, LOH of BRCA1 or BRCA2, and multiple breast cancer-associated CNAs in addition to one or more of the common CNAs in 1q, 10q or 16q. Notably, cells with intermediate levels of CNAs were not detected, arguing against a stepwise gradual accumulation of CNAs. Overall, our findings demonstrate that chromosomal alterations in normal breast epithelium including from high risk women, link early stage CNAs to clonally expanded tumor genomes .

  Study EGAS00001007716

• Bone marrow breakout lesions act as key sites for tumor-immune cell diversification and exhaustion in multiple myeloma

  Multiple myeloma is a disease characterized by the expansion of cancerous plasma cells in the bone marrow. The bone marrow microenvironment plays a pivotal role in supporting myeloma growth and modulating tumor immunity. As the disease progresses, myeloma cells may become independent of the bone marrow environment, leading to extramedullary disease associated with poor prognosis. However, the early processes associated with bone marrow independence and its implications for disease and immune control remain poorly understood. Here, we employed comprehensive single-cell and spatial mapping to identify the disruption of the cortical bone and subsequent expansion of myeloma cells in breakout lesions as a key event in multiple myeloma pathogenesis and tumor immunity. Breakout lesions harbor an adapted niche and a unique immune microenvironment conferred by sustained immune-tumor interactions. In particular, tumor-reactive T cells with a highly exhausted phenotype, alongside unique NK cell and macrophage subtypes, specifically expanded in breakout lesions. Spatially-resolved genomic, transcriptomic and cellular analyses uncovered extensive intra-lesion heterogeneity, suggesting a divergent co-evolution of genomic variation and locally confined T cell responses in distinct subregions. Ultra-high plex imaging revealed that tumor-reactive T cells, together with other immune-regulatory cells and dendritic cells, co-expand and exhaust upon sustained tumor interactions in locally confined immune islands. In contrast, immune-regulatory macrophage and NK cell subsets operate spatially separated within myeloma cell deserts. Jointly, our analyses uncover bone marrow breakout lesions as a hotspot for tumor-immune cell interactions, diversification and exhaustion, representing a key event in myeloma pathogenesis with significant therapeutic implications.

  Study EGAS50000000304

• Discovery and validation of an ancillary genomic test of malignancy for primary melanocytic tumors

  Background:Pathological diagnosis of melanocytic tumors can be difficult and prone to error. More accurate and reliable ancillary investigations than those currently available are required to improve diagnostic accuracy. Goals: (1) to generate a genomic model for distinguishing benign (nevi) from malignant (melanoma) melanocytic skin tumors, developed from a customized gene panel evaluated on a cohort of 250 nevi and melanomas. (2) To evaluate the performance of the genomic model in melanocytic tumors of uncertain diagnosis/behavior. Details: The cohort of 250 melanomas and naevi was split into a discovery set and a validation set. Our genomic model obtained an area under the curve of 95% in the discovery cohort and 96% for the validation cohort. Based on a Youden index threshold, the resulting classifier shows high specificity (>95% in the discovery and validation cohorts), with sensitivities >80% in the discovery and validation cohorts. Evaluation of the genomic model on a cohort of 110 melanocytic tumors identifies key driver mutations in agreement with the likely pathway of origin of the tumors. The probability of melanoma (p=0.02) from the genomic model was significantly associated with pathology-based estimates of malignancy for melanocytic tumors lacking pathway-defining genomic aberrations (conventional tumors). The genomic model identified 8/36 (22%) conventional tumors with unresolved malignancy to have a high probability of melanoma, with 7/8 of these cases showing negative/non-aberrant results with either fluorescent in-situ hybridization, PRAME and/or p16 immunostains. Furthermore, no evidence of malignant behavior based on the genomic model was seen in 4 conventional borderline tumors with more than 1 year follow-up. Conclusion: This genomic model can become a clinically useful ancillary tool that is highly specific for differentiating melanomas from nevi.

  Study EGAS50000000887

• Cholesterol homeostasis and lipid raft dynamics at the basis of tumor-induced immune dysfunction in chronic lymphocytic leukemia

  Autologous T-cell therapies show limited efficacy in chronic lymphocytic leukemia (CLL), where acquired immune dysfunction prevails. In CLL, disturbed mitochondrial metabolism has been linked to defective T-cell activation and proliferation. Recent research suggests that lipid metabolism regulates mitochondrial function and differentiation in T cells, yet its role in CLL remains unexplored. This comprehensive study compares T-cell lipid metabolism in CLL patients and healthy donors, revealing critical dependence on exogenous cholesterol for human T-cell expansion following TCR-mediated activation. Using multi-omics and functional assays, we found that T cells present in viably frozen samples of patients with CLL (CLL T cells) showed impaired adaptation to cholesterol deprivation and inadequate upregulation of key lipid metabolism transcription factors. CLL T cells exhibited altered lipid storage, with increased triacylglycerols and decreased cholesterol, and inefficient fatty acid oxidation (FAO). Functional consequences of reduced FAO in T cells were studied using samples from patients with inherent FAO disorders. Reduced FAO was associated with lower T-cell activation but did not affect proliferation. This implicates low cholesterol levels as a primary factor limiting T-cell proliferation in CLL. Consequently, CLL T cells displayed fewer and less clustered lipid rafts, potentially explaining the impaired immune synapse formation observed in these patients. Our findings highlight significant disruptions in lipid metabolism as drivers of functional deficiencies in CLL T cells, underscoring the pivotal role of cholesterol in T-cell proliferation. This study suggests that modulating cholesterol metabolism could enhance T-cell function in CLL, presenting novel immunotherapeutic approaches to improve outcome in this challenging disease.

  Study EGAS50000000933

• Bladder cancer subtyping study across 4 atezo clinical trials

  Checkpoint inhibitors targeting the PD-1/PD-L1 interaction have revolutionized cancer therapy across many indications. In urothelial carcinoma (UC), several PD-1/PD-L1 inhibitors have been approved for treatment in localized and metastatic disease. While these agents provide durable clinical benefit in a number of patients, the majority do not respond. A better understanding of the molecular mechanisms underlying response and resistance is needed to tailor existing treatments and develop new therapeutics in UC. We profiled tumors from 2803 patients from four late-stage randomized clinical trials evaluating the PD-L1 inhibitor atezolizumab in metastatic UC and muscle-invasive bladder cancer by RNA-seq, a targeted DNA panel, immunohistochemistry, and digital pathology. A machine learning algorithm identified four UC transcriptional subtypes, representing luminal desert, stromal, immune, and basal tumors. Overall survival benefit from atezolizumab over standard-of-care was observed in immune and basal tumors (45% of patients). Digital pathology showed enrichment of human interpretable features in specific molecular subtypes, and a self-supervised approach could classify molecular subtypes from H&E slides with high accuracy, suggesting that machine-based image interpretation could accelerate tumor molecular profiling. Somatic alterations showed limited capacity to identify responders. Finally, mechanisms of benefit with PD-L1 blockade appeared different in immune and basal tumors, with immune tumors benefiting from B/plasma cell enrichment, while basal tumors benefited from granulocyte enrichment. This study represents the largest integration of UC molecular and curated clinical data in randomized clinical trials to date, paving the way for prospective clinical studies tailoring treatment regimens to specific tumor molecular subtypes in UC and other indications.

  Study EGAS50000000497

• The phase II Neo-Pembro trial: neoadjuvant pembrolizumab in stage IV high-grade serous ovarian cancer

  While immune checkpoint inhibitors (ICIs) have revolutionized cancer treatment, their efficacy in high-grade serous ovarian cancer (HGSOC) remains limited. Nevertheless, some patients achieve lasting responses, emphasizing the need to understand how the tumor microenvironment and molecular characteristics influence ICI response. We hypothesized higher efficacy of ICI when used in the neoadjuvant setting and preceded by induction chemotherapy. The phase 2 Neo-Pembro study (NCT03126812) included 33 untreated patients with stage IV HGSOC who were planned for 6 cycles of carboplatin-paclitaxel and interval cytoreductive surgery. Pembrolizumab (pembro) was added from cycle two and continued for one year. The primary objective was to assess intratumoral immune activation, through multiplexed immunofluorescence and immune-related gene expression. Our data provide evidence of immune activation, characterized by an increase in CD3+, CD8+, CD8+/FOXP3+ ratio, and TNF-α and interferon-γ signaling, and was most evident in major pathologic responders, which were 9/33 patients (27%, 95%CI 14-46) at time of surgery. Pathologic response was strongly associated with overall survival. At a median follow-up of 52.8 months, 8/9 major responders were alive, including 6 patients without recurrence, and 4 of 24 minor responders were alive, including one patient without recurrence. PD-L1 expression and homologous recombination deficiency were predictive of major response and may serve as potential biomarkers, warranting further exploration. In a post-hoc comparison, major responders treated with pembro had significantly better overall survival than major responders in a historic cohort treated without pembro. Our data indicate that immune activation was stronger in the subset of major pathologic responders and these patients may derive long-term survival benefits from the addition of neoadjuvant pembro.

  Study EGAS50000000781

• Analysis of CD20 loss in patients treated with Mosunetuzumab

  CD20 is an established therapeutic target in B-cell malignancies. The CD20xCD3 bispecific antibody mosunetuzumab has significant efficacy in B-cell non-Hodgkin lymphomas (NHL). Since target antigen loss is a recognized mechanism of resistance, we evaluated the contribution of CD20 expression to clinical response in patients with relapsed and/or refractory NHL in the phase 1/2 GO29781 trial (ClinicalTrials.gov: NCT02500407) investigating mosunetuzumab monotherapy. CD20 was studied using immunohistochemistry (IHC), RNA sequencing, and whole-exome sequencing performed centrally in biopsies collected prior to treatment, or as paired samples at pre-dose and on-treatment, or at-progression. Prior to treatment, most patients exhibited a high proportion of tumor cells expressing CD20; however, the proportion of CD20+ tumor cells in 14/293 (6%) patients was &lt;10%. Analyses of paired biopsies from patients during treatment revealed that CD20 levels were maintained in 29/30 (97%), vs patients with samples at-progression, in which 11/32 (34%) showed CD20 loss. Reduced transcription or acquisition of truncating mutations explained most but not all cases of CD20 loss. In vitro modeling confirmed the effect of CD20 variants identified in clinical samples on reduction of CD20 expression and confirmed missense mutations in the extracellular domain that could block mosunetuzumab binding. This study expands knowledge about the occurrence of target antigen loss after anti-CD20 therapeutics to CD20-targeting bispecific antibodies and elucidates the underlying mechanisms for reduced CD20 expression at disease progression that may be generalizable to other anti-CD20 targeting agents. These results also establish the utility of readily-available IHC staining for CD20 as a tool to inform clinical decisions.

  Study EGAS50000000151

• Repeat expansions with small TTTCA insertions in MARCHF6 cause Familial Adult Myoclonus without Epilepsy

  Familial Adult Myoclonic Epilepsy (FAME) is a rare autosomal dominant disorder characterized by cortical tremulous myoclonus and generalized tonic clonic (GTCS) or myoclonic seizures. Intrafamilial variability regarding the prevalence of seizure has been reported. Despite apparent genetic heterogeneity, FAME is caused by the same TTTTA/TTTCA pentanucleotide repeat expansion in introns of seven different genes with no common function. Expansions containing only TTTTA repeats are benign, and only expansions including a TTTCA expansion insertion, which can be located 3' to the TTTTA expansion or in the middle of the TTTTA repeats, are pathogenic. We aimed to identify the underlying genetic cause of dominant cortical myoclonus without GTCS in two unrelated families. We used a combination of repeat-primed PCR, long-range PCR and nanopore sequencing to detect and characterize expansions at known FAME loci. We report a novel configuration of the FAME3 repeat expansion in MARCHF6 segregating with cortical myoclonus without epilepsy GTCS in both families. This expansion predominantly comprises elongated TTTTA repeats, ranging from 388 to 454 repeats, with the presence of only 5-11 TTTCA repeats at its 3' terminus (i.e., 5’-(TTTTA)~388-454(TTTCA)5-11-3’). This configuration exhibits greater stability upon meiotic transmission but somatic variability in TTTCA repeat number was detected at low level in blood. We observed an inverse correlation between the number of TTTCA repeats and the age at myoclonus onset. This study shows that as little as five TTTCA juxtaposed with expanded TTTTA expansions is sufficient to cause cortical myoclonus without GTCS. This finding has implications regarding possible pathophysiological mechanisms leading to FAME.

  Study EGAS50000000570

• Enhancer profiling identifies epigenetic markers of endocrine resistance and reveals therapeutic options for metastatic castration-resistant prostate cancer patients (LuCaP cell line data)

  Androgen Receptor (AR) signaling inhibitors, including enzalutamide, are treatment options for patients with metastatic castration-resistant prostate cancer (mCRPC), but resistance inevitably develops. Using metastatic samples from a prospective phase II clinical trial, we epigenetically profiled enhancer/promoter activities with H3K27ac chromatin immunoprecipitation followed by sequencing, before and after AR-targeted therapy. We identified a distinct subset of H3K27ac-differentially marked regions that associated with treatment responsiveness. These data were successfully validated in mCRPC patient-derived xenograft models (PDX). In silico analyses revealed HDAC3 as a critical factor that can drive resistance to hormonal interventions, which we validated in vitro. Critically, we identified the pan-HDAC inhibitor, vorinostat to be effective in decreasing tumor cell proliferation both in vitro and in vivo. Moreover, we uncovered evidence for HDAC3 working together with Glucocorticoid Receptor (GR) as a potential mechanism for this therapeutic effect. These findings demonstrate rationale for new therapeutic strategies including HDAC inhibitors to improve patient outcome in advanced stages of mCRPC. Project Goal: We identified H3K27ac regions associated with resistance in patient samples. Extensive analysis of these regions suggested that HDAC3 may play a role in resistance. To test this, we analyzed ChIP-seq and RNA-seq from cell lines generated from PDX LuCaP samples (recapitulating resistant and non-resistant) to validate findings in patient samples. ChIP-seq samples were used to verify the signals identified from specific factors, such as H3K27ac, AR, FOXA1, GR and HDAC3, their association with resistance and possible mechanisms of action. RNA-seq samples were analyzed to determine the effect of AR signaling in these samples in specific genes in the H3K27ac resistance-associated regions.

  Study EGAS50000000917

• Convergent evolution towards CD38 biallelic loss is a recurrent mechanism of resistance to anti-CD38 antibodies in multiple myeloma.

  Monoclonal antibodies directed against CD38 are a therapeutic mainstay in multiple myeloma (MM). While they have contributed to improved outcomes, most patients will still experience disease relapse. Despite their ubiquity, little is known about tumor-intrinsic mechanisms of resistance to these drugs. Antigen escape, whereby tumor cells alter expression of cell surface immunotherapy targets, has been proposed as a potential avenue, and has been seen with other immunotherapies in MM. Yet, it is unknown whether MM cells can develop permanent resistance to anti-CD38 antibodies by acquiring genomic events leading to biallelic disruption of the CD38 gene locus. Here, by using whole genome and whole exome sequencing data from 701 newly diagnosed patients, 72 patients at relapse with naivety to anti-CD38 antibodies, and 52 patients collected at relapse following anti-CD38 antibodies, we report biallelic inactivation of CD38 as a recurrent mechanism of resistance in 6% of patients, often via convergent evolution where distinct subclones are advantaged towards the same selective endpoint within the same patient. Monoallelic loss of the CD38 locus was seen in another 14% post-CD38-directed therapy. Three of the cases with biallelic events on CD38 presented with missense mutations and retained CD38 expression. Functional work demonstrated that two of these mutations decreased binding affinity and antibody-dependent cellular cytotoxicity of commercial antibodies Daratumumab and Isatuximab. However, the third mutation (R140G) provided relative resistance only towards one commercial product. Given that patients are often rechallenged with CD38 antibodies following disease progression, these data support the need for next generation sequencing to guide treatment selection.

  Study EGAS50000000709

• IDH-wildtype untreated human glioblastoma samples (GB-UK cohort), published in Noorani & Haughey et al 2025.

  Extrachromosomal DNA (ecDNA) oncogene amplification is associated with treatment resistance and shorter survival in cancer. Currently, the spatial dynamics of ecDNA, and their evolutionary impact, are poorly understood. Here, we investigate the spatiotemporal evolution of ecDNA by integrating computational modeling with samples from 94 treatment-naive human IDH wild-type glioblastoma patients. Random ecDNA segregation combined with ecDNA-conferred fitness advantages induce distinct spatial ecDNA copy number patterns which depend on ecDNA oncogenic makeup. EGFR-ecDNAs often reach high copy numbers and confer strong fitness advantages. In contrast, PDGFRA-ecDNAs reach lower copy numbers and confer weaker fitness advantages. EGFR structural variants occur exclusively on ecDNA, arise from and are intermixed with wild-type EGFR-ecDNAs. Computational modeling suggests wild-type and variant EGFR-ecDNAs often accumulate before clonal expansion. Using a genetically engineered mouse model, we confirm that the accumulation of oncogenic ecDNA prior to clonal expansion under strong positive selection is observable both in vivo and in vitro in mouse neural stem cells. Intermixed wild-type and variant ecDNAs remain prevalent at diagnosis, as less fit wild-type ecDNAs hitchhike on strongly selected variant EGFR-ecDNAs. A dN/dS analysis in the Genomics England cohort further suggests positive selection of variants on EGFR-ecDNAs but not PDGFRA-ecDNAs in glioblastoma. Overall, our results suggest ecDNA oncogenic makeup determines unique evolutionary trajectories in glioblastoma. New concepts like ecDNA tumor clonality, ecDNA heteroplasmy and multi-species ecDNAs require a refined evolutionary interpretation of genomic data in a large subset of glioblastoma patients and likely other solid cancers.

  Study EGAS00001008126

• RNAseq: Acquired non-permissive bone marrow microenvironment impairs hematopoietic stem cell proliferation and maintenance and B-cell development post-HSCT

  Defects in B-cell reconstitution upon hematopoietic stem-cell (HSC) transplantation (HSCT) are a common observation, yet the mechanism remains unexplained. The bone marrow (BM) stroma, including mesenchymal stromal cells (MSCs), guides HSC maintenance and B-lymphopoiesis by secreting crucial cytokines. We report acquired, permanent, selective and complete B-cell deficiency in the context of full donor-chimerism in a patient with X-linked lymphoproliferative disease and aimed to identify the contribution of the BM-microenvironment in disrupted B-cell reconstitution post-HSCT. We studied longitudinal BM samples from the patient and his identical twin, both of whom underwent HSCT with the same donor with opposite outcomes in B-cell reconstitution. In the index patient BM, we observed progressive loss of proliferation of HSCs and a selective block at the pre-BI cell stage. In vitro modeling studies showed limited survival of patient-HSCs and a relative accumulation of pre-B cells. Patient-derived MSCs failed to support survival and proliferation of HSCs and B-cell development of healthy-HSCs which was correlated with reduced CXCL12 levels. Using bulk RNA-sequencing of MSCs and in vitro functional studies, we showed global changes in the patients’ MSCs and a progressive loss of CXCL12 expression. Indeed, survival of patients HSCs improved supplementing in vitro development culture with CXCL12, suggesting a contribution of defective CXCL12 signaling to the phenotype. In summary, our data show that an acquired defect in the BM-stromal microenvironment and exhaustion of HSCs and committed progenitors may cause a permanent non-permissive state for normal B-cell development.

  Study EGAS50000001437

• ATAC-seq: Acquired non-permissive bone marrow microenvironment impairs hematopoietic stem cell proliferation and maintenance and B-cell development post-HSCT

  Defects in B-cell reconstitution upon hematopoietic stem-cell (HSC) transplantation (HSCT) are a common observation, yet the mechanism remains unexplained. The bone marrow (BM) stroma, including mesenchymal stromal cells (MSCs), guides HSC maintenance and B-lymphopoiesis by secreting crucial cytokines. We report acquired, permanent, selective and complete B-cell deficiency in the context of full donor-chimerism in a patient with X-linked lymphoproliferative disease and aimed to identify the contribution of the BM-microenvironment in disrupted B-cell reconstitution post-HSCT. We studied longitudinal BM samples from the patient and his identical twin, both of whom underwent HSCT with the same donor with opposite outcomes in B-cell reconstitution. In the index patient BM, we observed progressive loss of proliferation of HSCs and a selective block at the pre-BI cell stage. In vitro modeling studies showed limited survival of patient-HSCs and a relative accumulation of pre-B cells. Patient-derived MSCs failed to support survival and proliferation of HSCs and B-cell development of healthy-HSCs which was correlated with reduced CXCL12 levels. Using bulk RNA-sequencing of MSCs and in vitro functional studies, we showed global changes in the patients’ MSCs and a progressive loss of CXCL12 expression. Indeed, survival of patients HSCs improved supplementing in vitro development culture with CXCL12, suggesting a contribution of defective CXCL12 signaling to the phenotype. In summary, our data show that an acquired defect in the BM-stromal microenvironment and exhaustion of HSCs and committed progenitors may cause a permanent non-permissive state for normal B-cell development.

  Study EGAS50000001438

• Single-cell Transcriptome Profiling of Treatment-naïve and Post-treatment Colorectal Cancer: Insights into Putative Mechanisms of Chemoresistance

  Drug resistance remains a major clinical challenge in the treatment of colorectal cancer (CRC) with conventional chemotherapy. Investigation of cancer cells and tumor microenvironment (TME) following chemotherapy is critical for unraveling the mechanisms of resistance. In this study, we analyzed scRNA-Seq data from 56 CRCs including treatment-naïve tumors (25 microsatellite stable (MSS) tumors and 7 tumors with microsatellite instability (MSI)); and tumors treated with standard chemotherapy with the known response status (18 responders and 6 progressors). We observed in treatment-naïve MSS CRC that right-sided tumors are more immunogenic with higher numbers of B cells and CD8+ T memory effector cells than left-sided tumors. MSI CRC were even more immunogenic then MSS CRC and were characterized by elevated Tregs in their TME. In the post-treatment CRC there was a high prevalence of dendritic cells (DC) in the TME in the response group. The DC-derived signature was associated with better survival in a large CRC cohort from the TCGA. We also observed two fibroblast subtypes in CRC, one of which was enriched in progressors and was associated with poor survival in a CRC-TCGA cohort. Progressors also showed elevated exhausted CD8+ T memory cells and monocytes suggesting a pro-inflammatory TME. In tumor cells in progressor group we identified specific expression of chemo-protective markers MTRNR2L1 and CDX1; and their co-expression with stemness-related marker CD24. In summary, scRNA-Seq provides a valuable information for the discovery of prognostic markers, and reveals distinct features of CRC based on location, metastasis status, MSI status and treatment response status.

  Study EGAS50000000830

• Development and validation of multivariate predictors of primary endocrine resistance to tamoxifen and aromatase inhibitors in luminal breast cancer reveal drug-specific differences

  Background: Endocrine therapy is highly effective in blocking the estrogen receptor pathway in HR+/HER2- early breast cancer (EBC). However, up to 40% of patients experience relapse during or after adjuvant endocrine therapy. Here, we investigate molecular mechanisms associated with primary resistance to endocrine therapy and develop predictive models. Results: TP53 mutations were prominently associated with primary resistance to both tamoxifen (TAM) and aromatase inhibitors (AI), with AI non-responders exhibiting resistance in up to 32% of cases. Additionally, we identified distinct DNA methylation patterns in TAM and AI non-responders, with TAM non-responders showing global DNA methylation loss, associated with KRAS signaling, apical junctions and epithelial-mesenchymal transition (EMT). Conversely, we observed methylation gain in AI non-responders affecting developmental transcription factors, hypoxia and estrogen signaling. TAM or AI resistance was associated with increased methylation-inferred proportions of immune cells and decreased proportions of endothelial cells. Based on these findings and patient age, we developed the Predictive Endocrine ResistanCe Index (PERCI). PERCI stratified NR and R cases in both treatment groups and cohorts with high accuracy (ROC AUC TAM discovery 93.9%, validation 83%; AI discovery 98.6%, validation 76.9%). A simplified PERCI efficiently predicted progression-free survival in the TCGA-BRCA sub-cohort (Kaplan-Meier log-rank p-value = 0.03 between low and high PERCI groups). Conclusions: We identified genomic and epigenomic features associated with primary resistance to TMA and AI. By combining information on genomic alterations, patient age, differential methylation and tumor microenvironment (TME) composition, we developed PERCI TAM and PERCI AI as novel predictors of primary resistance, with potential additional prognostic value. Applying PERCI in a clinical setting may allow patient-specific drug selection to overcome resistance.

  Study EGAS50000001102

• Biallelic variants in the non-protein coding minor spliceosome components RNU6ATAC and RNU6ATAC cause syndromic monogenic autoimmune diabetes

  Non-protein coding genes are emerging as critical contributors to the aetiology of rare diseases, providing key insights into human biology and uncovering novel disease mechanisms. We identified 7 individuals from 4 families with early-onset diabetes (diagnosed <5 years) and immune dysregulatory features caused by biallelic variants in RNU6ATAC. RNU6ATAC encodes a small nuclear RNA (snRNA) acting as a catalytic component of the minor spliceosome, a protein-RNA complex mediating splicing of ~700 genes containing U12/minor-type introns. Variant screening of the other 64 minor spliceosome genes in 276 infants with diabetes identified 12 unrelated individuals with biallelic disease-causing variants in RNU4ATAC. Biallelic pathogenic RNU4ATAC variants are known to cause a variable spectrum of clinical features, which until now did not include diabetes. Clinically, 12/19 RNU6ATAC/RNU4ATAC patients had additional immune dysregulatory features, and 50% of patients tested were islet-autoantibody positive, strongly supporting an autoimmune aetiology for their diabetes. RNA-seq in 3 individuals with biallelic RNU6ATAC variants showed a pattern of intron retention in U12 intron-containing genes similar to that seen in RNU4ATAC-individuals (n=3), supporting a shared disease mechanism. Analysis of patients’ transcriptomic, methylation and immune data revealed impaired B cell development and maturation. We conclude that biallelic RNU6ATAC variants cause a novel syndrome of early-onset autoimmune diabetes and immune dysregulation. We further show that infancy-onset diabetes is a previously unrecognised feature of RNU4ATACopathy. Our work highlights the important role of two snRNAs critical to minor spliceosome function in immune system regulation, providing novel insights into the pathogenesis of autoimmune diabetes.

  Study EGAS50000001565

• Confronting historical legacies of biological anthropology in South Africa—Restitution, redress and community-centered science: The Sutherland Nine

  We describe a process of restitution of nine unethically acquired human skeletons to their families, together with attempts at redress. Between 1925–1927 C.E., the skeletonised remains of nine San or Khoekhoe people, eight of them known-in-life, were removed from their graves on the farm Kruisrivier, near Sutherland in the Northern Cape Province of South Africa. They were donated to the Anatomy Department at the University of Cape Town. This was done without the knowledge or permission of their families. The donor was a medical student who removed the remains from the labourers’ cemetery on his family farm. Nearly 100 years later, the remains are being returned to their community, accompanied by a range of community-driven interdisciplinary historical, archaeological and analytical (osteobiographic, craniofacial, ancient DNA, stable isotope) studies to document, as far as possible, their lives and deaths. The restitution process began by contacting families living in the same area with the same surnames as the deceased. The restitution and redress process prioritises the descendant families’ memories, wishes and desire to understand the situation, and learn more about their ancestors. The descendant families have described the process as helping them to reconnect with their ancestors. A richer appreciation of their ancestors’ lives, gained in part from scientific analyses, culminating with reburial, is hoped to aid the descendant families and wider community in [re-]connecting with their heritage and culture, and contribute to restorative justice, reconciliation and healing while confronting a traumatic historical moment. While these nine individuals were exhumed as specimens, they will be reburied as people.

  Study EGAS50000000971

• Repotrectinib in NTRK fusion–positive advanced solid tumors: a phase 1/2 trial

  Early-generation TRK tyrosine kinase inhibitors (TKIs) approved for treating NTRK fusion–positive (NTRK+) solid tumors provide clinical benefit; however, resistance emerges. Repotrectinib is a next-generation ROS1/TRK TKI with a compact macrocyclic structure designed to improve durability of response. TRIDENT-1 is a registrational phase 1/2 trial assessing repotrectinib, a next-generation ROS1/TRK TKI, in adults with advanced solid tumors, including NTRK+ disease. The primary endpoint was confirmed objective response; secondary endpoints included duration of response (DOR), progression-free survival (PFS), overall survival and safety. Median follow-up ranged between 21.3 months and 25.7 months. In the TKI-naive cohort (n = 51; 95% confidence interval (CI)), the response rate was 59% (44–72); the median DOR was not estimable (NE); and the median PFS was 30.3 months (9.0–NE). In the TKI-pretreated cohort (n = 69; 95% CI), the response rate was 48% (36–60); the median DOR was 9.8 months (7.4–13.0); and the median PFS was 7.4 months (3.9–9.7). Of 30 TKI-pretreated patients with NTRK solvent front mutations, 16 had a response (53%; 95% CI: 34–72). Intracranial responses were observed in two of three TKI-naive patients and in four of six TKI-pretreated patients with measurable intracranial disease at baseline. Among all treated patients (n = 565), the most common any-grade treatment-related adverse event (TRAE) was dizziness (57%); most TRAEs were low grade; and 4% discontinued repotrectinib due to a TRAE. Overall, repotrectinib demonstrated durable systemic and intracranial responses with generally low-grade adverse events in patients with NTRK+ solid tumors, including those with previous TRK TKI treatment and solvent front mutations. These results support the use repotrectinib to treat patients with NTRK+ solid tumors. ClinicalTrials.gov identifier: NCT03093116 .

  Study EGAS50000001572

• Distinct evolutionary trajectories of primary high grade serous ovarian cancers revealed through spatial mutational profiling

  High-grade serous ovarian cancer (HGSC) is characterized by poor outcome, often attributed to emergence of treatment-resistant sub-clones. We sought to measure the degree of genomic diversity within primary, untreated HGSC to examine the natural state of tumor evolution prior to therapy. We performed exome sequencing, copy number analysis, targeted amplicon deep sequencing and gene expression profiling on thirty-one spatially and temporally separated HGSC tumor specimens (six patients) including ovarian masses, distant metastases, and fallopian tube lesions. We found widespread intra-tumoral variation in mutation, copy number, and gene expression profiles, with key driver alterations in genes present in only a subset of samples (e.g. PIK3CA, CTNNB1, NF1). On average, only 51.5% of mutations were present in every sample of a given case (range: 10.2% to 91.4%), with TP53 as the only somatic mutation consistently present in all samples. Complex segmental aneuploidies, such as whole genome doubling, were present in a subset of samples from the same individual, with divergent copy number changes segregating independently of point mutation acquisition. Reconstruction of evolutionary histories showed one patient with mixed HGSC and endometrioid histology with common etiologic origin in the fallopian tube and subsequent selection of different driver mutations in the histologically distinct samples. In this patient, we observed mixed cell populations in the early fallopian tube lesion, indicating diversity arises at early stages of tumorigenesis. Our results reveal that HGSC exhibit highly individual evolutionary trajectories and diverse genomic tapestries prior to therapy, exposing an essential biological characteristic to inform future design of personalized therapeutic solutions and investigation of drug resistance mechanisms.

  Study EGAS00001000547

• PCGP Ph-like ALL

  Background: BCR-ABL1-like, or Ph-like acute lymphoblastic leukemia is characterized by a gene expression profile similar to BCR-ABL1 positive ALL, genetic alterations of lymphoid transcription factor genes, and poor outcome. Sequencing of small numbers of Ph-like ALL cases has identified genetic alterations activating kinase signaling suggesting Ph-like ALL may be amenable to treatment with tyrosine kinase inhibitors. However, the spectrum of genetic alterations in childhood and adult Ph-like ALL is incompletely understood. Methods: We performed genomic profiling of 1736 B-ALL cases and next-generation sequencing for 160 Ph-like cases. We examined the functional effects of chimeric fusion proteins in mouse cell lines. Results: The frequency of Ph-like ALL rose from 11% in children to 26% in young adults, and was associated with very poor outcome. Kinase-activating alterations were identified in 90% of Ph-like cases, most commonly fusions involving 10 kinase or cytokine receptor genes (ABL1, ABL2, CRLF2, CSF1R, EPOR, JAK2, NTRK3, PDGFRB, PTK2B and TYK2), and mutations involving FLT3, IL7R and SH2B3. Expression of ABL1, ABL2, CSF1R and JAK2 fusions resulted in cytokine-independent proliferation of cell lines and activation of pSTAT5. Cells expressing ABL1, ABL2, CSF1R and PDGFRB fusions were sensitive to dasatinib, and JAK2 fusions to ruxolitinib. Conclusions: Ph-like ALL is characterized by a diverse range of genomic alterations that converge on a limited number of kinase signaling pathways amenable to inhibition with currently available tyrosine kinase inhibitors. Trials identifying Ph-like ALL and testing the efficacy of tyrosine kinase inhibitor therapy are warranted

  Study EGAS00001000654

• Dynamics of genomic clones in breast cancer patient xenografts at single cell resolution

  Human cancers, including breast cancers, are comprised of clones differing in mutation content that evolve dynamically in space and time following principles of Darwinian evolution, underpinning important emergent features such as drug resistance and metastasis. Human breast cancer xenoengraftment is used as a means of capturing and studying tumour biology, and breast tumour xenografts are generally assumed to be reasonable models of the originating tumours. However the consequences and reproducibility of engraftment and propagation on the genomic clonal architecture of tumours has not been systematically examined at single cell resolution. Here we show by both deep genome and single cell sequencing methods, the clonal dynamics of initial engraftment and subsequent serial propagation of primary and metastatic human breast cancers in immunodeficient mice. In all cases examined, clonal selection on engraftment was observed in both primary and metastatic breast tumours, varying in degree from extreme selective engraftment of minor (&lt;5% of starting population) clones to moderate, polyclonal engraftment. Furthermore, ongoing clonal dynamics during serial passaging is a feature of tumours experiencing modest initial selection. Through single cell sequencing, we show that major mutation clusters estimated from tumour population sequencing relate predictably to the most abundant clonal genotypes, even in clonally complex and rapidly evolving cases. Finally, we show that similar clonal expansion patterns can emerge in independent grafts of the same starting tumour population, indicating that genomic aberrations can be reproducible determinants of evolutionary trajectories. Our results show that measurement of genomically defined clonal population dynamics will be highly informative for functional studies utilizing patient-derived breast cancer xenoengraftment.

  Study EGAS00001000952

• UK10K NEURO FSZ

  In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches. These Finnish schizophrenia samples have been collected from a population cohort using national registers. The entire sample collection consists of 2756 individuals from 458 families of whom 931 are diagnosed with schizophrenia spectrum disorder, each family having at least two affected siblings. 170 families originate from an internal isolate (Kuusamo) with a three-fold lifetime risk for the trait. The genealogy of the internal isolate is well documented and the individuals form a “megapedigree” reaching to the 17th Century. All diagnoses are based on DSM-IV and for a large fraction of cases there is cognitive data. For further details/descriptions with regard to this data set please contact Tiina Paunio (tiina.paunio@thl.fi)

  Study EGAS00001000118

• Molecular_diagnosis_of_albinism

  Albinism is genetically heterogeneous rare genetic condition affecting 1:17000 in the Western world (but more frequent in Africa) whose main feature is a profound visual impairment, characterised by foveal hypoplasia, abnormal chiasmatic connections, nystagmus and photofobia. All these features result in severly altered visual acuity (&lt;0,1), absent depth perception and poor night vision. People with albinism are primarily visually handicapped. In addition, for some types of albinism, the visual phenotype can be presented with partial or total hypopigmentation, hence resulting in a secondary phenotype which can lead to skin cancer if skin is not adequately protected. Recently a new syndrome has been described, FHONDA, with the same visual abnormalities of albinism but without pigment alteration. The traditional classification differentiates Oculoculatenous albinism (OCA), where hypopigmentation involves hair, skin and eyes versus Ocular Albinism (OA), where hypopigmentation only affects the eyes. These are non-sydrimic types of albinism. Some syndromic forms (Hermansky-Pudlak=HPS, Chediak-Higashi=CHS) affect cells beyond pigment cells, present in the lungs, immune system, platelets and intestines, resulting in more severe phenotypes that can be fatal. Mutations in at least 19 genes are assocaited with the corresponding types of albinism. Most hospitals will only diagnose the most frequent cases using traditional Sanger, MLPA approaches. Some will use CGH arrays. We aim to diagnose all cases of albinism through the Albinochip proposal, which combines a Sequenom first step of known mutations combined with subsequent NGS approaches. In some cases we fail to find a second mutation, these are good candidates for further full exome analyses. 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 EGAS00001002068

• Whole-exome Sequencing Combined with Functional Genomics Reveals Novel Candidate Driver Cancer Genes in Endometrial Cancer

  Endometrial cancer is the most common gynecological malignancy, with more than 280,000 cases occurring annually worldwide. Although previous studies have identified important common somatic mutations in endometrial cancer, they have primarily focused on a small set of known cancer genes and have thus provided a limited view of the molecular basis underlying this disease. Here we developed an integrated systems biology approach to identifying novel cancer genes contributing to endometrial tumorigenesis. We first performed whole-exome sequencing on 13 endometrial cancers and matched normal samples, systematically identifying somatic alterations with high precision and sensitivity. We then combined bioinformatics prioritization with high-throughput screening (including both shRNA-mediated knockdown and expression of wild type and mutant constructs) in a highly sensitive cell viability assay. Our results revealed 12 potential driver cancer genes including 10 tumor suppressor candidates (ARID1A, INHBA, KMO, TTLL5, GRM8, IGFBP3, AKTIP, PHKA2, TRPS1 and WNT11) and two oncogene candidates (ERBB3 and RPS6KC1). The results in the “sensor” cell line were recapitulated by siRNA-mediated knockdown in endometrial cancer cell lines. Focusing on ARID1A, we integrated mutation profiles with functional proteomics in 222 endometrial cancer samples, demonstrating that ARID1A mutations frequently co-occur with mutations in the phosphatidylinositol 3-kinase (PI3K) pathway and are associated with PI3K pathway activation. siRNA knockdown in endometrial cancer cell lines increased AKT phosphorylation supporting ARID1A as a novel regulator of PI3K pathway activity. Our study presents the first unbiased view of somatic coding mutations in endometrial cancer and provides functional evidence for diverse driver genes and mutations in this disease.

  Study EGAS00001000318

• Integrative genomic profiling of hepatocellular adenomas reveals recurrent FRK activating mutations and mutational processes of malignant transformation

  Hepatocellular adenomas (HCA) are benign liver tumors predominantly developed in women using oral contraception. Here, exome sequencing identified novel recurrent somatic FRK mutations that induce constitutive kinase activity, STAT3 activation and cell proliferation sensitive to Src inhibitors. We also found new recurrent mutations activating JAK1, gp130 or ß-catenin. Chromosome copy-number and methylation profiling revealed patterns that correlated with specific gene mutations and tumor phenotypes. Finally, integrative analysis of HCA transformed to hepatocellular carcinoma revealed ß-catenin mutation as an early alteration and TERT promoter mutations as associated with the last step of the adenoma-carcinoma transition. In conclusion, we identified the genomic diversity in benign hepatocyte proliferation, new therapeutic targets and the key genomic determinants of the adenoma-carcinoma transformation sequence. Significance: Malignant transformation of benign adenoma into malignant carcinoma is frequently observed in several epithelial tumor types. Adenoma-carcinoma transition is crucial for the patient prognosis, however, little is known on the molecular mechanisms involved. Here, we performed an integrated genomic analysis of hepatocellular adenomas (HCA). Among several novel gene mutations, we identified recurrent somatic mutation activating FRK, a Src-like kinase. FRK activating mutations induce STAT3 activation and cell proliferation that can be targeted by Src inhibitors. Focusing on malignant transformation of HCA in hepatocellular carcinoma, we identified CTNNB1 and TERT promoter mutation as early and late genomic events involved in adenoma-carcinoma transition. In conclusion, we propose new therapeutic targets and an identification of adenomas with the highest risk of malignant transformation. Highlights: - FRK is activated by somatic recurrent mutations. - FRK mutants promote cell proliferation sensitive to src inhibitors. - CTNNB1 activation is an early event in HCA transformation. - TERT promoter mutations are associated with adenoma-carcinoma transition.

  Study EGAS00001000679

• Systems biology of Colorectal Cancer

  Colorectal cancer (CRC) is one of the most common cancers in both males and females, and it is perhaps the best understood of all epithelial tumors in terms of its molecular origin. Yet, despite large amount of work that has concentrated on understanding of colon tumorigenesis, we still do not know the full complement of molecular lesions that are individually necessary – and together sufficient – to cause colorectal cancer. Neither do we understand why some specific mutations that are relatively rare in other tumors (e.g. loss of the APC tumor suppressor) are extremely common in colorectal cancer. We propose here to use the tools of systems biology to develop a quantitative and comprehensive model of colorectal tumorigenesis. The model will include identification of cell-type specific and oncogenic pathways that contribute to colon tumorigenesis, and explain in molecular detail how a genotype of an individual CRC leads to activation of downstream genes that drive uncontrolled cell growth. This model will subsequently be used to find novel therapeutic targets, to guide genetic screening to identify individuals with elevated risk for developing CRC, and to classify patients into molecular subgroups to select the treatment combination that is optimal for each patient (personalized medicine). The specific objectives of the SYSCOL project are: 1. Identify genetic markers for individual risk using genotyping and sequencing of constitutional DNA from sporadic and familial CRC cases and controls 2. Identify genes and regulatory elements that contribute to colorectal cancer cell growth 3. Use data from Aims 1-2 to develop a quantitative model for colorectal tumorigenesis 4. Apply the model for identification of high-risk individuals, for molecular classification of the disease, and for identification of novel molecular treatment targets

  Study EGAS00001000854

• The Genomic Landscape of Response to EGFR Blockade in Colorectal Cancer

  Colorectal cancer (CRC) is the third most common cancer world-wide with 1.2 million patients diagnosed yearly. In late stage CRC, the most commonly used targeted therapies are monoclonal antibodies cetuximab and panitumumab, which inactivate EGFR. Recent studies have identified alterations in KRAS and other genes as likely mechanisms of primary and secondary resistance to anti-EGFR antibody therapy. Despite these efforts, additional mechanisms of resistance to EGFR blockade are thought to be present in CRC and little is known about determinants of sensitivity to this therapy. To examine the effect of somatic genetic changes in CRC on response to anti-EGFR antibody therapy, we performed complete exome sequence and copy number analyses of 129 patient-derived tumorgrafts and targeted genomic analyses of 55 patient tumors, all of which were KRAS wild-type. We analyzed the response of tumors to anti-EGFR antibody blockade in tumorgraft models or in clinical settings. In addition to previously identified genes, we detected mutations in ERBB2, EGFR, FGFR1, PDGFRA, and MAP2K1 as potential mechanisms of primary resistance to this therapy. Novel alterations in the ectodomain of EGFR were identified in patients with acquired resistance to EGFR blockade. Amplifications and sequence changes in the tyrosine kinase receptor adaptor gene IRS2 were identified in tumors with increased sensitivity to anti-EGFR therapy. Therapeutic resistance to EGFR blockade could be overcome in tumorgraft models through combinatorial therapies targeting actionable genes. These analyses provide a systematic approach to evaluate response to targeted therapies in human cancer, highlight new mechanisms of responsiveness to anti-EGFR therapies, and provide new avenues for intervention in the management of CRC.

  Study EGAS00001001305

• MiR expression profiles of paired primary colorectal cancerand metastases by next-generation sequencing

  MicroRNAs (miRs) have been recognized as promising biomarkers. It is unknown to what extent tumor-derived miRs aredifferentially expressed between primary colorectal cancers (pCRCs) and metastatic lesions, and to what extent the expressionprofiles of tumor tissue differ from the surrounding normal tissue. Next-generation sequencing (NGS) of 220 fresh-frozen samples,including paired primary and metastatic tumor tissue and non-tumorous tissue from 38 patients, revealed expression of 2245known unique mature miRs and 515 novel candidate miRs. Unsupervised clustering of miR expression profiles of pCRC tissue withpaired metastases did not separate the two entities, whereas unsupervised clustering of miR expression profiles of pCRC withnormal colorectal mucosa demonstrated complete separation of the tumor samples from their paired normal mucosa. Twohundred and twenty-two miRs differentiated both pCRC and metastases from normal tissue samples (false discovery rate (FDR)o0.05). The highest expressed tumor-specific miRs were miR-21 and miR-92a, both previously described to be involved in CRC withpotential as circulating biomarker for early detection. Only eight miRs, 0.5% of the analysed miR transcriptome, were differentiallyexpressed between pCRC and the corresponding metastases (FDR o0.1), consisting of five known miRs (miR-320b, miR-320d,miR-3117, miR-1246 and miR-663b) and three novel candidate miRs (chr 1-2552-5p, chr 8-20656-5p and chr 10-25333-3p). Theseresults indicate that previously unrecognized candidate miRs expressed in advanced CRC were identified using NGS. In addition,miR expression profiles of pCRC and metastatic lesions are highly comparable and may be of similar predictive value for prognosisor response to treatment in patients with advanced CRC

  Study EGAS00001001127

• Genesis of Two Most Prevalent Variants Causing Combined Pituitary Hormone Deficiency in 21 Populations

  Two variants (c.[301_302delAG];[301_302delAG] and c.[150delA];[150delA]) in the PROP1 gene are the most common genetic causes of recessively inherited combined pituitary hormones deficiency (CPHD). Our objective was to analyze in detail the origin of the two most prevalent variants. In the multicentric study were included 237 patients with CPHD and their 15 relatives carrying c.[301_302delAG];[301_302delAG] or c.[150delA];[150delA] or c.[301_302delAG];[ 150delA]. They originated from 21 different countries worldwide. We genotyped 21 single nucleotide variant markers flanking the 9.6Mb region around the PROP1 gene that are not in mutual linkage disequilibrium in the general populations - a finding of a common haplotype would be indicative of ancestral origin of the variant. Haplotypes were reconstructed by Phase and Haploview software, and the variant age was estimated using an allelic association method. We demonstrated ancestral origin of both variants: c.[301_302delAG] was carried on 0.2Mb-long haplotype in a majority of European patients arising approximately 101 generations ago (confidence interval 90.1-116.4). Patients from the Iberian Peninsula displayed a different haplotype which was estimated to have emerged 23.3 (20.1-29.1) generations ago. Subsequently, the data indicated that both haplotypes were transmitted to Latin American patients approximately 13.8 (12.2-17.0) and 16.4 (14.4-20.1) generations ago, respectively. The c.[150delA] variant that was carried on a haplotype spanning about 0.3Mb was estimated to appear 43.7 (38.4-52.7) generations ago. We present strong evidence that the most frequent variants in the PROP1 gene are not a consequence of variant hot-spots as previously assumed, but are founder variants.

  Study EGAS00001001165

• 'Targeted High Throughput Sequencing in Clinical Cancer Settings: Formaldehyde fixed-paraffin embedded (FFPE) tumor tissues, input amount and tumor heterogeneity'

  "Background: Massively parallel sequencing technologies have brought an enormous increase in sequencing throughput. However, these technologies need to be further improved with regard to reproducibility and applicability to clinical samples and settings. Methods: Using identification of genetic variations in prostate cancer as an example we address three crucial challenges in the field of targeted re-sequencing: Small nucleotide variation (SNV) detection in samples of formalin-fixed paraffin embedded (FFPE) tissue material, minimal amount of input sample and sampling in view of tissue heterogeneity.Results: We show that FFPE tissue material can supplement for fresh frozen tissues for the detection of SNVs and that solution-based enrichment experiments can be accomplished with small amounts of DNA with only minimal effects on enrichment uniformity and data variance. Finally, we address the question whether the heterogeneity of a tumor is reflected by different genetic alterations, e.g. different foci of a tumor display different genomic patterns. We show that the tumor heterogeneity plays an important role for the detection of copy number variations, but is of minor importance for the detection of somatic variations.Conclusions: The application of high throughput sequencing technologies in cancer genomics opens up a new dimension for the identification of disease mechanisms. In particular the ability to use small amounts of FFPE samples available from surgical tumor resections and histopathological examinations facilitates the collection of precious tissue materials. However, care needs to be taken in regard to the locations of the biopsies taken, which can have an influence on the prediction of copy number variations. Bearing these technological challenges in mind will significantly improve many large-scale sequencing studies and will - on a long-range - result in a more reliable prediction of individual cancer therapies."

  Study EGAS00001000136

• Telomerase activation by genomic rearrangements in high-risk neuroblastoma

  Neuroblastoma is a malignant pediatric tumor of the sympathetic nervous system1. Roughly half of these tumors regress spontaneously or are cured by limited therapy. By contrast, high-risk neuroblastomas have an unfavorable clinical course despite intensive multimodal treatment, and their molecular basis has remained largely elusive2-4. We have performed whole-genome sequencing of 56 neuroblastomas (high-risk, n=39; low-risk, n=17) and discovered recurrent genomic rearrangements affecting a chromosomal region (5p15.22) proximal of the telomerase reverse transcriptase gene (TERT). These rearrangements occurred only in high-risk neuroblastomas (12/39, 31%) in mutually exclusive fashion with MYCN amplifications and ATRX mutations, which are known genetic events in this tumor type1,2,5. In an extended case series (n=217), TERT rearrangements defined a subgroup of high-risk tumors with particularly poor outcome. Despite the large diversity of these rearrangements, they all induced massive transcriptional upregulation of TERT. In the remaining high-risk tumors, TERT expression was also elevated in MYCN-amplified tumors, whereas alternative lengthening of telomeres was present in neuroblastomas without TERT or MYCN alterations, suggesting that telomere lengthening represents a central mechanism defining this subtype. The 5p15.22 rearrangements juxtapose the TERT coding sequence to strong enhancer elements, resulting in massive chromatin remodeling and DNA methylation of the affected region. Supporting a functional role of TERT, neuroblastoma cells bearing rearrangements or amplified MYCN exhibited both upregulated TERT expression and enzymatic telomerase activity. In summary, our findings show that remodeling of the genomic context abrogates transcriptional silencing of TERT in high-risk neuroblastoma and places telomerase activation in the center of transformation in a large fraction of these tumors.

  Study EGAS00001001308

• An integrative model of pathway convergence in genetically heterogeneous blast crisis chronic myeloid leukemia (CML)

  Targeted therapies against the BCR-ABL1 kinase have revolutionized treatment of chronic phase (CP) chronic myeloid leukemia (CML). In contrast, management of blast crisis (BC) CML remains challenging because BC cells acquire complex molecular alterations that confer stemness features to progenitor populations and resistance to BCR-ABL1 tyrosine kinase inhibitors. Comprehensive models of BC transformation have proved elusive because of the rarity and genetic heterogeneity of BC, but are important for developing biomarkers predicting BC progression and effective therapies. To better understand BC, we performed an integrated multiomics analysis of 74 CP and BC samples using whole-genome and exome sequencing, transcriptome and methylome profiling, and chromatin immunoprecipitation followed by high-throughput sequencing. Employing pathway-based analysis, we found the BC genome was significantly enriched for mutations affecting components of the polycomb repressive complex (PRC) pathway. While transcriptomically, BC progenitors were enriched and depleted for PRC1- and PRC2-related gene sets respectively. By integrating our data sets, we determined that BC progenitors undergo PRC-driven epigenetic reprogramming toward a convergent transcriptomic state. Specifically, PRC2 directs BC DNA hypermethylation, which in turn silences key genes involved in myeloid differentiation and tumor suppressor function via so-called epigenetic switching, whereas PRC1 represses an overlapping and distinct set of genes, including novel BC tumor suppressors. On the basis of these observations, we developed an integrated model of BC that facilitated the identification of combinatorial therapies capable of reversing BC reprogramming (decitabine+PRC1 inhibitors), novel PRC-silenced tumor suppressor genes (NR4A2), and gene expression signatures predictive of disease progression and drug resistance in CP.

  Study EGAS00001001751

• Tracing the origin of disseminated tumor cells in breast cancer using single-cell sequencing

  Background Single-cell micro-metastases of solid tumors often occur in the bone marrow. These disseminated tumor cells (DTCs) may resist therapy and lay dormant or progress to cause overt bone and visceral metastases. The molecular nature of DTCs remains elusive, as well as when and from where in the tumor they originate. Here, we apply single-cell sequencing to identify and trace the origin of DTCs in breast cancer. Results We sequence the genomes of 63 single cells isolated from six non-metastatic breast cancer patients. By comparing the cells DNA copy number aberration (CNA) landscapes with those of the primary tumors and lymph node metastasis, we establish that 53% of the single cells morphologically classified as tumor cells are DTCs disseminating from the observed tumor. The remaining cells represent either non-aberrant ‘normal’ cells or aberrant cells of unknown origin that have CNA landscapes discordant from the tumor. Further analyses suggest that the prevalence of aberrant cells of unknown origin is age-dependent, and that at least a subset is hematopoietic in origin. Evolutionary reconstruction analysis of bulk tumor and DTC genomes enables ordering of CNA events in molecular pseudo-time and traced the origin of the DTCs to either the main tumor clone, primary tumor subclones, or subclones in an axillary lymph node metastasis. Conclusions Single-cell sequencing of bone marrow epithelial-like cells, in parallel with intra-tumor genetic heterogeneity profiling from bulk DNA, is a powerful approach to identify and study DTCs, yielding insight into metastatic processes. A heterogeneous population of CNA-positive cells is present in the bone marrow of non-metastatic breast cancer patients, only part of which are derived from the observed tumor lineages.

  Study EGAS00001002102

• Targeting the DNA Repair Pathway in Ewing Sarcoma

  Ewing’s Sarcoma is a bone and soft tissue tumor that primarily affects adolescents and young adults. With current therapies, 70% of patients with localized disease survive but survival for metastatic and recurrent disease is poor. Whole genome sequencing of 19 Ewing’s Sarcoma tumors showed that STAG2 was mutated in 10% (2/19) of the tumors and STAG2 protein was absent in 14% (13/106) tumors by immunohistochemical staining. Previous studies have shown that glioblastoma cells lacking STAG2 are more sensitive to poly-ADP ribose polymerase (PARP) inhibitors. We found that Ewing’s Sarcoma cell lines are sensitive to PARP inhibitors irrespective of STAG2 protein expression. Ewing Sarcoma cell lines are defective in double strand DNA break repair. PARP inhibitor cytotoxicity in Ewing’s Sarcoma cells was potentiated 10-1,000 fold by DNA damaging agents (irinotecan and temozolomide). To extend these studies in vivo, we developed an orthotopic Ewing’s Sarcoma mouse model and performed pharmacokinetic and pharmacodynamic studies with three different PARP inhibitors (BMN-673, olaparib and veliparib) in clinical development for pediatric cancer. Those data were used to design preclinical phase I studies to identify tolerable drug combinations for pilot efficacy testing (preclinical phase II). Based on the results of the preclinical phase I/II data, we performed a double blind, randomized, placebo controlled preclinical phase III trial with 274 mice in 15 treatment groups. Irinotecan administered in a low-dose protracted schedule optimized for pediatric patients was an effective DNA damaging agent to combine with olaparib and BMN-673 and was better tolerated than combinations with temozolomide. Combining olaparib or BMN-673 with irinotecan and temozolomide gave complete and durable responses in over 80% of the mice.

  Study EGAS00001000839

• The Genetic Landscape of BCL2 Break Negative Follicular Lymphoma

  Background: 15% of high and 50% of low stage FL lack the t(14;18). Nevertheless, most of these cases express BCL2 by so far unknown mechanisms. Principals/Methodology: To find molecular mechanisms, triggering FL-pathogenesis in t(14;18)-negative FL, the exonic SNV-profiles of 33 FL with and 16 FL without t(14;18) were assessed by whole exome-sequencing and external SNV-data integration. The SNV-profile of the whole exome-sequenced FL was correlated with copy number, LOH and gene-expression data. Finally, the N-glycosylation status was assessed (LympTrack assay) in a separate validation cohort. Results: SNVs in FL with and without t(14;18) affected different genes, but accumulated in similar pathways, e.g. apoptosis, epigenetic processes and cell cycle/proliferation. Mutations in genes associated with immune-response were more prominently and those associated with NFkB and N-glycosylation selectively enriched in t(14;18)-negative FL. ~50% of mutated genes showed BCL2-asociation at high confidence in FL with and without t(14;18). Moreover, among the genes mutated in t(14;18)-negative FL 555 (~1/3) were affected by copy number alterations and/or copy-neutral LOH, 125 were differently expressed between FL with and without t(14;18) (p&lt;0.01) and 62 were differentially expressed and affected by copy number alterations and/or copy-neutral LOH. N-glycosylation was strongly reduced in t(14;18)-negative FL as compared to results reported for historical FL cohorts. Conclusion: These results suggest a stronger crosstalk with the microenvironment that may compensate for the lack of N-glycosylation and a diverse portfolio of SNVs that may contribute to the upregulation of BCL2 and FL pathogenesis in t(14;18)-negative FL.

  Study EGAS00001002164

• Molecular evolution of metastatic clear cell renal cell carcinoma progressing under tyrosine kinase inhibitor therapy (HIPO-045)

  Background:Acquired resistance to targeted drugs is largely due to intra-tumor molecular heterogeneity and clonal selection. Objective:The objective of this study was to determine molecular alterations and mediating resistant mechanisms and to suggest additional treatment options.Design, Setting, and Participants:Four patients were enrolled into the MORE (Molecular Renal Cancer Evolution) trial. We investigated the mutational spectra of metastatic clear cell renal cell carcinoma (ccRCC) at baseline and compared them to those upon clonal evolution at progression on systemic therapy with tyrosine kinase inhibitors.Outcome Measurements:DNA was isolated from tumor tissues after surgery, metastasis biopsies after progression on sunitinib or axitinib, and subjected to whole-exome sequencing (WES). In addition, the exomes of circulating tumor DNA (ctDNA) from two patients were sequenced at baseline and at progression.Results and Limitations:Our data provide evidence for clonal evolution and various acquired resistance mechanisms of ccRCC with subclonal mutations in FLT4, MTOR, ITGA3/5, SETD2, and VHL on metastatic progression. The limited number of shared mutations in tissue and plasma biopsies indicates that both sources provide complementary information, possibly due to mutations present in metastatic lesions which were not amenable to tissue sampling but are represented in the ctDNA. Conclusion:Acquired resistance to targeted therapies in metastatic ccRCC is due to intra-tumor heterogeneity and clonal evolution. Larger patient cohorts have to be studied to estimate the benefit of WES from biopsies and ctDNA for therapy decision in metastatic ccRCC.Patient Summary:We analyzed the molecular profiles of ccRCC patients in primary tissues and metastatic tissues upon progression on tyrosine kinase inhibitor therapy. Clinically relevant mutations upon progression were detected in the metastatic sites of all of four patients, suggesting that it might be possible to derive therapies based on sequencing the metastases.?HIPO-045

  Study EGAS00001001861

• Whole genome and transcriptome analysis of a sporadic and recurring parathyroid carcinoma

  Parathyroid carcinoma is an extremely rare endocrine malignancy with an estimated incidence of less than 1 per million population. Excessive secretion of parathyroid hormone, extremely high serum calcium level and the deleterious effects of hypercalcemia are the clinical manifestations of the disease. Up to 60% of patients develop multiple disease recurrences and although long-term survival is possible with palliative surgery, permanent remission is rarely achieved. Molecular drivers of sporadic parathyroid carcinoma have remained largely unknown. Previous studies, mostly based on familial cases of the disease, suggested a potential role for the tumor suppressor MEN1 and proto-oncogene RET in benign parathyroid tumorigenesis while the tumor suppressor HRPT2 and proto-oncogene CCND1 may also act as drivers in parathyroid cancer. Here, we report the complete genomic analysis of a sporadic and recurring parathyroid carcinoma. Mutational landscapes of the primary and recurrent tumor specimens were analyzed using high throughput sequencing technologies. Such molecular profiling allowed for identification of somatic mutations never previously identified in this malignancy. These included single nucleotide point mutations in well characterized cancer genes such as mTOR, MLL2, CDKN2C and PIK3CA. Comparison of acquired mutations in patient-matched primary and recurrent tumors revealed loss of PIK3CA activating mutation during the evolution of the tumor from the primary to the recurrence. Structural variations leading to gene fusions and regions of copy loss and gain were identified at a single-base resolution. Loss of short arm of chromosome 1 along with somatic mis-sense and truncating mutations in CDKN2C and THRAP3, respectively, provide new evidence for the potential role of these genes as tumor suppressors in parathyroid cancer. The key somatic mutations indentified in this study can serve as novel diagnostic markers as well as therapeutic targets.

  Study EGAS00001000484

• Intra-tumor heterogeneity and clonal evolution patterns towards platinum-resistant high-grade serous ovarian cancer

  Although many high-grade serous ovarian carcinoma (HGSOC) patients are initially sensitive to the standard-of-care platinum-based therapy, intra-tumor heterogeneity may fuel selection towards platinum-resistant clones over subsequent lines of therapy. We performed whole-exome coupled to ultra-deep re-sequencing and copy number alteration (CNA) profiling on 31 matched diagnosis and relapse tumors from HGSOC patients treated with platinum. Genetic heterogeneity within tumor pairs varied extremely, with some recurrent tumors being identical and others completely different from the primary tumor. All primary tumors showed evidence of homologous recombination (HR)-deficiency, but none of them carried reactivating mutations in the HR pathway at relapse, also not when becoming resistant to subsequent lines of platinum. We did also not identify other mutations, CNAs or pathways that were enriched in platinum-resistant clones. However, a platinum score consisting of 13 CNAs correlated significantly with platinum-free interval (PFI) after first-line therapy in multiple cohorts of &gt;450 primary HGSOCs. A relative increase of this score in the relapsed tumor correlated with a change in PFI during subsequent therapy. Clonal evolution under platinum therapy contributes to genetic heterogeneity within tumor pairs, implicating that genetic tests linked to targeted treatment options for recurrent HGSOC need to be performed on biopsies collected at relapse. Platinum resistance of an initially platinum-sensitive HR-deficient HGSOC at relapse is not caused by reactivation of the HR pathway, suggesting that these tumors are still responsive to PARP inhibitors, instead, several chromosomal aberrations contribute to a score predictive of disease outcome under platinum therapy, both in first and subsequent lines of platinum therapy.

  Study EGAS00001001244

• Genome wide association study on coronary heart disease in patients with familial hypercholesterolemia

  Mutations in the low-density lipoprotein receptor (LDLR) gene cause familial hypercholesterolemia (FH), a disorder characterized by coronary heart disease (CHD) at young age. We aimed to apply an extreme sampling method to enhance the statistical power to identify novel genetic risk variants for CHD in individuals with FH. We selected cases and controls with an extreme contrast in CHD risk from 17 000 FH patients from the Netherlands, whose functional LDLR mutation was unequivocally established. The genome-wide association (GWA) study was performed on 249 very young FH cases with CHD and 217 old FH controls without CHD (above 65 years for males and 70 years of age for females) using the Illumina HumanHap 550K chip. In the next stage, two independent samples (one from the Netherlands and one from Italy, Norway, Spain, and the United Kingdom) of FH patients were used as replication samples. In the initial GWA analysis we identified 29 independent single nucleotide polymorphisms (SNPs) with suggestive associations with premature CHD (p&lt;1×10-4). We examined the association of these SNPs with CHD risk in the replication samples. After Bonferroni correction, none of the SNPs either replicated or reached genome-wide significance after combining the discovery and replication samples. Therefore, we conclude that the genetics of CHD risk in FH is complex and even applying an ‘extreme genetics’ approach we did not identify new genetic risk variants. Most likely, this method is not as effective in leveraging effect size as anticipated, and may therefore not lead to significant gains in statistical power. (NB submitten and conditionally accepted at European Journal of Human Genetics-please use abstract only for internal review until submission complete)

  Study EGAS00001000734

• UK10K NEURO EDINBURGH

  In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches. This sample set consists of subjects with schizophrenia recruited from psychiatric in-patient and out-patient facilities in Scotland. All diagnoses are based on standard research procedures and family histories are available. Patients have IQ&gt;70 and the cohort includes the following groups: 100 cases with detailed clinical, cognitive and structural and functional neuroimaging phenotypes; 138 familial cases who are the probands of families where DNA has been collected from other affected members; 162 unrelated individuals. In most cases patients and their families may be re contacted to take part in further studies.For further information on this cohort please contact Andrew McIntosh (andrew.mcintosh@ed.ac.uk).

  Study EGAS00001000117

• UK10K NEURO IMGSAC

  In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches.The IMGSAC cohort is an international collection of families containing children ascertained for ASDs (autism spectrum disorders). The affected individuals are have been phenotyped, including using the ADI-R and ADOS instruments. Individuals with a past or current medical disorder of probable etiological significance or TSC have been excluded. Where possible, karyotyping has been performed on one affected individual per family to exclude Fragile X syndrome. Many of the samples have been genotyped, using the Affymetrix 10k and Illumina 1M platforms. All samples to be included in the current study are of UK origin.For further information on this cohort please contact Jeremy Parr (jeremy.parr@newcastle.ac.uk).

  Study EGAS00001000120

• Pharmacological and genomic profiling identifies NFκB-targeted treatment strategies for mantle cell lymphoma

  Mantle cell lymphoma (MCL) is an aggressive malignancy characterized by poor prognosis, underscoring the need for novel therapeutic strategies1. In this study, we employed large-scale pharmacological profiling of 16 pathway inhibitors across more than 100 hematological cell line models to identify novel therapeutic approaches. This screen identified a subset of MCL cell lines that were highly sensitive to the B-cell receptor (BCR) signaling inhibitors ibrutinib and sotrastaurin (STN), which are currently in clinical trials2, 3. Sensitive MCL models exhibited chronic activation of the CARD11-BCL10-MALT1 (CBM) complex and dependency on downstream classical NFκB signaling. In contrast, insensitive cell lines displayed activation of the alternative NFκB pathway. Massively parallel transcriptome sequencing identified mutations in negative regulators of the alternative NFκB pathway in ibrutinib-insensitive cell lines. Notably, sequencing of 165 MCL patient samples revealed that components of the alternative NFκB pathway, particularly TRAF2 and BIRC3, are mutated in roughly 15% of MCL patients. Functional analyses of patient-derived BIRC3 mutant alleles indicated that these are loss-of-function mutants that constitutively activate alternative NFκB signaling. While being associated with insensitivity to ibrutinib, genetic lesions in the alternative pathway conferred dependence on NIK, a kinase that mediates signaling downstream of the TRAF2/3-BIRC2/3 complex. Thus, NIK represents a novel therapeutic target for MCL treatment, particularly for tumors refractory to BCR pathway inhibitors. Collectively, these findings reveal a pattern of mutually exclusive activation of the CBM-NFκB or alternative NIK-NFκB pathways in MCL and provide critical insights into patient stratification strategies for NFκB-pathway targeted agents.

  Study EGAS00001000622

• Transcriptional_profiling_of_tauopathies_in_human_IPS_derived_neurons

  This project is a pilot study, in collaboration with Maria Grazia Spillantini and Mariangela Iovino (Cambridge Centre for Brain Repair), to investigate the utility of IPS-derived neurons for the study of neurodegenerative disorders. Our aim is to characterise the transcriptional consequences of tauopathies using neurons derived from differentiated IPSCs as a model system. We will use IPS cells derived from six individuals, four with known mutations in the tau protein, 2 without. RNA will be extracted at Day 0 and Day 65 of differentiation by which time the neuronal tauopathy is apparent. RNA will be extracted and the transcriptome of each line characterised using RNAseq. We will then search for genes that are differentially expressed between the transcriptomes of individuals with tau mutations versus those in controls. My lab will analyse the RNAseq data, comparing both affected and controls and both time-points, to establish candidate genes. Darren Logan’s lab, along with our collaborators, will experimentally verify and further investigate these genes in additional lines and animal models. From this analysis we will generate a list of candidate genes that are differentially expressed between cases and controls. This study will not only help us understand the molecular basis of tauopathies, but also identify gene candidates for biomarkers of neurodegenerative disease. It will serve as a proof of principle for future planned studies into generating and transcriptomically analysing an allelic series of Tau mutations in IPSCs with a controlled genetic background.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 EGAS00001000382

• Plasma DNA aberrations in systemic lupus erythematosus revealed by genomic and methylomic sequencing

  We performed a high-resolution analysis of the biological characteristics of plasma DNA in systemic lupus erythematosus (SLE) patients using massively parallel genomic and methylomic sequencing. A number of plasma DNA abnormalities were found. First, aberrations in measured genomic representations (MGRs) were identified in the plasma DNA of SLE patients. The extent of the aberrations in MGRs correlated with anti-double–stranded DNA (anti-dsDNA) antibody level. Second, the plasma DNA of active SLE patients exhibited skewed molecular size-distribution profiles with a significantly increased proportion of short DNA fragments. The extent of plasma DNA shortening in SLE patients correlated with the SLE disease activity index (SLEDAI) and anti-dsDNA antibody level. Third, the plasma DNA of active SLE patients showed decreased methylation densities. The extent of hypomethylation correlated with SLEDAI and anti-dsDNA antibody level. To explore the impact of anti-dsDNA antibody on plasma DNA in SLE, a column-based protein G capture approach was used to fractionate the IgG-bound and non–IgG-bound DNA in plasma. Compared with healthy individuals, SLE patients had higher concentrations of IgG-bound DNA in plasma. More IgG binding occurs at genomic locations showing increased MGRs. Furthermore, the IgG-bound plasma DNA was shorter in size and more hypomethylated than the non–IgG-bound plasma DNA. These observations have enhanced our understanding of the spectrum of plasma DNA aberrations in SLE and may provide new molecular markers for SLE. Our results also suggest that caution should be exercised when interpreting plasma DNA-based noninvasive prenatal testing and cancer testing conducted for SLE patients.

  Study EGAS00001000962

• Ischemic stroke in a Swedish case-control study.

  Background. Genetic risk scores (GRS), summing up the total effect of several single nucleotide polymorphisms (SNPs) in genes associated with either coronary risk or cardiovascular risk factors, have been tested for association with ischemic stroke with conflicting results. Recently an association was found between a GRS, based on 29 SNPs discovered by genome-wide association studies (GWAS) and hypertension. The aim of our study was to investigate the possible association of the same GRS with ischemic stroke on top of other “traditional risk factors”, also testing its potential improvement in indices of discrimination and reclassification, in a Swedish case-control study. Methods and results. Twenty-nine SNPs were genotyped in 3,677 stroke cases and 2,415 controls included in the Lund Stroke Register (LSR), the Malmö Diet and Cancer (MDC) study and the Sahlgrenska Academy Study on Ischemic Stroke (SAHLSIS). The analysis was conducted in the entire material, and separately for the three samples. After adjustment for hypertension, diabetes mellitus and smoking habits, the GRS was associated with ischemic stroke in the entire material [OR (95% CI) 1.086 (1.029-1.147) per SD increase in the GRS p=0.003] with similar trends in all three samples; LSR [1.050 (0.967-1.140)]; p=0.25], MDC [1.168 (1.060-1.288); p=0.002] and SAHLSIS [1.124 (0.997-1.267); P=0.055]. Measures of discrimination and reclassification improved marginally using the GRS. Conclusions. A blood pressure GRS is independently associated with ischemic stroke risk in three Swedish case-control studies, however, the effect size is low and adds marginally to prediction of stroke on top of traditional risk factors including hypertension.

  Study EGAS00001000936

• Prenatal lead exposure is associated with decreased cord blood DNA methylation of the glycoprotein VI gene involved in platelet activation and thrombus formation

  Early-life lead exposure impairs neurodevelopment and later exposure affects the cardiovascular system. Lead has been associated with reduced global 5-methylcytosine DNA methylation, suggesting that lead toxicity acts through epigenetic mechanisms. The objective of this study is to clarify how early life lead exposure alters DNA methylation of specific genes, using an epigenomic approach. We measured lead concentrations in urine (gestational week, GW, 8) and erythrocytes (GW 14), using ICP-MS, for 127 pregnant mothers recruited in the MINIMat food and supplementation cohort in rural Bangladesh. Cord blood DNA methylation was analyzed with the Infinium HumanMethylation450K BeadChip and top sites were validated by methylation-sensitive high-resolution melt curve analysis. Maternal urinary lead concentrations (divided into quartiles) showed significant (after adjustment for FDR) inverse associations with methylation at nine CpGs. Three of these sites were in the 5’ end, including the promoter, of glycoprotein IV (GP6); cg18355337 (q=0.029, β=-0.30), cg25818583 (q=0.041, β=-0.18) and cg23796967 (q=0.047, β=-0.17). The methylation in another CpG site in GP6 was close to significant (cg05374025, q=0.057, β=-0.23). The erythrocyte lead concentrations (divided into quartiles) were also inversely associated with CpG methylation in GP6, although this was not statistically significant after FDR-adjustments. Eight CpG sites in GP6 constituted a differentially methylated region in relation to urinary lead (p=0.005, q=0.48) and erythrocyte lead (p=0.007, q=0.46). In conclusion, we found that moderate prenatal lead exposure appear to epigenetically affect GP6, a key component of platelet aggregation and thrombus formation, suggesting a novel link between early lead exposure and cardiovascular disease later in life.

  Study EGAS00001001575

• UK10K NEURO UKSCZ

  In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches. These samples have been collected from throughout the UK and Ireland. The samples fall into two main categories, approximately 500 have a full diagnostic work up. A proportion of these are cases with a positive family history of schizophrenia, either collected as sib-pairs or from multiplex kindred's. The second group consist mainly of &gt;300 samples that have been systematically collected within South Wales and in addition to full diagnostic work up have undergone detailed cognitive testing. All samples have obtained a DSM IV diagnosis of schizophrenia or schizoaffective disorder.For further information on this cohort please contact Michael Owen (OwenMJ@cardiff.ac.uk).

  Study EGAS00001000123

• Induced Pluripotent Cells Derived from Differentiated Rod Photoreceptors Undergo Efficient Retinogenesis in Three-Dimensional Cultures

  "Embryonic stem cells (ESCs) and induced pluripotent stem cells (iPSCs) are used to study lineage-specific differentiation in culture. We developed a standardized quantitative protocol called STEM-RET to compare eye field specification, optic cup formation, and retinal differentiation across stem cell populations. We reprogrammed individual rod photoreceptors into iPSCs by inducible expression of Oct3/4, Klf4, Sox2, and Myc. We also developed ESC lines from rods by somatic cell nuclear transfer. Using the STEM-RET protocol, we compared the abilities of ESCs, fibroblast-derived iPSCs (f-iPSCs), rod photoreceptor-derived iPSCs (r-iPSCs), and rod photoreceptorderived ESCs to form retinae using STEM-RET. r-iPSCs were the most efficient at producing differentiated retina. Retinae derived from f-iPSCs lacked an inner nuclear layer (INL) and had a concomitant reduction in amacrine cells and other INL cell populations. The INL-specific LIM homeobox gene, Lhx9, was hypermethylated in fiPSCs and concomitantly downregulated in retinae derived from those stem cells. ChIP-seq analysis for H3K36me3, H3K4me1, H3K4me3, H3K9/14Ac, H3K27me3, H3K9me3, CTCF and 5 hydroxymethyl cytosine showed that the major difference between the r-iPSCs and the f-iPSCs were in the genes marked by the CTCF insulator protein. Importantly, CTCF binding enrichment at rod photoreceptor specific genes was at least 10-fold higher in r-iPSCs relative to f-iPSCs. Taken together, our data suggest that both inhibitory and permissive epigenetic marks important for retinal development can be retained in iPSCs and such marks can be exploited to select the most favorable stem cell population to study retinal development or produce photoreceptor precursors for cell transplantation."

  Study EGAS00001001288

• Targeting FGFR1 for treatment of soft-tissue sarcoma (H021)

  Purpose: Altered FGFR1 signaling has emerged as therapeutic target in various epithelial malignancies. In contrast, the role of FGFR1 in soft-tissue sarcoma (STS) has not been established. Prompted by the detection and subsequent therapeutic inhibition of amplified FGFR1 in a patient with metastatic leiomyosarcoma, we investigated the oncogenic properties and potential as drug target of FGFR1 in STS.Experimental Design: The frequency of FGFR1 amplification and overexpression, as assessed by fluorescence in situ hybridization, microarray-based comparative genomic hybridization and mRNA expression profiling, SNP array profiling, and RNA sequencing, was determined in three independent patient cohorts. The sensitivity of STS cell lines with or without FGFR1 alterations to genetic and pharmacologic FGFR1 inhibition and the signaling pathways engaged by FGFR1 were investigated using short-term viability assays, long-term colony formation assays, and biochemical analysis.Results: Increased FGFR1 copy number was detected in 74 of 190 (38.9%; Cohort 1), 13 of 79 (16.5%; Cohort 2), and 80 of 254 (31.5%; Cohort 3) patients. FGFR1 overexpression occurred in 16 of 79 (20.2%, Cohort 2) and 39 of 254 (15.4%; Cohort 3) patients. Targeting of FGFR1 by short hairpin RNA-mediated knockdown and different small-molecule inhibitors (PD173074, AZD4547, BGJ398) revealed that the requirement for FGFR1 signaling in STS cells is primarily dictated by FGFR1 expression levels, and identified the MAPK-ERK1/2 axis as critical FGFR1 effector pathway.Conclusion: These data identify FGFR1 as driver gene in multiple STS subtypes and support FGFR1 inhibition, guided by patient selection according to FGFR1 expression and monitoring of MAPK-ERK1/2 signaling, as therapeutic option in this challenging group of diseases.H021

  Study EGAS00001001844

• Whole-exome sequencing of human pancreatic cancers and characterization of genomic instability caused by MLH1 haploinsufficiency and complete deficiency

  Whole-exome sequencing (Exome-Seq) has been successfully applied in several recent studies. We here performed Exome-Seq of 15 pancreatic tumor-normal pairs. We captured 162,073 exons of 16,954 genes and sequenced the targeted regions to a mean coverage of 56 folds. This study identified a total of 1,517 somatic mutations and validated 934 of them by transcriptome sequencing. We detected recurrent mutations in 56 genes. Among them, 41 have not been described. The mutation rates varied widely among tumors. The diversity of mutation rates was significantly correlated with the distinct MLH1 copy-number status. Exome-Seq revealed intensive genomic instability in a tumor with MLH1 homozygous deletion, indicated by a dramatically elevated mutation rate of somatic substitutions, small insertions/deletions (indels), as well as indels presented at microsatellites. Notably, we detected MLH1 allelic loss in 4 tumors. MLH1 expression was decreased by nearly half in these cases. While they were negative in the conventional microsatellite instability assay, they showed a 10.5-fold increase in the mutation rate of somatic indels, e.g. truncating indels in TP53 and TGFBR2, indicating MLH1 haploinsufficiency in correction of DNA indel errors. We further analyzed the exomes of 15 renal cell carcinomas (RCCs) and confirmed MLH1 haploinsufficiency. We observed a significant increase in the mutation rate of somatic indels in affected cases and identified recurrent truncating indels in several cancer genes such as VHL, PBRM1 and JARID1C. Together, our data suggest that MLH1 hemizygous deletion, through increasing the mutation rate of somatic indels, could drive the development and progression of sporadic cancers.

  Study EGAS00001000149

• Malignant pheochromocytomas/paragangliomas harbor mutations in transport and cell adhesion genes

  One out of ten patients with pheochromocytoma (PCC) and paraganglioma (PGL) develop malignant disease. Today there are no reliable pathological methods to predict malignancy at the time of diagnosis. Tumors harboring mutations in the succinate dehydrogenase subunit B (SDHB) gene often metastasize but the sequential genetic events resulting in malignant progression are not fully understood. The aim of this study was to identify somatic mutations that contribute to the malignant transformation of PCC/PGL. We performed pair-wise (tumor-normal) whole-exome sequencing to analyze the somatic mutational landscape in five malignant and four benign primary PCC/sympathetic PGL (sPGL), including two biological replicates from each specimen. In total, 225 unique somatic mutations were identified in 215 genes, with an average mutation rate of 0.54 mutations/megabase. Malignant tumors had a significantly higher number of mutations compared to benign tumors (p&lt;0.001). Three novel genes were identified as recurrently mutated; MYCN, MYO5B and VCL, and mutations in these genes were exclusively found in malignant sPGL tumors. Mutations in the MYO5B gene could be verified in two publicly available data sets. A gene ontology analysis of mutated genes showed enrichment of cellular functions related to cytoskeletal protein binding, myosin complex, and motor activity, many of which had functions in Rab and Rac/Rho GTPase pathways. Conclusively, we have identified recurrent mutations in genes related to intracellular transport and cell adhesion, and we have confirmed recurrent mutations of MYO5B in PCC/PGL cases with malignant potential. Our study suggests that deregulated Rab and Rac/Rho pathways may be important in PCC/PGL tumorigenesis.

  Study EGAS00001001601

• Characterization of a human iPSC-derived endocrine pancreas model

  Directed differentiation of stem cells offers a scalable solution to the need for human cell models recapitulating islet biology and T2D pathogenesis. We profiled mRNA expression at six stages of an induced pluripotent stem cell (iPSC) model of endocrine pancreas development from two donors, and characterized the distinct transcriptomic profiles associated with each stage. Established regulators of endodermal lineage commitment, such as SOX17 (log2 fold change [FC] compared to iPSCs=14.2, p-value=4.9x10-5) and the pancreatic agenesis gene GATA6 (log2 FC=12.1, p-value=8.6x10-5), showed transcriptional variation consistent with their known developmental roles. However, these analyses highlighted many other genes with stage-specific expression patterns, some of which may be novel drivers or markers of islet development. For example, the leptin receptor gene, LEPR, was most highly expressed in published data from in vivo-matured cells compared to our endocrine pancreas-like cells (log2 FC=5.5, p-value=2.0x10-12), suggesting a role for the leptin pathway in the maturation process. Endocrine pancreas-like cells showed significant stage-selective expression of adult islet genes, including INS, ABCC8, and GLP1R, and enrichment of relevant GO-terms (e.g. "insulin secretion"; odds ratio=4.2, p-value=1.9x10-3): however, principal component analysis indicated that in vitro-differentiated cells were more immature than adult islets. Integration of the stage-specific expression information with genetic data from T2D genome-wide association studies revealed that 46 of 82 T2D-associated loci harbor genes present in at least one developmental stage, facilitating refinement of potential effector transcripts. Together, these data show that expression profiling in an iPSC islet development model can further understanding of islet biology and T2D pathogenesis.

  Study EGAS00001001803

• Spiradenocarcinoma

  • Spiradenocarcinoma is a rare cutaneous sweat gland adnexal cancer with potential for aggressive behaviour. They are classified histologically into low- and high-grade tumours, with morphologically low-grade tumours thought to behave more favourably. However, limited information is available, with only 18 published cases. • We have collected morphologically low-grade spiroadenocarcinomas (one with a lung metastasis) and high-grade spiroadenocarcinomas, as well as some spiradenomas (benign lesions), cylindromas (another type of malignant cutaneous sweat gland adnexal tumour) and hybrid spiradenoma-cylindromas. • H&amp;E-stained sections were reviewed, follow-up was obtained, and immunohistochemistry for Ki-67, p53 and, MYB has been performed. The tumours were solitary, measuring 0.8-7?cm (median: 2.7?cm), with a predilection for the head and neck of elderly patients (median age: 72 years; range 53-92) without gender bias. Histologically, the tumours were multinodular and located in deep dermis and subcutis. A pre-existing spiradenoma was present in all cases. The malignant component was characterized by expansile growth with loss of the dual cell population, up to moderate cytological atypia and increased mitotic activity (median: 10/10 HPF; range 1-28). Additional findings included squamoid differentiation (n=9), necrosis (n=7), and ulceration (n=5). P53 expression was variable and no significant differences were noted in the benign compared with the malignant parts of the tumours. In contrast, in the malignant components the Ki-67 proliferative index was slightly increased, and MYB expression was lost. Follow-up (median: 67 months; range: 13-132) available for 16 patients (84%) revealed a local recurrence rate of 19% but no metastases or disease-related mortality. • Here we wish to exome sequence these cases to define the first genomic landscape for this malignancy.

  Study EGAS00001001799

• panALL exome sequencing, data set2, 700 samples

  panALL exome sequencing, data set2, 700 samples

  Dataset EGAD00001008268

• 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

• 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. &lt;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

• 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 &gt;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

• 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

• The University of Hong Kong Gastric Cancer Organoids RHO Signaling Study

  Objective: Cell-cell (CC) and cell-matrix (CM) adhesions are essential for epithelial cell survival, yet dissociation-induced apoptosis is frequently circumvented in malignant cells. Design: We explored CC and CM dependence in 58 gastric cancer (GC) organoids by withdrawing either ROCK inhibitor, matrix, or both, to evaluate their tumorigenic potential in terms of apoptosis resistance, correlation with oncogenic driver mutations and clinical behaviour. We performed mechanistic studies to determine the role of diffuse-type GC drivers: ARHGAP fusions, RHOA, and CDH1, in modulating CC (CCi) or CM (CMi) adhesion independence. Results: 97% of the tumour organoids were CMi, 66% were CCi and 52% were resistant to double withdrawal (CCi/CMi), while normal organoids were neither CMi nor CCi. Clinically, the CCi/CMi phenotype was associated with an infiltrative tumour edge and advanced tumour stage. Moreover, the CCi/CMi transcriptome signature was associated with poor patient survival when applied to 3 public GC datasets. CCi/CMi and CCi phenotypes were enriched in diffuse-type GC organoids, especially in those with oncogenic driver perturbation of RHO signalling via RHOA mutation or ARHGAP fusions. Inducible knockout of ARHGAP fusions in CCi/CMi tumour organoids led to re-sensitization to CC/CM dissociation-induced apoptosis, upregulation of focal adhesion and tight junction genes, partial reversion to a more normal cystic phenotype, and inhibited xenograft formation. Normal gastric organoids engineered with CDH1 or RHOA mutations became CMi or CCi, respectively. Conclusions: The CCi/CMi phenotype has a critical role in malignant transformation and tumour progression, offering new mechanistic information on RHO-ROCK pathway inhibition that contributes to GC pathogenicity.

  Study EGAS00001006252

• 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

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

  B-cell prolymphocytic leukemia (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&lt;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 EGAS00001003274

• 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

• 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

• 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

• Transcriptomics-driven classification of ALAL

  Acute leukemia of ambiguous lineage (ALAL) is a rare, poor-prognosis subtype of acute leukemia that cannot be assigned to a single hematopoietic lineage. Although ALAL patients are typically treated with acute myeloid leukemia (AML) or acute lymphoblastic leukemia (ALL) regimens, optimal treatment choice is hindered by their lineage ambiguity. Therefore, we investigated the added value of transcriptomics for improving lineage assignment, currently based mainly on surface markers. First, we used an in-house pipeline to comprehensively detect genetic lesions in RNA-sequencing data (n=23) with a sensitivity &gt;90%, as confirmed by targeted DNA sequencing. Second, we compared ALAL gene expression profiles (GEPs) with representative AML (n=145), B-ALL (n=223) and T-ALL (n=85) cases. In a principal component analysis, ALALs did not form a separate group, but clustered with AML, B-ALL or T-ALL. Accordingly, a machine learning classifier trained with GEPs of acute leukemias segregated 20/23 ALALs into myeloid-, B- or T-lymphoid leukemia. These 20 cases harbored genetic abnormalities consistent with the classifier-assigned leukemia. Furthermore, ALAL GEPs were deconvoluted with single-cell transcriptional profiles of normal hematopoietic cells using CIBERSORTx, revealing enrichment for signatures of lineages corresponding to the leukemic type predicted by our algorithm. The classifier was validated on an external ALAL cohort (n=24), assigning 75% of the patients to a lineage matching their immunophenotypic and methylation profiles. In conclusion, integrative analysis of RNA-sequencing data can accurately classify most ALAL cases. The pipeline and classifier developed in our study are valuable tools to improve ALAL diagnosis and guide therapeutic decisions.

  Study EGAS00001007967

• 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

• 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 (&lt;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

• 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

• Genomic and epigenomic characterization of juvenile myelomonocytic leukemia (JMML)

  Juvenile myelomonocytic leukemia (JMML) is an aggressive myeloproliferative disorder of early childhood. While some cases show spontaneous remission, allogeneic hematopoietic stem cell transplantation (HSCT) remains the only curative treatment option for the majority of patients, however, the 5-year event-free survival reaches only about 50%. Hyperactive RAS signaling is assumed to be the main driving event in JMML. It is caused by genetic alterations in CBL, KRAS, NF1, NRAS, or PTPN11 in about 90% of patients. So far, there is no clear understanding of how RAS pathway mutations relate to the heterogeneous disease biology and variable clinical outcome seen in JMML patients. As a consequence, established clinical and genetic markers fail to fully represent the observed disease heterogeneity. We hypothesized that DNA methylation profiling, either alone or in combination with genetic alterations, might provide a molecular basis for disease classification. Genome wide DNA methylation analysis using the HumanMethylation450 Bead Chip array was performed in a discovery cohort of 20 patients (19 JMML, 1 Noonan syndrome). The clinical relevance of our findings was assessed in an unselected sample set consisting of 147 consecutive patients with JMML (n=129) or Noonan syndrome (n=18) associated myeloproliferative disorder (n=18) registered in the EWOG-MDS 1998 &amp; 2006 trials. Data integration was performed in a subset of patients with available exome sequencing (n=50) and expression profiling (n=15) data. This integrated approach identified three JMML subgroups characterized by distinct clinical and biological features. This analysis further suggested a molecular mechanism by which additional genetic events, presumably further activating the RAS-RAF-MEK-ERK pathway, mediate DNA hypermethylation via up-regulation of DNMTs in more aggressive JMML cases.

  Study EGAS00001002511

• Multi-site tumor sampling highlights molecular intra-tumor heterogeneity in malignant pleural mesothelioma

  Background:Malignant pleural mesothelioma (MPM) is a heterogeneous cancer. Better knowledge of molecular and cellular intra-tumor heterogeneity throughout the thoracic cavity is required to develop efficient therapies. This study focuses on molecular intra-tumor heterogeneity using the largest series to date in MPM and is the first to report on the multi-omics profiling of a substantial series of multi-site tumor samples. Methods:Intra-tumor heterogeneity was investigated in 16 patients from whom biopsies were taken at distinct anatomical sites. The paired biopsies collected from apex, side wall, costo-diaphragmatic or highest metabolic sites as well as 5 derived cell lines were screened using targeted sequencing. Whole exome sequencing, RNA sequencing, and DNA methylation were performed on a subset of the cohort for deep characterization. Molecular classification, recently defined histo-molecular gradients, and cell populations of the tumor microenvironment were assessed. Results:Sequencing analysis identified heterogeneous variants notably in NF2, a key tumor suppressor gene of mesothelial carcinogenesis. Subclonal tumor populations were shared among paired biopsies, suggesting a polyclonal dissemination of the tumor. Transcriptome analysis highlighted dysregulation of cell adhesion and extracellular matrix pathways, linked to changes in histo-molecular gradient proportions between anatomic sites. Methylome analysis revealed contribution of epigenetic mechanisms in two patients. Finally, significant changes in the expression of immune mediators and genes related to immunological synapse, as well as differential infiltration of immune populations in the tumor environment were observed and led to a switch from a hot to a cold immune profile in three patients. Conclusions:This comprehensive analysis reveals patient-dependent spatial intra-tumor heterogeneity at the genetic, transcriptomic and epigenetic levels and in the immune landscape of the tumor microenvironment. Results support

  Study EGAS00001005328

• Integrated genomic analyses reveal molecular correlates of clinical response and resistance to atezolizumab in combination with bevacizumab in advanced hepatocellular carcinoma

  Atezolizumab (anti–programmed death-ligand 1 [anti-PD-L1]) and bevacizumab (anti-vascular endothelial growth factor [anti-VEGF]) combination therapy has become the new standard of care in patients with unresectable hepatocellular carcinoma. However, potential predictive biomarkers and mechanisms of response and resistance remain less understood. We report integrated molecular analyses of tumor samples from 358 patients enrolled in the GO30140 phase 1b or IMbrave150 phase 3 trials who were treated with atezolizumab combined with bevacizumab, atezolizumab alone, or sorafenib. Pre-existing immunity assessed by high expression of CD274 (PD-L1 mRNA), T-effector signature, and intratumoral CD8+ T cell density were associated with response and better clinical outcomes with the combination treatment. Less clinical benefit was associated with high regulatory T cell (Treg) to effector T-cell ratio and high expression of oncofetal genes such as GPC3 and AFP. Improved outcomes from the combination vs atezolizumab alone appeared to be associated with high expression of KDR (VEGFR2), Treg, and myeloid inflammation signatures, indicating anti-VEGF may contribute by targeting these axes. These findings were further validated by transcriptome analysis of paired pre- and post-treatment biopsies, in situ analysis by multiplex IHC and digital pathology, and in vivo study of an immunogenic hepatocellular carcinoma mouse model. Our findings highlight that anti-VEGF therapy might augment antitumor immunity and enhance anti–PD-L1 immunotherapy by targeting VEGF-mediated angiogenesis, Treg proliferation, and myeloid cell inflammation. Furthermore, our study identified candidate biomarkers for predicting response and resistance to anti–PD-L1 and anti-VEGF combination therapy.

  Study EGAS00001005503

• The clinical utility of genomics in childhood cancer extends beyond targetable mutations

  Survival of children with relapsed cancer remains unacceptably low. Cancer genomic profiling may identify new therapeutic options, enhance diagnostic accuracy, and improve understanding of tumors’ origins and continued evolution. We deeply sequenced 864 cancer-associated genes and the complete genomes and transcriptomes for 300 pediatric and adolescent/young adult (AYA) patients with poor prognosis or rare tumors. Using integrative somatic-germline analyses, we assessed the clinical utility of cancer genomics in this population. We analyzed tumor mutational signatures and burden to interrogate the relevance of germline variants in DNA repair genes and the impact of therapeutic exposures on tumor evolution. We also evaluated the frequency of changes to tumor drivers and therapeutic targets at relapse in those with serial samples. Integrative analysis yielded clinically actionable variants in a majority of patients. Improved diagnostic accuracy led to modified management in a subset. Therapeutically targetable variants were of unanticipated timing and type, with many of these derived from the germline. Enrichment in corroborating mutational signature 3 (‘BRCAness’) in patients with germline homologous recombination defects suggests they are drivers of pediatric cancers and potential targets for PARP inhibition. We show that mutation burden, which may in part be driven by prior therapeutic exposures, was significantly elevated in a minority of patients. Over one-third of patients with sequential samples analyzed showed a change in therapeutically-targetable drivers. Comprehensive somatic-germline cancer genomic profiling is useful at multiple points in the care trajectory for pediatric/AYA patients, supporting its integration into early clinical management. Re-biopsy at the time of relapse should be considered, and specific attention to mutational burden is indicated. Defective DNA repair genes in the germline may represent important targetable drivers.

  Study EGAS00001006034